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 investigate the effects of dissipation on the quantum dynamics of many-body systems at quantum transitions, especially considering those of the first order. This issue is studied within the paradigmatic one-dimensional quantum Ising model. We analyze the out-of-equilibrium dynamics arising from quenches of the Hamiltonian parameters and dissipative mechanisms modeled by a Lindblad master equation, with either local or global spin operators acting as dissipative operators. Analogously to what happens at continuous quantum transitions, we observe a regime where the system develops a nontrivial dynamic scaling behavior, which is realized when the dissipation parameter
u (globally controlling the decay rate of the dissipation within the Lindblad framework) scales as the energy difference $\mathrm{\ensuremath{\Delta}}$ of the lowest levels of the Hamiltonian, i.e., $u\ensuremath{\sim}\mathrm{\ensuremath{\Delta}}$. However, unlike continuous quantum transitions where $\mathrm{\ensuremath{\Delta}}$ is power-law suppressed, at first-order quantum transitions $\mathrm{\ensuremath{\Delta}}$ is exponentially suppressed with increasing the system size (provided the boundary conditions do not favor any particular phase).
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.
By means of the discrete truncated Wigner approximation, we study dynamical phase transitions arising in the steady state of transverse-field Ising models after a quantum quench. Starting from a fully polarized ferromagnetic initial condition, these transitions separate a phase with nonvanishing magnetization along the ordering direction from a disordered symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $\ensuremath{\propto}1/{r}^{\ensuremath{\alpha}}$ decaying algebraically as a function of distance
r and a two-dimensional system with short-range nearest-neighbor interactions. In the former case, we identify dynamical phase transitions for $\ensuremath{\alpha}\ensuremath{\lesssim}2$ and we extract the critical exponents from a data collapse of the steady-state magnetization for up to 1200 lattice sites. We find identical exponents for $\ensuremath{\alpha}\ensuremath{\lesssim}0.5$, suggesting that the dynamical transitions in this regime fall into the same universality class as the nonergodic mean-field limit. The two-dimensional Ising model is believed to be thermalizing, which we also confirm using exact diagonalization for small system sizes. Thus, the dynamical transition is expected to correspond to the thermal phase transition, which is consistent with our data upon comparing to equilibrium quantum Monte Carlo simulations. We further test the accuracy of the discrete truncated Wigner approximation by comparing against numerically exact methods such as exact diagonalization, tensor network, as well as artificial neural network states and we find good quantitative agreement on the accessible time scales. Finally, our work provides an additional contribution to the understanding of the range and the limitations of qualitative and quantitative applicability of the discrete truncated Wigner approximation.
We address the quantum dynamics of a system composed of a qubit globally coupled to a many-body system characterized by short-range interactions. We employ a dynamic finite-size scaling framework to investigate the out-of-equilibrium dynamics arising from the sudden variation (turning on) of the interaction between the qubit and the many-body system, in particular when the latter is in proximity to a quantum first-order or continuous-phase transition. Although the approach is quite general, we consider
d-dimensional quantum Ising spin models in the presence of transverse and longitudinal fields as paradigmatic quantum many-body systems. To characterize the out-of-equilibrium dynamics, we focus on a number of quantum-information-oriented properties of the model. Namely, we concentrate on the decoherence features of the qubit, the energy interchange between the qubit and the many-body system during the out-of-equilibrium dynamics, and the work distribution associated with the quench. The scaling behaviors predicted by the dynamic finite-size scaling theory are verified through extensive numerical computations for the one-dimensional Ising model, which reveal a fast convergence to the expected asymptotic behavior with increasing system size.
We study the critical behavior of the nonequilibrium dynamics and of the steady states emerging from the competition between coherent and dissipative dynamics close to quantum phase transitions. The latter is induced by the coupling of the system with a Markovian bath, such that the evolution of the system's density matrix can be effectively described by a Lindblad master equation. We devise general scaling behaviors for the out-of-equilibrium evolution and the stationary states emerging in the large-time limit for generic initial conditions in terms of the parameters of the Hamiltonian providing the coherent driving and those associated with the dissipative interactions with the environment. Our framework is supported by numerical results for the dynamics of a one-dimensional lattice fermion gas undergoing a quantum Ising transition in the presence of dissipative mechanisms which include local pumping and decay of particles.