Dynamical Critical Exponent Research Articles

Disorder-Induced Quantum Griffiths Singularity Revealed in an Artificial 2D Superconducting System.

Disorder‐induced Griffiths singularity of quantum phase transition (QPT) is a crucial issue in 2D superconductors (2DSC). In a superconducting system, the strength of disorder is found to be associated with the vortex pinning energy, which is closely related to the quantum Griffiths singularity; however, a direct study to elucidate the role of vortex pinning energy on the quantum Griffiths singularity in 2DSC remains to be undertaken. Here, an artificial 2DSC system is designed by randomly depositing superconducting nanoislands on 2Delectron gas (2DEG). Quantum Griffiths singularity is present in a graphene/Pb‐islands‐array hybrid, where the superconducting behavior transits to weakly localized metallic behavior induced by the vertical magnetic field and exhibits critical behavior with a diverging dynamical critical exponent approaching zero temperature. Compared to the study of graphene/Sn‐islands‐array hybrid where the sharp QPT is observed, the vortex pinning energy acquired from the Arrhenius plot analysis is greater in graphene/Pb‐islands‐array hybrid, which may contribute to the presence of the quantum Griffiths singularity. This work may provide a comprehensive interpretation of the QPT in 2DSC.

Bilocal quantum criticality

We consider 2+1 dimensional conformal gauge theories coupled to additional degrees of freedom which induce a spatially local but long-range in time $1/(\tau-\tau')^2$ interaction between gauge-neutral local operators. Such theories have been argued to describe the hole-doped cuprates near optimal doping. We focus on a SU(2) gauge theory with $N_h$ flavors of adjoint Higgs fields undergoing a quantum transition between Higgs and confining phases: the $1/(\tau-\tau')^2$ interaction arises from a spectator large Fermi surface of electrons. The large $N_h$ expansion leads to an effective action containing fields which are bilocal in time but local in space. We find a strongly-coupled fixed point at order $1/N_h$, with dynamic critical exponent $z > 1$. We show that the entropy preserves hyperscaling, but nevertheless leads to a linear in temperature specific heat with a co-efficient which has a finite enhancement near the quantum critical point.

Antiferromagnetism, spin-glass state, H–T phase diagram, and inverse magnetocaloric effect in Co2RuO4

Static and dynamic magnetic properties of normal spinel Co2RuO4 = (Co2+) are reported based on our investigations of the temperature (T), magnetic field (H) and frequency (f) dependence of the ac-magnetic susceptibilities and dc-magnetization (M) covering the temperature range T = 2 K–400 K and H up to 90 kOe. These investigations show that Co2RuO4 exhibits an antiferromagnetic (AFM) transition at TN ∼ 15.2 K, along with a spin-glass state at slightly lower temperature (TSG) near 14.2 K. It is argued that TN is mainly governed by the ordering of the spins of Co2+ ions occupying the A-site, whereas the exchange interaction between the Co2+ ions on the A-site and randomly distributed Ru3+ on the B-site triggers the spin-glass phase, Co3+ ions on the B-site being in the low-spin non-magnetic state. Analysis of measurements of M (H, T) for T < TN are used to construct the H–T phase diagram showing that TSG shifts to lower T varying as H2/3.2 expected for spin-glass state whereas TN is nearly H-independent. For T > TN, analysis of the paramagnetic susceptibility (χ) vs. T data are fit to the modified Curie–Weiss law, χ = χ0 + C/(T + θ), with χ0 = 0.0015 emu mol−1Oe−1 yielding θ = 53 K and C = 2.16 emu-K mol−1Oe−1, the later yielding an effective magnetic moment μeff = 4.16 μB comparable to the expected value of μeff = 4.24 μB per Co2RuO4. Using TN, θ and high temperature series for χ, dominant exchange constant J1/kB ∼ 6 K between the Co2+ on the A-sites is estimated. Analysis of the ac magnetic susceptibilities near TSG yields the dynamical critical exponent zν = 5.2 and microscopic spin relaxation time τ0 ∼ 1.16 × 10−10 sec characteristic of cluster spin-glasses and the observed time-dependence of M(t) is supportive of the spin-glass state. Large M–H loop asymmetry at low temperatures with giant exchange bias effect (HEB ∼ 1.8 kOe) and coercivity (HC ∼ 7 kOe) for a field cooled sample further support the mixed magnetic phase nature of this interesting spinel. The negative magnetocaloric effect observed below TN is interpreted to be due to the AFM and SG ordering. It is argued that the observed change from positive MCE (magnetocaloric effect) for T > TN to inverse MCE for T < TN observed in Co2RuO4 (and reported previously in other systems also) is related to the change in sign of (∂M/∂T) vs. T data.

Super slowing down in the bond-diluted Ising model.

In models in statistical physics, the dynamics often slows down tremendously near the critical point. Usually, the correlation time τ at the critical point increases with system size L in power-law fashion: τ∼L^{z}, which defines the critical dynamical exponent z. We show that this also holds for the two-dimensional bond-diluted Ising model in the regime p>p_{c}, where p is the parameter denoting the bond concentration, but with a dynamical critical exponent z(p) which shows a strong p dependence. Moreover, we show numerically that z(p), as obtained from the autocorrelation of the total magnetization, diverges when the percolation threshold p_{c}=1/2 is approached: z(p)-z(1)∼(p-p_{c})^{-2}. We refer to this observed extremely fast increase of the correlation time with size as super slowing down. Independent measurement data from the mean-square deviation of the total magnetization, which exhibits anomalous diffusion at the critical point, support this result.

Nonequilibrium Criticality in Quench Dynamics of Long-Range Spin Models.

Long-range interacting spin systems are ubiquitous in physics and exhibit a variety of ground-state disorder-to-order phase transitions. We consider a prototype of infinite-range interacting models known as the Lipkin-Meshkov-Glick model describing the collective interaction of N spins and investigate the dynamical properties of fluctuations and correlations after a sudden quench of the Hamiltonian. Specifically, we focus on critical quenches, where the initial state and/or the postquench Hamiltonian are critical. Depending on the type of quench, we identify three distinct behaviors where both the short-time dynamics and the stationary state at long times are effectively thermal, quantum, and genuinely nonequilibrium, characterized by distinct universality classes and static and dynamical critical exponents. These behaviors can be identified by an infrared effective temperature that is finite, zero, and infinite (the latter scaling with the system size as N^{1/3}), respectively. The quench dynamics is studied through a combination of exact numerics and analytical calculations utilizing the nonequilibrium Keldysh field theory. Our results are amenable to realization in experiments with trapped-ion experiments where long-range interactions naturally arise.

On the universality of the scaling relations during sol-gel transition

The evolution of viscoelastic properties near the sol-gel transition is studied by performing oscillatory rheological measurements on two different types of systems: a colloidal dispersion and a thermoresponsive polymer solution under isothermal and non-isothermal conditions. While undergoing sol-gel transition, both the systems pass through a critical point. An approach to the critical point is characterized in terms of divergence of zero-shear viscosity and the subsequent appearance of the low-frequency modulus. In the vicinity of the critical gel state, both the viscosity and the modulus show a power-law dependence on the relative distance from the critical point. Interestingly, the longest relaxation time has been observed to diverge symmetrically on both sides of the critical point and also shows a power-law dependence on relative distance from the critical gel point. The critical (power-law) exponents of the zero-shear viscosity and modulus are observed to be related to the exponents of the longest relaxation time by the hyperscaling laws. The dynamic critical exponent has also been calculated from the growth of the dynamic moduli. Remarkably, the critical relaxation exponent and dynamic critical exponent predicted from the scaling laws precisely agree with the experimental values from the isothermal as well as non-isothermal experiments. The associated critical exponents show remarkable internal consistency and advocate the universality of scaling relations for different kinds of systems undergoing the sol-gel transition.

Bose-glass phase of a one-dimensional disordered Bose fluid: Metastable states, quantum tunneling, and droplets.

We study a one-dimensional disordered Bose fluid using bosonization, the replica method, and a nonperturbative functional renormalization-group approach. We find that the Bose-glass phase is described by a fully attractive strong-disorder fixed point characterized by a singular disorder correlator whose functional dependence assumes a cuspy form that is related to the existence of metastable states. At nonzero momentum scale k, quantum tunneling between the ground state and low-lying metastable states leads to a rounding of the cusp singularity into a quantum boundary layer (QBL). The width of the QBL depends on an effective Luttinger parameter K_{k}∼k^{θ} that vanishes with an exponent θ=z-1 related to the dynamical critical exponent z. The QBL encodes the existence of rare "superfluid" regions, controls the low-energy dynamics, and yields a (dissipative) conductivity vanishing as ω^{2} in the low-frequency limit. These results reveal the glassy properties (pinning, "shocks," or static avalanches) of the Bose-glass phase and can be understood within the "droplet" picture put forward for the description of glassy (classical) systems.

Classical critical dynamics in quadratically driven Kerr resonators

Driven-dissipative kerr lattices with two-photon driving are experimentally relevant systems known to exhibit a symmetry-breaking phase transition, which belongs to the universality class of the thermal Ising model for the parameter regime studied here. In this work, we perform finite-size scaling of this system as it is quenched to the transition and the dynamical critical exponent is found to be compatible with $z\approx2.18$ corresponding with metropolis dynamics in classical simulations. Furthermore, we show that the Liouvillian gap scales with the same exponent, similar to scaling of the Hamiltonian gap at quantum phase transitions in closed systems.

The negativity contour: a quasi-local measure of entanglement for mixed states

In this paper, we study the entanglement structure of mixed states in quantum many-body systems using the negativity contour, a local measure of entanglement that determines which real-space degrees of freedom in a subregion are contributing to the logarithmic negativity and with what magnitude. We construct an explicit contour function for Gaussian states using the fermionic partial-transpose. We generalize this contour function to generic many-body systems using a natural combination of derivatives of the logarithmic negativity. Though the latter negativity contour function is not strictly positive for all quantum systems, it is simple to compute and produces reasonable and interesting results. In particular, it rigorously satisfies the positivity condition for all holographic states and those obeying the quasi-particle picture. We apply this formalism to quantum field theories with a Fermi surface, contrasting the entanglement structure of Fermi liquids and holographic (hyperscale violating) non-Fermi liquids. The analysis of non-Fermi liquids show anomalous temperature dependence of the negativity depending on the dynamical critical exponent. We further compute the negativity contour following a quantum quench and discuss how this may clarify certain aspects of thermalization.

Hierarchy of energy scales in an O(3) symmetric antiferromagnetic quantum critical metal: A Monte Carlo study

We present numerically exact results from sign-problem free quantum Monte Carlo simulations for a spin-fermion model near an $O(3)$ symmetric antiferromagnetic (AFM) quantum critical point. We find a hierarchy of energy scales that emerges near the quantum critical point. At high energy scales, there is a broad regime characterized by Landau-damped order parameter dynamics with dynamical critical exponent $z=2$, while the fermionic excitations remain coherent. The quantum critical magnetic fluctuations are well described by Hertz-Millis theory, except for a $T^{-2}$ divergence of the static AFM susceptibility. This regime persists down to a lower energy scale, where the fermions become overdamped and concomitantly, a transition into a $d-$wave superconducting state occurs. These findings resemble earlier results for a spin-fermion model with easy-plane AFM fluctuations of an $O(2)$ SDW order parameter, despite noticeable differences in the perturbative structure of the two theories. In the $O(3)$ case, perturbative corrections to the spin-fermion vertex are expected to dominate at an additional energy scale, below which the $z=2$ behavior breaks down, leading to a novel $z=1$ fixed point with emergent local nesting at the hot spots [Schlief et al., PRX 7, 021010 (2017)]. Motivated by this prediction, we also consider a variant of the model where the hot spots are nearly locally nested. Within the available temperature range in our study ($T\ge E_F/200$), we find substantial deviations from the $z=2$ Hertz-Millis behavior, but no evidence for the predicted $z=1$ criticality.

Formation and critical dynamics of topological defects in Lifshitz holography

We examine the formation and critical dynamics of topological defects via Kibble-Zurek mechanism in a (2+1)-dimensional quantum critical point, which is conjectured to dual to a Lifshitz geometry. Quantized magnetic fluxoids are spontaneously generated and trapped in the cores of order parameter vortices, which is a feature of type-II superconductors. Time evolution of the average condensate is found to lag behind the instantaneous equilibrium value of the order parameter, a typical phenomenon in non- equilibrium dynamics. Scalings of vortex number density and the “freeze-out” time match the predictions from Kibble-Zurek mechanism. From these scalings, the dynamic and static critical exponents in the boundary field theory are found, at least at finite temperature, to be irrespective of the Lifshitz exponent in the bulk.

Anisotropic critical behavior of current-driven skyrmion dynamics in chiral magnets with disorder

The dynamic pinning effects are significant in manipulating skymions in chiral magnetic materials with quenched disorder. Through numerical simulations of the non-stationary current-driven dynamics of skyrmions with the Landau–Lifshitz–Gilbert equation, the critical current, static and dynamic critical exponents of the depenning phase transition are accurately determined for both adiabatic and non-adiabatic spin-transfer torques and with different strengths of disorder, based on the dynamic scaling behavior far from stationary. We find that the threshold current is insensitive to a small non-adiabatic coefficient of the spin-transfer torque, but dramatically reduced for a large one. The critical exponents indicate that the critical dynamic behavior is robust for different spin-transfer torques in the perpendicular component of the Hall motion, while exhibits a weak universality class in the direction of the driving current. The anisotropic behavior around the depinning phase transition provides a quantitative analysis of the drive-dependent skyrmion Hall effect in experiments. Further, the theoretical analysis using the Thiele’s approach is presented, and the critical current and the static exponents support the simulation results.

Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model.

We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.

Study of the phase diagram of the Kitaev-Hubbard chain

We present a detailed study of the phase diagram of the Kitaev-Hubbard chain, that is the Kitaev chain in the presence of a nearest-neighbour density-density interaction, using both analytical techniques as well as DMRG. In the case of a moderate attractive interaction, the model has the same phases as the non-interacting chain, a trivial and a topological phase. For repulsive interactions, the phase diagram is more interesting. Apart from the previously observed topological, incommensurate and charge density wave phases, we identify the `excited state charge density wave' phase. In this phase, the ground state resembles an excited state of an ordinary charge density phase, but is lower in energy due to the frustrated nature of the model. We find that the dynamical critical exponent takes the value $z\simeq 1.8$. Interestingly, this phase only appears for even system sizes, and is sensitive to the chemical potential on the edges of the chain. For the topological phase, we present an argument that excludes the presence of a strong zero mode for a large part of the topological phase. For the remaining region, we study the time dependence of the edge magnetization (using the bosonic incarnation of the model). These results further expand the region where a strong zero mode does not occur.

Universal Relaxation Dynamics of Sphere Packings below Jamming.

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Hot electron effects and electric field scaling near the metal-insulator transition in multilayer MoS2

The layered transition metal dichalcogenides have emerged as valuable platforms to study the challenging problem of metal-insulator transition in two dimensions. It was demonstrated that multilayer $\mathrm{Mo}{\mathrm{S}}_{2}$ exhibits clearly distinctive metallic and insulating behaviors in conductivity in response to both temperature and the electric field. Here, we report on the scaling analyses of conductivity for the electric field in addition to the temperature, which is performed with the consideration of electron-electron interactions for multilayer $\mathrm{Mo}{\mathrm{S}}_{2}$. Based on the analysis of hot electron effects in the electric field, we find that scaling for the electric field is relevant for the metallic phase in the high-field regime, enabling one to extract the dynamical critical exponent $z$ close to 1. This result supports that the metal-insulator transition in multilayer $\mathrm{Mo}{\mathrm{S}}_{2}$ is a true quantum critical phenomenon, in which strong interactions induce the transition.

Cluster Monte Carlo method for the 3D Ashkin–Teller model

This paper presents a cluster Monte Carlo method suitable for both, the first- and the second-order phase transitions in the 3D Ashkin–Teller (AT) model. A cluster algorithm is necessary to verify correctness of the results obtained so far using Metropolis algorithm exhibiting significant critical slowing down. Moreover, metastable states have recently been investigated in this model which can affect the results, especially obtained using the Metropolis algorithm. Our Wolff type algorithm is described and its dynamic critical behavior is demonstrated. Our computer experiments exploit the properties of Binder and Challa cumulants and additionally the one proposed by Lee and Kosterlitz, the last two adapted by us to give clear results for the AT model. The energy distribution histogram method is also independently applied for the first time for the 3D AT model using the Wolff type algorithm. For validation of the previous results and of our algorithm, it is demonstrated that the results of our computations along the line between Ising and Potts points, which are rescaled to their thermodynamic limits, are consistent with those obtained using the Metropolis algorithm. It is also shown that the presented cluster algorithm of the Wolff type significantly reduces the problem of critical slowing down for the 3D AT model, and the dynamic critical exponent reaches values close to zero. As the best strategy, it is suggested to use the cluster algorithm in the critical region and the Metropolis one beyond.

Renormalization group crossover in the critical dynamics of field theories with mode coupling terms

Motivated by the collective behavior of biological swarms, we study the critical dynamics of field theories with coupling between order parameter and conjugate momentum in the presence of dissipation. Under a fixed-network approximation, we perform a dynamical renormalization group calculation at one loop in the near-critical disordered region, and we show that the violation of momentum conservation generates a crossover between an unstable fixed point, characterized by a dynamic critical exponent z=d/2, and a 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. The crossover is regulated by a conservation length scale R_{0}, given by the ratio between the transport coefficient and the effective friction, which is larger as the dissipation is smaller: Beyond R_{0}, the stable fixed point dominates, while at shorter distances dynamics is ruled by the unstable fixed point and critical exponent, a behavior 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 slowdown 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 toward 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.

Dynamical Renormalization Group Approach to the Collective Behavior of Swarms.

We study the critical behavior of a model with nondissipative couplings aimed at describing the collective behavior of natural swarms, using the dynamical renormalization group under a fixed-network approximation. At one loop, we find a crossover between an unstable fixed point, characterized by a dynamical critical exponent z=d/2, and a stable fixed point with z=2, a result we confirm through numerical simulations. The crossover is regulated by a length scale given by the ratio between the transport coefficient and the effective friction, so that in finite-size biological systems with low dissipation, dynamics is ruled by the unstable fixed point. In three dimensions this mechanism gives z=3/2, a value significantly closer to the experimental window, 1.0≤z≤1.3, than the value z≈2 numerically found in fully dissipative models, either at or off equilibrium. This result indicates that nondissipative dynamical couplings are necessary to develop a theory of natural swarms fully consistent with experiments.

Computation of the dynamic critical exponent of the three-dimensional Heisenberg model.

Working in and out of equilibrium and using state-of-the-art techniques we have computed the dynamic critical exponent of the three-dimensional Heisenberg model. By computing the integrated autocorrelation time at equilibrium, for lattice sizes L≤64, we have obtained z=2.033(5). In the out-of-equilibrium regime we have run very large lattices (L≤250) obtaining z=2.04(2) from the growth of the correlation length. We compare our values with that previously computed at equilibrium with relatively small lattices (L≤24), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, η and ν, in the out-of-equilibrium regime.