Source author record

Alessandro Torcini

Alessandro Torcini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2020arXiv

Theta-nested gamma oscillations in next generation neural mass models

Theta-nested gamma oscillations have been reported in many areas of the brain and are believed to represent a fundamental mechanism to transfer information across spatial and temporal scales. In a series of recent experiments in vitro it has been possible to replicate with an optogenetic theta frequency stimulation several features of cross-frequency coupling among theta and gamma rhythms observed in behaving animals. In order to reproduce the main findings of these experiments we have considered a new class of neural mass models able to reproduce exactly the macroscopic dynamics of spiking neural networks. In this framework, we have examined two set-ups able to support collective gamma oscillations: the pyramidal interneuronal network gamma (PING) and the interneuronal network gamma (ING). In both set-ups we observe the emergence of theta-nested gamma oscillations by driving the system with a sinusoidal theta-forcing in proximity of a Hopf bifurcation. These mixed rhythms display always phase amplitude coupling. However 2 different types of nested oscillations can be identified: one characterized by a perfect phase locking between theta and gamma rhythms, corresponding to an overall periodic behaviour; another one where the locking is imperfect and the dynamics is quasi-periodic or even chaotic. From our analysis it emerges that the locked states are more frequent in the ING set-up. In agreement with the experiments, we find theta-nested gamma oscillations for forcing frequencies in the range [1:10] Hz, whose amplitudes grow proportionally to the forcing one and which are clearly modulated by the theta phase. At variance with experimental findings, the gamma-power peak does not shift to higher frequencies by increasing the theta frequency. This effect can be obtained, in or model, only by incrementing, at the same time, also the noise or the forcing amplitude.

preprint2019arXiv

Coexistence of fast and slow gamma oscillations in one population of inhibitory spiking neurons

Oscillations are a hallmark of neural population activity in various brain regions with a spectrum covering a wide range of frequencies. Within this spectrum gamma oscillations have received particular attention due to their ubiquitous nature and to their correlation with higher brain functions. Recently, it has been reported that gamma oscillations in the hippocampus of behaving rodents are segregated in two distinct frequency bands: slow and fast. These two gamma rhythms correspond to dfferent states of the network, but their origin has been not yet clarified. Here, we show theoretically and numerically that a single inhibitory population can give rise to coexisting slow and fast gamma rhythms corresponding to collective oscillations of a balanced spiking network. The slow and fast gamma rhythms are generated via two different mechanisms: the fast one being driven by the coordinated tonic neural firing and the slow one by endogenous fluctuations due to irregular neural activity. We show that almost instantaneous stimulations can switch the collective gamma oscillations from slow to fast and vice versa. Furthermore, to make a closer contact with the experimental observations, we consider the modulation of the gamma rhythms induced by a slower (theta) rhythm driving the network dynamics. In this context, depending on the strength of the forcing, we observe phase-amplitude and phase-phase coupling between the fast and slow gamma oscillations and the theta forcing. Phase-phase coupling reveals different theta-phases preferences for the two coexisting gamma rhythms.

preprint2019arXiv

Cross frequency coupling in next generation inhibitory neural mass models

Coupling among neural rhythms is one of the most important mechanisms at the basis of cognitive processes in the brain. In this study we consider a neural mass model, rigorously obtained from the microscopic dynamics of an inhibitory spiking network with exponential synapses, able to autonomously generate collective oscillations (COs). These oscillations emerge via a super-critical Hopf bifurcation, and their frequencies are controlled by the synaptic time scale, the synaptic coupling and the excitability of the neural population. Furthermore, we show that two inhibitory populations in a master-slave configuration with different synaptic time scales can display various collective dynamical regimes: namely, damped oscillations towards a stable focus, periodic and quasi-periodic oscillations, and chaos. Finally, when bidirectionally coupled the two inhibitory populations can exhibit different types of theta-gamma cross-frequency couplings (CFCs): namely, phase-phase and phase-amplitude CFC. The coupling between theta and gamma COs is enhanced in presence of a external theta forcing, reminiscent of the type of modulation induced in Hippocampal and Cortex circuits via optogenetic drive.

preprint2019arXiv

Quantitative and qualitative analysis of asynchronous neural activity

The activity of a sparse network of leaky integrate-and-fire neurons is carefully revisited with reference to a regime of a bona-fide asynchronous dynamics. The study is preceded by a finite-size scaling analysis, carried out to identify a setup where collective synchronization is negligible. The comparison between quenched and annealed networks reveals the emergence of substantial differences when the coupling strength is increased, via a scenario somehow reminiscent of a phase transition. For sufficiently strong synaptic coupling, quenched networks exhibit a highly bursting neural activity, well reproduced by a self-consistent approach, based on the assumption that the input synaptic current is the superposition of independent renewal processes. The distribution of interspike intervals turns out to be relatively long-tailed; a crucial feature required for the self-sustainment of the bursting activity in a regime where neurons operate on average (much) below threshold. A semi-quantitative analogy with Ornstein-Uhlenbeck processes helps validating this interpretation. Finally, an alternative explanation in terms of Poisson processes is offered under the additional assumption of mutual correlations among excitatory and inhibitory spikes.

preprint2015arXiv

Cell assembly dynamics of sparsely-connected inhibitory networks: a simple model for the collective activity of striatal projection neurons

Striatal projection neurons form a sparsely-connected inhibitory network, and this arrangement may be essential for the appropriate temporal organization of behavior. Here we show that a simplified, sparse inhibitory network of Leaky-Integrate-and-Fire neurons can reproduce some key features of striatal population activity, as observed in brain slices [Carrillo-Reid et al., J. Neurophysiology 99 (2008) 1435{1450]. In particular we develop a new metric to determine the conditions under which sparse inhibitory networks form anti-correlated cell assemblies with time-varying activity of individual cells. We found that under these conditions the network displays an input-specific sequence of cell assembly switching, that effectively discriminates similar inputs. Our results support the proposal [Ponzi and Wickens, PLoS Comp Biol 9 (2013) e1002954] that GABAergic connections between striatal projection neurons allow stimulus-selective, temporally-extended sequential activation of cell assemblies. Furthermore, we help to show how altered intrastriatal GABAergic signaling may produce aberrant network-level information processing in disorders such as Parkinson's and Huntington's diseases.

preprint2015arXiv

Dynamics of fully coupled rotators with unimodal and bimodal frequency distribution

We analyze the synchronization transition of a globally coupled network of N phase oscillators with inertia (rotators) whose natural frequencies are unimodally or bimodally distributed. In the unimodal case, the system exhibits a discontinuous hysteretic transition from an incoherent to a partially synchronized (PS) state. For sufficiently large inertia, the system reveals the coexistence of a PS state and of a standing wave (SW) solution. In the bimodal case, the hysteretic synchronization transition involves several states. Namely, the system becomes coherent passing through traveling waves (TWs), SWs and finally arriving to a PS regime. The transition to the PS state from the SW occurs always at the same coupling, independently of the system size, while its value increases linearly with the inertia. On the other hand the critical coupling required to observe TWs and SWs increases with N suggesting that in the thermodynamic limit the transition from incoherence to PS will occur without any intermediate states. Finally a linear stability analysis reveals that the system is hysteretic not only at the level of macroscopic indicators, but also microscopically as verified by measuring the maximal Lyapunov exponent.

preprint2015arXiv

Dynamics versus energetics in phase separation

Phase separation may be driven by the minimization of a suitable free energy ${\cal F}$. This is the case, e.g., for diblock copolimer melts, where ${\cal F}$ is minimized by a steady periodic pattern whose wavelength $λ_{GS}$ depends on the segregation strength $α^{-1}$ and it is know since long time that in one spatial dimension $λ_{GS} \simeq α^{-1/3}$. Here we study in details the dynamics of the system in 1D for different initial conditions and by varying $α$ by five orders of magnitude. We find that, depending on the initial state, the final configuration may have a wavelength $λ_{D}$ with $λ_{min}(α)<λ_{D}<λ_{max}(α)$, where $λ_{min} \approx \ln (1/α)$ and $λ_{max}\approx α^{-1/2}$. In particular, if the initial state is homogeneous, the system exhibits a logarithmic coarsening process which arrests whenever $λ_{D}\approxλ_{min}$.

preprint2015arXiv

Intermittent chaotic chimeras for coupled rotators

Two symmetrically coupled populations of N oscillators with inertia $m$ display chaotic solutions with broken symmetry similar to experimental observations with mechanical pendula. In particular, we report the first evidence of intermittent chaotic chimeras, where one population is synchronized and the other jumps erratically between laminar and turbulent phases. These states have finite life-times diverging as a power-law with N and m. Lyapunov analyses reveal chaotic properties in quantitative agreement with theoretical predictions for globally coupled dissipative systems.

preprint2014arXiv

Clique of functional hubs orchestrates population bursts in developmentally regulated neural networks

It has recently been discovered that single neuron stimulation can impact network dynamics in immature and adult neuronal circuits. Here we report a novel mechanism which can explain in neuronal circuits, at an early stage of development, the peculiar role played by a few specific neurons in promoting/arresting the population activity. For this purpose, we consider a standard neuronal network model, with short-term synaptic plasticity, whose population activity is characterized by bursting behavior. The addition of developmentally inspired constraints and correlations in the distribution of the neuronal connectivities and excitabilities leads to the emergence of functional hub neurons, whose stimulation/deletion is critical for the network activity. Functional hubs form a clique, where a precise sequential activation of the neurons is essential to ignite collective events without any need for a specific topological architecture. Unsupervised time-lagged firings of supra-threshold cells, in connection with coordinated entrainments of near-threshold neurons, are the key ingredients to orchestrate

preprint2014arXiv

Hysteretic transitions in the Kuramoto model with inertia

We report finite size numerical investigations and mean field analysis of a Kuramoto model with inertia for fully coupled and diluted systems. In particular, we examine for a Gaussian distribution of the frequencies the transition from incoherence to coherence for increasingly large system size and inertia. For sufficiently large inertia the transition is hysteretic and within the hysteretic region clusters of locked oscillators of various sizes and different levels of synchronization coexist. A modification of the mean field theory developed by Tanaka, Lichtenberg, and Oishi [Physica D, 100 (1997) 279] allows to derive the synchronization profile associated to each of these clusters. We have also investigated numerically the limits of existence of the coherent and of the incoherent solutions. The minimal coupling required to observe the coherent state is largely independent of the system size and it saturates to a constant value already for moderately large inertia values. The incoherent state is observable up to a critical coupling whose value saturates for large inertia and for finite system sizes, while in the thermodinamic limit this critical value diverges proportionally to the mass. By increasing the inertia the transition becomes more complex, and the synchronization occurs via the emergence of clusters of whirling oscillators. The presence of these groups of coherently drifting oscillators induces oscillations in the order parameter. We have shown that the transition remains hysteretic even for randomly diluted networks up to a level of connectivity corresponding to few links per oscillator. Finally, an application to the Italian high-voltage power grid is reported, which reveals the emergence of quasi-periodic oscillations in the order parameter due to the simultaneous presence of many competing whirling clusters.

preprint2014arXiv

Sisyphus Effect in Pulse Coupled Excitatory Neural Networks with Spike-Timing Dependent Plasticity

The collective dynamics of excitatory pulse coupled neural networks with spike timing dependent plasticity (STDP) is studied. Depending on the model parameters stationary states characterized by High or Low Synchronization can be observed. In particular, at the transition between these two regimes, persistent irregular low frequency oscillations between strongly and weakly synchronized states are observable, which can be identified as infraslow oscillations with frequencies 0.02 - 0.03 Hz. Their emergence can be explained in terms of the Sisyphus Effect, a mechanism caused by a continuous feedback between the evolution of the coherent population activity and of the average synaptic weight. Due to this effect, the synaptic weights have oscillating equilibrium values, which prevents the neuronal population from relaxing into a stationary macroscopic state.

preprint2014arXiv

Spatio-temporal dynamics induced by competing instabilities in two asymmetrically coupled nonlinear evolution equations

Pattern formation often occurs in spatially extended physical, biological and chemical systems due to an instability of the homogeneous steady state. The type of the instability usually prescribes the resulting spatio-temporal patterns and their characteristic length scales. However, patterns resulting from the simultaneous occurrence of instabilities cannot be expected to be simple superposition of the patterns associated with the considered instabilities. To address this issue we design two simple models composed by two asymmetrically coupled equations of non-conserved (Swift-Hohenberg equations) or conserved (Cahn-Hilliard equations) order parameters with different characteristic wave lengths. The patterns arising in these systems range from coexisting static patterns of different wavelengths to traveling waves. A linear stability analysis allows to derive a two parameter phase diagram for the studied models, in particular revealing for the Swift-Hohenberg equations a co-dimension two bifurcation point of Turing and wave instability and a region of coexistence of stationary and traveling patterns. The nonlinear dynamics of the coupled evolution equations is investigated by performing accurate numerical simulations. These reveal more complex patterns, ranging from traveling waves with embedded Turing patterns domains to spatio-temporal chaos, and a wide hysteretic region, where waves or Turing patterns coexist. For the coupled Cahn-Hilliard equations the presence of an weak coupling is sufficient to arrest the coarsening process and to lead to the emergence of purely periodic patterns. The final states are characterized by domains with a characteristic length, which diverges logarithmically with the coupling amplitude.

preprint2014arXiv

Stable chaos in fluctuation driven neural circuits

We study the dynamical stability of pulse coupled networks of leaky integrate-and-fire neurons against infinitesimal and finite perturbations. In particular, we compare current versus fluctuations driven networks, the former (latter) is realized by considering purely excitatory (inhibitory) sparse neural circuits. In the excitatory case the instabilities of the system can be completely captured by an usual linear stability (Lyapunov) analysis, on the other hand the inhibitory networks can display the coexistence of linear and nonlinear instabilities. The nonlinear effects are associated to finite amplitude instabilities, which have been characterized in terms of suitable indicators. For inhibitory coupling one observes a transition from chaotic to non chaotic dynamics by decreasing the pulse width. For sufficiently fast synapses the system, despite showing an erratic evolution, is linearly stable, thus representing a prototypical example of Stable Chaos.

preprint2013arXiv

Emergence of slow collective oscillations in neural networks with spike timing dependent plasticity

The collective dynamics of excitatory pulse coupled neurons with spike timing dependent plasticity (STDP) is studied. The introduction of STDP induces persistent irregular oscillations between strongly and weakly synchronized states, reminiscent of brain activity during slow-wave sleep. We explain the oscillations by a mechanism, the Sisyphus Effect, caused by a continuous feedback between the synaptic adjustments and the coherence in the neural firing. Due to this effect, the synaptic weights have oscillating equilibrium values, and this prevents the system from relaxing into a stationary macroscopic state.

preprint2013arXiv

Linear stability in networks of pulse-coupled neurons

In a first step towards the comprehension of neural activity, one should focus on the stability of the various dynamical states. Even the characterization of idealized regimes, such as a perfectly periodic spiking activity, reveals unexpected difficulties. In this paper we discuss a general approach to linear stability of pulse-coupled neural networks for generic phase-response curves and post-synaptic response functions. In particular, we present: (i) a mean-field approach developed under the hypothesis of an infinite network and small synaptic conductances; (ii) a "microscopic" approach which applies to finite but large networks. As a result, we find that no matter how large is a neural network, its response to most of the perturbations depends on the system size. There exists, however, also a second class of perturbations, whose evolution typically covers an increasingly wide range of time scales. The analysis of perfectly regular, asynchronous, states reveals that their stability depends crucially on the smoothness of both the phase-response curve and the transmitted post-synaptic pulse. The general validity of this scenarion is confirmed by numerical simulations of systems that are not amenable to a perturbative approach.

preprint2012arXiv

Chaos in the Hamiltonian mean field model

We study the dynamical properties of the canonical ordered phase of the Hamiltonian mean-field (HMF) model, in which $N$ particles, globally-coupled via pairwise attractive interactions, form a rotating cluster. Using a combination of numerical and analytical arguments, we first show that the largest Lyapunov exponent remains strictly positive in the infinite-size limit, converging to its asymptotic value with $1/\ln N$ corrections. We then elucidate the scaling laws ruling the behavior of this asymptotic value in the critical region separating the ordered, clustered phase and the disordered phase present at high energy densities. We also show that the full spectrum of Lyapunov exponents consists of a bulk component converging to the (zero) value taken by a test oscillator forced by the mean field, plus subextensive bands of ${\cal O}(\ln N)$ exponents taking finite values. We finally investigate the robustness of these results by studying a "2D" extension of the HMF model where each particle is endowed with 4 degrees of freedom, thus allowing the emergence of chaos at the level of single particle. Altogether, these results illustrate the subtle effects of global (or long-range) coupling, and the importance of the order in which the infinite-time and infinite-size limits are taken: for an infinite-size HMF system represented by the Vlasov equation, no chaos is present, while chaos exists and subsists for any finite system size.

preprint2012arXiv

Collective dynamics in sparse networks

The microscopic and macroscopic dynamics of random networks is investigated in the strong-dilution limit (i.e. for sparse networks). By simulating chaotic maps, Stuart-Landau oscillators, and leaky integrate-and-fire neurons, we show that a finite connectivity (of the order of a few tens) is able to sustain a nontrivial collective dynamics even in the thermodynamic limit. Although the network structure implies a non-additive dynamics, the microscopic evolution is extensive (i.e. the number of active degrees of freedom is proportional to the number of network elements).

preprint2012arXiv

Convective Lyapunov Spectra

We generalize the concept of convective (or velocity-dependent) Lyapunov exponent $Λ(v)$ to an entire spectrum $Λ(v,n)$. Our results are supported by the consistency between the outcome of the chronotopic approach [{\it S. Lepri et al. J. Stat. Phys., 82 5/6 (1996) 1429}] and a more direct method. There exists a critical integrated density $n=n_c$, beyond which the convective exponent exhibits a discontinuous dependence on the velocity, which originates from the appearance of multiple branches. This phenomenon can be traced back to a change of concavity of the so-called {\it temporal} Lyapunov spectrum for $n>n_c$, which is therefore a dynamical invariant.

preprint2012arXiv

Splay states in finite pulse-coupled networks of excitable neurons

The emergence and stability of splay states is studied in fully coupled finite networks of N excitable quadratic integrate-and-fire neurons, connected via synapses modeled as pulses of finite amplitude and duration. For such synapses, by introducing two distinct types of synaptic events (pulse emission and termination), we were able to write down an exact event-driven map for the system and to evaluate the splay state solutions. For M overlapping post synaptic potentials the linear stability analysis of the splay state should take in account, besides the actual values of the membrane potentials, also the firing times associated to the M previous pulse emissions. As a matter of fact, it was possible, by introducing M complementary variables, to rephrase the evolution of the network as an event-driven map and to derive an analytic expression for the Floquet spectrum. We find that, independently of M, the splay state is marginally stable with N-2 neutral directions. Furthermore, we have identified a family of periodic solutions surrounding the splay state and sharing the same neutral stability directions. In the limit of $δ$-pulses, it is still possible to derive an event-driven formulation for the dynamics, however the number of neutrally stable directions, associated to the splay state, becomes N. Finally, we prove a link between the results for our system and a previous theory [Watanabe and Strogatz, Physica D, 74 (1994), pp. 197- 253] developed for networks of phase oscillators with sinusoidal coupling.

preprint2012arXiv

Stability of the splay state in networks of pulse-coupled neurons

We analytically investigate the stability of {\it splay states} in networks of $N$ pulse-coupled phase-like models of neurons. By developing a perturbative technique, we find that, in the limit of large $N$, the Floquet spectrum scales as $1/N^2$ for generic discontinuous velocity fields. Moreover, the stability of the so-called short-wavelength component is determined by the sign of the jump at the discontinuity. Altogether, the form of the spectrum depends on the pulse shape but is independent of the velocity field.

preprint2011arXiv

Coherent periodic activity in excitatory Erdos-Renyi neural networks:The role of network connectivity

We consider an excitatory random network of leaky integrate-and-fire pulse coupled neurons. The neurons are connected as in a directed Erdös-Renyi graph with average connectivity $<k>$ scaling as a power law with the number of neurons in the network. The scaling is controlled by a parameter $γ$, which allows to pass from massively connected to sparse networks and therefore to modify the topology of the system. At a macroscopic level we observe two distinct dynamical phases: an Asynchronous State (AS) corresponding to a desynchronized dynamics of the neurons and a Partial Synchronization (PS) regime associated with a coherent periodic activity of the network. At low connectivity the system is in an AS, while PS emerges above a certain critical average connectivity $<k>_c$. For sufficiently large networks, $<k>_c$ saturates to a constant value suggesting that a minimal average connectivity is sufficient to observe coherent activity in systems of any size irrespectively of the kind of considered network: sparse or massively connected. However, this value depends on the nature of the synapses: reliable or unreliable. For unreliable synapses the critical value required to observe the onset of macroscopic behaviors is noticeably smaller than for reliable synaptic transmission. Due to the disorder present in the system, for finite number of neurons we have inhomogeneities in the neuronal behaviors, inducing a weak form of chaos, which vanishes in the thermodynamic limit. In such a limit the disordered systems exhibit regular (non chaotic) dynamics and their properties correspond to that of a homogeneous fully connected network for any $γ$-value. Apart for the peculiar exception of sparse networks, which remain intrinsically inhomogeneous at any system size.

preprint2011arXiv

Extensive and Sub-Extensive Chaos in Globally-Coupled Dynamical Systems

Using a combination of analytical and numerical techniques, we show that chaos in globally-coupled identical dynamical systems, be they dissipative or Hamiltonian, is both extensive and sub-extensive: their spectrum of Lyapunov exponents is asymptotically flat (thus extensive) at the value $λ_0$ given by a single unit forced by the mean-field, but sandwiched between sub-extensive bands containing typically $\mathcal{O}(\log N)$ exponents whose values vary as $λ\simeq λ_\infty + c/\log N$ with $λ_\infty \neq λ_0$.

preprint2010arXiv

Collective chaos in pulse-coupled neural networks

We study the dynamics of two symmetrically coupled populations of identical leaky integrate-and-fire neurons characterized by an excitatory coupling. Upon varying the coupling strength, we find symmetry-breaking transitions that lead to the onset of various chimera states as well as to a new regime, where the two populations are characterized by a different degree of synchronization. Symmetric collective states of increasing dynamical complexity are also observed. The computation of the the finite-amplitude Lyapunov exponent allows us to establish the chaoticity of the (collective) dynamics in a finite region of the phase plane. The further numerical study of the standard Lyapunov spectrum reveals the presence of several positive exponents, indicating that the microscopic dynamics is high-dimensional.

preprint2010arXiv

Collective oscillations in disordered neural networks

We investigate the onset of collective oscillations in a network of pulse-coupled leaky-integrate-and-fire neurons in the presence of quenched and annealed disorder. We find that the disorder induces a weak form of chaos that is analogous to that arising in the Kuramoto model for a finite number N of oscillators [O.V. Popovych at al., Phys. Rev. E 71} 065201(R) (2005)]. In fact, the maximum Lyapunov exponent turns out to scale to zero for N going to infinite, with an exponent that is different for the two types of disorder. In the thermodynamic limit, the random-network dynamics reduces to that of a fully homogenous system with a suitably scaled coupling strength. Moreover, we show that the Lyapunov spectrum of the periodically collective state scales to zero as 1/N^2, analogously to the scaling found for the `splay state'.

preprint2009arXiv

Stable chaos

Stable chaos is a generalization of the chaotic behaviour exhibited by cellular automata to continuous-variable systems and it owes its name to an underlying irregular and yet linearly stable dynamics. In this review we discuss analogies and differences with the usual deterministic chaos and introduce several tools for its characterization. Some examples of transitions from ordered behavior to stable chaos are also analyzed to further clarify the underlying dynamical properties. Finally, two models are specifically discussed: the diatomic hard-point gas chain and a network of globally coupled neurons.

preprint2009arXiv

Unfolding times for proteins in a force clamp

The escape process from the native valley for proteins subjected to a constant stretching force is examined using a model for a Beta-barrel. For a wide range of forces, the unfolding dynamics can be treated as one-dimensional diffusion, parametrized in terms of the end-to-end distance. In particular, the escape times can be evaluated as first passage times for a Brownian particle moving on the protein free-energy landscape, using the Smoluchowski equation. At strong forces, the unfolding process can be viewed as a diffusive drift away from the native state, while at weak forces thermal activation is the relevant mechanism. An escape-time analysis within this approach reveals a crossover from an exponential to an inverse Gaussian escape-time distribution upon passing from weak to strong forces. Moreover, a single expression valid at weak and strong forces can be devised both for the average unfolding time as well as for the corresponding variance. The analysis offers a possible explanation of recent experimental findings for ddFLN4 and ubiquitin.

preprint1998arXiv

Analytical Estimation of the Maximal lyapunov Exponent in Oscillator Chains

An analytical expression for the maximal Lyapunov exponent $λ_1$ in generalized Fermi-Pasta-Ulam oscillator chains is obtained. The derivation is based on the calculation of modulational instability growth rates for some unstable periodic orbits. The result is compared with numerical simulations and the agreement is good over a wide range of energy densities $ε$. At very high energy density the power law scaling of $λ_1$ with $ε$ can be also obtained by simple dimensional arguments, assuming that the system is ruled by a single time scale. Finally, we argue that for repulsive and hard core potentials in one dimension $λ_1 \sim \sqrtε$ at large $ε$.

preprint1998arXiv

Anomalous diffusion as a signature of collapsing phase in two dimensional self-gravitating systems

A two dimensional self-gravitating Hamiltonian model made by $N$ fully-coupled classical particles exhibits a transition from a collapsing phase (CP) at low energy to a homogeneous phase (HP) at high energy. From a dynamical point of view, the two phases are characterized by two distinct single-particle motions : namely, superdiffusive in the CP and ballistic in the HP. Anomalous diffusion is observed up to a time $τ$ that increases linearly with $N$. Therefore, the finite particle number acts like a white noise source for the system, inhibiting anomalous transport at longer times.

preprint1998arXiv

Equilibrium and dynamical properties of two dimensional self-gravitating systems

A system of N classical particles in a 2D periodic cell interacting via long-range attractive potential is studied. For low energy density $U$ a collapsed phase is identified, while in the high energy limit the particles are homogeneously distributed. A phase transition from the collapsed to the homogeneous state occurs at critical energy U_c. A theoretical analysis within the canonical ensemble identifies such a transition as first order. But microcanonical simulations reveal a negative specific heat regime near $U_c$. The dynamical behaviour of the system is affected by this transition : below U_c anomalous diffusion is observed, while for U > U_c the motion of the particles is almost ballistic. In the collapsed phase, finite $N$-effects act like a noise source of variance O(1/N), that restores normal diffusion on a time scale diverging with N. As a consequence, the asymptotic diffusion coefficient will also diverge algebraically with N and superdiffusion will be observable at any time in the limit N \to \infty. A Lyapunov analysis reveals that for U > U_c the maximal exponent λdecreases proportionally to N^{-1/3} and vanishes in the mean-field limit. For sufficiently small energy, in spite of a clear non ergodicity of the system, a common scaling law λ\propto U^{1/2} is observed for any initial conditions.

preprint1997arXiv

Localization and Equipartition of Energy in the beta-FPU Chain : Chaotic Breathers

The evolution towards equipartition in the $β$-FPU chain is studied considering as initial condition the highest frequency mode. Above an analytically derived energy threshold, this zone-boundary mode is shown to be modulationally unstable and to give rise to a striking localization process. The spontaneously created excitations have strong similarity with moving exact breathers solutions. But they have a finite lifetime and their dynamics is chaotic. These chaotic breathers are able to collect very efficiently the energy in the chain. Therefore their size grows in time and they can transport a very large quantity of energy. These features can be explained analyzing the dynamics of perturbed exact breathers of the FPU chain. In particular, a close connection between the Lyapunov spectrum of the chaotic breathers and the Floquet spectrum of the exact ones has been found. The emergence of chaotic breathers is convincingly explained by the absorption of high frequency phonons whereas a breather's metastability is for the first time identified. The lifetime of the chaotic breather is related to the time necessary for the system to reach equipartition. The equipartition time turns out to be dependent on the system energy density $ε$ only. Moreover, such time diverges as $ε^{-2}$ in the limit $ε\to 0$ and vanishes as $ε^{-1/4}$ for $ε\to \infty$.