-
Discrete Laplacian thermostat for flocks and swarms: the fully conserved Inertial Spin Model
Authors:
Andrea Cavagna,
Javier Cristín,
Irene Giardina,
Tomas S. Grigera,
Mario Veca
Abstract:
Experiments on bird flocks and midge swarms reveal that these natural systems are well described by an active theory in which conservation laws play a crucial role. By building a symplectic structure that couples the particles' velocities to the generator of their internal rotations (spin), the Inertial Spin Model (ISM) reinstates a second-order temporal dynamics that captures many phenomenologica…
▽ More
Experiments on bird flocks and midge swarms reveal that these natural systems are well described by an active theory in which conservation laws play a crucial role. By building a symplectic structure that couples the particles' velocities to the generator of their internal rotations (spin), the Inertial Spin Model (ISM) reinstates a second-order temporal dynamics that captures many phenomenological traits of flocks and swarms. The reversible structure of the ISM predicts that the total spin is a constant of motion, the central conservation law responsible for all the novel dynamical features of the model. However, fluctuations and dissipation introduced in the original model to make it relax, violate the spin conservation law, so that the ISM aligns with the biophysical phenomenology only within finite-size regimes, beyond which the overdamped dynamics characteristic of the Vicsek model takes over. Here, we introduce a novel version of the ISM, in which the irreversible terms needed to relax the dynamics strictly respect the conservation of the spin. We perform a numerical investigation of the fully conservative model, exploring both the fixed-network case, which belongs to the equilibrium class of Model G, and the active case, characterized by self-propulsion of the agents and an out-of-equilibrium reshuffling of the underlying interaction network. Our simulations not only capture the correct spin wave phenomenology of the ordered phase, but they also yield dynamical critical exponents in the near-ordering phase that agree very well with the theoretical predictions.
△ Less
Submitted 8 October, 2024; v1 submitted 12 March, 2024;
originally announced March 2024.
-
Mean-field solution of the neural dynamics in a Greenberg-Hastings model with excitatory and inhibitory units
Authors:
Joaquin Almeira,
Tomas S. Grigera,
Daniel A. Martin,
Dante R. Chialvo,
Sergio A. Cannas
Abstract:
We present a mean field solution of the dynamics of a Greenberg-Hastings neural network with both excitatory and inhibitory units. We analyse the dynamical phase transitions that appear in the stationary state as the model parameters are varied. Analytical solutions are compared with numerical simulations of the microscopic model defined on a fully connected network. We found that the stationary s…
▽ More
We present a mean field solution of the dynamics of a Greenberg-Hastings neural network with both excitatory and inhibitory units. We analyse the dynamical phase transitions that appear in the stationary state as the model parameters are varied. Analytical solutions are compared with numerical simulations of the microscopic model defined on a fully connected network. We found that the stationary state of this system exhibits a first order dynamical phase transition (with the associated hysteresis) when the fraction of inhibitory units $f< f_t \leq 1/2$, even for a finite system. In finite systems, when $f > f_t$ the first order transition is replaced by a pseudo critical one, namely a continuous crossover between regions of low and high activity that resembles the finite size behaviour of a continuous phase transition order parameter. However, in the thermodynamic limit, we found that $f_t\to 1/2$ and the activity for $f\geq f_t$ becomes negligible for any value of $T>0$, while the first order transition for $f<f_t$ remains.
△ Less
Submitted 29 December, 2023;
originally announced December 2023.
-
The Inertial Spin Model of flocking with position-dependent forces
Authors:
Sebastián Carruitero,
Alejo Costa Duran,
Giulia Pisegna,
Mauricio B. Sturla,
Tomás S. Grigera
Abstract:
We consider the ISM of flocking and swarming. The model has been introduced to explain certain dynamic features of swarming (second sound, a lower than expected dynamic critical exponent) while preserving the mechanism for onset of order provided by the Vicsek model. The ISM has only been formulated with an imitation (''ferromagnetic'') interaction between velocities. Here we show how to add posit…
▽ More
We consider the ISM of flocking and swarming. The model has been introduced to explain certain dynamic features of swarming (second sound, a lower than expected dynamic critical exponent) while preserving the mechanism for onset of order provided by the Vicsek model. The ISM has only been formulated with an imitation (''ferromagnetic'') interaction between velocities. Here we show how to add position-dependent forces in the model, which allows to consider effects such as cohesion, excluded volume, confinement and perturbation with external position-dependent field. We study numerically a one-particle case.
△ Less
Submitted 20 December, 2023; v1 submitted 11 December, 2023;
originally announced December 2023.
-
Tricritical behavior in a neural model with excitatory and inhibitory units
Authors:
Joaquin Almeira,
Tomas S. Grigera,
Dante R. Chialvo,
Sergio A. Cannas
Abstract:
While the support for the relevance of critical dynamics to brain function is increasing, there is much less agreement on the exact nature of the advocated critical point. Thus, a considerable number of theoretical efforts are currently concentrated on which mechanisms and what type/s of transition can be exhibited by neuronal networks models. In that direction, the present work describes the effe…
▽ More
While the support for the relevance of critical dynamics to brain function is increasing, there is much less agreement on the exact nature of the advocated critical point. Thus, a considerable number of theoretical efforts are currently concentrated on which mechanisms and what type/s of transition can be exhibited by neuronal networks models. In that direction, the present work describes the effect of incorporating a fraction of inhibitory neurons on the collective dynamics. As we show, this results in the appearence of a tricritical point for highly connected networks and non-zero fraction of inhibitory neurons. We discuss the relation of the present results with relevant experimental evidence.
△ Less
Submitted 5 July, 2022;
originally announced July 2022.
-
Scale free density and correlations fluctuations in the dynamics of large microbial ecosystems
Authors:
Nahuel Zamponi,
Tomas S. Grigera,
Ewa Gudowska-Nowak,
Dante R. Chialvo
Abstract:
Microorganisms self-organize in very large communities exhibiting complex fluctuations. Despite recent advances, still the mechanism by which these systems are able to exhibit large variability at the one hand and dynamical robustness on the other, is not fully explained. With that motivation, here we analyze three aspects of the dynamics of the microbiota and plankton: the density fluctuations, t…
▽ More
Microorganisms self-organize in very large communities exhibiting complex fluctuations. Despite recent advances, still the mechanism by which these systems are able to exhibit large variability at the one hand and dynamical robustness on the other, is not fully explained. With that motivation, here we analyze three aspects of the dynamics of the microbiota and plankton: the density fluctuations, the correlation structure and the avalanching dynamics. In all communities under study we find that the results exhibits scale-free density fluctuations, anomalous variance' scaling, scale-free abundance correlations and stationary scale-free avalanching dynamics. These behaviors, typical in systems exhibiting critical dynamics, suggest criticality as a potential mechanism to explain both the robustness and (paradoxical) high irregularity of processes observed in very large microbial communities.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
Scale-free correlations in the dynamics of a small (N ~ 10000) cortical network
Authors:
Sabrina Camargo,
Daniel A. Martin,
Eyisto J. Aguilar Trejo,
Aylen de Florian,
Maciej A. Nowak,
Sergio A. Cannas,
Tomas S. Grigera,
Dante R. Chialvo
Abstract:
The advent of novel opto-genetics technology allows the recording of brain activity with a resolution never seen before. The characterisation of these very large data sets offers new challenges as well as unique theory-testing opportunities. Here we discuss whether the spatial and temporal correlation of the collective activity of thousands of neurons are tangled as predicted by the theory of crit…
▽ More
The advent of novel opto-genetics technology allows the recording of brain activity with a resolution never seen before. The characterisation of these very large data sets offers new challenges as well as unique theory-testing opportunities. Here we discuss whether the spatial and temporal correlation of the collective activity of thousands of neurons are tangled as predicted by the theory of critical phenomena. The analysis shows that both, the correlation length $ξ$ and the correlation time $τ$ scale as predicted as a function of the system size. With some peculiarities that we discuss, the analysis uncovers new evidence consistent with the view that the large scale brain cortical dynamics corresponds to critical phenomena.
△ Less
Submitted 13 September, 2023; v1 submitted 15 June, 2022;
originally announced June 2022.
-
Finite-size correlation behavior near a critical point: a simple metric for monitoring the state of a neural network
Authors:
Eyisto J. Aguilar Trejo,
Daniel A. Martin,
Dulara De Zoysa,
Zac Bowen,
Tomas S. Grigera,
Sergio A. Cannas,
Wolfgang Losert,
Dante R. Chialvo
Abstract:
In this article, a correlation metric $κ_C$ is proposed for the inference of the dynamical state of neuronal networks. $κ_C$ is computed from the scaling of the correlation length with the size of the observation region, which shows qualitatively different behavior near and away from the critical point of a continuous phase transition. The implementation is first studied on a neuronal network mode…
▽ More
In this article, a correlation metric $κ_C$ is proposed for the inference of the dynamical state of neuronal networks. $κ_C$ is computed from the scaling of the correlation length with the size of the observation region, which shows qualitatively different behavior near and away from the critical point of a continuous phase transition. The implementation is first studied on a neuronal network model, where the results of this new metric coincide with those obtained from neuronal avalanche analysis, thus well characterizing the critical state of the network. The approach is further tested with brain optogenetic recordings in behaving mice from a publicly available database. Potential applications and limitations for its use with currently available optical imaging techniques are discussed.
△ Less
Submitted 23 January, 2023; v1 submitted 23 May, 2022;
originally announced May 2022.
-
Short-time dynamics in active systems: the Vicsek model
Authors:
M. Leticia Rubio Puzzo,
Ernesto S. Loscar,
Andres De Virgiliis,
Tomas S. Grigera
Abstract:
We study the short-time dynamics (STD) of the Vicsek model with vector noise. The study of STD has proved to be very useful in the determination of the critical point, critical exponents, and spinodal points in equilibrium phase transitions. Here we aim to test its applicability in active systems. We find that, despite the essential non-equilibrium characteristics of the VM (absence of detailed ba…
▽ More
We study the short-time dynamics (STD) of the Vicsek model with vector noise. The study of STD has proved to be very useful in the determination of the critical point, critical exponents, and spinodal points in equilibrium phase transitions. Here we aim to test its applicability in active systems. We find that, despite the essential non-equilibrium characteristics of the VM (absence of detailed balance, activity), the STD presents qualitatively the same phenomenology as in equilibrium systems. From the STD one can distinguish whether the transition is continuous or discontinuous (which we have checked also computing the Binder cumulant). When the transition is continuous, one can determine the critical point and the critical exponents.
△ Less
Submitted 9 May, 2022;
originally announced May 2022.
-
On the neural network flow of spin configurations
Authors:
Santiago Acevedo,
Carlos A. Lamas,
Alejo Costa Duran,
Mauricio B. Sturla,
Tomás S. Grigera
Abstract:
We study the so-called neural network flow of spin configurations in the 2-d Ising ferromagnet. This flow is generated by successive reconstructions of spin configurations, obtained by an artificial neural network like a restricted Boltzmann machine or an autoencoder. It was reported recently that this flow may have a fixed point at the critical temperature of the system, and even allow the comput…
▽ More
We study the so-called neural network flow of spin configurations in the 2-d Ising ferromagnet. This flow is generated by successive reconstructions of spin configurations, obtained by an artificial neural network like a restricted Boltzmann machine or an autoencoder. It was reported recently that this flow may have a fixed point at the critical temperature of the system, and even allow the computation of critical exponents. Here we focus on the flow produced by a fully-connected autoencoder, and we refute the claim that this flow converges to the critical point of the system by directly measuring physical observables, and showing that the flow strongly depends on the network hyperparameters. We explore the network metric, the reconstruction error, and we relate it to the so called intrinsic dimension of data, to shed light on the origin and properties of the flow.
△ Less
Submitted 29 June, 2022; v1 submitted 13 March, 2022;
originally announced March 2022.
-
Renormalization group study of marginal ferromagnetism
Authors:
Andrea Cavagna,
Antonio Culla,
Tomás S. Grigera
Abstract:
When studying the collective motion of biological groups a useful theoretical framework is that of ferromagnetic systems, in which the alignment interactions are a surrogate of the effective imitation among the individuals. In this context, the experimental discovery of scale-free correlations of speed fluctuations in starling flocks poses a challenge to the common statistical physics wisdom, as i…
▽ More
When studying the collective motion of biological groups a useful theoretical framework is that of ferromagnetic systems, in which the alignment interactions are a surrogate of the effective imitation among the individuals. In this context, the experimental discovery of scale-free correlations of speed fluctuations in starling flocks poses a challenge to the common statistical physics wisdom, as in the ordered phase of standard ferromagnetic models with $\mathrm{O}(n)$ symmetry, the modulus of the order parameter has finite correlation length. To make sense of this anomaly a novel ferromagnetic theory has been proposed, where the bare confining potential has zero second derivative (i.e.\ it is marginal) along the modulus of the order parameter. The marginal model exhibits a zero-temperature critical point, where the modulus correlation length diverges, hence allowing to boost both correlation and collective order by simply reducing the temperature. Here, we derive an effective field theory describing the marginal model close to the $T=0$ critical point and calculate the renormalization group equations at one loop within a momentum shell approach. We discover a non-trivial scenario, as the cubic and quartic vertices do not vanish in the infrared limit, while the coupling constants effectively regulating the exponents $ν$ and $η$ have upper critical dimension $d_c=2$, so that in three dimensions the critical exponents acquire their free values, $ν=1/2$ and $η=0$. This theoretical scenario is verified by a Monte Carlo study of the modulus susceptibility in three dimensions, where the standard finite-size scaling relations have to be adapted to the case of $d>d_c$. The numerical data fully confirm our theoretical results.
△ Less
Submitted 12 January, 2023; v1 submitted 9 February, 2022;
originally announced February 2022.
-
Apparently similar neuronal dynamics may lead to different collective repertoire
Authors:
Margarita M. Sánchez Díaz,
Eyisto J. Aguilar Trejo,
Daniel A. Martin,
Sergio A. Cannas,
Tomás S. Grigera,
Dante R. Chialvo
Abstract:
This report is concerned with the relevance of the microscopic rules, that implement individual neuronal activation, in determining the collective dynamics, under variations of the network topology. To fix ideas we study the dynamics of two cellular automaton models, commonly used, rather in-distinctively, as the building blocks of large scale neuronal networks. One model, due to Greenberg \& Hast…
▽ More
This report is concerned with the relevance of the microscopic rules, that implement individual neuronal activation, in determining the collective dynamics, under variations of the network topology. To fix ideas we study the dynamics of two cellular automaton models, commonly used, rather in-distinctively, as the building blocks of large scale neuronal networks. One model, due to Greenberg \& Hastings, (GH) can be described by evolution equations mimicking an integrate-and-fire process, while the other model, due to Kinouchi \& Copelli, (KC) represents an abstract branching process, where a single active neuron activates a given number of postsynaptic neurons according to a prescribed "activity" branching ratio. Despite the apparent similarity between the local neuronal dynamics of the two models, it is shown that they exhibit very different collective dynamics as a function of the network topology. The GH model shows qualitatively different dynamical regimes as the network topology is varied, including transients to a ground (inactive) state, continuous and discontinuous dynamical phase transitions. In contrast, the KC model only exhibits a continuous phase transition, independently of the network topology. These results highlight the importance of paying attention to the microscopic rules chosen to model the inter-neuronal interactions in large scale numerical simulations, in particular when the network topology is far from a mean field description. One such case is the extensive work being done in the context of the Human Connectome, where a wide variety of types of models are being used to understand the brain collective dynamics.
△ Less
Submitted 3 September, 2021;
originally announced September 2021.
-
Natural Swarms in $\bf 3.99$ Dimensions
Authors:
Andrea Cavagna,
Luca Di Carlo,
Irene Giardina,
Tomas S. Grigera,
Stefania Melillo,
Leonardo Parisi,
Giulia Pisegna,
Mattia Scandolo
Abstract:
The dynamical critical exponent $z$ of natural swarms of insects is calculated using the renormalization group to order $ε= 4-d$. A novel fixed point emerges, where both activity and inertia are relevant. In three dimensions the critical exponent at the new fixed point is $z = 1.35$, in agreement with both experiments ($1.37 \pm 0.11$) and numerical simulations ($1.35 \pm 0.04$).
The dynamical critical exponent $z$ of natural swarms of insects is calculated using the renormalization group to order $ε= 4-d$. A novel fixed point emerges, where both activity and inertia are relevant. In three dimensions the critical exponent at the new fixed point is $z = 1.35$, in agreement with both experiments ($1.37 \pm 0.11$) and numerical simulations ($1.35 \pm 0.04$).
△ Less
Submitted 19 December, 2022; v1 submitted 9 July, 2021;
originally announced July 2021.
-
Correlation functions as a tool to study collective behaviour phenomena in biological systems
Authors:
Tomás S. Grigera
Abstract:
Much of interesting complex biological behaviour arises from collective properties. Important information about collective behaviour lies in the time and space structure of fluctuations around average properties, and two-point correlation functions are a fundamental tool to study these fluctuations. We give a self-contained presentation of definitions and techniques for computation of correlation…
▽ More
Much of interesting complex biological behaviour arises from collective properties. Important information about collective behaviour lies in the time and space structure of fluctuations around average properties, and two-point correlation functions are a fundamental tool to study these fluctuations. We give a self-contained presentation of definitions and techniques for computation of correlation functions aimed at providing students and researchers outside the field of statistical physics a practical guide to calculating correlation functions from experimental and simulation data. We discuss some properties of correlations in critical systems, and the effect of finite system size, which is particularly relevant for most biological experimental systems. Finally we apply these to the case of the dynamical transition in a simple neuronal model,
△ Less
Submitted 9 July, 2021;
originally announced July 2021.
-
Dynamical renormalization group for mode-coupling field theories with solenoidal constraint
Authors:
Andrea Cavagna,
Luca di Carlo,
Irene Giardina,
Tomas S. Grigera,
Giulia Pisegna,
Mattia Scandolo
Abstract:
The recent inflow of empirical data about the collective behaviour of strongly correlated biological systems has brought field theory and the renormalization group into the biophysical arena. Experiments on bird flocks and insect swarms show that social forces act on the particles' velocity through the generator of its rotations, namely the spin, indicating that mode-coupling field theories are ne…
▽ More
The recent inflow of empirical data about the collective behaviour of strongly correlated biological systems has brought field theory and the renormalization group into the biophysical arena. Experiments on bird flocks and insect swarms show that social forces act on the particles' velocity through the generator of its rotations, namely the spin, indicating that mode-coupling field theories are necessary to reproduce the correct dynamical behaviour. Unfortunately, a theory for three coupled fields - density, velocity and spin - has a prohibitive degree of intricacy. A simplifying path consists in getting rid of density fluctuations by studying incompressible systems. This requires imposing a solenoidal constraint on the primary field, an unsolved problem even for equilibrium mode-coupling theories. Here, we perform an equilibrium dynamic renormalization group analysis of a mode-coupling field theory subject to a solenoidal constraint; using the classification of Halperin and Hohenberg, we can dub this case as a solenoidal Model G. We demonstrate that the constraint produces a new vertex that mixes static and dynamical coupling constants, and that this vertex is essential to grant the closure of the renormalization group structure and the consistency of dynamics with statics. Interestingly, although the solenoidal constraint leads to a modification of the static universality class, we find that it does not change the dynamical universality class, a result that seems to represent an exception to the general rule that dynamical universality classes are narrower than static ones. Our results constitute a solid stepping stone in the admittedly large chasm towards developing an off-equilibrium mode-coupling theory of biological groups.
△ Less
Submitted 19 March, 2021;
originally announced March 2021.
-
Marginal speed confinement resolves the conflict between correlation and control in natural flocks of birds
Authors:
Andrea Cavagna,
Antonio Culla,
Xiao Feng,
Irene Giardina,
Tomás S. Grigera,
Willow Kion-Crosby,
Stefania Melillo,
Giulia Pisegna,
Lorena Postiglione,
Pablo Villegas
Abstract:
Speed fluctuations of individual birds in natural flocks are moderate, due to the aerodynamic and biomechanical constraints of flight. Yet the spatial correlations of such fluctuations are scale-free, namely they have a range as wide as the entire group, a property linked to the capacity of the system to collectively respond to external perturbations. Scale-free correlations and moderate fluctuati…
▽ More
Speed fluctuations of individual birds in natural flocks are moderate, due to the aerodynamic and biomechanical constraints of flight. Yet the spatial correlations of such fluctuations are scale-free, namely they have a range as wide as the entire group, a property linked to the capacity of the system to collectively respond to external perturbations. Scale-free correlations and moderate fluctuations set conflicting constraints on the mechanism controlling the speed of each agent, as the factors boosting correlation amplify fluctuations, and vice versa. Here, using a statistical field theory approach, we suggest that a marginal speed confinement that ignores small deviations from the natural reference value while ferociously suppressing larger speed fluctuations, is able to reconcile scale-free correlations with biologically acceptable group's speed. We validate our theoretical predictions by comparing them with field experimental data on starling flocks with group sizes spanning an unprecedented interval of over two orders of magnitude.
△ Less
Submitted 11 May, 2022; v1 submitted 24 January, 2021;
originally announced January 2021.
-
Equilibrium to off-equilibrium crossover in homogeneous active matter
Authors:
Andrea Cavagna,
Luca Di Carlo,
Irene Giardina,
Tomás S. Grigera,
Giulia Pisegna
Abstract:
We study the crossover between equilibrium and off-equilibrium dynamical universality classes in the Vicsek model near its ordering transition. Starting from the incompressible hydrodynamic theory of Chen et al \cite{chen2015critical}, we show that increasing the activity leads to a renormalization group (RG) crossover between the equilibrium ferromagnetic fixed point, with dynamical critical expo…
▽ More
We study the crossover between equilibrium and off-equilibrium dynamical universality classes in the Vicsek model near its ordering transition. Starting from the incompressible hydrodynamic theory of Chen et al \cite{chen2015critical}, we show that increasing the activity leads to a renormalization group (RG) crossover between the equilibrium ferromagnetic fixed point, with dynamical critical exponent $z = 2$, and the off-equilibrium active fixed point, with $z = 1.7$ (in $d=3$). We run simulations of the classic Vicsek model in the near-ordering regime and find that critical slowing down indeed changes with activity, displaying two exponents that are in remarkable agreement with the RG prediction. The equilibrium-to-off-equilibrium crossover is ruled by a characteristic length scale beyond which active dynamics takes over. Such length scale is smaller the larger the activity, suggesting the existence of a general trade-off between activity and system's size in determining the dynamical universality class of active matter.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Joint assessment of density correlations and fluctuations for analysing spatial tree patterns
Authors:
P. Villegas,
A. Cavagna,
M. Cencini,
H. Fort,
T. S. Grigera
Abstract:
Inferring the processes underlying the emergence of observed patterns is a key challenge in theoretical ecology. Much effort has been made in the past decades to collect extensive and detailed information about the spatial distribution of tropical rainforests, as demonstrated, e.g., in the 50 ha tropical forest plot on Barro Colorado Island, Panama. These kind of plots have been crucial to shed li…
▽ More
Inferring the processes underlying the emergence of observed patterns is a key challenge in theoretical ecology. Much effort has been made in the past decades to collect extensive and detailed information about the spatial distribution of tropical rainforests, as demonstrated, e.g., in the 50 ha tropical forest plot on Barro Colorado Island, Panama. These kind of plots have been crucial to shed light on diverse qualitative features, emerging both at the single-species or the community level, like the spatial aggregation or clustering at short scales. Here, we build on the progress made in the study of the density correlation functions applied to biological systems, focusing on the importance of accurately defining the borders of the set of trees, and removing the induced biases. We also pinpoint the importance of combining the study of correlations with the scale dependence of fluctuations in density, which are linked to the well known empirical Taylor's power law. Density correlations and fluctuations, in conjunction, provide an unique opportunity to interpret the behaviors and possibly to allow comparisons between data and models. We also study such quantities in models of spatial patterns and, in particular, we find that a spatially explicit neutral model generates patterns with many qualitative features in common with the empirical ones.
△ Less
Submitted 2 December, 2020; v1 submitted 31 August, 2020;
originally announced August 2020.
-
Box-scaling as a proxy of finite-size correlations
Authors:
Daniel A. Martin,
Tiago L. Ribeiro,
Sergio A. Cannas,
Tomas S. Grigera,
Dietmar Plenz,
Dante R. Chialvo
Abstract:
The scaling of correlations as a function of system size provides important hints to understand critical phenomena on a variety of systems. Its study in biological systems offers two challenges: usually they are not of infinite size, and in the majority of cases sizes can not be varied at will. Here we discuss how finite-size scaling can be approximated in an experimental system of fixed and relat…
▽ More
The scaling of correlations as a function of system size provides important hints to understand critical phenomena on a variety of systems. Its study in biological systems offers two challenges: usually they are not of infinite size, and in the majority of cases sizes can not be varied at will. Here we discuss how finite-size scaling can be approximated in an experimental system of fixed and relatively small size by computing correlations inside of a reduced field of view of various sizes (i.e., "box-scaling"). Numerical simulations of a neuronal network are used to verify such approximation, as well as the ferromagnetic 2D Ising model. The numerical results support the validity of the heuristic approach, which should be useful to characterize relevant aspects of critical phenomena in biological systems.
△ Less
Submitted 16 July, 2020;
originally announced July 2020.
-
Everything you wish to know about correlations but are afraid to ask
Authors:
Tomas S. Grigera
Abstract:
We discuss the various definitions of time correlation functions and how to estimate them from experimental or simulation data. We start with the various definitions, both in real and in Fourier space, and explain how to extract from them a characteristic time scale. We then study how to estimate the correlation functions, i.e.\ how to obtain a good approximation to them from a sample of data obta…
▽ More
We discuss the various definitions of time correlation functions and how to estimate them from experimental or simulation data. We start with the various definitions, both in real and in Fourier space, and explain how to extract from them a characteristic time scale. We then study how to estimate the correlation functions, i.e.\ how to obtain a good approximation to them from a sample of data obtained experimentally. Finally we discuss some practical details that arise in the actual computation of these estimates, and we describe some relevant algorithms.
△ Less
Submitted 5 February, 2020;
originally announced February 2020.
-
Controlling a complex system near its critical point via temporal correlations
Authors:
Dante R. Chialvo,
Sergio A. Cannas,
Dietmar Plenz,
Tomas S. Grigera
Abstract:
A wide variety of complex systems exhibit large fluctuations both in space and time that often can be attributed to the presence of some kind of critical phenomena. Under such critical scenario it is well known that the properties of the correlation functions in space and time are two sides of the same coin. Here we test wether systems exhibiting a phase transition could self-tune to its critical…
▽ More
A wide variety of complex systems exhibit large fluctuations both in space and time that often can be attributed to the presence of some kind of critical phenomena. Under such critical scenario it is well known that the properties of the correlation functions in space and time are two sides of the same coin. Here we test wether systems exhibiting a phase transition could self-tune to its critical point taking advantage of such correlation properties. We describe results in three models: the 2D Ising ferromagnetic model, the 3D Vicsek flocking model and a small-world neuronal network model. We illustrate how the feedback of the autocorrelation function of the order parameter fluctuations is able to shift the system towards its critical point. Since the results rely on universal properties they are expected to be relevant to a variety of other settings.
△ Less
Submitted 28 May, 2019;
originally announced May 2019.
-
Renormalization group crossover in the critical dynamics of field theories with mode coupling terms
Authors:
Andrea Cavagna,
Luca Di Carlo,
Irene Giardina,
Luca Grandinetti,
Tomas S. Grigera,
Giulia Pisegna
Abstract:
Motivated by the collective behaviour of biological swarms, we study the critical dynamics of field theories with coupling between order parameter and conjugate momentum in the presence of dissipation. By performing a dynamical renormalization group calculation at one loop, we show that the violation of momentum conservation generates a crossover between a conservative yet IR-unstable fixed point,…
▽ More
Motivated by the collective behaviour of biological swarms, we study the critical dynamics of field theories with coupling between order parameter and conjugate momentum in the presence of dissipation. By performing a dynamical renormalization group calculation at one loop, we show that the violation of momentum conservation generates a crossover between a conservative yet IR-unstable fixed point, characterized by a dynamic critical exponent $z=d/2$, and a dissipative IR-stable fixed point with $z=2$. Interestingly, the two fixed points have different upper critical dimensions. The interplay between these two fixed points gives rise to a crossover in the critical dynamics of the system, characterized by a crossover exponent $κ=4/d$. Such crossover is regulated by a conservation length scale, $\mathcal R_0$, which is larger the smaller the dissipation: beyond $\mathcal R_0$ the dissipative fixed point dominates, while at shorter distances dynamics is ruled by the conservative fixed point and critical exponent, a behaviour which is all the more relevant in finite-size systems with weak dissipation. We run numerical simulations in three dimensions and find a crossover between the exponents $z=3/2$ and $z=2$ in the critical slowing down of the system, confirming the renormalization group results. From the biophysical point of view, our calculation indicates that in finite-size biological groups mode-coupling terms in the equation of motion can significantly change the dynamical critical exponents even in the presence of dissipation, a step towards reconciling theory with experiments in natural swarms. Moreover, our result provides the scale within which fully conservative Bose-Einstein condensation is a good approximation in systems with weak symmetry-breaking terms violating number conservation, as quantum magnets or photon gases.
△ Less
Submitted 3 May, 2019;
originally announced May 2019.
-
Dynamical renormalization group approach to the collective behaviour of swarms
Authors:
Andrea Cavagna,
Luca Di Carlo,
Irene Giardina,
Luca Grandinetti,
Tomas S. Grigera,
Giulia Pisegna
Abstract:
We study the critical behaviour of a model with non-dissipative couplings aimed at describing the collective behaviour of natural swarms, using the dynamical renormalization group. At one loop, we find a crossover between a conservative yet unstable fixed point, characterized by a dynamical critical exponent $z=d/2$, and a dissipative stable fixed point with $z=2$, a result we confirm through nume…
▽ More
We study the critical behaviour of a model with non-dissipative couplings aimed at describing the collective behaviour of natural swarms, using the dynamical renormalization group. At one loop, we find a crossover between a conservative yet unstable fixed point, characterized by a dynamical critical exponent $z=d/2$, and a dissipative stable fixed point with $z=2$, a result we confirm through numerical simulations. The crossover is regulated by a conservation length scale that is larger the smaller the effective friction, so that in finite-size biological systems with low dissipation, dynamics is ruled by the conservative fixed point. In three dimensions this mechanism gives $z=3/2$, a value significantly closer to the experimental result $z\approx 1$ than the value $z\approx 2$ found in fully dissipative models, either at or off equilibrium. This result indicates that non-dissipative dynamical couplings are necessary to develop a theory of natural swarms fully consistent with experiments
△ Less
Submitted 3 May, 2019;
originally announced May 2019.
-
Low-temperature marginal ferromagnetism explains anomalous scale-free correlations in natural flocks
Authors:
Andrea Cavagna,
Antonio Culla,
Luca Di Carlo,
Irene Giardina,
Tomas S. Grigera
Abstract:
We introduce a new ferromagnetic model capable of reproducing one of the most intriguing properties of collective behaviour in starling flocks, namely the fact that strong collective order of the system coexists with scale-free correlations of the modulus of the microscopic degrees of freedom, that is the birds' speeds. The key idea of the new theory is that the single-particle potential needed to…
▽ More
We introduce a new ferromagnetic model capable of reproducing one of the most intriguing properties of collective behaviour in starling flocks, namely the fact that strong collective order of the system coexists with scale-free correlations of the modulus of the microscopic degrees of freedom, that is the birds' speeds. The key idea of the new theory is that the single-particle potential needed to bound the modulus of the microscopic degrees of freedom around a finite value, is marginal, that is has zero curvature. We study the model by using mean-field approximation and Monte Carlo simulations in three dimensions, complemented by finite-size scaling analysis. While at the standard critical temperature, $T_c$, the properties of the marginal model are exactly the same as a normal ferromagnet with continuous symmetry-breaking, our results show that a novel zero-temperature critical point emerges, so that in its deeply ordered phase the marginal model develops divergent susceptibility and correlation length of the modulus of the microscopic degrees of freedom, in complete analogy with experimental data on natural flocks of starlings.
△ Less
Submitted 17 September, 2019; v1 submitted 18 December, 2018;
originally announced December 2018.
-
Self-propelled Vicsek particles at low speed and low density
Authors:
M. Leticia Rubio Puzzo,
Andres De Virgiliis,
Tomas S. Grigera
Abstract:
We study through numerical simulation the Vicsek model for very low speeds and densities. We consider scalar noise in 2-d and 3-d, and vector noise in 3-d. We focus on the behavior of the critical noise with density and speed, trying to clarify seemingly contradictory earlier results. We find that, for scalar noise, the critical noise is a power law both in density and speed, but although we confi…
▽ More
We study through numerical simulation the Vicsek model for very low speeds and densities. We consider scalar noise in 2-d and 3-d, and vector noise in 3-d. We focus on the behavior of the critical noise with density and speed, trying to clarify seemingly contradictory earlier results. We find that, for scalar noise, the critical noise is a power law both in density and speed, but although we confirm the density exponent in 2-d, we find a speed exponent different from earlier reports (we consider lower speeds than previous studies). On the other hand, for the vector noise case we find that the dependence of the critical noise cannot be separated as a product of power laws in speed and density. Finally, we study the dependence of the relaxation time with speed and find the same power law in 2-d and 3-d, with and exponent that depends on whether the noise is above or below the critical value.
△ Less
Submitted 4 October, 2018;
originally announced October 2018.
-
Propagating speed waves in flocks: a mathematical model
Authors:
Andrea Cavagna,
Daniele Conti,
Irene Giardina,
Tomas S. Grigera
Abstract:
Efficient collective response to external perturbations is one of the most striking abilities of a biological system. Signal propagation through the group is an important condition for the imple- mentation of such a response. Information transfer has been experimentally observed in the turning mechanism of birds flocks. In this context it is well-known also the existence of density waves: birds un…
▽ More
Efficient collective response to external perturbations is one of the most striking abilities of a biological system. Signal propagation through the group is an important condition for the imple- mentation of such a response. Information transfer has been experimentally observed in the turning mechanism of birds flocks. In this context it is well-known also the existence of density waves: birds under predation, attempting to escape, give rise to self-organized density waves that propagates linearly on the flock. Most aspects of this phenomenon are still not fully captured by theoretical models. In this work we present a new model for the propagation of the speed (the modulus of the velocity) fluctuations inside a flock, which is the simplest way to reproduce the observed density waves. We have studied the full solution of the model in d = 1 and we found that there is a line in the parameter space along which the system relaxes in the fastest way with no oscillation after a signal has passed. This is the critical damping condition. By analyzing the parameters plane we show that critical damping represents an attractor for a steepest descent dynamics of the return time of the system. Finally we propose a method to test the validity of the model through through future experiments.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Dynamic scaling in natural swarms
Authors:
A. Cavagna,
D. Conti,
C. Creato,
L. Del Castello,
I. Giardina,
T. S. Grigera,
S. Melillo,
L. Parisi,
M. Viale
Abstract:
Collective behaviour in biological systems pitches us against theoretical challenges way beyond the borders of ordinary statistical physics. The lack of concepts like scaling and renormalization is particularly grievous, as it forces us to negotiate with scores of details whose relevance is often hard to assess. In an attempt to improve on this situation, we present here experimental evidence of t…
▽ More
Collective behaviour in biological systems pitches us against theoretical challenges way beyond the borders of ordinary statistical physics. The lack of concepts like scaling and renormalization is particularly grievous, as it forces us to negotiate with scores of details whose relevance is often hard to assess. In an attempt to improve on this situation, we present here experimental evidence of the emergence of dynamic scaling laws in natural swarms. We find that spatio-temporal correlation functions in different swarms can be rescaled by using a single characteristic time, which grows with the correlation length with a dynamical critical exponent z~1. We run simulations of a model of self-propelled particles in its swarming phase and find z~2, suggesting that natural swarms belong to a novel dynamic universality class. This conclusion is strengthened by experimental evidence of non-exponential relaxation and paramagnetic spin-wave remnants, indicating that previously overlooked inertial effects are needed to describe swarm dynamics. The absence of a purely relaxational regime suggests that natural swarms are subject to a near-critical censorship of hydrodynamics.
△ Less
Submitted 24 November, 2016;
originally announced November 2016.
-
Structure and dynamics of the t154 lattice glass
Authors:
Alejandro Seif,
Tomas S. Grigera
Abstract:
We revisit the $t154$ variant of the Biroli-Mezard lattice glass, complementing previous studies by studying statics and dynamics under periodic boundary conditions as well as systems confined in cavities with amorphous boundaries. We compute the point-to-se correlation and relaxation times under the different boundary conditions. Results point to a scenario with dynamics ruled by structural corre…
▽ More
We revisit the $t154$ variant of the Biroli-Mezard lattice glass, complementing previous studies by studying statics and dynamics under periodic boundary conditions as well as systems confined in cavities with amorphous boundaries. We compute the point-to-se correlation and relaxation times under the different boundary conditions. Results point to a scenario with dynamics ruled by structural correlations.
△ Less
Submitted 21 November, 2016;
originally announced November 2016.
-
Response to "Comment on Static correlations functions and domain walls in glass-forming liquids: The case of a sandwich geometry" [J. Chem. Phys. 144, 227101 (2016)]
Authors:
Giacomo Gradenigo,
Roberto Trozzo,
Andrea Cavagna,
Tomas S. Grigera
Abstract:
The point-to-set correlation function has proved to be a very valuable tool to probe structural correlations in disordered systems, but more than that, its detailed behavior has been used to try to draw information on the mechanisms leading to glassy behavior in supercooled liquids. For this reason it is of primary importance to discern which of those details are peculiar to glassy systems, and wh…
▽ More
The point-to-set correlation function has proved to be a very valuable tool to probe structural correlations in disordered systems, but more than that, its detailed behavior has been used to try to draw information on the mechanisms leading to glassy behavior in supercooled liquids. For this reason it is of primary importance to discern which of those details are peculiar to glassy systems, and which are general features of confinement. Within the present response we provide an answer to the concerns raised in [J. Chem. Phys. 144, 227101 (2016)].
△ Less
Submitted 24 June, 2016;
originally announced June 2016.
-
Spatio-temporal correlations in models of collective motion ruled by different dynamical laws
Authors:
Andrea Cavagna,
Daniele Conti,
Irene Giardina,
Tomas S. Grigera,
Stefania Melillo,
Massimiliano Viale
Abstract:
Information transfer is an essential factor in determining the robustness of collective behaviour in biological systems with distributed control. The most direct way to study the information transfer mechanisms is to experimentally detect the propagation across the system of a signal triggered by some perturbation. However, for field experiments this method is inefficient, as the possibilities of…
▽ More
Information transfer is an essential factor in determining the robustness of collective behaviour in biological systems with distributed control. The most direct way to study the information transfer mechanisms is to experimentally detect the propagation across the system of a signal triggered by some perturbation. However, for field experiments this method is inefficient, as the possibilities of the observer to perturb the group are limited and empirical observations must rely on rare natural perturbations. An alternative way is to use spatio-temporal correlations to assess the information transfer mechanism directly from the spontaneous fluctuations of the system, without the need to have an actual propagating signal on record. We test the approach on ground truth data provided by numerical simulations in three dimensions of two models of collective behaviour characterized by very different dynamical equations and information transfer mechanisms: the classic Vicsek model, describing an overdamped noninertial dynamics and the inertial spin model, characterized by an un- derdamped inertial dynamics. By using dynamical finite size scaling, we show that spatio-temporal correlations are able to distinguish unambiguously the diffusive information transfer mechanism of the Vicsek model from the linear mechanism of the inertial spin model.
△ Less
Submitted 12 December, 2016; v1 submitted 31 May, 2016;
originally announced May 2016.
-
An intermediate state between the kagome-ice and the fully polarized state in Dy$_2$Ti$_2$O$_7$
Authors:
S. A. Grigera,
R. Borzi,
D. G. Slobinsky,
A. S. Gibbs,
R. Higashinaka,
Y. Maeno,
T. S. Grigera
Abstract:
Dy$_2$Ti$_2$O$_7$ is at present the cleanest example of a spin-ice material. Previous theoretical and experimental work on the first-order transition between the kagome-ice and the fully polarized state has been taken as a validation for the dipolar spin-ice model. Here we investigate in further depth this phase transition using ac-susceptibility and dc-magnetization, and compare this results with…
▽ More
Dy$_2$Ti$_2$O$_7$ is at present the cleanest example of a spin-ice material. Previous theoretical and experimental work on the first-order transition between the kagome-ice and the fully polarized state has been taken as a validation for the dipolar spin-ice model. Here we investigate in further depth this phase transition using ac-susceptibility and dc-magnetization, and compare this results with Monte-Carlo simulations and previous magnetization and specific heat measurements. We find signatures of an intermediate state between the kagome-ice and full polarization. This signatures are absent in current theoretical models used to describe spin-ice materials.
△ Less
Submitted 20 October, 2015;
originally announced October 2015.
-
Ageing and crystallization in a lattice glass model
Authors:
Alejandro Seif,
Ernesto S. Loscar,
Tomas S. Grigera
Abstract:
We have studied a the 3-$d$ lattice glass of Pica Ciamarra, Tarzia, de Candia and Coniglio [Phys.\ Rev.\ E. 67, 057105 (2013)], which has been shown to reproduce several features of the structural glass phenomenology, such as the cage effect, exponential increase of relaxation times and ageing. We show, using short-time dynamics, that the metastability limit is above the estimated Kauzmann tempera…
▽ More
We have studied a the 3-$d$ lattice glass of Pica Ciamarra, Tarzia, de Candia and Coniglio [Phys.\ Rev.\ E. 67, 057105 (2013)], which has been shown to reproduce several features of the structural glass phenomenology, such as the cage effect, exponential increase of relaxation times and ageing. We show, using short-time dynamics, that the metastability limit is above the estimated Kauzmann temperature. We also find that in the region where the metastable liquid exists the aging exponent is lower that 0.5, indicating that equilibrium is reached relatively quickly. We conclude that the usefulness of this model to study the deeply supercooled regime is rather limited.
△ Less
Submitted 15 December, 2014;
originally announced December 2014.
-
Silent Flocks
Authors:
Andrea Cavagna,
Irene Giardina,
Tomas S. Grigera,
Asja Jelic,
Dov Levine,
Sriram Ramaswamy,
Massimiliano Viale
Abstract:
Experiments find coherent information transfer through biological groups on length and time scales distinctly below those on which asymptotically correct hydrodynamic theories apply. We present here a new continuum theory of collective motion coupling the velocity and density fields of Toner and Tu to the inertial spin field recently introduced to describe information propagation in natural flocks…
▽ More
Experiments find coherent information transfer through biological groups on length and time scales distinctly below those on which asymptotically correct hydrodynamic theories apply. We present here a new continuum theory of collective motion coupling the velocity and density fields of Toner and Tu to the inertial spin field recently introduced to describe information propagation in natural flocks of birds. The long-wavelength limit of the new equations reproduces Toner-Tu theory, while at shorter wavelengths (or, equivalently, smaller damping), spin fluctuations dominate over density fluctuations and second sound propagation of the kind observed in real flocks emerges. We study the dispersion relation of the new theory and find that when the speed of second sound is large, a gap sharply separates first from second sound modes. This gap implies the existence of `silent' flocks, namely medium-sized systems across which neither first nor second sound can propagate.
△ Less
Submitted 10 October, 2014;
originally announced October 2014.
-
Specific heat anomaly in a supercooled liquid with amorphous boundary conditions
Authors:
Daniel A. Martin,
Andrea Cavagna,
Tomas S. Grigera
Abstract:
We study the specific heat of a model supercooled liquid confined in a spherical cavity with amorphous boundary conditions. We find the equilibrium specific heat has a cavity-size-dependent peak as a function of temperature. The cavity allows us to perform a finite-size scaling (FSS) analysis, which indicates that the peak persists at a finite temperature in the thermodynamic limit. We attempt to…
▽ More
We study the specific heat of a model supercooled liquid confined in a spherical cavity with amorphous boundary conditions. We find the equilibrium specific heat has a cavity-size-dependent peak as a function of temperature. The cavity allows us to perform a finite-size scaling (FSS) analysis, which indicates that the peak persists at a finite temperature in the thermodynamic limit. We attempt to collapse the data onto a FSS curve according to different theoretical scenarios, obtaining reasonable results in two cases: a "not-so-simple" liquid with nonstandard values of the exponents α and ν, and random first-order theory, with two different length scales.
△ Less
Submitted 16 August, 2015; v1 submitted 23 July, 2014;
originally announced July 2014.
-
Superfluid transport of information in turning flocks of starlings
Authors:
Alessandro Attanasi,
Andrea Cavagna,
Lorenzo Del Castello,
Irene Giardina,
Tomas S. Grigera,
Asja Jelić,
Stefania Melillo,
Leonardo Parisi,
Oliver Pohl,
Edward Shen,
Massimiliano Viale
Abstract:
Collective decision-making in biological systems requires all individuals in the group to go through a behavioural change of state. During this transition, the efficiency of information transport is a key factor to prevent cohesion loss and preserve robustness. The precise mechanism by which natural groups achieve such efficiency, though, is currently not fully understood. Here, we present an expe…
▽ More
Collective decision-making in biological systems requires all individuals in the group to go through a behavioural change of state. During this transition, the efficiency of information transport is a key factor to prevent cohesion loss and preserve robustness. The precise mechanism by which natural groups achieve such efficiency, though, is currently not fully understood. Here, we present an experimental study of starling flocks performing collective turns in the field. We find that the information to change direction propagates across the flock linearly in time with negligible attenuation, hence keeping group decoherence to a minimum. This result contrasts with current theories of collective motion, which predict a slower and dissipative transport of directional information. We propose a novel theory whose cornerstone is the existence of a conserved spin current generated by the gauge symmetry of the system. The theory turns out to be mathematically identical to that of superfluid transport in liquid helium and it explains the dissipationless propagating mode observed in turning flocks. Superfluidity also provides a quantitative expression for the speed of propagation of the information, according to which transport must be swifter the stronger the group's orientational order. This prediction is verified by the data. We argue that the link between strong order and efficient decision-making required by superfluidity may be the adaptive drive for the high degree of behavioural polarization observed in many living groups. The mathematical equivalence between superfluid liquids and turning flocks is a compelling demonstration of the far-reaching consequences of symmetry and conservation laws across different natural systems.
△ Less
Submitted 28 March, 2013;
originally announced March 2013.
-
Static correlations functions and domain walls in glass-forming liquids: the case of a sandwich geometry
Authors:
G. Gradenigo,
R. Trozzo,
A. Cavagna,
T. S. Grigera,
P. Verrocchio
Abstract:
The problem of measuring nontrivial static correlations in deeply supercooled liquids made recently some progress thanks to the introduction of amorphous boundary conditions, in which a set of free particles is subject to the effect of a different set of particles frozen into their (low temperature) equilibrium positions. In this way, one can study the crossover from nonergodic to ergodic phase, a…
▽ More
The problem of measuring nontrivial static correlations in deeply supercooled liquids made recently some progress thanks to the introduction of amorphous boundary conditions, in which a set of free particles is subject to the effect of a different set of particles frozen into their (low temperature) equilibrium positions. In this way, one can study the crossover from nonergodic to ergodic phase, as the size of the free region grows and the effect of the confinement fades. Such crossover defines the so-called point-to-set correlation length, which has been measured in a spherical geometry, or cavity. Here, we make further progress in the study ofcorrelations under amorphous boundary conditions by analyzing the equilibrium properties of a glass-forming liquid, confined in a planar ("sandwich") geometry. The mobile particles are subject to amorphous boundary conditions with the particles in the surrounding walls frozen into their low temperature equilibrium configurations. Compared to the cavity, the sandwich geometry has three main advantages: i) the width of the sandwich is decoupled from its longitudinal size, making the thermodynamic limit possible; ii) for very large width, the behaviour off a single wall can be studied; iii) we can use "anti-parallel" boundary conditions to force a domain wall and measure its excess energy. Our results confirm that amorphous boundary conditions are indeed a very useful new tool inthe study of static properties of glass-forming liquids, but also raise some warning about the fact that not all correlation functions that can be calculated in this framework give the same qualitative results.
△ Less
Submitted 26 September, 2012;
originally announced September 2012.
-
glsim: A general library for numerical simulation
Authors:
Tomas S. Grigera
Abstract:
We describe glsim, a C++ library designed to provide routines to perform basic housekeeping tasks common to a very wide range of simulation programs, such as reading simulation parameters or reading and writing self-describing binary files with simulation data. The design also provides a framework to add features to the library while preserving its structure and interfaces.
We describe glsim, a C++ library designed to provide routines to perform basic housekeeping tasks common to a very wide range of simulation programs, such as reading simulation parameters or reading and writing self-describing binary files with simulation data. The design also provides a framework to add features to the library while preserving its structure and interfaces.
△ Less
Submitted 11 May, 2011;
originally announced May 2011.
-
On the high-density expansion for Euclidean Random Matrices
Authors:
T. S. Grigera,
V. Martin-Mayor,
G. Parisi,
P. Urbani,
P. Verrocchio
Abstract:
Diagrammatic techniques to compute perturbatively the spectral properties of Euclidean Random Matrices in the high-density regime are introduced and discussed in detail. Such techniques are developed in two alternative and very different formulations of the mathematical problem and are shown to give identical results up to second order in the perturbative expansion. One method, based on writing th…
▽ More
Diagrammatic techniques to compute perturbatively the spectral properties of Euclidean Random Matrices in the high-density regime are introduced and discussed in detail. Such techniques are developed in two alternative and very different formulations of the mathematical problem and are shown to give identical results up to second order in the perturbative expansion. One method, based on writing the so-called resolvent function as a Taylor series, allows to group the diagrams in a small number of topological classes, providing a simple way to determine the infrared (small momenta) behavior of the theory up to third order, which is of interest for the comparison with experiments. The other method, which reformulates the problem as a field theory, can instead be used to study the infrared behaviour at any perturbative order.
△ Less
Submitted 11 November, 2010;
originally announced November 2010.
-
Numerical simulations of liquids with amorphous boundary conditions
Authors:
A. Cavanga,
T. S. Grigera,
P. Verrocchio
Abstract:
It has recently become clear that simulations under amorphpous boundary conditions (ABCs) can provide valuable information on the dynamics and thermodynamics of disordered systems with no obvious ordered parameter. In particular, they allow to detect a correlation length that is not measurable with standard correlation functions. Here we explain what exactly is meant by ABCs, discuss their relatio…
▽ More
It has recently become clear that simulations under amorphpous boundary conditions (ABCs) can provide valuable information on the dynamics and thermodynamics of disordered systems with no obvious ordered parameter. In particular, they allow to detect a correlation length that is not measurable with standard correlation functions. Here we explain what exactly is meant by ABCs, discuss their relation with point-to-set correlations and briefly describe some recent results obtained with this technique.
△ Less
Submitted 1 September, 2010;
originally announced September 2010.
-
Dynamic relaxation of a liquid cavity under amorphous boundary conditions
Authors:
Andrea Cavagna,
Tomas S. Grigera,
Paolo Verrocchio
Abstract:
The growth of cooperatively rearranging regions was invoked long ago by Adam and Gibbs to explain the slowing down of glass-forming liquids. The lack of knowledge about the nature of the growing order, though, complicates the definition of an appropriate correlation function. One option is the point-to-set correlation function, which measures the spatial span of the influence of amorphous boundary…
▽ More
The growth of cooperatively rearranging regions was invoked long ago by Adam and Gibbs to explain the slowing down of glass-forming liquids. The lack of knowledge about the nature of the growing order, though, complicates the definition of an appropriate correlation function. One option is the point-to-set correlation function, which measures the spatial span of the influence of amorphous boundary conditions on a confined system. By using a swap Monte Carlo algorithm we measure the equilibration time of a liquid droplet bounded by amorphous boundary conditions in a model glass-former at low temperature, and we show that the cavity relaxation time increases with the size of the droplet, saturating to the bulk value when the droplet outgrows the point-to-set correlation length. This fact supports the idea that the point-to-set correlation length is the natural size of the cooperatively rearranging regions. On the other hand, the cavity relaxation time computed by a standard, nonswap dynamics, has the opposite behavior, showing a very steep increase when the cavity size is decreased. We try to reconcile this difference by discussing the possible hybridization between MCT and activated processes, and by introducing a new kind of amorphous boundary conditions, inspired by the concept of frozen external state as an alternative to the commonly used frozen external configuration.
△ Less
Submitted 28 July, 2011; v1 submitted 18 June, 2010;
originally announced June 2010.
-
A phase-separation perspective on dynamic heterogeneities in glass-forming liquids
Authors:
C. Cammarota,
A. Cavagna,
I. Giardina,
G. Gradenigo,
T. S. Grigera,
G. Parisi,
P. Verrocchio
Abstract:
We study dynamic heterogeneities in a model glass-former whose overlap with a reference configuration is constrained to a fixed value. The system phase-separates into regions of small and large overlap, so that dynamical correlations remain strong even for asymptotic times. We calculate an appropriate thermodynamic potential and find evidence of a Maxwell's construction consistent with a spinoda…
▽ More
We study dynamic heterogeneities in a model glass-former whose overlap with a reference configuration is constrained to a fixed value. The system phase-separates into regions of small and large overlap, so that dynamical correlations remain strong even for asymptotic times. We calculate an appropriate thermodynamic potential and find evidence of a Maxwell's construction consistent with a spinodal decomposition of two phases. Our results suggest that dynamic heterogeneities are the expression of an ephemeral phase-separating regime ruled by a finite surface tension.
△ Less
Submitted 14 January, 2010;
originally announced January 2010.
-
Numerical determination of the exponents controlling the relationship between time, length and temperature in glass-forming liquids
Authors:
Chiara Cammarota,
Andrea Cavagna,
Giacomo Gradenigo,
Tomas S. Grigera,
Paolo Verrocchio
Abstract:
There is a certain consensus that the very fast growth of the relaxation time $τ$ occurring in glass-forming liquids on lowering the temperature must be due to the thermally activated rearrangement of correlated regions of growing size. Even though measuring the size of these regions has defied scientists for a while, there is indeed recent evidence of a growing correlation length $ξ$ in glass-f…
▽ More
There is a certain consensus that the very fast growth of the relaxation time $τ$ occurring in glass-forming liquids on lowering the temperature must be due to the thermally activated rearrangement of correlated regions of growing size. Even though measuring the size of these regions has defied scientists for a while, there is indeed recent evidence of a growing correlation length $ξ$ in glass-formers. If we use Arrhenius law and make the mild assumption that the free-energy barrier to rearrangement scales as some power $ψ$ of the size of the correlated regions, we obtain a relationship between time and length, $T\logτ\sim ξ^ψ$. According to both the Adam-Gibbs and the Random First Order theory the correlation length grows as $ξ\sim (T-T_k)^{-1/(d-θ)}$, even though the two theories disagree on the value of $θ$. Therefore, the super-Arrhenius growth of the relaxation time with the temperature is regulated by the two exponents $ψ$ and $θ$ through the relationship $T\logτ\sim (T-T_k)^{-ψ/(d-θ)}$. Despite a few theoretical speculations, up to now there has been no experimental determination of these two exponents. Here we measure them numerically in a model glass-former, finding $ψ=1$ and $θ=2$. Surprisingly, even though the values we found disagree with most previous theoretical suggestions, they give back the well-known VFT law for the relaxation time, $T\logτ\sim (T-T_k)^{-1}$.
△ Less
Submitted 19 November, 2009; v1 submitted 21 June, 2009;
originally announced June 2009.
-
Surface tension fluctuations and a new spinodal point in glass-forming liquids
Authors:
C. Cammarota,
A. Cavagna,
G. Gradenigo,
T. S. Grigera,
P. Verrocchio
Abstract:
The dramatic slowdown of glass-forming liquids has been variously linked to increasing dynamic and static correlation lengths. Yet, empirical evidence is insufficient to decide among competing theories. The random first order theory (RFOT) links the dynamic slowdown to the growth of amorphous static order, whose range depends on a balance between configurational entropy and surface tension. This…
▽ More
The dramatic slowdown of glass-forming liquids has been variously linked to increasing dynamic and static correlation lengths. Yet, empirical evidence is insufficient to decide among competing theories. The random first order theory (RFOT) links the dynamic slowdown to the growth of amorphous static order, whose range depends on a balance between configurational entropy and surface tension. This last quantity is expected to vanish when the temperature surpasses a spinodal point beyond which there are no metastable states. Here we measure for the first time the surface tension in a model glass-former, and find that it vanishes at the energy separating minima from saddles, demonstrating the existence of a spinodal point for amorphous metastable order. Moreover, the fluctuations of surface tension become smaller for lower temperatures, in quantitative agreement with recent theoretical speculation that spatial correlations in glassy systems relax nonexponentially because of the narrowing of the surface tension distribution.
△ Less
Submitted 9 April, 2009;
originally announced April 2009.
-
Non-equilibrium Characterization of Spinodal Points using Short Time Dynamics
Authors:
Ernesto S. Loscar,
Ezequiel E. Ferrero,
Tomas S. Grigera,
Sergio A. Cannas
Abstract:
Though intuitively appealing, the concept of spinodal is rigourously defined only in systems with infinite range interactions (mean field systems). In short-range systems, a pseudo-spinodal can be defined by extrapolation of metastable measurements, but the point itself is not reachable because it lies beyond the metastability limit. In this work we show that a sensible definition of spinodal po…
▽ More
Though intuitively appealing, the concept of spinodal is rigourously defined only in systems with infinite range interactions (mean field systems). In short-range systems, a pseudo-spinodal can be defined by extrapolation of metastable measurements, but the point itself is not reachable because it lies beyond the metastability limit. In this work we show that a sensible definition of spinodal points can be obtained through the short time dynamical behavior of the system deep inside the metastable phase, by looking for a point where the system shows critical behavior. We show that spinodal points obtained by this method agree both with the thermodynamical spinodal point in mean field systems and with the pseudo-spinodal point obtained by extrapolation of meta-equilibrium behavior in short range systems. With this definition, a practical determination can be achieved without regard for equilibration issues.
△ Less
Submitted 5 August, 2009; v1 submitted 17 November, 2008;
originally announced November 2008.
-
Thermodynamic signature of growing amorphous order in glass-forming liquids
Authors:
G. Biroli,
J. -P. Bouchaud,
A. Cavagna,
T. S. Grigera,
P. Verrocchio
Abstract:
Although several theories relate the steep slowdown of glass formers to increasing spatial correlations of some sort, standard static correlation functions show no evidence for this. We present results that reveal for the first time a qualitative thermodynamic difference between the high temperature and deeply supercooled equilibrium glass-forming liquid: the influence of boundary conditions pro…
▽ More
Although several theories relate the steep slowdown of glass formers to increasing spatial correlations of some sort, standard static correlation functions show no evidence for this. We present results that reveal for the first time a qualitative thermodynamic difference between the high temperature and deeply supercooled equilibrium glass-forming liquid: the influence of boundary conditions propagates into the bulk over larger and larger lengthscales upon cooling, and, as this static correlation length grows, the influence decays nonexponentially. Increasingly long-range susceptibility to boundary conditions is expected within the random firt-order theory (RFOT) of the glass transition, but a quantitative account of our numerical results requires a generalization of RFOT where the surface tension between states fluctuates.
△ Less
Submitted 28 May, 2008;
originally announced May 2008.
-
Mosaic multi-state scenario vs. one-state description of supercooled liquids
Authors:
Andrea Cavagna,
Tomas S. Grigera,
Paolo Verrocchio
Abstract:
According to the mosaic scenario, relaxation in supercooled liquids is ruled by two competing mechanisms: surface tension, opposing the creation of local excitations, and entropy, providing the drive to the configurational rearrangement of a given region. We test this scenario through numerical simulations well below the Mode Coupling temperature. For an equilibrated configuration, we freeze all…
▽ More
According to the mosaic scenario, relaxation in supercooled liquids is ruled by two competing mechanisms: surface tension, opposing the creation of local excitations, and entropy, providing the drive to the configurational rearrangement of a given region. We test this scenario through numerical simulations well below the Mode Coupling temperature. For an equilibrated configuration, we freeze all the particles outside a sphere and study the thermodynamics of this sphere. The frozen environment acts as a pinning field. Measuring the overlap between the unpinned and pinned equilibrium configurations of the sphere, we can see whether it has switched to a different state. We do not find any clear evidence of the mosaic scenario. Rather, our results seem compatible with the existence of a single (liquid) state. However, we find evidence of a growing static correlation length, apparently unrelated to the mosaic one.
△ Less
Submitted 7 May, 2007; v1 submitted 31 July, 2006;
originally announced July 2006.
-
Geometrical properties of the potential energy of the soft-sphere binary mixture
Authors:
Tomas S. Grigera
Abstract:
We report a detailed study of the stationary points (zero-force points) of the potential energy surface (PES) of a model structural glassformer. We compare stationary points found with two different algorithms (eigenvector following and square gradient minimization), and show that the mapping between instantaneous configuration and stationary points defined by those algorithms is as different as…
▽ More
We report a detailed study of the stationary points (zero-force points) of the potential energy surface (PES) of a model structural glassformer. We compare stationary points found with two different algorithms (eigenvector following and square gradient minimization), and show that the mapping between instantaneous configuration and stationary points defined by those algorithms is as different as to strongly influence the instability index K vs. temperature plot, which relevance in analyzing the liquid dynamics is thus questioned. On the other hand, the plot of K vs. energy is much less sensitive to the algorithm employed, showing that the energy is the good variable to discuss geometric properties of the PES. We find new evidence of a geometric transition between a minima-dominated phase and a saddle-point-dominated one. We analyze the distances between instantaneous configurations and stationary points, and find that above the glass transition, the system is closer to saddle points than to minima.
△ Less
Submitted 12 September, 2005;
originally announced September 2005.
-
Localization threshold of Instantaneous Normal Modes from level-spacing statistics
Authors:
Stefano Ciliberti,
Tomas S. Grigera
Abstract:
We study the statistics of level-spacing of Instantaneous Normal Modes in a supercooled liquid. A detailed analysis allows to determine the mobility edge separating extended and localized modes in the negative tail of the density of states. We find that at temperature below the mode coupling temperature only a very small fraction of negative eigenmodes are localized.
We study the statistics of level-spacing of Instantaneous Normal Modes in a supercooled liquid. A detailed analysis allows to determine the mobility edge separating extended and localized modes in the negative tail of the density of states. We find that at temperature below the mode coupling temperature only a very small fraction of negative eigenmodes are localized.
△ Less
Submitted 20 April, 2004;
originally announced April 2004.
-
Anderson Localization in Euclidean Random Matrices
Authors:
S. Ciliberti,
T. S. Grigera,
V. Martin-Mayor,
G. Parisi,
P. Verrocchio
Abstract:
We study spectra and localization properties of Euclidean random matrices. The problem is approximately mapped onto that of a matrix defined on a random graph. We introduce a powerful method to find the density of states and the localization threshold. We solve numerically an exact equation for the probability distribution function of the diagonal element of the the resolvent matrix, with a popu…
▽ More
We study spectra and localization properties of Euclidean random matrices. The problem is approximately mapped onto that of a matrix defined on a random graph. We introduce a powerful method to find the density of states and the localization threshold. We solve numerically an exact equation for the probability distribution function of the diagonal element of the the resolvent matrix, with a population dynamics algorithm, and we show how this can be used to find the localization threshold. An application of the method in the context of the Instantaneous Normal Modes of a liquid system is given.
△ Less
Submitted 3 March, 2004;
originally announced March 2004.
-
The Boson peak and the phonons in glasses
Authors:
S. Ciliberti,
T. S. Grigera,
V. Martin-Mayor,
G. Parisi,
P. Verrocchio
Abstract:
Despite the presence of topological disorder, phonons seem to exist also in glasses at very high frequencies (THz) and they remarkably persist into the supercooled liquid. A universal feature of such a systems is the Boson peak, an excess of states over the standard Debye contribution at the vibrational density of states. Exploiting the euclidean random matrix theory of vibrations in amorphous s…
▽ More
Despite the presence of topological disorder, phonons seem to exist also in glasses at very high frequencies (THz) and they remarkably persist into the supercooled liquid. A universal feature of such a systems is the Boson peak, an excess of states over the standard Debye contribution at the vibrational density of states. Exploiting the euclidean random matrix theory of vibrations in amorphous systems we show that this peak is the signature of a phase transition in the space of the stationary points of the energy, from a minima-dominated phase (with phonons) at low energy to a saddle-point dominated phase (without phonons). The theoretical predictions are checked by means of numeric simulations.
△ Less
Submitted 2 December, 2003;
originally announced December 2003.
-
Asymptotic aging in structural glasses
Authors:
T. S. Grigera,
V. Martin-Mayor,
G. Parisi,
P. Verrocchio
Abstract:
Using a non-local Monte Carlo algorithm, we study the aging of a fragile glass, being able to follow it up to equilibibrium down to 0.89 Tmc (Tmc is the Mode-Coupling temperature) and up to unprecedentedly large waiting times at lower temperatures. We show that the fluctuation-dissipation ratio is independent of the dynamics chosen and is compatible with a phase transition, and that the scaling…
▽ More
Using a non-local Monte Carlo algorithm, we study the aging of a fragile glass, being able to follow it up to equilibibrium down to 0.89 Tmc (Tmc is the Mode-Coupling temperature) and up to unprecedentedly large waiting times at lower temperatures. We show that the fluctuation-dissipation ratio is independent of the dynamics chosen and is compatible with a phase transition, and that the scaling behaviour of the aging part of the correlation supports the full-aging scenario.
△ Less
Submitted 13 April, 2004; v1 submitted 25 July, 2003;
originally announced July 2003.