No Arabic abstract
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 symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $propto 1/r^{alpha}$ decaying algebraically as a function of distance $r$ and a two-dimensional system with short-range nearest-neighbour interactions. In the former case we identify dynamical phase transitions for $alpha 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 $alpha 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.
In recent years, dynamical phase transitions and out-of-equilibrium criticality have been at the forefront of ultracold gases and condensed matter research. Whereas universality and scaling are established topics in equilibrium quantum many-body physics, out-of-equilibrium extensions of such concepts still leave much to be desired. Using exact diagonalization and the time-dependent variational principle in uniform martrix product states, we calculate the time evolution of the local order parameter and Loschmidt return rate in transverse-field Ising chains with antiferromagnetic power law-decaying interactions, and map out the corresponding rich dynamical phase diagram. textit{Anomalous} cusps in the return rate, which are ubiquitous at small quenches within the ordered phase in the case of ferromagnetic long-range interactions, are absent within the accessible timescales of our simulations. We attribute this to much weaker domain-wall binding in the antiferromagnetic case. For quenches across the quantum critical point, textit{regular} cusps appear in the return rate and connect to the local order parameter changing sign, indicating the concurrence of two major concepts of dynamical phase transitions. Our results consolidate conclusions of previous works that a necessary condition for the appearance of anomalous cusps in the return rate after quenches within the ordered phase is for topologically trivial local spin flips to be the energetically dominant excitations in the spectrum of the quench Hamiltonian. Our findings are readily accessible in modern trapped-ion setups, and we outline the associated experimental considerations.
In recent years, dynamical quantum phase transitions (DQPTs) have emerged as a useful theoretical concept to characterize nonequilibrium states of quantum matter. DQPTs are marked by singular behavior in an textit{effective free energy} $lambda(t)$, which, however, is a global measure, making its experimental or theoretical detection challenging in general. We introduce two local measures for the detection of DQPTs with the advantage of requiring fewer resources than the full effective free energy. The first, called the textit{real-local} effective free energy $lambda_M(t)$, is defined in real space and is therefore suitable for systems where locally resolved measurements are directly accessible such as in quantum-simulator experiments involving Rydberg atoms or trapped ions. We test $lambda_M(t)$ in Ising chains with nearest-neighbor and power-law interactions, and find that this measure allows extraction of the universal critical behavior of DQPTs. The second measure we introduce is the textit{momentum-local} effective free energy $lambda_k(t)$, which is targeted at systems where momentum-resolved quantities are more naturally accessible, such as through time-of-flight measurements in ultracold atoms. We benchmark $lambda_k(t)$ for the Kitaev chain, a paradigmatic system for topological quantum matter, in the presence of weak interactions. Our introduced local measures for effective free energies can further facilitate the detection of DQPTs in modern quantum-simulator experiments.
Using a Fermi-Bose mixture of ultra-cold atoms in an optical lattice, we construct a quantum simulator for a U(1) gauge theory coupled to fermionic matter. The construction is based on quantum links which realize continuous gauge symmetry with discrete quantum variables. At low energies, quantum link models with staggered fermions emerge from a Hubbard-type model which can be quantum simulated. This allows us to investigate string breaking as well as the real-time evolution after a quench in gauge theories, which are inaccessible to classical simulation methods.
The discrete truncated Wigner approximation (DTWA) is a powerful tool for analyzing dynamics of quantum spin systems. Since the DTWA includes the leading-order quantum corrections to a mean-field approximation, it is naturally expected that the DTWA becomes more accurate when the range of interactions of the system increases. However, quantitative corroboration of this expectation is still lacking mainly because it is generally difficult in a large system to evaluate a timescale on which the DTWA is quantitatively valid. In order to investigate how the validity timescale depends on the interaction range, we analyze dynamics of quantum spin models with a step function type interaction subjected to a sudden quench of a magnetic field by means of both DTWA and its extension including the second-order correction, which is derived from the Bogoliubov-Born-Green-Kirkwood-Yvon equation. We also develop a formulation for calculating the second-order Renyi entropy within the framework of the DTWA. By comparing the time evolution of the Renyi entropy computed by the DTWA with that by the extension including the correction, we find that both in the one- and two-dimensional systems the validity timescale increases algebraically with the range of the step function type interaction.
We investigate two solvable models for Bose-Einstein condensates and extract physical information by studying the structure of the solutions of their Bethe ansatz equations. A careful observation of these solutions for the ground state of both models, as we vary some parameters of the Hamiltonian, suggests a connection between the behavior of the roots of the Bethe ansatz equations and the physical behavior of the models. Then, by the use of standard techniques for approaching quantum phase transition - gap, entanglement and fidelity - we find that the change in the scenery in the roots of the Bethe ansatz equations is directly related to a quantum phase transition, thus providing an alternative method for its detection.