Dynamics of quenching temperature is studied in pure and random Ising chains. Using the Kibble-Zurek argument, we obtain for the pure Ising model that the density of kinks after quenching decays as 1/\sqrtτ with the quench rate of temperature 1/τfor large τ. For the random Ising model, we show that decay rates of the density of kinks and the residual energy are 1/\lnτand 1/(\lnτ)^2 for large τrespectively. Analytic results for the random Ising model are confirmed by the Monte-Carlo simulation. Our results reveal a clear difference between classical and quantum quenches in the random Ising chain.
Dynamics of quenching temperature is studied in pure and random Ising chains. Using the Kibble-Zurek argument, we obtain for the pure Ising model that the density of kinks after quenching decays as 1/\sqrtτ with the quench rate of temperature 1/τfor large τ. For the random Ising model, we show that decay rates of the density of kinks and the residual energy are 1/\lnτand 1/(\lnτ)^2 for large τrespectively. Analytic results for the random Ising model are confirmed by the Monte-Carlo simulation. Our results reveal a clear difference between classical and quantum quenches in the random Ising chain.
We investigate the time evolution of the density of kinks in the spin-1/2 quantum Ising spin chain after a sudden quench in the transverse field strength, and find that it relaxes to a value which depends on the initial and the final values of the transverse field, with an oscillating power-law decay. We provide analytical estimates of the long-time behavior and of the asymptotic value reached after complete relaxation, and discuss the role of quantum criticality in the quench dynamics. We show that, for a dynamics at the critical point, the residual density of kinks after the quench can be described by equilibrium statistical mechanics at a finite temperature dictated by the energy of the state after the quench. On the other hand, outside of criticality it does not exhibit thermalization.
We consider a linear quench from the paramagnetic to ferromagnetic phase in the quantum Ising chain interacting with a static spin environment. Both decoherence from the environment and nonadiabaticity of the evolution near a critical point excite the system from the final ferromagnetic ground state. For weak decoherence and relatively fast quenches the excitation energy, proportional to the number of kinks in the final state, decays like an inverse square root of a quench time, but slow transitions or strong decoherence makes it decay in a much slower logarithmic way. We also find that fidelity between the final ferromagnetic ground state and a final state after a quench decays exponentially with a size of a chain, with a decay rate proportional to average density of excited kinks and a proportionality factor evolving from 1.3 for weak decoherence and fast quenches to approximately 1 for slow transitions or strong decoherence. Simultaneously, correlations between kinks randomly distributed along the chain evolve from a near-crystalline antibunching to a Poissonian distribution of kinks in a number of isolated Anderson localization centers randomly scattered along the chain.
We present here our study of the adiabatic quantum dynamics of a random Ising chain across its quantum critical point. The model investigated is an Ising chain in a transverse field with disorder present both in the exchange coupling and in the transverse field. The transverse field term is proportional to a function $\ensuremath{\Gamma}(t)$ which, as in the Kibble-Zurek mechanism, is linearly reduced to zero in time with a rate ${\ensuremath{\tau}}^{\ensuremath{-}1}$, $\ensuremath{\Gamma}(t)=\ensuremath{-}t∕\ensuremath{\tau}$, starting at $t=\ensuremath{-}\ensuremath{\infty}$ from the quantum disordered phase $(\ensuremath{\Gamma}=\ensuremath{\infty})$ and ending at t=0 in the classical ferromagnetic phase $(\ensuremath{\Gamma}=0)$. We first analyze the distribution of the gaps, occurring at the critical point ${\ensuremath{\Gamma}}_{c}=1$, which are relevant for breaking the adiabaticity of the dynamics. We then present extensive numerical simulations for the residual energy Eres and density of defects ${\ensuremath{\rho}}_{k}$ at the end of the annealing, as a function of the annealing inverse rate $\ensuremath{\tau}$. Both the average ${E}_{\mathrm{res}}(\ensuremath{\tau})$ and ${\ensuremath{\rho}}_{k}(\ensuremath{\tau})$ are found to behave logarithmically for large $\ensuremath{\tau}$, but with different exponents, $[{E}_{\mathrm{res}}(\ensuremath{\tau})∕L{]}_{\mathrm{av}}\ensuremath{\sim}1∕{\mathrm{ln}}^{\ensuremath{\zeta}}(\ensuremath{\tau})$ with $\ensuremath{\zeta}\ensuremath{\approx}3.4$, and $[{\ensuremath{\rho}}_{k}(\ensuremath{\tau}){]}_{\mathrm{av}}\ensuremath{\sim}1∕{\mathrm{ln}}^{2}(\ensuremath{\tau})$. We propose a mechanism for $1∕{\mathrm{ln}}^{2}\ensuremath{\tau}$ behavior of $[{\ensuremath{\rho}}_{k}{]}_{\mathrm{av}}$ based on the Landau-Zener tunneling theory and on a Fisher's-type real-space renormalization group analysis of the relevant gaps. The model proposed shows therefore a paradigmatic example of how an adiabatic quantum computation can become very slow when disorder is at play, even in the absence of any source of frustration.
We consider the annealing dynamics of a one-dimensional Ising ferromagnet induced by a temperature quench in finite time. In the limit of slow cooling, the asymptotic two-point correlator is analytically found under Glauber dynamics, and the distribution of the number of kinks in the final state is shown to be consistent with a Poissonian distribution. The mean kink number, the variance, and the third centered moment take the same value and obey a universal power-law scaling with the quench time in which the temperature is varied. The universal power-law scaling of cumulants is corroborated by numerical simulations based on Glauber dynamics for moderate cooling times away from the asymptotic limit, when the kink-number distribution takes a binomial form. We analyze the relation of these results to physics beyond the Kibble-Zurek mechanism for critical dynamics, using the kink-number distribution to assess adiabaticity and its breakdown. We consider linear, nonlinear, and exponential cooling schedules, among which the last provides the most efficient shortcuts to cooling in a given quench time. The nonthermal behavior of the final state is established by considering the trace norm distance to a canonical Gibbs state.
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.
Quantum Ising model in one dimension is an exactly solvable example of a quantum phase transition. We investigate its behavior during a quench caused by a gradual turning off of the transverse bias field. The system is then driven at a fixed rate characterized by the quench time ${\ensuremath{\tau}}_{Q}$ across the critical point from a paramagnetic to ferromagnetic phase. In agreement with Kibble-Zurek mechanism (which recognizes that evolution is approximately adiabatic far away, but becomes approximately impulse sufficiently near the critical point), quantum state of the system after the transition exhibits a characteristic correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ proportional to the square root of the quench time ${\ensuremath{\tau}}_{Q}$: $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}=\sqrt{{\ensuremath{\tau}}_{Q}}$. The inverse of this correlation length is known to determine average density of defects (e.g., kinks) after the transition. In this paper, we show that this same $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ controls the entropy of entanglement, e.g., entropy of a block of L spins that are entangled with the rest of the system after the transition from the paramagnetic ground state induced by the quench. For large L, this entropy saturates at $\frac{1}{6}\phantom{\rule{0.2em}{0ex}}{\mathrm{log}}_{2}\phantom{\rule{0.2em}{0ex}}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, as might have been expected from the Kibble-Zurek mechanism. Close to the critical point, the entropy saturates when the block size $L\ensuremath{\approx}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but---in the subsequent evolution in the ferromagnetic phase---a somewhat larger length scale $l=\sqrt{{\ensuremath{\tau}}_{Q}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}\phantom{\rule{0.2em}{0ex}}{\ensuremath{\tau}}_{Q}$ develops as a result of a dephasing process that can be regarded as a quantum analog of phase ordering, and the entropy saturates when $L\ensuremath{\approx}l$. We also study the spin-spin correlation using both analytic methods and real time simulations with the Vidal algorithm. We find that at an instant when quench is crossing the critical point, ferromagnetic correlations decay exponentially with the dynamical correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but (as for entropy of entanglement) in the following evolution length scale l gradually develops. The correlation function becomes oscillatory at distances less than this scale. However, both the wavelength and the correlation length of these oscillations are still determined by $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$. We also derive probability distribution for the number of kinks in a finite spin chain after the transition.