We consider the finite-time quench dynamics in the quantum transverse field Ising model which exhibits a second order phase transition from a paramagnetic to a ferromagnetic phase, as the transverse magnetic field is decreased. These dynamics have been thoroughly investigated in previous studies when the critical point is crossed during the quench; here, we quench the system from deep in the paramagnetic phase to just above the critical field so that the system remains in the gapped phase throughout the quench duration. On linearly quenching the infinitely large system, we find that the behavior of mean longitudinal defect density and mean transverse magnetization at the end of the quench falls into three distinct scaling regimes as the quench time is increased. For sufficiently small quench times, these observables remain roughly constant, but for larger quench times, a crossover occurs from the Kibble-Zurek scaling law to the quadratic quench rate law when the Kibble-Zurek time is of the order of relaxation time at the final quench field. These results are shown analytically using power series and uniform asymptotic expansions of the exact solution of the model, and also compared with an adiabatic perturbation theory in the third regime. We find that the above mentioned scaling regimes hold for quenches within the ferromagnetic phase also, and provide a general scaling argument for crossover from the Kibble-Zurek regime to an adiabatic regime for slow quenches within a gapped phase.
We consider the finite-time quench dynamics in the quantum transverse field Ising model which exhibits a second order phase transition from a paramagnetic to a ferromagnetic phase, as the transverse magnetic field is decreased. These dynamics have been thoroughly investigated in previous studies when the critical point is crossed during the quench; here, we quench the system from deep in the paramagnetic phase to just above the critical field so that the system remains in the gapped phase throughout the quench duration. On linearly quenching the infinitely large system, we find that the behavior of mean longitudinal defect density and mean transverse magnetization at the end of the quench falls into three distinct scaling regimes as the quench time is increased. For sufficiently small quench times, these observables remain roughly constant, but for larger quench times, a crossover occurs from the Kibble-Zurek scaling law to the quadratic quench rate law when the Kibble-Zurek time is of the order of relaxation time at the final quench field. These results are shown analytically using power series and uniform asymptotic expansions of the exact solution of the model, and also compared with an adiabatic perturbation theory in the third regime. We find that the above mentioned scaling regimes hold for quenches within the ferromagnetic phase also, and provide a general scaling argument for crossover from the Kibble-Zurek regime to an adiabatic regime for slow quenches within a gapped phase.
According to the Kibble-Zurek mechanism, there is a universal power-law relationship between the defect density and the quench rate during a slow linear quench through a critical point. It is generally accepted that a fast quench results in a deviation from the Kibble-Zurek scaling law and leads to the formation of a saturated plateau in the defect density. By adjusting the quench rate from slow to very fast limits, we observe the varying quench dynamics and identify a presaturated regime that lies between the saturated and Kibble-Zurek regimes. This significant result is first elucidated through the adiabatic-impulse approximation, then verified by a rigorous analysis on the transverse Ising chain as well. As we approach the turning point from the saturated to presaturated regimes, we notice a change in scaling laws and, with an increase in the initial transverse field, a shrinking of the saturated regime until it disappears. During another turning point from the Kibble-Zurek to presaturated regimes, we observe an attenuation of the dephasing effect and a change in the behavior of the kink-kink correlation function from a Gaussian decay to an exponential decay. Finally, the coherent many-body oscillation after quench exhibits different behaviors in the three regimes and shows a significant change of scaling behavior between the S and PS regimes.
When a quantum system exhibiting a second order phase transition is quenched across the critical point in large but finite time, the dynamics are not adiabatic in the critical region and the Kibble-Zurek (KZ) mechanism provides a framework to determine local observables such as the mean defect density. However, to find higher-point functions, one has to go beyond the KZ paradigm asshown in recent works on one-dimensional transverse field Ising model (TFIM) following a linear quench. It has been found that (i) besides the KZ scale, the quench dynamics depend on another length scale that arises due to the finite phase difference between the low energy modes, and (ii) contrary to the expectations based on the KZ mechanism, in general, the correlation functions do not decay exponentially with distance. Motivated by these results for the linear quench, we are interested in understanding if these properties are universal, and consider the 1D TFIM when the transverse field varies algebraically in the vicinity of the critical field. We focus on the equal-time,longitudinal kink-kink correlation function at the end of the quench from the paramagnetic to the ferromagnetic phase, and find that (i) the correlator depends only on the KZ length for superlinear quenches, otherwise an additional dephasing length is required to describe it, and (ii) the dephased correlator decays as a compressed exponential with an exponent that changes continuously with the quench exponent. Our results are obtained using an adiabatic perturbation theory, analytical arguments and exact numerical integration of the relevant equations.
Slow quenches of the magnetic field across the paramagnetic-ferromagnetic phase transition of spin systems produce heat. In systems with short-range interactions the heat exhibits universal power-law scaling as a function of the quench rate, known as Kibble-Zurek scaling. In this work we analyze slow quenches of the magnetic field in the Lipkin-Meshkov-Glick (LMG) model, which describes fully connected quantum spins. We analytically determine the quantum contribution to the residual heat as a function of the quench rate δ by means of a Holstein-Primakoff expansion about the mean-field value. Unlike in the case of short-range interactions, scaling laws in the LMG model are only found for a ramp starting or ending at the critical point. If instead the ramp is symmetric, as in the typical Kibble-Zurek scenario, then the number of excitations exhibits a crossover behavior as a function of δ and tends to a constant in the thermodynamic limit. Previous, and seemingly contradictory, theoretical studies are identified as specific limits of this dynamics. Our results can be tested on several experimental platforms, including quantum gases and trapped ions.
Quantum quenches in continuum field theory across critical points are known to display different scaling behaviours in different regimes of the quench rate. We extend these results to integrable lattice models such as the transverse field Ising model on a one-dimensional chain and the Kitaev model on a two-dimensional honeycomb lattice using a nonlinear quench protocol which allows for exact analytical solutions of the dynamics. Our quench protocol starts with a finite mass gap at early times and crosses a critical point or a critical region, and we study the behaviour of one point functions of the quenched operator at the critical point or in the critical region as a function of the quench rate. For quench rates slow compared to the initial mass gap, we find the expected Kibble-Zurek scaling. In contrast, for rates fast compared to the mass gap, but slow compared to the inverse lattice spacing, we find scaling behaviour similar to smooth fast continuum quenches. For quench rates of the same order of the lattice scale, the one point function saturates as a function of the rate, approaching the results of an abrupt quench. The presence of an extended critical surface in the Kitaev model leads to a variety of scaling exponents depending on the starting point and on the time where the operator is measured. We discuss the role of the amplitude of the quench in determining the extent of the slow (Kibble-Zurek) and fast quench regimes, and the onset of the saturation.
We study the crossing of the quantum phase transition in the transverse-field Ising model after modulating the magnetic field at an arbitrary rate, exploring the critical dynamics from the slow to the sudden-quench regime. We do so by analyzing the defect density, the complete kink-number distribution, and its cumulants upon completion of a linearized quench. Our analysis relies on the diagonalization of the model using the standard Jordan-Wigner and Fourier transformations, along with the exact solution of the time evolution in each mode in terms of parabolic cylinder functions. The free-fermion nature of the problem dictates that the kink-number distribution is associated with independent and distinguishable Bernoulli variables, each with a success probability pk. We employ a combination of convergent and asymptotic series expansions to characterize pk without restrictions on the driving rate. When the latter is approximated by the Landau-Zener formula, the kink density is described by the Kibble-Zurek mechanism, and higher-order cumulants follow a universal power-law behavior, as recently predicted theoretically and verified in the laboratory. By contrast, for moderate and sudden-driving protocols, the cumulants exhibit a nonuniversal behavior that is not monotonic on the driving rate and changes qualitatively with the cumulant order. The kink-number statistics remain sub-Poissonian for any driving rate, as revealed by an analysis of the cumulant ratios that vary nonmonotonically from the sudden- to the slow-driving regime. Due to the determination of pk for an arbitrary rate, our study provides a complete analytical understanding of kink-number statistics from the slow-driving regime to the sudden-quench limit.
Abstract Quantum quenches in continuum field theory across critical points are known to display different scaling behaviours in different regimes of the quench rate. We extend these results to integrable lattice models such as the transverse field Ising model on a one-dimensional chain and the Kitaev model on a two-dimensional honeycomb lattice using a nonlinear quench protocol which allows for exact analytical solutions of the dynamics. Our quench protocol starts with a finite mass gap at early times and crosses a critical point or a critical region, and we study the behaviour of one point functions of the quenched operator at the critical point or in the critical region as a function of the quench rate. For quench rates slow compared to the initial mass gap, we find the expected Kibble-Zurek scaling. In contrast, for rates fast compared to the mass gap, but slow compared to the inverse lattice spacing, we find scaling behaviour similar to smooth fast continuum quenches. For quench rates of the same order of the lattice scale, the one point function saturates as a function of the rate, approaching the results of an abrupt quench. The presence of an extended critical surface in the Kitaev model leads to a variety of scaling exponents depending on the starting point and on the time where the operator is measured. We discuss the role of the amplitude of the quench in determining the extent of the slow (Kibble-Zurek) and fast quench regimes, and the onset of the saturation.