Mapping finite-temperature dynamical phase diagrams of quantum many-body models is a necessary step towards establishing a framework of far-from-equilibrium quantum many-body universality. However, this is quite difficult due, in part, to the severe challenges in representing the volume-law entanglement that is generated under nonequilibrium dynamics at finite temperatures. Here, we address these challenges with an efficient equilibrium quantum Monte Carlo (QMC) framework for computing the finite-temperature dynamical phase diagram. Our method uses energy conservation and the self-thermalizing properties of ergodic quantum systems to determine observables at late times after a quantum quench. We use this technique to chart the dynamical phase diagram of the
2+1D quantum Ising model generated by quenches of the transverse field in initial thermal states. Our approach allows us to track the evolution of dynamical phases as a function of both the initial temperature and transverse field. Surprisingly, we identify quenches in the ordered phase that cool the system as well as an interval of initial temperatures where it is possible to quench from the paramagnetic (PM) to ferromagnetic (FM) phases. Our method gives access to dynamical properties without explicitly simulating unitary time evolution, and is immediately applicable to other lattice geometries and interacting many-body systems. Finally, we propose a quantum simulation experiment on state-of-the-art digital quantum hardware to directly probe the predicted dynamical phases and their real-time formation.
Quantum entanglement is fragile to thermal fluctuations, which raises the question whether finite temperature phase transitions support long-range entanglement similar to their zero temperature counterparts. Here we use quantum Monte Carlo simulations to study the third Renyi negativity, a generalization of entanglement negativity, as a proxy of mixed-state entanglement in the 2D transverse field Ising model across its finite temperature phase transition. We find that the area-law coefficient of the Renyi negativity is singular across the transition, while its subleading constant is zero within the statistical error. This indicates that the entanglement is short-range at the critical point despite a divergent correlation length. Renyi negativity in several exactly solvable models also shows qualitative similarities to that in the 2D transverse field Ising model.
Entanglement is a key physical resource in quantum information that enables basic quantum tasks, such as quantum teleportation and quantum key distribution, which are classically impossible. In this work, we investigate the entanglement theory of the 2D transverse field Ising model (TFIM) under the influence of a longitudinal magnetic field at zero temperature. For the fixed transverse field, the entanglement of the model initially increases to a maximum with the strength of the longitudinal magnetic field and then exponentially decreases to zero. Importantly, the maximal entanglement decreases with increasing transverse field, leading to a fully separable state in the limit of an infinite transverse magnetic field.