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.
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.
The Kibble-Zurek mechanism (KZM) and finite-time scaling (FTS) provide a foundational framework for driven critical dynamics, yet their predictive power has been largely confined to local observables. Here, we establish a universal finite-time scaling theory for the nonequilibrium dynamics of quantum entanglement. Using unbiased quantum Monte Carlo simulations, we investigate the corner entanglement entropy of (2+1)-dimensional interacting Dirac fermions driven from ordered phases toward a quantum critical point. We find that the corner entanglement accurately obeys a universal driven scaling governed by the driving rate and system size, persisting whether the initial ordered state is fully gapped or hosts gapless Goldstone modes. Crucially, this dynamical entanglement exhibits a logarithmic dependence on the driving rate, from which the universal corner coefficient of the underlying conformal field theory can be robustly extracted far from equilibrium. These results generalize the KZM from local observables to the intrinsic nonlocal quantum information measures, offering a practical blueprint for characterizing quantum criticality and entanglement on programmable quantum simulators.
The entanglement entropy is a unique probe to reveal universal features of strongly interacting many-body systems. In two or more dimensions these features are subtle, and detecting them numerically requires extreme precision, a notoriously difficult task. This is especially challenging in models of interacting fermions, where many such universal features have yet to be observed. In this Letter we tackle this challenge by introducing a new method to compute the Rényi entanglement entropy in auxiliary-field quantum Monte Carlo simulations, where we treat the entangling region itself as a stochastic variable. We demonstrate the efficiency of this method by extracting, for the first time, universal subleading logarithmic terms in a two-dimensional model of interacting fermions, focusing on the half-filled honeycomb Hubbard model at T=0. We detect the universal corner contribution due to gapless fermions throughout the Dirac semi-metal phase and at the Gross-Neveu-Yukawa critical point, where the latter shows a pronounced enhancement depending on the type of entangling cut. Finally, we observe the universal Goldstone mode contribution in the antiferromagnetic Mott insulating phase.
Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.
Entanglement entropy has become an important tool for identifying phases of matter and phase transitions in quantum systems using quantum Monte Carlo methods. While such approaches have achieved great success in systems composed of interacting spins/bosons, their application to interacting fermionic systems remains limited due to substantially higher computational cost. Here we show a more efficient Monte Carlo algorithm for fermionic systems that enables systematic exploration of entanglement entropy across a wide range of physical parameters. Based on the incremental technique along physical parameters, the method significantly reduces computational effort while maintaining accuracy. Using this approach, we study a two-dimensional square lattice Hubbard model, revealing the phase diagram that includes the Fermi surface and Goldstone modes. We further apply the method to the Gross-Neveu criticality and find that the scaling of the entanglement entropy follows a universal form quantified by the critical exponent ν. Further investigation reveals that the leading coefficient of the entanglement entropy decreases monotonically, rather than developing a local maximum as O(N) transition point. Entanglement entropy is a powerful tool for identifying quantum phases and phase transitions but calculating it in interacting fermionic systems remains challenging. Here, the authors introduce an efficient Monte Carlo method that reduces computational cost and reveals universal scaling across Gross-Neveu fermionic phase transitions.
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.
Entanglement scaling provides a powerful probe of universal properties at quantum critical points. In two dimensions, contributions originating from a corner-shaped bipartition exhibit a universal scaling, which is determined by the underlying conformal field theory. We develop a method to extract this corner entanglement entropy from projected entangled pair states directly in the thermodynamic limit. When applied to models at a quantum critical point, we show that the corner contribution exhibits scaling with the effective correlation length, in agreement with the hypothesis of finite-entanglement scaling. Our results for the corner coefficients are consistent with other methods, demonstrating the efficiency of our method for diagnosing strongly-correlated quantum critical points in two dimensions.