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.

Averaging, slaving and balance dynamics in a simple atmospheric model

We report numerical results from a study of balance dynamics using a simple model of atmospheric motion that is designed to help address the question of why balance dynamics is so stable. The non-autonomous Hamiltonian model has a chaotic slow degree of freedom (representing vortical modes) coupled to one or two linear fast oscillators (representing inertia-gravity waves). The system is said to be balanced when the fast and slow degrees of freedom are separated. We find adiabatic invariants that drift slowly in time. This drift is consistent with a random-walk behaviour at a speed which qualitatively scales, even for modest time scale separations, as the upper bound given by Neishtadt’s and Nekhoroshev’s theorems. Moreover, a similar type of scaling is observed for solutions obtained using a singular perturbation (‘slaving’) technique in resonant cases where Nekhoroshev’s theorem does not apply. We present evidence that the smaller Lyapunov exponents of the system scale exponentially as well. The results suggest that the observed stability of nearly-slow motion is a consequence of the approximate adiabatic invariance of the fast motion.

Bubbling transition to spatial mode excitation in an extended dynamical system

We investigated the transition to spatio-temporal chaos in spatially extended nonlinear dynamical systems possessing an invariant subspace with a low-dimensional attractor. When the latter is chaotic and the subspace is transversely stable we have a spatially homogeneous state only. The onset of spatio-temporal chaos, i.e. the excitation of spatially inhomogeneous modes, occur through the loss of transversal stability of some unstable periodic orbit embedded in the chaotic attractor lying in the invariant subspace. This is a bubbling transition, since there is a switching between spatially homogeneous and nonhomogeneous states with statistical properties of on-off intermittency. Hence the onset of spatio-temporal chaos depends critically both on the existence of a chaotic attractor in the invariant subspace and its being transversely stable or unstable. (C) 2008 Elsevier B.V. All rights reserved.

Spatiotemporal Pattern Formation in the Oscillatory Electro-Oxidation of Sulfide on a Platinum Disk

Spatiotemporal pattern formation in the electrocatalytic oxidation of sulfide on a platinum disk is investigated using electrochemical methods and a charge-coupled device (CCD) camera simultaneously. The system is characterized by different oscillatory regions spread over a wide potential range. An additional series resistor and a large electrode area facilitate observation of multiple regions of kinetic instabilities along the current/potential curve. Spatiotemporal patterns on the working electrode, such as fronts, pulses, spirals, twinkling eyes, labyrinthine stripes, and alternating synchronized deposition and dissolution, are observed at different operating conditions of series resistance and sweep rate.

Time-periodic perturbation of a Lienard equation with an unbounded homoclinic loop

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Global dynamics of the Rikitake system

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Continuously deformable topological structure

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Using a one-dimensional lattice applied to the thermodynamic study of DNA

In this work it is analyzed a one-dimensional lattice which is composed by mass-spring systems with one additional Rosen-Morse potential on site. This kind of lattice is used to study thermodynamic properties of DNA, especially its thermal denaturation. on the context of this work, the Rosen-Morse potential simulates hydrogen bonds between double strands of the molecule. From the graphic of the average stretching of base pairs versus temperature it is possible to observe the thermal denaturation of the system. This result shows that it is possible to obtain phase transition with an asymmetric potential without an infinite barrier.

Stability of breathers in simple mechanical models for DNA

This work studies through the Floquet theory the stability of breathers generated by the anti-continuous limit. We used the Peyrard-Bishop model for DNA and two kinds of nonlinear potential: the Morse potential and a potential with a hump. The comparison of their stability was done in function of the coupling parameter. We also investigate the dynamic behaviour of the system in stable and unstable regions. Qualitatively, the dynamic of mobile breathers resembles DNA.

Regularization and singular perturbation techniques for non-smooth systems

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Theory of small aspect ratio waves in deep water

In the limit of small values of the aspect ratio parameter (or wave steepness) which measures the amplitude of a surface wave in units of its wave-length, a model equation is derived from the Euler system in infinite depth (deep water) without potential flow assumption. The resulting equation is shown to sustain periodic waves which on the one side tend to the proper linear limit at small amplitudes, on the other side possess a threshold amplitude where wave crest peaking is achieved. An explicit expression of the crest angle at wave breaking is found in terms of the wave velocity. By numerical simulations, stable soliton-like solutions (experiencing elastic interactions) propagate in a given velocities range on the edge of which they tend to the peakon solution. (c) 2005 Elsevier B.V. All rights reserved.

A mathematical model for wave propagation in elastic tubes with inhomogeneities: Application to blood waves propagation

We study the propagation of waves in an elastic tube filled with an inviscid fluid. We consider the case of inhomogeneity whose mechanical and geometrical properties vary in space. We deduce a system of equations of the Boussinesq type as describing the wave propagation in the tube. Numerical simulations of these equations show that inhomogeneities prevent separation of right-going from left-going waves. Then reflected and transmitted coefficients are obtained in the case of localized constriction and localized rigidity. Next we focus on wavetrains incident on various types of anomalous regions. We show that the existence of anomalous regions modifies the wavetrain patterns. (c) 2007 Elsevier B.V. All rights reserved.

Evolution equation for short surface waves on water of finite depth

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

Gap solitons in a model of a superfluid fermion gas in optical lattices

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