Quantum computers promise a highly efficient approach to investigate quantum phase transitions, which describe abrupt changes between different ground states of many-body systems. At quantum critical points, the divergent correlation length and entanglement entropy render the ground state preparation difficult. In this work, we explore the imaginary-time evolution for probing the universal critical behavior as the universal information of the ground state can be extracted in the early time relaxation process. We propose a systematic and scalable scheme to probe the universal behaviors via imaginary-time critical dynamics on quantum computers and demonstrate the validness of our approach by both numerical simulation and quantum hardware experiments. With the full form of the universal scaling function in terms of imaginary time, system size, and circuit depth, we successfully probe the universality by scaling analysis of the critical dynamics at an early time and with shallower quantum circuit depth. Equipped with quantum error mitigation, we also confirm the expected scaling behavior from experimental results on a superconducting quantum processor which stands as the first experimental demonstration of universal imaginary-time quantum critical dynamics.
We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at
(h/J)c=1.00003(8) in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length
L=32−−256 and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by
3.04451(7) in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size
L=24--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.
Characterizing universal entanglement features in higher-dimensional quantum matter is a central goal of quantum information science and condensed matter physics. While the subleading corner terms in two-dimensional quantum systems encapsulate essential universal information of the underlying conformal field theory, our understanding of these features remains remarkably limited compared to their one-dimensional counterparts. We address this challenge by investigating the entanglement dynamics of fermionic systems along the imaginary-time evolution. We uncover a pioneering non-equilibrium scaling law where the corner entanglement entropy grows linearly with the logarithm of imaginary time, dictated solely by the universality class of the quantum critical point. Through unbiased Quantum Monte Carlo simulations, we verify this scaling in the interacting Gross-Neveu-Yukawa model, demonstrating that universal data can be accurately recovered from the early stages of relaxation. Our findings significantly circumvent the computational bottlenecks inherent in reaching full equilibrium convergence. This work establishes a direct link between the fundamental theory of non-equilibrium critical phenomena and the high-precision determination of universal entanglement properties on both classical and quantum platforms, paving the way for probing the rich entanglement structure of quantum critical systems.
While imaginary-time evolution has long served as a standard paradigm for ground-state preparation in numerical simulations and quantum devices, its intrinsic dynamical properties has been largely overlooked. Here, we investigate the short-imaginary-time critical dynamics in quantum systems with boundaries. A universal scaling theory is developed and verified in the two-dimensional quantum Ising model, uncovering rich dynamic critical behaviors dictated by boundary universality classes. For ordered initial states, the boundary order parameter
Ms decays with imaginary time
τ as
Ms∝τ−β1/νz, where
β1 denotes the boundary order parameter exponent, and
ν and
z correspond to the correlation length exponent and the dynamic exponent, respectively. For disordered initial states, the autocorrelation of the boundary order parameter is governed by a novel critical exponent
θ1, which is closely related to the critical initial slip behavior of
Ms characterized by the corresponding exponent
θ1′. In contrast to its positive bulk counterpart, the boundary initial-slip exponent
θ1′ is negative for the ordinary transition while remaining positive for the special transition. Although the static universality classes of
d-dimensional quantum phase transitions generally coincide with those of
(d+1)-dimensional classical phase transitions, we show that
θ1 does not follow this conventional quantum-classical mapping. We further discuss the implications of our results for more exotic forms of boundary criticality. Our findings provide new physical insights into boundary critical dynamics and offer a novel route for probing exotic boundary critical behaviors in quantum many-body systems.
We investigate the imaginary-time relaxation critical dynamics of the two-dimensional transverse-field Ising model using infinite projected entangled pair states (iPEPS) with the full-update strategy. Simulating directly in the thermodynamic limit, we explore the relaxation process near the critical point with two types of initial states: a fully polarized state and a product state with a small magnetization. For the fully polarized state, the magnetization shows a power law scaling
M∝τ−β/(νz) in the imaginary-time evolution, from which both the critical point and critical exponent can be determined with high accuracy. For the nearly paramagnetic state, the relaxation process exhibits a behavior of
M∝τθ with
θ=0.1958 being the critical initial-slip exponent, which is in good agreement with that obtained from the dynamic scaling of the self-correlation in quantum Monte Carlo method. These universal features emerge well before the system converges to the ground state, demonstrating the efficiency of imaginary-time evolution for probing quantum criticality. Our results demonstrate that iPEPS can serve as a robust and scalable method for studying dynamical critical phenomena in two-dimensional quantum many-body systems.
Computing the ground-state properties of quantum many-body systems is a promising application of near-term quantum hardware with a potential impact in many fields. The conventional algorithm quantum phase estimation uses deep circuits and requires fault-tolerant technologies. Many quantum simulation algorithms developed recently work in an inexact and variational manner to exploit shallow circuits. In this work, we combine quantum Monte Carlo with quantum computing and propose an algorithm for simulating the imaginary-time evolution and solving the ground-state problem. By sampling the real-time evolution operator with a random evolution time according to a modified Cauchy-Lorentz distribution, we can compute the expected value of an observable in imaginary-time evolution. Our algorithm approaches the exact solution given a circuit depth increasing polylogarithmically with the desired accuracy. Compared with quantum phase estimation, the Trotter step number, i.e. the circuit depth, can be thousands of times smaller to achieve the same accuracy in the ground-state energy. We verify the resilience to Trotterisation errors caused by the finite circuit depth in the numerical simulation of various models. The results show that Monte Carlo quantum simulation is promising even without a fully fault-tolerant quantum computer.
The scaling of the entanglement entropy at a quantum critical point allows us to extract universal properties of the state, e.g., the central charge of a conformal field theory. With the rapid improvement of noisy intermediate-scale quantum (NISQ) devices, these quantum computers present themselves as a powerful tool to study critical many-body systems. We use finite-depth quantum circuits suitable for NISQ devices as a variational ansatz to represent ground states of critical, infinite systems. We find universal finite-depth scaling relations for these circuits and verify them numerically at two different critical points, i.e., the critical Ising model with an additional symmetry-preserving term and the critical XXZ model.