Mean field model for spatially extended systems in the presence of multiplicative noise

We present a mean field model that describes the effect of multiplicative noise in spatially extended systems. The model can be solved analytically. For the case of the phi4 potential it predicts that the phase transition is shifted. This conclusion is supported by numerical simulations of this model in two dimensions.

Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

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.

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 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.

Mean-field model of interaction between bright vortex solitons in Bose-Einstein condensates

Using the explicit numerical solution of the axially symmetric Gross-Pitaevskii equation we study the dynamics of interaction among vortex solitons in a rotating matter-wave bright soliton train in a radially trapped and axially free Bose-Einstein condensate to understand certain features of the experiment by Strecker et al (2002 Nature 417 150). In a soliton train, solitons of opposite phase (phase δ = π) repel and stay apart without changing shape; solitons with δ = 0 attract, interact and coalesce, but eventually come out; solitons with a general δ usually repel but interact inelastically by exchanging matter. We study this and suggest future experiments with vortex solitons.

Mean-field model of jet formation in a collapsing Bose-Einstein condensate

We perform a systematic numerical study, based on the time-dependent Gross-Pitaevskii equation, of jet formation in collapsing and exploding Bose-Einstein condensates as in the experiment by Donley et al (2001 Nature 412 295). In the actual experiment, via a Feshbach resonance, the scattering length of atomic interaction was suddenly changed from positive to negative on a pre-formed condensate. Consequently, the condensate collapsed and ejected atoms via explosion. On disruption of the collapse by suddenly changing the scattering length to zero, a radial jet of atoms was formed in the experiment. We present a satisfactory account of jet formation under the experimental conditions and also make predictions beyond experimental conditions which can be verified in future experiments.

Experimental validation of a mean field model of mineralized collagen fiber arrays at two levels of hierarchy

In the course of this study, stiffness of a fibril array of mineralized collagen fibrils modeled with a mean field method was validated experimentally at site-matched two levels of tissue hierarchy using mineralized turkey leg tendons (MTLT). The applied modeling approaches allowed to model the properties of this unidirectional tissue from nanoscale (mineralized collagen fibrils) to macroscale (mineralized tendon). At the microlevel, the indentation moduli obtained with a mean field homogenization scheme were compared to the experimental ones obtained with microindentation. At the macrolevel, the macroscopic stiffness predicted with micro finite element (μFE) models was compared to the experimental stiffness measured with uniaxial tensile tests. Elastic properties of the elements in μFE models were injected from the mean field model or two-directional microindentations. Quantitatively, the indentation moduli can be properly predicted with the mean-field models. Local stiffness trends within specific tissue morphologies are very weak, suggesting additional factors responsible for the stiffness variations. At macrolevel, the μFE models underestimate the macroscopic stiffness, as compared to tensile tests, but the correlations are strong.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

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.

Origin of the neutron skin thickness of 208Pb in nuclear mean-field models

We study whether the neutron skin thickness Δrnp of 208Pb originates from the bulk or from the surface of the nucleon density distributions, according to the mean-field models of nuclear structure, and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to Δrnp arises from an extended sharp radius of neutrons, whereas the surface contribution arises from different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy is stiff, as typical of relativistic models, predict a bulk contribution in Δrnp of 208Pb about twice as large as the surface contribution. In contrast, models with a soft symmetry energy like common nonrelativistic models predict that Δrnp of 208Pb is divided similarly into bulk and surface parts. Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note that the linear correlation of Δrnp of 208Pb with the density derivative of the nuclear symmetry energy arises from the bulk part of Δrnp. We also note that most models predict a mixed-type (between halo and skin) neutron distribution for 208Pb. Although the halo-type limit is actually found in the models with a supersoft symmetry energy, the skin-type limit is not supported by any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of the 208Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating asymmetry in terms of the parameters that characterize the shape of the 208Pb nucleon densities.

Connecting mean field models of neural activity to EEG and fMRI data

Progress in functional neuroimaging of the brain increasingly relies on the integration of data from complementary imaging modalities in order to improve spatiotemporal resolution and interpretability. However, the usefulness of merely statistical combinations is limited, since neural signal sources differ between modalities and are related non-trivially. We demonstrate here that a mean field model of brain activity can simultaneously predict EEG and fMRI BOLD with proper signal generation and expression. Simulations are shown using a realistic head model based on structural MRI, which includes both dense short-range background connectivity and long-range specific connectivity between brain regions. The distribution of modeled neural masses is comparable to the spatial resolution of fMRI BOLD, and the temporal resolution of the modeled dynamics, importantly including activity conduction, matches the fastest known EEG phenomena. The creation of a cortical mean field model with anatomically sound geometry, extensive connectivity, and proper signal expression is an important first step towards the model-based integration of multimodal neuroimages.

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.

Study of pK Values and Effective Dielectric Constants of Ionizable Residues in Pentapeptides and in Staphylococcal Nuclease (SNase) Using a Mean-Field Approach

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