We investigate the quantum dynamics generated by quantum quenches (QQs) of the Hamiltonian parameters in many-body systems, focusing on protocols that cross first-order and continuous quantum transitions, both in finite-size systems and in the thermodynamic limit. As a paradigmatic example, we consider the quantum Ising chain in the presence of homogeneous transverse (
g) and longitudinal (
h) magnetic fields. This model exhibits a continuous quantum transition (CQT) at
g=gc and
h=0, and first-order quantum transitions (FOQTs) driven by
h along the line
h=0 (
g<gc). In the integrable limit
h=0, the system can be mapped onto a quadratic fermionic theory; however, any nonvanishing longitudinal field breaks integrability and the spectrum of the resulting Hamiltonian is generally expected to enter a chaotic regime. We analyze QQs in which the longitudinal field is suddenly changed from a negative value
hi<0 to a positive value
hf>0. We focus on values of
hf such that the spectrum of the post-QQ Hamiltonian
H(g,hf) lies in the chaotic regime, where thermalization may emerge at asymptotically long times. We study the out-of-equilibrium dynamics for different values of
g, finding qualitatively distinct behaviors for
g>gc (where the chain is in the disordered phase), for
g=gc (QQ across the CQT), and for
g<gc (QQ across the FOQT line).
We consider the non-equilibrium dynamics after a sudden quench of the magnetic field in the transverse field Ising chain starting from excited states of the pre-quench Hamiltonian. We prove that stationary values of local correlation functions can be described by the generalized Gibbs ensemble. Then we study the full time evolution of the transverse magnetization by means of stationary phase methods. The equal-time two-point longitudinal correlation function is analytically derived for a particular class of excited states for quenches within the ferromagnetic phase, and studied numerically in general. The full time dependence of the entanglement entropy of a block of spins is also obtained analytically for the same class of states and for arbitrary quenches.
Abstract We show that the symmetry-breaking gap of the quantum Ising model in the transverse field can be extracted from free evolution of the longitudinal magnetization taking place after a gradual quench of the magnetic field. We perform for this purpose numerical simulations of the Ising chains with either periodic or open boundaries. We also study the condition for adiabaticity of evolution of the longitudinal magnetization finding excellent agreement between our simulations and the prediction based on the Kibble–Zurek theory of non-equilibrium phase transitions. Our results should be relevant for ongoing cold atom and ion experiments targeting either equilibrium or dynamical aspects of quantum phase transitions. They could be also useful for benchmarking D-Wave machines.
We study the out-of-equilibrium dynamics of one-dimensional quantum Ising models in a transverse field
g, driven by a time-dependent longitudinal field
h across their {\em magnetic} first-order quantum transition at
h=0, for sufficiently small values of
∣g∣. We consider nearest-neighbor Ising chains of size
L with periodic boundary conditions. We focus on the out-of-equilibrium behavior arising from Kibble-Zurek protocols, in which
h is varied linearly in time with time scale
ts, i.e.,
h(t)=t/ts. The system starts from the ground state at
hi≡h(ti)<0, where the longitudinal magnetization
M is negative. Then it evolves unitarily up to positive values of
h(t), where
M(t) becomes eventually positive. We identify several scaling regimes characterized by a nontrivial interplay between the size
L and the time scale
ts, which can be observed when the system is close to one of the many avoided level crossings that occur for
h≥0. In the
L→∞ limit, all these crossings approach
h=0+, making the study of the thermodynamic limit, defined as the limit
L→∞ keeping
t and
ts constant, problematic. We study such limit numerically, by first determining the large-
L quantum evolution at fixed
ts, and then analyzing its behavior with increasing
ts. Our analysis shows that the system switches from the initial state with
M<0 to a positively magnetized state at
h=h∗(ts)>0, where
h∗(ts) decreases with increasing
ts, apparently as
h∗∼1/lnts. This suggests the existence of a scaling behavior in terms of the rescaled time
Ω=tlnts/ts. The numerical results also show that the system converges to a nontrivial stationary state in the large-
t limit, characterized by an energy significantly larger than that of the corresponding homogeneously magnetized ground state.
Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of the interesting cases of disordered hyperuniformity are provided by complex many-body systems such as liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. Such spin-chain systems are ideal models for exploring the effects of quantum mechanics on hyperuniformity. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform an extensive analysis of these systems using density matrix renormalization group simulations to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. Even for a spin chain of length 30, we see discontinuous changes in properties like the ``$\ensuremath{\tau}$ order metric'' of the ground state, the measure of hyperuniformity, and the second cumulant of the total magnetization along the
x-direction, all suggestive of first-order QPTs. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of a transverse magnetic field may lead remarkably to an increase in the order of the ground state as measured by the ``$\ensuremath{\tau}$ order metric,'' even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping onto a system of spinless fermions, showing that under the pairwise-interaction approximation and a mean-field treatment, there can be at most two quantum phase transitions. Based on these numerical and theoretical explorations, we conjecture that for spin chains with long-range pair interactions that have convergent cosine transforms, there can be a maximum of two quantum phase transitions at zero temperature.
Recent experimental achievements in controlling ultracold gases in optical lattices open a new perspective on quantum many-body physics. In these experimental setups, it is possible to study coherent time evolution of isolated quantum systems. These dynamics reveal new physics beyond the low-energy properties that are usually relevant in solid-state many-body systems. In this paper, we study the time evolution of antiferromagnetic order in the Heisenberg chain after a sudden change of the anisotropy parameter, using various numerical and analytical methods. As a generic result, we find that the order parameter, which can show oscillatory or non-oscillatory dynamics, decays exponentially except for the effectively non-interacting case of the XX limit. For weakly ordered initial states, we also find evidence for an algebraic correction to the exponential law. The study is based on numerical simulations using a numerical matrix product method for infinite system sizes (iMPS), for which we provide a detailed description and an error analysis. Additionally, we investigate in detail the exactly solvable XX limit. These results are compared to approximative analytical approaches including an effective description by the XZ model as well as by mean-field, Luttinger-liquid and sine-Gordon theories. The comparison reveals which aspects of non-equilibrium dynamics can, as in equilibrium, be described by low-energy theories and which are the novel phenomena specific to quantum quench dynamics. The relevance of the energetically high part of the spectrum is illustrated by means of a full numerical diagonalization of the Hamiltonian.