Search

We analyse biased ensembles of trajectories for the random-field Ising model on a fully-connected lattice, which is described exactly by mean-field theory. By coupling the activity of the system to a dynamical biasing field, we find a range of dynamical phase transitions, including spontaneous symmetry breaking into ordered states. For weak bias, the phase behaviour is controlled by extrema of the free energy, which may be local minima or saddle points. For large bias, the system tends to states of extremal activity, which may differ strongly from free energy minima. We discuss connections of these results to random first-order transition theory of glasses, which motivates an extension of the analysis to random-field Ising models where the dynamical activity is not symmetric under magnetisation reversal., 26 pages, 10 figures

We perform numerical simulations of a long-range spherical spin glass with two and three body interaction terms. We study the gradient descent dynamics and the inherent structures found after a quench from initial conditions, well thermalized at temperature $T_{in}$. In large systems, the dynamics strictly agrees with the integration of the mean-field dynamical equations. In particular, we confirm the existence of an onset initial temperature, within the liquid phase, below which the energy of the inherent structures undoubtedly depends on $T_{in}$. This behavior is in contrast with that of pure models, where there is a 'threshold energy' that attracts all the initial configurations in the liquid. Our results strengthen the analogy between mean-field spin glass models and supercooled liquids., 15 pages, 8 figures

A mean-field kinetic model suggests that the relaxation dynamics of wormlike micellar networks is a long and complex process due to the problem of reducing the number of free end-caps (or dangling ends) while also reaching an equilibrium level of branching after an earlier overgrowth. The model is validated against mesoscopic molecular dynamics simulations and is based on kinetic equations accounting for scission and synthesis processes of blobs of surfactants. A long relaxation time scale is reached both with thermal quenches and small perturbations of the system. The scaling of this relaxation time is exponential with the free energy of an end cap and with the branching free energy. We argue that the subtle end-recombination dynamics might yield effects that are difficult to detect in rheology experiments, with possible underestimates of the typical time scales of viscoelastic fluids., Comment: 10 pages, 6 figures

We introduce a kinetically constrained spin model with a local softness parameter, such that spin flips can violate the kinetic constraint with an (annealed) site-dependent rate. We show that adding MC swap moves to this model can dramatically accelerate structural relaxation. We discuss the connection of this observation with the fact that swap moves are also able to accelerate relaxation in structural glasses. We analyse the rates of relaxation in the model. We also show that the extent of dynamical heterogeneity is strongly suppressed by the swap moves., Comment: 20 pages, accepted version

We study the evolution of the maximum energy $E_\max(t)$ reached between time $0$ and time $t$ in the dynamics of simple models with glassy energy landscapes, in instant quenches from infinite temperature to a target temperature $T$. Through a detailed description of the activated dynamics, we are able to describe the evolution of $E_\max(t)$ from short times, through the aging regime, until after equilibrium is reached, thus providing a detailed description of the long-time dynamics. Finally, we compare our findings with numerical simulations of the $p$-spin glass and show how the maximum energy record can be used to identify the threshold energy in this model., 22 pages, 5 figures; revised presentation

Colloidal suspensions have long been studied as a model for atomic and molecular systems, due to the ability to fluorescently label and individually track each particle, yielding particle-resolved structural information. This allows various local order parameters to be probed that are otherwise inaccessible for a comparable molecular system. For phase transitions such as crystallisation, appropriate order parameters which emphasise 6-fold symmetry are a natural choice, but for vitrification the choice of order parameter is less clear cut. Previous work has highlighted the importance of icosahedral local structure as the glass transition is approached. However, counting icosahedra or related motifs is not a continuous order parameter in the same way as, for example, the bond-orientational order parameters $Q_{6}$ and $W_6$. In this work we investigate the suitability of using pentagonal bipyramid membership, a structure which can be assembled into larger, five-fold symmetric structures, as a finer order parameter to investigate the glass transition. We explore various structural and dynamic properties and show that this new approach produces many of the same findings as simple icosahedral membership, but we also find that large instantaneous displacements are often correlated with significant changes in pentagonal bipyramid membership, and unlike the population of defective icosahedra, the pentagonal bypyramid membership and spindle number do not saturate for any measured volume fraction, but continues to increase., Comment: accepted by JStat Mech: Theory and Experiment 2019

We have introduced a variational method to improve the computation of integrated correlation times in the Parallel Tempering Dynamics, obtaining a better estimate (a lower bound, at least) of the exponential correlation time. Using this determination of the correlation times, we revisited the problem of the characterization of the chaos in temperature in finite dimensional spin glasses by means of the study of correlations between different chaos indicators computed in the static and the correlation times of the Parallel Tempering dynamics. The sample-distribution of the characteristic time for the Parallel Tempering dynamics turns out to be fat-tailed and it obeys finite-size scaling., Comment: New appendix added. Version to appear in JSTAT

We study the out-of-equilibrium aging dynamics of the Random Energy Model (REM) ruled by a single spin-flip Metropolis dynamics. We focus on the dynamical evolution taking place on time-scales diverging with the system size. Our aim is to show to what extent the activated dynamics displayed by the REM can be described in terms of an effective trap model. We identify two time regimes: the first one corresponds to the process of escaping from a basin in the energy landscape and to the subsequent exploration of high energy configurations, whereas the second one corresponds to the evolution from a deep basin to the other. By combining numerical simulations with analytical arguments we show why the trap model description does not hold in the former but becomes exact in the second., 26 pages, 12 figures. Version accepted for publication in JSTAT

In recent years, studies of long-range interacting (LRI) systems have taken centre stage in the arena of statistical mechanics and dynamical system studies, due to new theoretical developments involving tools from as diverse a field as kinetic theory, non-equilibrium statistical mechanics, and large deviation theory, but also due to new and exciting experimental realizations of LRI systems. In this invited contribution, we discuss the general features of long-range interactions, emphasizing in particular the main physical phenomenon of non-additivity, which leads to a plethora of distinct effects, both thermodynamic and dynamic, that are not observed with short-range interactions: Ensemble inequivalence, slow relaxation, broken ergodicity. We also discuss several physical systems with long-range interactions: mean-field spin systems, self-gravitating systems, Euler equations in two dimensions, Coulomb systems, one-component electron plasma, dipolar systems, free-electron lasers, atoms trapped in optical cavities., Invited contribution to a special issue on the occasion of Abdus Salam's 90th birthday; v2: Refs. added

We study numerically the coarsening kinetics of a two-dimensional ferromagnetic system with aleatory bond dilution. We show that interfaces between domains of opposite magnetisation are fractal on every lengthscale, but with different properties at short or long distances. Specifically, on lengthscales larger than the typical domains' size the topology is that of critical random percolation, similarly to what observed in clean systems or models with different kinds of quenched disorder. On smaller lengthscales a dilution dependent fractal dimension emerges. The Hausdorff dimension increases with increasing dilution $d$ up to the value $4/3$ expected at the bond percolation threshold $d=1/2$. We discuss how such different geometries develop on different lengthscales during the phase-ordering process and how their simultaneous presence determines the scaling properties of observable quantities., Comment: 15 pages, 9 figures

We investigate the relaxation dynamics of magnetic vortex lines in type-II superconductors following rapid changes of the external driving current by means of an elastic line model simulated with Langevin molecular dynamics. A system of flux vortices in a sample with randomly distributed point-like defects is subjected to an external current of appropriate strength for a sufficient period of time so as to be in a moving non-equilibrium steady state. The current is then instantaneously lowered to a value that pertains to either the moving or pinned regime. The ensuing relaxation of the flux lines is studied via one-time observables such as their mean velocity and radius of gyration. We have in addition measured the two-time flux line height autocorrelation function to investigate dynamical scaling and aging behavior in the system, which in particular emerge after quenches into the glassy pinned state., 16 pages, 5 figures

We investigate the static and dynamic properties of a weakly polydisperse hard sphere system in the deeply supercooled state, i.e. at densities higher than that corresponding to the mode-coupling transition. The structural analysis reveals the emergence of icosahedral locally favoured structures, previously only found in trace quantities. We present a new approach to probe the shape of dynamically heterogeneous regions, which is insensitive to particle packing effects that can hamper such analysis. Our results indicate that the shape of the dynamically heterogeneous regions changes only weakly and moreover hint that the often-used four-point correlation length may exhibit a growth in excess of that which our method identifies. The growth of the size of the dynamically heterogeneous regions appears instead to be in line with the one of structural and dynamic propensity correlations., Comment: under review at J. Stat Mech: Theory and Experiment

We consider the long-ranged Ising spin-glass with random couplings decaying as a power-law of the distance, in the region of parameters where the spin-glass phase exists with a positive droplet exponent. For the Metropolis single-spin-flip dynamics near zero temperature, we construct via real-space renormalization the full hierarchy of relaxation times of the master equation for any given realization of the random couplings. We then analyze the probability distribution of dynamical barriers as a function of the spatial scale. This real-space renormalization procedure represents a simple explicit example of the droplet scaling theory, where the convergence towards local equilibrium on larger and larger scales is governed by a strong hierarchy of activated dynamical processes, with valleys within valleys., v2=final version

Behavior of dynamical process of complex systems is investigated. Specifically we analyse two types of ideal complex systems. For analysing the ideal complex systems, we define the response functions describing the internal states to an external force. The internal states are obtained as a relaxation process showing a “power law” distribution, such as scale free behaviors observed in actual measurements. By introducing a hybrid system, the logarithmic time, and double logarithmic time, we show how the “slow relaxation” (SR) process and “super slow relaxation” (SSR) process occur. Regarding the irregular variations of the internal states as an activation process, we calculate the response function to the external force. The behaviors are classified into “power”, “exponential”, and “stretched exponential” type. Finally we construct a fractional differential equation (FDE) describing the time evolution of these complex systems. In our theory, the exponent of the FDE or that of the power law distribution is expressed in terms of the parameters characterizing the structure of the system.

Under the influence of an external field many systems exhibit slow relaxations processes. In that sense, we study the Physics of this behavior/process in the early stages of its evolution on bidimensional Ising systems using Monte Carlo simulations at a broad range of temperature and spins up concentration. The Monte Carlo simulations were done up to second neighbors interactions shown a more stable dynamic regime than first neighbors interactions, when the systems energy is considered. The behavior of relative energy/spin and the relative perimeter (the line between spin up and spin down regions) are monitored against Monte Carlo steps. Also, our findings show that the stretching index ( a 2) exhibit values below the borderline (1) in a broad range of concentration below 50% and temperature from 0.5 ≤ T ≤ 4.5 (in units of J / kB , where J is the ferromagnetic constant and kB is the Boltzmann constant).

Trap models have been initially proposed as toy models for dynamical relaxation in extremely simplified rough potential energy landscapes. Their importance has considerably grown recently thanks to the discovery that the trap like aging mechanism is directly controlling the out-of-equilibrium relaxation processes of more sophisticated spin models, that are considered as the solvable counterpart of real disordered systems. Establishing on a firmer ground the connection between these spin model out-of-equilibrium behavior and the trap like aging mechanism would shed new light on the properties, still largely mysterious, of the activated out-of-equilibrium dynamics of disordered systems. In this work we discuss numerical evidences of emergent trap-like aging behavior in a variety of disordered models. Our numerical results are backed by analytic derivations and heuristic discussions. Such exploration reveals some of the tricks needed to analyze the trap behavior in spite of the occurrence of secondary processes, of the existence of dynamical correlations and of finite system's size effects., Comment: 25 pages, 15 figures

The analytical description of the dynamics in models with discrete variables (e.g. Ising spins) is a notoriously difficult problem, that can be tackled only under some approximation. Recently a novel variational approach to solve the stationary dynamical regime has been introduced by Pelizzola [Eur. Phys. J. B, 86 (2013) 120], where simple closed equations are derived under mean-field approximations based on the cluster variational method. Here we propose to use the same approximation based on the cluster variational method also for the non-stationary regime, which has not been considered up to now within this framework. We check the validity of this approximation in describing the non-stationary dynamical regime of several Ising models defined on Erdos-R\'enyi random graphs: we study ferromagnetic models with symmetric and partially asymmetric couplings, models with random fields and also spin glass models. A comparison with the actual Glauber dynamics, solved numerically, shows that one of the two studied approximations (the so-called 'diamond' approximation) provides very accurate results in all the systems studied. Only for the spin glass models we find some small discrepancies in the very low temperature phase, probably due to the existence of a large number of metastable states. Given the simplicity of the equations to be solved, we believe the diamond approximation should be considered as the 'minimal standard' in the description of the non-stationary regime of Ising-like models: any new method pretending to provide a better approximate description to the dynamics of Ising-like models should perform at least as good as the diamond approximation., Comment: 21 pages, 12 figures