Fast Quenches Research Articles

Dynamical evolution of spinodal decomposition in holographic superfluids

We study the nonlinear dynamical evolution of spinodal decomposition in a first-order superfluid phase transition using a simple holographic model in the probe limit. We first confirm the linear stability analysis based on quasinormal modes and verify the existence of a critical length scale related to a gradient instability — negative speed of sound squared — of the superfluid sound mode, which is a consequence of a negative thermodynamic charge susceptibility. We present a comparison between our case and the standard Cahn-Hilliard equation for spinodal instability, in which a critical length scale can be also derived based on a diffusive instability. We then perform several numerical tests which include the nonlinear time evolution directly from an unstable state and fast quenches from a stable to an unstable state in the spinodal region. Our numerical results provide a real time description of spinodal decomposition and phase separation in one and two spatial dimensions. We reveal the existence of four different stages in the dynamical evolution, and characterize their main properties. Finally, we investigate the strength of dynamical heterogeneity using the spatial variance of the local chemical potential and we correlate the latter to other features of the dynamical evolution.

Meson excitation time as a probe of holographic critical point

We study the time evolution of expectation value of Wilson loop as a non-local observable in a strongly coupled field theory with a critical point at finite temperature and nonzero chemical potential, which is dual to an asymptotically AdS charged black hole via gauge/gravity duality. Due to inject of energy into the plasma, the temperature and a chemical potential increase to finite values and the plasma experiences an out-of-equilibrium process. By defining meson excitation time t_{ex} as a time at which the meson falls into the final excited state, we investigate the behavior of t_{ex} near the critical point as the system evolves towards the critical point. We observe that by increasing the interquark distance the dynamical critical exponent increases smoothly. Also, we obtain for slow quenches different values of the dynamical critical exponent, although for fast quenches our result for the dynamical critical exponent is in agreement with the one that is reported for studying the quasi-normal modes. Consequently, this indicates that in this model for fast quenches and small values of interquark distances the gauge invariant Wilson loop is a good non-local observable to probe the critical point.

Role of boundary conditions in the full counting statistics of topological defects after crossing a continuous phase transition

In a scenario of spontaneous symmetry breaking in finite time, topological defects are generated at a density that scale with the driving time according to the Kibble-Zurek mechanism (KZM). Signatures of universality beyond the KZM have recently been unveiled: The number distribution of topological defects has been shown to follow a binomial distribution, in which all cumulants inherit the universal power-law scaling with the quench rate, with cumulant rations being constant. In this work, we analyze the role of boundary conditions in the statistics of topological defects. In particular, we consider a lattice system with nearest-neighbor interactions subject to soft anti-periodic, open, and periodic boundary conditions implemented by an energy penalty term. We show that for fast and moderate quenches, the cumulants of the kink number distribution present a universal scaling with the quench rate that is independent of the boundary conditions except by an additive term, that becomes prominent in the limit of slow quenches, leading to the breaking of power-law behavior. We test our theoretical predictions with a one-dimensional scalar theory on a lattice.

History-dependent phase transition character

We consider history-dependent behavior in domain-type configurations in orientational order that are formed in configurations reached via continuous symmetry-breaking phase transitions. In equilibrium, these systems exhibit in absence of impurities a spatially homogeneous order. We focus on cases where domains are formed via (i) Kibble-Zurek mechanism in fast enough quenches or by (ii) Kibble mechanism in strongly supercooled phases. In both cases, domains could be arrested due to pinned topological defects that are formed at domain walls. In systems exhibiting polar or quadrupolar order, point and line defects (disclinations) dominate, respectively. In particular, the disclinations could form complex entangled structures and are more efficient in stabilizing domains. Domain patterns formed by fast quenches could be arrested by impurities imposing a strong enough random-field type disorder, as suggested by the Imry-Ma theorem. On the other hand, domains formed in supercooled systems could be also formed if large enough energy barriers arresting domains are established due to large enough systems’ stiffness. The resulting effective interactions in established domain-type patterns could be described by random matrices. The resulting eigenvectors reveal expected structural excitations formed in such structures. The most important role is commonly played by the random matrix largest eigenvector. Qualitatively different behavior is expected if this eigenvector exhibits a localized or extended character. In the former case, one expects a gradual, non-critical-type transition into a glass-type structure. However, in the latter case, a critical-like phase behavior could be observed.Graphical abstract

Yielding in an amorphous solid subject to constant stress at finite temperatures.

Understanding the nature of the yield transition is a long-standing problem in the physics of amorphous solids. Here we use molecular dynamics simulations to study the response of amorphous solids to constant stresses at finite temperatures. We compare amorphous solids that are prepared using fast and slow quenches and show that for thermal systems, the steady-state velocity exhibits a continuous transition from very slow creep to a finite strain rate as a function of the stress. This behavior is observed for both well-annealed and poorly annealed systems. However, the transient dynamics is different in the latter and involves overcoming an energy barrier. Due to the different simulation protocol, the strain rate as a function of stress and temperature follows a scaling relation that is different from the ones that are shown for systems where the strain is controlled. Collapsing the data using this scaling relation allows us to calculate critical exponents for the dynamics close to yield, including an exponent for thermal rounding. We also demonstrate that strain slips due to avalanche events above yield follow standard scaling relations and we extract critical exponents that are comparable to the ones obtained in previous studies that performed simulations of both molecular dynamics and elastoplastic models using strain-rate control.

Topological thermalization via vortex formation in ultrafast quenches.

We investigate the thermalization of a two-component scalar field across a second-order phase transition under extremely fast quenches. We find that vortices start developing once the thermal bath sets the control parameter to its final value in the nonsymmetric phase. Specifically, we find that vortices emerge as the fluctuating field relaxes and departs macroscopically from its symmetric configuration. The density of primordial vortices at the relaxation time is a decreasing function of the final temperature of the quench. Subsequently, vortices and antivortices annihilate at a rate that eventually determines the total thermalization time. This rate decreases if the theory contains a discrete anisotropy term, which otherwise leaves the primordial vortex density unaffected. Our results thus establish a link between the topological processes involved in the vortex dynamics and the delay in the thermalization of the system.

Observation of soft glassy behavior in a magnetic colloid exposed to an external magnetic field.

We provide the first experimental evidence for soft glassy behavior in a sterically stabilized magnetic colloid (ferrofluid) of relatively low volume fraction (φ = 0.037) when a uniform magnetic field is applied at a sufficiently high rate (fast quench). Fast magnetic-field quenches favor structural arrest of field-induced aggregates, owing to insufficient time to settle into lower energy states, thereby pushing the system to a frustrated metastable configuration like a repulsive glass. Brownian dynamics simulations are used to show that the polydisperse ferrofluid (as in experiments) forms thick ropes aligned along the field direction, while a monodisperse ferrofluid does not. The simulations show that there is practically no ordering of the thin, monodisperse chains, while the thick, polydisperse ropes show positional ordering with a typical center-center separation between the particles in different ropes of about 0.39μm. As a consequence of structural arrest, the ferrofluid exhibits aging with broken time-translational invariance, a hallmark of glassy dynamics. The superposition of strain and creep compliance curves obtained from rheological measurements at different waiting times in the effective time domain corroborates the soft glassy behavior when exposed to a magnetic field applied at a fast ramp rate.

Universal scaling in fast quenches near Lifshitz-like fixed points

We study critical dynamics through time evolution of quantum field theories driven to a Lifshitz-like fixed point, with z>1, under relevant deformations. The deformations we consider are fast smooth quantum quenches, namely when the quench scale δt−z is large compared to the deformation scale. We show that in holographic models the response of the system merely depends on the scaling dimension of the quenched operator as δλ⋅δtd−2Δ+z−1, where δλ is the deformation amplitude. This scaling behavior is enhanced logarithmically in certain cases. We also study free Lifshitz scalar theory deformed by mass operator and show that the universal scaling of the response completely matches with holographic analysis. We argue that this scaling behavior is universal for any relevant deformation around Lifshitz-like UV fixed points.

Quantum quench in non-relativistic fermionic field theory: harmonic traps and 2d string theory

We investigate a class of exactly solvable quantum quench protocols with a finite quench rate in systems of one dimensional non-relativistic fermions in external harmonic oscillator or inverted harmonic oscillator potentials, with time dependent masses and frequencies. These hamiltonians arise, respectively, in harmonic traps, and the c = 1 Matrix Model description of two dimensional string theory with time dependent string coupling. We show how the dynamics is determined by a single function of time which satisfies a generalized Ermakov-Pinney equation. The quench protocols we consider asymptote to constant masses and frequencies at early times, and cross or approach a gapless potential. In a right side up harmonic oscillator potential we determine the scaling behavior of the one point function and the entanglement entropy of a subregion by obtaining analytic approximations to the exact answers. The results are consistent with Kibble-Zurek scaling for slow quenches and with perturbation calculations for fast quenches. For cis-critical quench protocols the entanglement entropy oscillates at late times around its initial value. For end-critical protocols the entanglement entropy monotonically goes to zero inversely with time, reflecting the spread of fermions over the entire line. For the inverted harmonic oscillator potential, the dual collective field description is a scalar field in a time dependent metric and dilaton background.

Momentum-space entanglement after smooth quenches

We compute the total amount of entanglement produced between momentum modes at late times after a smooth mass quench in free bosonic and fermionic quantum field theories. The entanglement and Rényi entropies are obtained in closed form as a function of the parameters characterizing the quench protocol. For bosons, we show that the entanglement production is more significant for light modes and for fast quenches. In particular, infinitely slow or adiabatic quenches do not produce any entanglement. Depending on the quench profile, the decrease as a function of the quench rate delta t can be either monotonic or oscillating. In the fermionic case the situation is more subtle and there is a critical value for the quench amplitude above which this behavior is changed and the entropies become peaked at intermediate values of momentum and of the quench rate. We also show that the results agree with the predictions of a Generalized Gibbs Ensemble and obtain explicitly its parameters in terms of the quench data.

3D two-temperature magnetohydrodynamic modeling of fast thermal quenches due to injected impurities in tokamaks

An integrated model for the ionization, radiation, and advection of impurities in the extended-magnetohydrodynamic code M3D-C1 is described. This implementation makes use of the KPRAD model, which calculates bremsstrahlung radiation and impurity ionization, recombination, and radiation rates using a model in which the density of each charge state is advanced separately. The integrated model presented here allows the independent evolution of electron and ion temperatures, which is necessary to accurately model cases where the electron temperature drops more quickly than the electron–ion thermal equilibration time. This model is used to simulate the disruption of a model NSTX discharge caused by the introduction of argon impurities, using physically realistic resistivity. Despite well-mixed impurities, contraction of the current channel is found to lead to magnetohydrodynamic instabilities that result in stochastization of the magnetic field, a fast thermal quench, and localized parallel electric fields that can exceed the axisymmetric values by a factor of five for brief periods.

Dynamically probing strongly-coupled field theories with critical point

The dependence of the rescaled equilibration time on different parameters of the field theories with a holographic dual has been investigated in this paper. We consider field theories with nonzero chemical potential at finite temperature which are dual to asymptotically AdS charged black holes. We examine a dynamical probe scalar operator where its dynamics is due to a time-dependent source, quantum quench, or out-of-equilibrium initial condition in field theory with fixed or varying temperature and chemical potential. We observe that the behavior of the scalar operator equilibration time with respect to temperature or chemical potential can not be predicted merely by field theory parameters and depends on how fast the energy is injected into the system. It is shown that in field theories with critical point the rescaled equilibration time is shorter for thermodynamically stable systems. We also observe that the rescaled equilibration time as one approaches the critical point enhances and acquires an infinite slope though its value remains finite. We show that for fast quenches, even though the system is far from equilibrium, the dynamical critical exponent is the same as the one reported for quasi-normal modes in the same background. However for slow quenches the dynamical critical exponent picks up a different value.

On fast quenches and spinning correlators

We study global quantum quenches in a continuous field theoretic system with UV fixed point. Assuming that the characteristic inverse time scale of the smooth quench is much larger than all scales inherent to the system except for the UV-cutoff, we derive the universal scaling behavior of the two-point correlation functions associated with Dirac fields and spin-1 currents. We argue that in certain regimes our results can be recovered using the technique of operator product expansion.

Interferometric sensitivity and entanglement by scanning through quantum phase transitions in spinor Bose-Einstein condensates

Recent experiments have demonstrated the generation of entanglement by quasi-adiabatically driving through quantum phase transitions of a ferromagnetic spin-1 Bose-Einstein condensate in the presence of a tunable quadratic Zeeman shift. We analyze, in terms of the Fisher information, the interferometric value of the entanglement accessible by this approach. In addition to the Twin-Fock phase studied experimentally, we unveil a second regime, in the broken axisymmetry phase, which provides Heisenberg scaling of the quantum Fisher information and can be reached on shorter time scales. We identify optimal unitary transformations and an experimentally feasible optimal measurement prescription that maximize the interferometric sensitivity. We further ascertain that the Fisher information is robust with respect to non-adiabaticity and measurement noise. Finally, we show that the quasi-adiabatic entanglement preparation schemes admit higher sensitivities than dynamical methods based on fast quenches.

Universality of fast quenches from the conformal perturbation theory

We consider global quantum quenches, a protocol when a continuous field theoretic system in the ground state is driven by a homogeneous time-dependent external interaction. When the typical inverse time scale of the interaction is much larger than all relevant scales except for the UV-cutoff the system’s response exhibits universal scaling behavior. We provide both qualitative and quantitative explanations of this universality and argue that physics of the response during and shortly after the quench is governed by the conformal perturbation theory around the UV fixed point. We proceed to calculate the response of one and two-point correlation functions confirming and generalizing universal scalings found previously. Finally, we discuss late time behavior after the quench and argue that all local quantities will equilibrate to their thermal values specified by an excess energy acquired by the system during the quench.

Observing Dissipative Topological Defects with Coupled Lasers.

Topological defects have been observed and studied in a wide range of systems, such as cosmology, spin systems, cold atoms, and optics, as they are quenched across a phase transition into an ordered state. These defects limit the coherence of the system and its ability to approach a fully ordered state, so revealing their origin and control is becoming an increasingly important field of research. We observe dissipative topological defects in a one-dimensional ring of phased-locked lasers, and show how their formation is related to the Kibble-Zurek mechanism and is governed in a universal manner by two competing time scales. The ratio between these two time scales depends on the system parameters, and thus offers the possibility of enabling the system to dissipate to a fully ordered, defect-free state that can be exploited for solving hard computational problems in various fields.

Simulations of Cold Electroweak Baryogenesis: finding the optimal quench time

We revisit the numerical computation of the baryon asymmetry from Cold Electroweak Baryogenesis given the physical Higgs mass. We investigate the dependence of the asymmetry on the speed at which electroweak symmetry breaking takes place. The maximum asymmetry does not occur for arbitrarily fast quenches, but at quench times of about τq ≃ 16mH− 1, with no asymmetry created for quenches slower than τq > 30mH− 1. Curiously, we also find that the overall sign of the asymmetry depends on the quench time.

Creation and counting of defects in a temperature-quenched Bose-Einstein condensate

We study the spontaneous formation of defects in the order parameter of a trapped ultracold bosonic gas while crossing the critical temperature for Bose-Einstein Condensation (BEC) at different rates. The system has the shape of an elongated ellipsoid, whose transverse width can be varied to explore dimensionality effects. For slow enough temperature quenches we find a power-law scaling of the average defect number with the quench rate, as predicted by the Kibble-Zurek mechanism. A breakdown of such a scaling is found for fast quenches, leading to a saturation of the average defect number. We suggest an explanation for this saturation in terms of the mutual interactions among defects.

Irreversibility and the Arrow of Time in a Quenched Quantum System.

Irreversibility is one of the most intriguing concepts in physics. While microscopic physical laws are perfectly reversible, macroscopic average behavior has a preferred direction of time. According to the second law of thermodynamics, this arrow of time is associated with a positive mean entropy production. Using a nuclear magnetic resonance setup, we measure the nonequilibrium entropy produced in an isolated spin-1/2 system following fast quenches of an external magnetic field. We experimentally demonstrate that it is equal to the entropic distance, expressed by the Kullback-Leibler divergence, between a microscopic process and its time reversal. Our result addresses the concept of irreversibility from a microscopic quantum standpoint.

Smooth and fast versus instantaneous quenches in quantum field theory

We examine in detail the relationship between smooth fast quantum quenches, characterized by a time scale $\delta t$, and {\em instantaneous quenches}, within the framework of exactly solvable mass quenches in free scalar field theory. Our earlier studies \cite{dgm1,dgm2} highlighted that the two protocols remain distinct in the limit $\delta t \rightarrow 0$ because of the relation of the quench rate to the UV cut-off, i.e., $1/\delta t\ll\Lambda$ always holds in the fast smooth quenches while $1/\delta t\sim\Lambda$ for instantaneous quenches. Here we study UV finite quantities like correlators at finite spatial distances and the excess energy produced above the final ground state energy. We show that at late times and large distances (compared to the quench time scale) the smooth quench correlator approaches that for the instantaneous quench. At early times, we find that for small spatial separation and small $\delta t$, the correlator scales universally with $\delta t$, exactly as in the scaling of renormalized one point functions found in earlier work. At larger separation, the dependence on $\delta t$ drops out. The excess energy density is finite (for finite $m\delta t$) and scales in a universal fashion for all $d$. However, the scaling behaviour produces a divergent result in the limit $m\delta t \rightarrow 0$ for $d\ge4$, just as in an instantaneous quench, where it is UV divergent for $d \geq 4$. We argue that similar results hold for arbitrary interacting theories: the excess energy density produced is expected to diverge for scaling dimensions $\Delta > d/2$.