Dynamic Universality Class Research Articles

Non-equilibrium dynamics of the open quantum O(n)-model with non-Markovian noise: exact results

The collective and purely relaxational dynamics of quantum many-body systems after a quench at temperature T = 0, from a disordered state to various phases is studied through the exact solution of the quantum Langevin equation of the spherical and the O(n)-model in the limit n → ∞. The stationary state of the quantum dynamics is shown to be a non-equilibrium state. The quantum spherical and the quantum O(n)-model for n → ∞ are in the same dynamical universality class. The long-time behaviour of single-time and two-time correlation and response functions is analysed and the universal exponents which characterise quantum coarsening and quantum ageing are derived. The importance of the non-Markovian long-time memory of the quantum noise is elucidated by comparing it with an effective Markovian noise having the same scaling behaviour and with the case of non-equilibrium classical dynamics.

Dynamical Renormalization Group for Mode-Coupling Field Theories with Solenoidal Constraint

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.

Fluctuations of the critical Casimir force.

The critical Casimir force (CCF) arises from confining fluctuations in a critical fluid and thus it is a fluctuating quantity itself. While the mean CCF is universal, its (static) variance has previously been found to depend on the microscopic details of the system which effectively set a large-momentum cutoff in the underlying field theory, rendering it potentially large. This raises the question how the properties of the force variance are reflected in experimentally observable quantities, such as the thickness of a wetting film or the position of a suspended colloidal particle. Here, based on a rigorous definition of the instantaneous force, we analyze static and dynamic correlations of the CCF for a conserved fluid in film geometry for various boundary conditions within the Gaussian approximation. We find that the dynamic correlation function of the CCF is independent of the momentum cutoff and decays algebraically in time. Within the Gaussian approximation, the associated exponent depends only on the dynamic universality class but not on the boundary conditions. We furthermore consider a fluid film, the thickness of which can fluctuate under the influence of the time-dependent CCF. The latter gives rise to an effective non-Markovian noise in the equation of motion of the film boundary and induces a distinct contribution to the position variance. Within the approximations used here, at short times, this contribution grows algebraically in time whereas, at long times, it saturates and contributes to the steady-state variance of the film thickness.

Effect of attractions on the hard-sphere dynamic universality class.

The fundamental understanding of the dynamic and transport properties of liquids is crucial for the better processing of most materials. The usefulness of this understanding increases when it involves general scaling rules, such as the concept of the hard-sphere dynamic universality class, which provides a unifying scaling of the dynamics of soft-sphere repulsive systems. A relevant question is how far this concept extends to systems that also involve attractive interactions. To answer this question, in this work we performed systematic molecular and Brownian dynamics simulations with the Lennard-Jones system in a wide range of temperatures and densities and verify the extent to which its static and dynamic properties map onto those of the hard-sphere system. We determine that in most of the fluid regime, the Lennard-Jones liquid exhibits the same dynamic equivalence with the hard-sphere system as most purely repulsive fluids, thus establishing the degree of its inclusion in the hard-sphere dynamic universality class.

Equilibrium to off-equilibrium crossover in homogeneous active matter

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.

Symmetries and effect of time dimension in non-equilibrium quantum matter

Non-equilibrium quantum many-body systems have attracted considerable attention in the past decades. The scope of the research of this kind of novel system involves interdisciplinary research of condensed matter, atomic and molecular physics, quantum optics, quantum information and quantum computation, as well as the non-equilibrium statistical physics. The non-equilibrium phenomena emerging from the aforementioned quantum systems can exhibit rich and universal behaviors, which have far from being well understood due to the novelties and complexities of these systems, and hence the quantum many-body physics becomes the research highlight. At the same time, with the rapid development of quantum techniques, the understanding of these complex systems is of important practical significance due to their potential applications in quantum computation and quantum manipulation. In this paper, we show our recent progress of non-equilibrium quantum many-body systems. We focus on the novel phenomena closely related to the temporary symmetry breaking, including the exotic quantum matter, quasi-particles as well as the dynamical universality classes in non-equilibrium quantum many-body systems.

Spectral functions and dynamic critical behavior of relativistic Z2 theories

We investigate the dynamic critical behavior of a relativistic scalar field theory with Z2 symmetry by calculating spectral functions of the order parameter at zero and non-vanishing momenta from first-principles classical-statistical lattice simulations in real-time. We find that at temperatures above the critical point (T>Tc), the spectral functions are well described by relativistic quasi-particle peaks. Close to the transition temperature (T∼Tc), we observe strong infrared contributions building up. In the ordered phase at low temperatures (T<Tc), in addition to the quasi-particle peak, we observe a soft mode with a dispersion relation indicative of collective excitations. Investigating the spectral functions close to Tc, we demonstrate that the behavior in the vicinity of the critical point is controlled by dynamic scaling functions and the dynamic critical exponent z, which we determine from our simulations. By considering the equations of motion for a closed system and a system coupled to a heat bath, we extract the dynamic critical behavior for two different dynamic universality classes (Models A & C) in two and three spatial dimensions.

Too Close to Integrable: Crossover from Normal to Anomalous Heat Diffusion.

Energy transport in one-dimensional chains of particles with three conservation laws is generically anomalous and belongs to the Kardar-Parisi-Zhang dynamical universality class. Surprisingly, some examples where an apparent normal heat diffusion is found over a large range of length scales were reported. We propose a novel physical explanation of these intriguing observations. We develop a scaling analysis that explains how this may happen in the vicinity of an integrable limit, such as, but not only, the famous Toda model. In this limit, heat transport is mostly supplied by quasiparticles with a very large mean free path ℓ. Upon increasing the system size L, three different regimes can be observed: a ballistic one, an intermediate diffusive range, and, eventually, the crossover to the anomalous (hydrodynamic) regime. Our theoretical considerations are supported by numerical simulations of a gas of diatomic hard-point particles for almost equal masses and of a weakly perturbed Toda chain. Finally, we discuss the case of the perturbed harmonic chain, which exhibits a yet different scenario.

Kardar-Parisi-Zhang universality from soft gauge modes

The emergence of superdiffusive spin dynamics in integrable classical and quantum magnets is well established by now, but there is no generally valid theoretical explanation for this phenomenon. A fundamental difficulty is that the hydrodynamic fluctuations of conserved quasiparticle modes are purely diffusive. We argue that in isotropic integrable magnets, a complete hydrodynamic description must include soft "gauge" degrees of freedom, that arise from spontaneous breaking of the Bethe pseudovacuum symmetry. We show that the coarse-grained time evolution of these modes lies in the Kardar-Parisi-Zhang universality class of dynamics.

Epidemic outbreaks on random Voronoi–Delaunay triangulations

We study epidemic outbreaks on random Delaunay triangulations by applying the Asynchronous SIR (susceptible–infected–removed) dynamics coupled to two-dimensional Voronoi–Delaunay triangulations. In order to investigate the critical behavior of the model, we obtain the cluster size distribution by using Newman–Ziff algorithm, allowing to simulate random inhomogeneous lattices and measure any desired observable related to percolation. We numerically calculate the order parameter, defined as the wrapping cluster density, the mean cluster size, and Binder cumulant ratio defined for percolation in order to estimate the epidemic threshold. Our findings suggest that the system falls into two-dimensional dynamic percolation universality class and the quenched random disorder is irrelevant, in agreement with results for classical percolation.

Kardar-Parisi-Zhang universality of the Nagel-Schreckenberg model.

Dynamical universality classes are distinguished by their dynamical exponent z and unique scaling functions encoding space-time asymmetry for, e.g., slow-relaxation modes or the distribution of time-integrated currents. So far the universality class of the Nagel-Schreckenberg (NaSch) model, which is a paradigmatic model for traffic flow on highways, was not known. Only the special case v_{max}=1, where the model corresponds to the totally asymmetric simple exclusion process, is known to belong to the superdiffusive Kardar-Parisi-Zhang (KPZ) class with z=3/2. In this paper, we show that the NaSch model also belongs to the KPZ class for general maximum velocities v_{max}>1. Using nonlinear fluctuating hydrodynamics theory we calculate the nonuniversal coefficients, fixing the exact asymptotic solutions for the dynamical structure function and the distribution of time-integrated currents. The results of large-scale Monte Carlo simulations match the exact asymptotic KPZ solutions without any fitting parameter left. Additionally, we find that nonuniversal early-time effects or the choice of initial conditions might have a strong impact on the numerical determination of the dynamical exponent and therefore lead to inconclusive results. We also show that the universality class is not changed by extending the model to a two-lane NaSch model with lane-changing rules.

Dynamics of coupled modes for sliding particles on a fluctuating landscape.

The recently developed formalism of nonlinear fluctuating hydrodynamics (NLFH) has been instrumental in unraveling many new dynamical universality classes in coupled driven systems with multiple conserved quantities. In principle, this formalism requires knowledge of the exact expression of locally conserved current in terms of local density of the conserved components. However, for most nonequilibrium systems an exact expression is not available and it is important to know what happens to the predictions of NLFH in these cases. We address this question here in a system with coupled time evolution of sliding particles on a fluctuating energy landscape. In the disordered phase this system shows short-ranged correlations, the exact form of which is not known, and so the exact expression for current cannot be obtained. We use approximate expressions based on mean-field theory and corrections to it, to test the prediction of NLFH using numerical simulations. In this process we also discover important finite size effects and show how they affect the predictions of NLFH. We find that our system is rich enough to show a large variety of universality classes. From our analytics and simulations we have been able to find parameter values which lead to diffusive, Karder-Parisi-Zhang (KPZ), 5/3 Lévy, and modified KPZ universality classes. Interestingly, the scaling function in the modified KPZ case turns out to be close to the Prähofer-Spohn function, which is known to describe usual KPZ scaling. Our analytics also predict the golden mean and the 3/2 Lévy universality classes within our model but our simulations could not verify this, perhaps due to strong finite size effects.

Measurement-Induced Phase Transitions in the Dynamics of Entanglement

We define dynamical universality classes for many-body systems whose unitary evolution is punctuated by projective measurements. In cases where such measurements occur randomly at a finite rate $p$ for each degree of freedom, we show that the system has two dynamical phases: `entangling' and `disentangling'. The former occurs for $p$ smaller than a critical rate $p_c$, and is characterized by volume-law entanglement in the steady-state and `ballistic' entanglement growth after a quench. By contrast, for $p > p_c$ the system can sustain only area-law entanglement. At $p = p_c$ the steady state is scale-invariant and, in 1+1D, the entanglement grows logarithmically after a quench. To obtain a simple heuristic picture for the entangling-disentangling transition, we first construct a toy model that describes the zeroth R\'{e}nyi entropy in discrete time. We solve this model exactly by mapping it to an optimization problem in classical percolation. The generic entangling-disentangling transition can be diagnosed using the von Neumann entropy and higher R\'{e}nyi entropies, and it shares many qualitative features with the toy problem. We study the generic transition numerically in quantum spin chains, and show that the phenomenology of the two phases is similar to that of the toy model, but with distinct `quantum' critical exponents, which we calculate numerically in $1+1$D. We examine two different cases for the unitary dynamics: Floquet dynamics for a nonintegrable Ising model, and random circuit dynamics. We obtain compatible universal properties in each case, indicating that the entangling-disentangling phase transition is generic for projectively measured many-body systems. We discuss the significance of this transition for numerical calculations of quantum observables in many-body systems.

Does the chiral magnetic effect change the dynamic universality class in QCD?

In QCD matter under an external magnetic field, the chiral magnetic effect (CME) leads to the collective gapless mode called the chiral magnetic wave (CMW). Since dynamic universality class generally depends on low-energy gapless modes, it is nontrivial whether the CME and the resulting CMW change that of the second-order chiral phase transition in QCD. To address this question, we study the critical dynamics near the chiral phase transition in massless two-flavor QCD under an external magnetic field. By performing the dynamic renormalization-group analysis within the ϵ expansion, we find that the presence of the CME changes the dynamic universality class to that of model A. We also show that the transport coefficient of the CME is not renormalized by the critical fluctuations of the order parameter.

Numerical simulations of critical dynamics in anisotropic magnetic films with the stochastic Landau-Lifshitz-Gilbert equation.

With the stochastic Landau-Lifshitz-Gilbert (sLLG) equation, critical dynamic behaviors far from equilibrium or stationary around the order-disorder and pinning-depinning phase transitions in anisotropic magnetic films are investigated. From the dynamic relaxation with and without an external field, the Curie temperature and critical exponents of the order-disorder phase transition are accurately determined. For the pinning-depinning phase transition induced by quenched disorder, the nonstationary creep motion of domain wall activated by finite temperatures is simulated, and the thermal rounding exponent is extracted. The results show that the dynamic universality class of the sLLG equation is different from those of the Monte Carlo dynamics and quenched Edwards-Wilkinson equation, and it may lead to alternative understanding of experiments. The dynamic approach shows its great efficiency for the sLLG equation.

Aging dynamics in quenched noisy long-range quantum Ising models

We consider the $d$-dimensional transverse-field Ising model with power-law interactions $J/r^{d+\sigma}$ in the presence of a noisy longitudinal field with zero average. We study the longitudinal-magnetization dynamics of an initial paramagnetic state after a sudden switch-on of both the interactions and the noisy field. While the system eventually relaxes to an infinite-temperature state with vanishing magnetization correlations, we find that two-time correlation functions show aging at intermediate times. Moreover, for times shorter than the inverse noise strength $\kappa$ and distances longer than $a(J/\kappa)^{2/\sigma}$ with $a$ being the lattice spacing, we find a critical scaling regime of correlation and response functions consistent with the model A dynamical universality class with an initial-slip exponent $\theta=1$ and dynamical critical exponent $z=\sigma/2$. We obtain our results analytically by deriving an effective action for the magnetization field including the noise in a non-perturbative way. The above scaling regime is governed by a non-equilibrium fixed point dominated by the noise fluctuations.

Weak ergodicity breaking from quantum many-body scars

Certain wave functions of non-interacting quantum chaotic systems can exhibit "scars" in the fabric of their real-space density profile. Quantum scarred wave functions concentrate in the vicinity of unstable periodic classical trajectories. We introduce the notion of many-body quantum scars which reflect the existence of a subset of special many-body eigenstates concentrated in certain parts of the Hilbert space. We demonstrate the existence of scars in the Fibonacci chain -- the one- dimensional model with a constrained local Hilbert space realized in the 51 Rydberg atom quantum simulator [H. Bernien et al., arXiv:1707.04344]. The quantum scarred eigenstates are embedded throughout the thermalizing many-body spectrum, but surprisingly lead to direct experimental signatures such as robust oscillations following a quench from a charge-density wave state found in experiment. We develop a model based on a single particle hopping on the Hilbert space graph, which quantitatively captures the scarred wave functions up to large systems of L = 32 atoms. Our results suggest that scarred many-body bands give rise to a new universality class of quantum dynamics, which opens up opportunities for creating and manipulating novel states with long-lived coherence in systems that are now amenable to experimental study.

Unconventional Universality Class of One-Dimensional Isolated Coarsening Dynamics in a Spinor Bose Gas.

By studying the coarsening dynamics of a one-dimensional spin-1 Bose-Hubbard model in a superfluid regime, we analytically find an unconventional universal dynamical scaling for the growth of the spin correlation length, which is characterized by the exponential integral unlike the conventional power law or simple logarithmic behavior, and numerically confirmed with the truncated Wigner approximation.

Dynamic universality class of model H with frustrated diffusion: ε expansion

We study a variation of the dynamic universality class of model H in a spatial dimension of $d=4-\epsilon$, by frustrating charge diffusion and momentum density fluctuations along $d_T=1$ or $d_T=2$ dimensions, while keeping the same dynamics of model H in the other $d_L=d-d_T$ dimensions. The case of $d_T=2$ describes the QCD critical point in a background magnetic field. We find that these models belong to a different dynamical universality class due to extended conservation laws compared to the model H, although the static universality class remains the same as the 3-dimensional Ising model. We compute the dynamic critical exponents of these models in first order of $\epsilon$-expansion to find that $x_\lambda\approx 0.847\,\epsilon$, $x_{\bar\eta}\approx 0.153\,\epsilon$, and $z=4-x_\lambda\approx 3.15$ when $\epsilon=1$ and $d_T=2$. For $d_T=1$ the results are numerically similar to the model H values: $z\approx 3.08$.

Dynamical universality classes of simple growth and lattice gas models

Large scale, dynamical simulations have been performed for the two dimensional octahedron model, describing the Kardar–Parisi–Zhang (KPZ) for nonlinear, or the Edwards–Wilkinson class for linear surface growth. The autocorrelation functions of the heights and the dimer lattice gas variables are determined with high precision. Parallel random-sequential (RS) and two-sub-lattice stochastic dynamics (SCA) have been compared. The latter causes a constant correlation in the long time limit, but after subtracting it one can find the same height functions as in case of RS. On the other hand the ordered update alters the dynamics of the lattice gas variables, by increasing (decreasing) the memory effects for nonlinear (linear) models with respect to RS. Additionally, we support the KPZ ansatz and the Kallabis–Krug conjecture in dimensions and provide a precise growth exponent value . We show the emergence of finite size corrections, which occur long before the steady state roughness is reached.