Search

The notorious fermion sign problem, arising from fermion statistics, constitutes one of the main obstacles of deciphering quantum many-body systems by numerical approach. The progress in overcoming sign problem will definitely lead to a great leap in various areas of modern physics. Here, by deviating from the conventional cognition that nonequilibrium studies should be more complicated than equilibrium cases, we propose an innovative framework based on nonequilibrium critical dynamics to preempt sign problem and investigate quantum critical point in fermionic model through numerically exact quantum Monte Carlo (QMC) simulation. By virtue of universal scaling theory of imaginary-time relaxation dynamics, we demonstrate that accurate critical point and critical exponents can be obtained in the short-time stage, in which the sign problem is not severe such that the QMC is accessible. After confirming the effectiveness of the method in two typical interacting fermionic models featuring Dirac quantum critical point (QCP), we for the first time reveal the quantum phase diagram in the Hubbard model hosting $\rm SU(3)$-symmetric Dirac fermions, and find that the QCP between Dirac semi-metal and a $\lambda_8$-antiferromagnetic phase belongs to a new universality class different from the previously known Gross-Neveu transitions., Comment: 7+13 pages

Various exotic phenomena emerge in non-equilibrium quantum many-body systems. The Mpemba effect, denoting the situation where a hot system freezes faster than the colder one, is a counterintuitive non-equilibrium phenomenon that has attracted enduring interest for more than half a century. In this Letter, we report a novel phenomenon of the Mpemba effect in the imaginary-time relaxation dynamics in quantum many-body systems, dubbed as imaginary-time Mpemba effect (ITME). Through numerically exact quantum Monte-Carlo (QMC) simulation, we unambiguously demonstrate that in different classes of interacting quantum models, the initial states with higher energy are relaxed faster than lower-energy initial states in the process of imaginary-time relaxation. The emergence of ITME is intimately associated with the low-energy excitations in quantum many-body systems. More crucially, since imaginary-time dynamics is broadly applied in numerical simulation on the quantum many-body ground states, the discovery of ITME potentially provides a new pathway to expedite the quantum many-body computation, particularly for QMC involving the sign problem., Comment: 4.5+8 pages, 4+6 figures

The nonequilibrium dynamics of quantum many-body systems have attracted growing attention due to various intriguing phenomena absent in equilibrium physics. One famous example is the quantum Mpemba effect, where the subsystem symmetry is restored faster under a symmetric quench from a more asymmetric initial state. The quantum Mpemba effect has been extensively studied in integrable and chaotic systems. In this Letter, we investigate symmetry restoration and quantum Mpemba effect in many-body localized systems with various initial states. We reveal that the symmetry can still be fully restored in many-body localization phases without approaching thermal equilibrium. Furthermore, we demonstrate that the presence of the quantum Mpemba effect is universal for any initial tilted product state, contrasting to the cases in the chaotic systems where the presence of the quantum Mpemba effect relies on the choice of initial states. We also provide a theoretical analysis of symmetry restoration and quantum Mpemba effects with the help of the effective model for many-body localization. This Letter not only sheds light on extending the quantum Mpemba effect to more non-equilibrium settings but also contributes to a deeper understanding of the many-body localization., Comment: 23 pages (including supplemental materials), 11 figures

Entanglement asymmetry, which serves as a diagnostic tool for symmetry breaking and a proxy for thermalization, has recently been proposed and studied in the context of symmetry restoration for quantum many-body systems undergoing a quench. In this Letter, we investigate symmetry restoration in various symmetric random quantum circuits, particularly focusing on the U(1) symmetry case. In contrast to non-symmetric random circuits where the U(1) symmetry of a small subsystem can always be restored at late times, we reveal that symmetry restoration can fail in U(1)-symmetric circuits for certain weak symmetry-broken initial states in finite-size systems. In the early-time dynamics, we observe an intriguing quantum Mpemba effect implying that symmetry is restored faster when the initial state is more asymmetric. Furthermore, we also investigate the entanglement asymmetry dynamics for SU(2) and $Z_{2}$ symmetric circuits and identify the presence and absence of the quantum Mpemba effect for the corresponding symmetries, respectively. A unified understanding of these results is provided through the lens of quantum thermalization with conserved charges., Comment: 4.5 pages, 5 figures, and Supplemental Material

While the zeros of complex partition functions, such as Lee-Yang zeros and Fisher zeros, have been pivotal in characterizing temperature-driven phase transitions, extending this concept to zero temperature remains an open question. In this work, we propose a solution to this issue by calculating the imaginary-temperature zeros (ITZs), which are defined as the roots of the imaginary-temperature partition function. We illustrate the analytical properties of ITZs in the transverse-field Ising chain, showing that the ITZs' distribution can distinguish between various phases and signify the critical exponents. Universal singular behaviors manifest in such quantities as the edge density of ITZs and the magnetization, with the scaling exponents remarkably differing from those in Lee-Yang theory. We further illuminate the consistency between ITZs and the zeros of the spectral form factor, which offers a practical path for the experimental detection of ITZs., Comment: 10 pages, 10 figures

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 on universal imaginary-time quantum critical dynamics., Comment: 5 pages, 4 figures and supplemental materials

Emergent symmetry is one of the characteristic phenomena in deconfined quantum critical point (DQCP). As its nonequilibrium generalization, the dual dynamic scaling was recently discovered in the nonequilibrium imaginary-time relaxation dynamics in the DQCP of the $J$-$Q_3$ model. In this work, we study the nonequilibrium imaginary-time relaxation dynamics in the $J$-$Q_2$ model, which also hosts a DQCP belonging to the same equilibrium universality class. We not only verify the universality of the dual dynamic scaling at the critical point, but also investigate the breakdown and the vestige of the dual dynamic scaling when the tuning parameter is away from the critical point. We also discuss its possible experimental realizations in devices of quantum computers., Comment: 18 pages, 11 figures

The conventional Kibble-Zurek mechanism (KZM) describes the driven critical dynamics in the Landau-Ginzburg-Wilson (LGW) spontaneous symmetry-breaking phase transitions. However, whether the KZM is still applicable in the deconfined quantum criticality, which is beyond the LGW paradigm, has not been explored. In this paper, we study the driven critical dynamics near a one-dimensional incarnation of deconfined quantum critical point between a ferromagnetic (FM) phase and a valance-bond-solid (VBS) phase. By investigating the dependence of the density of the topological defects on the driving rate, we verify the KZM in this Landau-forbidden critical point. Moreover, we find that both the FM and the VBS order parameters satisfy the finite-time scaling in the whole driven process. The effects of the emergent symmetry in the nonequilibrium dynamics are also studied., Comment: 5+3 pages, 4+2 figures

Fluctuations can change the phase transition properties drastically. An example is the fermion-induced quantum critical point (FIQCP), in which fluctuations of the massless Dirac fermions turn a putative Landau-de Gennes first-order phase transition (FOPT) with a cubic boson interaction into a continuous one. However, for the Landau-Devonshire theory, which characterizes another very large class of FOPTs, its fate under the coupling with extra fluctuations has not been explored. Here, we discover a new type of FIQCP, in which the Dirac fermion fluctuations round the boson Landau-Devonshire FOPT into a continuous phase transition. By using the functional renormalization group analyses, we determine the condition for the appearance of this FIQCP. Moreover, we point out that the present FIQCP can be a supersymmetric critical point. We finally show that the low-temperature phase diagram can provide distinct experimental evidences to detect this FIQCP., Comment: 9 pages, 7 figures, Accepted in PRB

We study the driven critical dynamics of the quantum link model, whose Hamiltonian describes the one-dimensional $U(1)$ lattice gauge theory. We find that combined topological defects emerge after the quench and they consist of both gauge field and matter field excitations. Furthermore, the ratio of gauge field and matter field excitation is $1/2$ due to the constraint of the Gauss' law. We show that the scaling of these combined topological defects satisfies the usual Kibble-Zurek mechanism. We verify that both the electric flux and the entanglement entropy satisfy the finite-time scaling theory in the whole driven process. Possible experimental realizations are discussed., Comment: 8.1 pages, 7 figures

We study the behavior of the out-of-time-ordered correlator (OTOC) in a non-Hermitian quantum Ising system. We show that the OTOC can diagnose not only the ground state exceptional point, which hosts the Yang-Lee edge singularity, but also the \textit{dynamical} exceptional point at the excited state. We find that the evolution of the OTOC in the parity-time symmetric phase can be divided into two stages: in the short-time stage, the OTOC oscillates periodically, and when the parameter is near the ground state exceptional point, this oscillation behavior can be described by both the scaling theory of the $(0+1)$D Yang-Lee edge singularity and the scaling theory of the $(1+1)$D quantum Ising model; while in the long-time stage the OTOC increases exponentially, controlled by the dynamical exceptional point. Possible experimental realizations are then discussed., Comment: 5+2 pages, 6+4 figures

We study the prethermal dynamics of the Gross-Neveu-Yukawa quantum field theory, suddenly quenched in the vicinity of a quantum critical point. We find that the universal prethermal dynamics is controlled by two fixed points depending on the size of the quench. Besides the usual equilibrium chiral Ising fixed point for a shallow quench, a dynamical chiral Ising fixed point is identified for a deep quench. Intriguingly, the latter is a non-thermal fixed point without any equilibrium counterpart due to the participation of gapless fermionic fields. We also find that in the scaling regime controlled by the equilibrium fixed point, the initial slip exponent is rendered negative if there are enough flavors of fermions, thus providing a unique signature of fermionic prethermal dynamics. We then explore the temporal crossover between the universal scaling regimes governed by the two universality classes. Possible experimental realizations are also discussed., Comment: 4.2 pages + supplemental materials, 2 figures

The conventional Kibble-Zurek scaling describes the scaling behavior in the driven dynamics across a single critical region. In this paper, we study the driven dynamics across an overlapping critical region, in which a critical region (Region-A) is overlaid by another critical region (Region-B). We develop a hybridized Kibble-Zurek scaling (HKZS) to characterize the scaling behavior in the driven process. According to the HKZS, the driven dynamics in the overlapping region can be described by the critical theories for both Region-A and Region-B simultaneously. This results in a constraint on the scaling function in the overlapping critical region. We take the quantum Ising chain in an imaginary longitudinal-field as an example. In this model, the critical region of the Yang-Lee edge singularity and the critical region of the ferromagnetic-paramagnetic phase transition point overlap with each other. We numerically confirm the HKZS by simulating the driven dynamics in this overlapping region. The HKZSs in other models are also discussed., Comment: 6.1 pages, 6 figures, added more references

We study the driven dynamics across the critical points of the Yang-Lee edge singularities (YLESes) in a finite-size quantum Ising chain with an imaginary symmetry-breaking field. In contrast to the conventional classical or quantum phase transitions, these phase transitions are induced by tuning the strength of the dissipation in a non-Hermitian system and can occur even at finite size. For conventional phase transitions, universal behaviors in driven dynamics across critical points are usually described by the Kibble-Zurek mechanism, which states that the scaling in dynamics is dictated by the critical exponents associated with one critical point and topological defects will emerge after the quench. While the mechanism leading to topological defects breaks down in the YLES, we find that for small lattice size, the driven dynamics can still be described by the Kibble-Zurek scaling with the exponents determined by the $(0+1)$-dimensional YLES. For medium finite size, however, the driven dynamics can be described by the Kibble-Zurek scaling with two sets of critical exponents determined by both the $(0+1)$-dimensional and the $(1+1)$-dimensional YLESes., Comment: 4.5+3 pages, 4+8 figures; The Experimental feasibility and the discussion on the KZM in Parity-time symmetry-breaking phase transition have been added

We study the driven critical dynamics with an equilibrium initial state near a quantum critical point. In contrast to the original Kibble-Zurek mechanism, which describes the driven dynamics starting from an adiabatic stage that is far from the critical point, the initial adiabacity is broken in this scenario. As a result, the scaling behavior cannot be described by the original Kibble-Zurek scaling. In this work we propose a scaling theory, which includes the initial parameters as additional scaling variables, to characterize the scaling behavior. In particular, this scaling theory can be used to describe the driven scaling behavior starting from a finite-temperature equilibrium state near a quantum critical point. We numerically confirm the scaling theory by simulating the real-time dynamics of the one-dimensional quantum Ising model at both zero and finite temperatures., Comment: 6 pages, 5 figures

We study the scaling behavior of the relaxation dynamics to thermal equilibrium when a quantum system is near the quantum critical point. In particular, we investigate systems whose relaxation dynamics is described by a Lindblad master equation. We find that the universal scaling behavior not only exhibits in the equilibrium stage at the long-time limit, but also manifests itself in the non-equilibrium relaxation process. While the critical behavior is dictated by the low-lying energy levels of the Hamiltonian, the dissipative part in the Lindblad equation also plays important roles in two aspects: First, the dissipative part makes the high energy levels decay fast after which the universal behavior controlled by the low-lying modes emerges. Second, the dissipation rate gives rise to a time scale that affects the scaling behavior. We confirm our theory by solving the Lindblad equation for the one-dimensional transverse-field Ising model., Comment: 5 pages, 4 figures

Starting from the Liouville-von Neumann equation, under a weak coupling limit we derive the Lindblad master equation for the one-dimensional quantum Ising model in a Markov approximation and a rotating wave approximation. We also prove that the steady solution of the Lindblad equation is the canonical distribution independent of the dissipation rate.