Search

We investigate operator dynamics and entanglement growth in dual-unitary circuits, a class of locally scrambled quantum systems that enables efficient simulation beyond the exponential complexity of the Hilbert space. By mapping the operator evolution to a classical Markov process,we perform Monte Carlo simulations to access the time evolution of local operator density and entanglement with polynomial computational cost. Our results reveal that the operator density converges exponentially to a steady-state value, with analytical bounds that match our simulations. Additionally, we observe a volume-law scaling of operator entanglement across different subregions,and identify a critical transition from maximal to sub-maximal entanglement growth, governed by the circuit's gate parameter. This transition, confirmed by both mean-field theory and Monte Carlo simulations, provides new insights into operator entanglement dynamics in quantum many-body systems. Our work offers a scalable computational framework for studying long-time operator evolution and entanglement, paving the way for deeper exploration of quantum information dynamics., Comment: 12 pages,12 figures

We propose the convex-roof extension of quantum conditional mutual information ("co(QCMI)") as a diagnostic of long-range entanglement in a mixed state. We focus primarily on topological states subjected to local decoherence, and employ the Levin-Wen scheme to define co(QCMI), so that for a pure state, co(QCMI) equals topological entanglement entropy (TEE). By construction, co(QCMI) is zero if and only if a mixed state can be decomposed as a convex sum of pure states with zero TEE. We show that co(QCMI) is non-increasing with increasing decoherence when Kraus operators are proportional to the product of onsite unitaries. This implies that unlike a pure state transition between a topologically trivial and a non-trivial phase, the long-range entanglement at a decoherence-induced topological phase transition as quantified by co(QCMI) is less than or equal to that in the proximate topological phase. For the 2d toric code decohered by onsite bit/phase-flip noise, we show that co(QCMI) is non-zero below the error-recovery threshold and zero above it. Relatedly, the decohered state cannot be written as a convex sum of short-range entangled pure states below the threshold. We conjecture and provide evidence that in this example, co(QCMI) equals TEE of a recently introduced pure state. In particular, we develop a tensor-assisted Monte Carlo (TMC) computation method to efficiently evaluate the R\'enyi TEE for the aforementioned pure state and provide non-trivial consistency checks for our conjecture. We use TMC to also calculate the universal scaling dimension of the anyon-condensation order parameter at this transition., Comment: 17 pages main text, 3 pages of appendices, 7 figures

The subleading corner logarithmic corrections in entanglement entropy (EE) are crucial for revealing universal characteristics of the quantum critical points (QCPs), but they are challenging to detect. Motivated by recent developments in the stable computation of EE in (2+1)D quantum many-body systems, we have developed a new method for directly measuring the corner contribution in EE with less computational cost. The cornerstone of our approach is to measure the subtracted corner entanglement entropy (SCEE) defined as the difference between the EEs of subregions with the same boundary length for smooth and cornered boundaries during the sign-problem free quantum Monte Carlo simulation. Our improved method inherently eliminates not only the area law term of EE but also the subleading log-corrections arising from Goldstone modes, leaving the universal corner contribution as the leading term of SCEE with greatly improved data quality. Utilizing this advanced approach, we calculate the SCEE of the bilayer Heisenberg model on both square and honeycomb lattices across their (2+1)D O(3) QCPs with different opening angles on entanglement boundary, and obtain the accurate values of the corresponding universal corner log-coefficients. These findings will encourage further theoretical investigations to access controlled universal information for interacting CFTs at (2+1)D., Comment: 6+6 pages, 4+6 figures

We employ a protocol, dubbed entanglement microscopy, to reveal the multipartite entanglement encoded in the full reduced density matrix of microscopic subregion both in spin and fermionic many-body systems. We exemplify our method by studying the phase diagram near quantum critical points (QCP) in 2 spatial dimensions: the transverse field Ising model and a Gross-Neveu-Yukawa transition of Dirac fermions. Our main results are: i) the Ising QCP exhibits short-range entanglement with a finite sudden death of the LN both in space and temperature; ii) the Gross-Neveu QCP has a power-law decaying fermionic LN consistent with conformal field theory (CFT) exponents; iii) going beyond bipartite entanglement, we find no detectable 3-party entanglement with our two witnesses in a large parameter window near the Ising QCP in 2d, in contrast to 1d. We further establish the singular scaling of general multipartite entanglement measures at criticality, and present an explicit analysis in the tripartite case. We also analytically obtain the large-temperature power-law scaling of the fermionic LN for general interacting systems. Entanglement microscopy opens a rich window into quantum matter, with countless systems waiting to be explored., Comment: 10+10 pages, 4+12 figures

We systematically investigate the finite size scaling behavior of the R\'enyi entanglement entropy (EE) of several representative 2d quantum many-body systems between a subregion and its complement, with smooth boundaries as well as boundaries with corners. In order to reveal the subleading correction, we investigate the quantity ``subtracted EE" $S^s(l) = S(2l) - 2S(l)$ for each model, which is designed to cancel out the leading perimeter law. We find that $\mathbf{(1)}$ for a spin-1/2 model on a 2d square lattice whose ground state is the Neel order, the coefficient of the logarithmic correction to the perimeter law is consistent with the prediction based on the Goldstone modes; $\mathbf{(2)}$ for the $(2+1)d$ O(3) Wilson-Fisher quantum critical point (QCP), realized with the bilayer antiferromagnetic Heisenberg model, a logarithmic subleading correction exists when there is sharp corner of the subregion, but for subregion with a smooth boundary our data suggests the absence of the logarithmic correction to the best of our efforts; $\mathbf{(3)}$ for the $(2+1)d$ SU(2) J-Q$_2$ and J-Q$_3$ model for the deconfined quantum critical point (DQCP), we find a logarithmic correction for the EE even with smooth boundary.

Based on the recently developed resummation-based quantum Monte Carlo method for the SU($N$) spin and loop-gas models, we develop a new algorithm, dubbed ResumEE, to compute the entanglement entropy (EE) with greatly enhanced efficiency. Our ResumEE exponentially speeds up the computation of the exponentially small value of the $\langle e^{-S^{(2)}}\rangle$, where $S^{(2)}$ is the 2nd order R\'enyi EE, such that the $S^{(2)}$ for a generic 2D quantum SU($N$) spin models can be readily computed with high accuracy. We benchmark our algorithm with the previously proposed estimators of $S^{(2)}$ on 1D and 2D SU($2$) Heisenberg spin systems to reveal its superior performance and then use it to detect the entanglement scaling data of the N\'eel-to-VBS transition on 2D SU($N$) Heisenberg model with continuously varying $N$. Our ResumEE algorithm is efficient for precisely evaluating the entanglement entropy of SU($N$) spin models with continuous $N$ and reliable access to the conformal field theory data for the highly entangled quantum matter., Comment: 11 pages, 9 figures

In recent years a consensus has gradually been reached that the previously proposed deconfined quantum critical point (DQCP) for spin-1/2 systems, an archetypal example of quantum phase transition beyond the classic Landau's paradigm, actually does not correspond to a true unitary conformal field theory (CFT). In this work we carefully investigate another type of quantum phase transition supposedly beyond the similar classic paradigm, the so called ``symmetric mass generation" (SMG) transition proposed in recent years. We employ the sharp diagnosis including the scaling of disorder operator and R\'enyi entanglement entropy in large-scale lattice model quantum Monte Carlo simulations. Our results strongly suggest that the SMG transition is indeed an unconventional quantum phase transition and it should correspond to a true $(2+1)d$ unitary CFT., Comment: 16 pages, 12 figures

Over the past two decades, the enigma of the deconfined quantum critical point (DQCP) has attracted broad attention across the condensed matter, quantum field theory, and high-energy physics communities, as it is expected to offer a new paradigm in theory, experiment, and numerical simulations that goes beyond the Landau-Ginzburg-Wilson framework of symmetry breaking and phase transitions. However, the nature of DQCP has been controversial. For instance, in the square-lattice spin-1/2 $J$-$Q$ model, believed to realize the DQCP between N\'eel and valence bond solid states, conflicting results, such as first-order versus continuous transition, and critical exponents incompatible with conformal bootstrap bounds, have been reported. The enigma of DQCP is exemplified in its anomalous logarithmic subleading contribution in its entanglement entropy (EE), which was discussed in recent studies. In the current work, we demonstrate that similar anomalous logarithmic behavior persists in a class of models analogous to the DQCP. We systematically study the quantum EE of square-lattice SU($N$) DQCP spin models. Based on large-scale quantum Monte Carlo computation of the EE, we show that for a series of $N$ smaller than a critical value, the anomalous logarithmic behavior always exists in the EE, which implies that the previously determined DQCPs in these models do not belong to conformal fixed points. In contrast, when $N\ge N_c$ with a finite $N_c$ that we evaluate to lie between $7$ and $8$, the DQCPs are consistent with conformal fixed points that can be understood within the Abelian Higgs field theory with $N$ complex components., Comment: The revised version focuses on the anomalous logarithmic correction to the EE arising from the smooth boundary, rather than corners. And the critical $N_c$ is determined based on the anomalous EE signal from smooth boundary

The study of quantum criticality and entanglement in systems with long-range (LR) interactions is still in its early stages, with many open questions remaining. In this work, we investigate critical exponents and scaling of entanglement entropies (EE) in the LR bilayer Heisenberg model using large-scale quantum Monte Carlo (QMC) simulations and the recently developed nonequilibrium increment algorithm for measuring EE. By applying modified (standard) finite-size scaling (FSS) above (below) the upper critical dimension and field theory analysis, we obtain precise critical exponents in three regimes: the LR Gaussian regime with a Gaussian fixed point, the short-range (SR) regime with Wilson-Fisher (WF) exponents, and a LR non-Gaussian regime where the critical exponents vary continuously from LR Gaussian to SR values. We compute the R\'enyi EE both along the critical line and in the N\'eel phase and observe that as the LR interaction is enhanced, the area-law contribution in EE gradually vanishes both at quantum critical points (QCPs) and in the N\'eel phase. The log-correction in EE arising from sharp corners at the QCPs also decays to zero as LR interaction grows, whereas the log-correction for N\'eel states, caused by the interplay of Goldstone modes and restoration of the symmetry in a finite system, is enhanced as LR interaction becomes stronger. We also discuss relevant experimental settings to detect these nontrivial properties in critical behavior and entanglement information for quantum many-body systems with LR interactions., Comment: 5pages, 4 figures

The Mermin-Wagner theorem states that spontaneous continuous symmetry breaking is prohibited in systems with short-range interactions at spatial dimension $D\le 2$. For long-range interactions with a power-law form ($1/r^{\alpha}$), the theorem further forbids ferromagnetic or antiferromagnetic order at finite temperature when $\alpha\ge 2D$. However, the situation for $\alpha \in (2,4)$ at $D=2$ is not covered by the theorem. To address this, we conduct large-scale quantum Monte Carlo simulations and field theoretical analysis. Our findings show spontaneous breaking of $SU(2)$ symmetry in the ferromagnetic Heisenberg model with $1/r^{\alpha}$-form long-range interactions at $D=2$. We determine critical exponents through finite-size analysis for $\alpha<3$ (above the upper critical dimension with Gaussian fixed point) and $3\le\alpha<4$ (below the upper critical dimension with non-Gaussian fixed point). These results reveal new critical behaviors in 2D long-range Heisenberg models, encouraging further experimental studies of quantum materials with long-range interactions beyond the Mermin-Wagner theorem's scope.

Employing large-scale quantum Monte Carlo simulations, we systematically compute the energy spectra of the 2D spin-1/2 Heisenberg model with long-range interactions. With the $1/r^{\alpha}$ ferromagnetic and staggered antiferromagnetic interactions, we find the explicit range in $\alpha$ for {\color{black} the short-range Goldstone-type (gapless), anomalous Goldstone-type (gapless) and Higgs-type (gapped) spectra. Accompanied by the spin wave analysis, our numerical results vividly reveal how the long-range interactions alter the usual linear and quadratic magnon dispersions in 2D quantum magnets and give rise to anomalous dynamical exponents. Moreover, we find explicit case where the gapped excitation exists even when the Hamiltonian is extensive. This work provides the first set of unbiased dynamical data} of long-range quantum many-body systems and suggests that many universally accepted low-energy customs for short-range systems need to be substantially modified for long-range ones which are of immediate relevance to the ongoing experimental efforts from quantum simulators to 2D quantum moir\'e materials., Comment: 5 pages,3 figures

We propose a physical picture based on the wormhole effect of the path-integral formulation to explain the mechanism of entanglement spectrum (ES), such that, our picture not only explains the topological state with bulk-edge correspondence of the energy spectrum and ES (the Li and Haldane conjecture), but is generically applicable to other systems independent of their topological properties. We point out it is ultimately the relative strength of bulk energy gap (multiplied with inverse temperature $\beta=1/T$) with respect to the edge energy gap that determines the behavior of the low-lying ES of the system. Depending on the circumstances, the ES can resemble the energy spectrum of the virtual edge, but can also represent that of the virtual bulk. We design models both in 1D and 2D to successfully demonstrate the bulk-like low-lying ES at finite temperatures, in addition to the edge-like case conjectured by Li and Haldane at zero temperature. Our results support the generality of viewing the ES as the wormhole effect in the path integral and the different temperature-dependence for the edge and bulk of ES., Comment: 5 pages, 3 figures