Dynamical complexity in a mean-field model of human EEG

A recently proposed mean-field theory of mammalian cortex rhythmogenesis describes the salient features of electrical activity in the cerebral macrocolumn, with the use of inhibitory and excitatory neuronal populations (Liley et al 2002). This model is capable of producing a range of important human EEG (electroencephalogram) features such as the alpha rhythm, the 40 Hz activity thought to be associated with conscious awareness (Bojak & Liley 2007) and the changes in EEG spectral power associated with general anesthetic effect (Bojak & Liley 2005). From the point of view of nonlinear dynamics, the model entails a vast parameter space within which multistability, pseudoperiodic regimes, various routes to chaos, fat fractals and rich bifurcation scenarios occur for physiologically relevant parameter values (van Veen & Liley 2006). The origin and the character of this complex behaviour, and its relevance for EEG activity will be illustrated. The existence of short-lived unstable brain states will also be discussed in terms of the available theoretical and experimental results. A perspective on future analysis will conclude the presentation.

Electrorhythmogenesis and anaesthesia in a physiological mean field theory

We solve eight partial-differential, two-dimensional, nonlinear mean field equations, which describe the dynamics of large populations of cortical neurons. Linearized versions of these equations have been used to generate the strong resonances observed in the human EEG, in particular the α-rhythm (8–), with physiologically plausible parameters. We extend these results here by numerically solving the full equations on a cortex of realistic size, which receives appropriately “colored” noise as extra-cortical input. A brief summary of the numerical methods is provided. As an outlook to future applications, we explain how the effects of GABA-enhancing general anaesthetics can be simulated and present first results.

Relating pseudospin and spin symmetries through charge conjugation and chiral transformations: the case of the relativistic harmonic oscillator

We solve the generalized relativistic harmonic oscillator in 1+1 dimensions, i.e., including a linear pseudoscalar potential and quadratic scalar and vector potentials which have equal or opposite signs. We consider positive and negative quadratic potentials and discuss in detail their bound-state solutions for fermions and antifermions. The main features of these bound states are the same as the ones of the generalized three-dimensional relativistic harmonic oscillator bound states. The solutions found for zero pseudoscalar potential are related to the spin and pseudospin symmetry of the Dirac equation in 3+1 dimensions. We show how the charge conjugation and gamma(5) chiral transformations relate the several spectra obtained and find that for massless particles the spin and pseudospin symmetry-related problems have the same spectrum but different spinor solutions. Finally, we establish a relation of the solutions found with single-particle states of nuclei described by relativistic mean-field theories with scalar, vector, and isoscalar tensor interactions and discuss the conditions in which one may have both nucleon and antinucleon bound states.

pK determination. A mean field, Poisson-Boltzmann approach

A model describing dissociation of monoprotonic acid and a method for the determination of its pK value are presented. The model is based on a mean field approximation. The Poisson-Boltzmann equation, adopting spherical symmetry, is numerically solved, and the solution of its linearized form is written. By use of the pH values of a dilution experiment of galacturonic acid as the entry data, the proposed method allowed estimation of the value of pK = 3.25 at a temperature of 25 degrees C. Values for the complex dimensions and dissociation degree are calculated using experimental pH values for solution concentration values ranging from 0.1 to 60 mM. The present analysis leads to the conclusion that the Poisson-Boltzmann equation or its linear form is equally suited for the description of such systems.

Mean-field model for Josephson oscillation in a Bose-Einstein condensate on an one-dimensional optical trap

Using the axially-symmetric time-dependent Gross-Pitaevskii equation we study the phase coherence in a repulsive Bose-Einstein condensate (BEC) trapped by a harmonic and an one-dimensional optical lattice potential to describe the experiment by Cataliotti et al. on atomic Josephson oscillation [Science 293, 843 (2001)]. The phase coherence is maintained after the BEC is set into oscillation by a small displacement of the magnetic trap along the optical lattice. The phase coherence in the presence of oscillating neutral current across an array of Josephson junctions manifests in an interference pattern formed upon free expansion of the BEC. The numerical response of the system to a large displacement of the magnetic trap is a classical transition from a coherent superfluid to an insulator regime and a subsequent destruction of the interference pattern in agreement With the more recent experiment by Cataliotti et al. [New J. Phys. 5, 71 (2003)].

Mean-field description of a dynamical collapse of a fermionic condensate in a trapped boson-fermion mixture

We suggest a time-dependent dynamical mean-field-hydrodynamic model for the collapse of a trapped boson-fermion condensate and perform numerical simulation based on it to understand some aspects of the experiment by Modugno et al. [Science 297, 2240 (2002)] on the collapse of the fermionic condensate in the K-40-Rb-87 mixture. We show that the mean-field model explains the formation of a stationary boson-fermion condensate at zero temperature with relative sizes compatible with experiment. This model is also found to yield a faithful representation of the collapse dynamics in qualitative agreement with experiment. In particular we consider the collapse of the fermionic condensate associated with (a) an increase of the number of bosonic atoms as in the experiment and (b) an increase of the attractive boson-fermion interaction using a Feshbach resonance. Suggestion for experiments of fermionic collapse using a Feshbach resonance is made.

Limits of validity for a semiclassical mean-field two-fluid model for Bose-Einstein condensation thermodynamics

We reinvestigate the Bose-Einstein condensation (BEC) thermodynamics of a weakly interacting dilute Bose gas under the action of a trap using a semi-classical two-fluid mean-field model in order to find the domain of applicability of the model. Such a model is expected to break down once the condition of diluteness and weak interaction is violated. We find that this breakdown happens for values of coupling and density near the present experimental scenario of BEG. With the increase of the interaction coupling and density the model may lead to unphysical results for thermodynamic observables. (C) 2000 Published by Elsevier B.V. B.V, All rights reserved.

Mean-field model for the interference of matter-waves from a three-dimensional optical trap

Using the mean-field time-dependent Gross-Pitaevskii equation we study the formation of a repulsive Bose-Einstein condensate on a combined optical and harmonic traps in two and three dimensions and subsequent generation of the interference pattern upon the removal of the combined traps as in the experiment by, Greiner et al. [Nature (London 415 (2002) 39]. For optical traps of moderate strength, interference pattern of 27 (9) prominent bright spots is found to be formed in three. (two) dimensions on a cubic (square) lattice in agreement with experiment. Similar interference pattern can also be formed upon removal of the optical lattice trap only. The pattern so formed can oscillate for a long time in the harmonic trap which can be observed experimentally. (C) 2003 Elsevier B.V. B.V. All rights reserved.

Dynamical mean-field study of strongly interacting Bose-Einstein condensate

The experimental results of Rb-85 Bose-Einstein condensates are analyzed within the mean-field approximation with time-dependent two-body interaction and dissipation due to three-body recombination. We found that the magnitude of the dissipation is consistent with the three-body theory for longer rise times. However, for shorter rise times, it occurs an enhancement of this parameter, consistent with a coherent dimer formation. (C) 2004 Elsevier B.V. All rights reserved.

Mean-field equations for cigar- and disc-shaped Bose and Fermi superfluids

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

QUARK MEAN-FIELD THEORY AND CONSISTENCY WITH NUCLEAR-MATTER

1/N(c) expansion in QCD (with N(c) the number of colors) suggests using a potential from meson sector (e.g., Richardson) for baryons. For light quarks a sigma-field has to be introduced to ensure chiral symmetry breaking (chi-SB). It is found that nuclear matter properties can be used to pin down the chi-SB modeling. All masses, M(N), m-sigma, m-omega, are found to scale with density. The equations are solved self-consistently.

Self-bound droplet of Bose and Fermi atoms in one dimension: Collective properties in mean-field and Tonks-Girardeau regimes

We investigate a dilute mixture of bosons and spin-polarized fermions in one dimension. With an attractive Bose-Fermi scattering length the ground state is a self-bound droplet, i.e., a Bose-Fermi bright soliton where the Bose and Fermi clouds are superimposed. We find that the quantum fluctuations stabilize the Bose-Fermi soliton such that the one-dimensional bright soliton exists for any finite attractive Bose-Fermi scattering length. We study density profile and collective excitations of the atomic bright soliton showing that they depend on the bosonic regime involved: mean-field or Tonks-Girardeau.

QCD phase transitions from relativistic hadron models

The models of translationally invariant infinite nuclear matter in the relativistic mean field models are very interesting and simple, since the nucleon can connect only to a constant vector and scalar meson field. Can one connect these to the complicated phase transitions of QCD? For an affirmative answer to this question, one must consider models where the coupling contstants to the scalar and vector fields depend on density in a nonlinear way, since as such the models are not explicitly chirally invariant. Once this is ensured, indeed one can derive a quark condensate indirectly from the energy density of nuclear matter which goes to zero at large density and temperature. The change to zero condensate indicates a smooth phase transition. © Springer-Verlag 1996.

Relativistic RPA for isobaric analogue and Gamow-Teller resonances in closed shell nuclei

We perform a self-consistent relativistic RPA calculation for the isobaric analogue and Gamow-Teller resonances based on relativistic mean field theory results for the ground states of 48Ca, 90Zr and 208Pb. We use the parameter set NL1 for the σ, ω and ρ mesons, and experimental values for the pion and nucleon. An extra parameter, related to the intensity of the contact term in the pion-exchange interaction, is crucial to reproduce the latter resonances. © 1998 Published by Elsevier Science B.V. All rights reserved.

Mean-field description of collapsing and exploding Bose-Einstein condensates

Numerical simulations based on the time-dependent mean-field Gross-Pitaevskii equation was performed to explain the dynamics of collapsing and exploding Bose-Einstein condensates (BEC) of 85Rb atoms. The atomic interaction was manipulated by an external magnetic field via a Feshbach resonance. On changing the scattering length of atomic interaction from a positive to a large negative value, the condensate collapsed and ejected atoms via explosion.