-
Vestigial Gapless Boson Density Wave Emerging between $ν= 1/2$ Fractional Chern Insulator and Finite-Momentum Supersolid
Authors:
Hongyu Lu,
Han-Qing Wu,
Bin-Bin Chen,
Zi Yang Meng
Abstract:
The roton-triggered charge-density-wave (CDW)is widely studied in fractional quantum Hall (FQH) and fractional Chern insulator (FCI) systems, and there also exist field theoretical and numerical realizations of continuous transition from FCI to superfluid (SF). However, the theory and numerical explorations of the transition between FCI and supersolid (SS) are still lacking. In this work, we study…
▽ More
The roton-triggered charge-density-wave (CDW)is widely studied in fractional quantum Hall (FQH) and fractional Chern insulator (FCI) systems, and there also exist field theoretical and numerical realizations of continuous transition from FCI to superfluid (SF). However, the theory and numerical explorations of the transition between FCI and supersolid (SS) are still lacking. In this work, we study the topological flat-band lattice models with $ν$ = 1/2 hard-core bosons, where the previous studies have discovered the existence of FCI states and possible direct FCI-SS transitions. While the FCI is robust, we find the direct FCI-SS transition is absent, and there exist more intriguing scenarios. In the case of checkerboard lattice, we find an intermediate gapless CDW state without SF, sandwiched between FCI and SS. This novel state is triggered by the roton instability in FCI and it further continuously brings about the intertwined finite-momentum SF fluctuation when the CDW order is strong enough, eventually transiting into an unconventional finite-momentum SS state. The intermediate gapless CDW state is a vestige from the SS state, since the increasing quantum fluctuation melts only the Larkin-Ovchinnikov-type SF order in SS but its (secondary) product -- the CDW order -- survives. On honeycomb lattice, we find no evidence of SS, but discover an interesting sequence of FCI-Solid I-Solid II transitions, with both solids incompressible. Moreover, in contrast to previous single-roton condensation, this sequence of FCI-Solid I-Solid II transitions is triggered by the softening of multi-roton modes in FCI. Considering the intertwined wave vectors of the CDW orders, Solid I is a vestige of Solid II. Our work provides new horizon not only for the quantum phase transitions in FCI but also for the intertwined orders and gapless states in bosonic systems, which will inspire future studies.
△ Less
Submitted 13 August, 2024;
originally announced August 2024.
-
An analog of topological entanglement entropy for mixed states
Authors:
Ting-Tung Wang,
Menghan Song,
Zi Yang Meng,
Tarun Grover
Abstract:
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…
▽ More
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ényi 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.
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
Continuous Transition between Bosonic Fractional Chern Insulator and Superfluid
Authors:
Hongyu Lu,
Han-Qing Wu,
Bin-Bin Chen,
Zi Yang Meng
Abstract:
The properties of fractional Chern insulator (FCI) phase and the phase transitions between FCI and Mott insulators (MI) in bosonic systems are well studied. The continuous transitions between FCI and superfluid (SF), however, despite the inspiring field theoretical predictions, have not been directly verified. The existing numerical results of FCI-SF transition are either indirect or clearly first…
▽ More
The properties of fractional Chern insulator (FCI) phase and the phase transitions between FCI and Mott insulators (MI) in bosonic systems are well studied. The continuous transitions between FCI and superfluid (SF), however, despite the inspiring field theoretical predictions, have not been directly verified. The existing numerical results of FCI-SF transition are either indirect or clearly first-order. Here, by simply tuning the bandwidth of Haldane honeycomb lattice model, we find direct transitions from a bosonic FCI at $ν$ = 1/2 filling of a flat Chern band to two SF states with bosons condensed at momenta $M$ and $Γ$, respectively. While the FCI-SF($M$) transition is first-order, the FCI-SF($Γ$) transition is found continuous. Through finite size criticality analysis, the obtained critical exponents $β\approx$0.35(5) and $ν\approx$0.62(12) are both compatible with those of the 3D XY universality class and more exotic beyond-Landau ones. Our work thence presents a direct numerical evidence of a continuous FCI-SF transition between topological ordered phase and spontaneous continuous symmetry-breaking phase, and further indicates the zero-field bosonic FCI might be realized from a SF state by gradually flattening the dispersion of Chern band, through the (quasi)adiabatic preparation in ultracold atom systems.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Extracting Universal Corner Entanglement Entropy during the Quantum Monte Carlo Simulation
Authors:
Yuan Da Liao,
Menghan Song,
Jiarui Zhao,
Zi Yang Meng
Abstract:
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 co…
▽ More
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.
△ Less
Submitted 28 August, 2024; v1 submitted 22 April, 2024;
originally announced April 2024.
-
Entanglement Microscopy: Tomography and Entanglement Measures via Quantum Monte Carlo
Authors:
Ting-Tung Wang,
Menghan Song,
Liuke Lyu,
William Witczak-Krempa,
Zi Yang Meng
Abstract:
We develop a protocol, dubbed entanglement microscopy, to obtain the full reduced density matrix associated with subregions in quantum Monte Carlo simulations for bosonic and fermionic manybody systems. Our microscopy allows to perform quantum state tomography, and thus gives access to true entanglement measures, such as the logarithmic negativity (LN). We exemplify our method by studying the phas…
▽ More
We develop a protocol, dubbed entanglement microscopy, to obtain the full reduced density matrix associated with subregions in quantum Monte Carlo simulations for bosonic and fermionic manybody systems. Our microscopy allows to perform quantum state tomography, and thus gives access to true entanglement measures, such as the logarithmic negativity (LN). 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 in a large parameter window near the Ising QCP in 2d, in contrast to 1d. We also analytically obtain the large-temperature power law scaling of the fermionic LN for general interacting systems. Our approach allows one to perform quantum state tomography locally in a way that is analogous to atomic-scale imaging with a scanning tunneling microscope. Controlled entanglement microscopy opens a new window into quantum matter, with countless systems waiting to be explored.
△ Less
Submitted 6 May, 2024; v1 submitted 22 February, 2024;
originally announced February 2024.
-
An integral algorithm of exponential observables for interacting fermions in quantum Monte Carlo simulation
Authors:
Xu Zhang,
Gaopei Pan,
Bin-Bin Chen,
Kai Sun,
Zi Yang Meng
Abstract:
Exponential observables, formulated as $\log \langle e^{\hat{X}}\rangle$ where $\hat{X}$ is an extensive quantity, play a critical role in study of quantum many-body systems, examples of which include the free-energy and entanglement entropy. Given that $e^{X}$ becomes exponentially large (or small) in the thermodynamic limit, accurately computing the expectation value of this exponential quantity…
▽ More
Exponential observables, formulated as $\log \langle e^{\hat{X}}\rangle$ where $\hat{X}$ is an extensive quantity, play a critical role in study of quantum many-body systems, examples of which include the free-energy and entanglement entropy. Given that $e^{X}$ becomes exponentially large (or small) in the thermodynamic limit, accurately computing the expectation value of this exponential quantity presents a significant challenge. In this Letter, we propose a comprehensive algorithm for quantifying these observables in interacting fermion systems, utilizing the determinant quantum Monte Carlo (DQMC) method. We have applied this novel algorithm to the 2D half-filled Hubbard model. At the strong coupling limit, our method showcases a significant accuracy improvement compared to conventional methods that are derived from the internal energy. We also illustrate that this novel approach delivers highly efficient and precise measurements of the nth Rényi entanglement entropy. Even more noteworthy is that this improvement comes without incurring increases in computational complexity. This algorithm effectively suppresses exponential fluctuations and can be easily generalized to other models.
△ Less
Submitted 27 May, 2024; v1 submitted 6 November, 2023;
originally announced November 2023.
-
Resummation-based Quantum Monte Carlo for Entanglement Entropy Computation
Authors:
Menghan Song,
Ting-Tung Wang,
Zi Yang Meng
Abstract:
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ényi…
▽ More
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ényi 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éel-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.
△ Less
Submitted 2 September, 2024; v1 submitted 2 October, 2023;
originally announced October 2023.
-
Relevant long-range interaction of the entanglement Hamiltonian emerges from a short-range gapped system
Authors:
Chuhao Li,
Rui-Zhen Huang,
Yi-Ming Ding,
Zi Yang Meng,
Yan-Cheng Wang,
Zheng Yan
Abstract:
Beyond the Li-Haldane-Poilblanc conjecture, we find the entanglement Hamiltonian (EH) is actually not closely similar to the original Hamiltonian on the virtual edge. Unexpectedly, the EH has some relevant long-range interacting terms which hugely affect the physics. Without loss of generality, we study a spin-1/2 Heisenberg bilayer to obtain the entanglement information between the two layers thr…
▽ More
Beyond the Li-Haldane-Poilblanc conjecture, we find the entanglement Hamiltonian (EH) is actually not closely similar to the original Hamiltonian on the virtual edge. Unexpectedly, the EH has some relevant long-range interacting terms which hugely affect the physics. Without loss of generality, we study a spin-1/2 Heisenberg bilayer to obtain the entanglement information between the two layers through our newly developed quantum Monte Carlo scheme, which can simulate large-scale EH. Although the entanglement spectrum carrying the Goldstone mode seems like a Heisenberg model on a single layer, which is consistent with Li-Haldane-Poilblanc conjecture, we demonstrate that there actually exists a finite-temperature phase transition of the EH. The results violate the Mermin-Wagner theorem, which means there should be relevant long-range terms in the EH. It reveals that the Li-Haldane-Poilblanc conjecture ignores necessary corrections for the EH which may lead totally different physics.
△ Less
Submitted 8 June, 2024; v1 submitted 27 September, 2023;
originally announced September 2023.
-
Cubic* criticality emerging from a quantum loop model on triangular lattice
Authors:
Xiaoxue Ran,
Zheng Yan,
Yan-Cheng Wang,
Junchen Rong,
Yang Qi,
Zi Yang Meng
Abstract:
Quantum loop and dimer models are archetypal examples of correlated systems with local constraints. Obtaining generic solutions for these models is difficult due to the lack of controlled methods to solve them in the thermodynamic limit. Nevertheless, these solutions are of immediate relevance to both statistical and quantum field theories, as well as the rapidly growing experiments in Rydberg ato…
▽ More
Quantum loop and dimer models are archetypal examples of correlated systems with local constraints. Obtaining generic solutions for these models is difficult due to the lack of controlled methods to solve them in the thermodynamic limit. Nevertheless, these solutions are of immediate relevance to both statistical and quantum field theories, as well as the rapidly growing experiments in Rydberg atom arrays and quantum moiré materials, where the interplay between correlation and local constraints gives rise to a plethora of novel phenomena. In a recent work [X. Ran, Z. Yan, Y.-C. Wang, et al, arXiv:2205.04472 (2022)], it was found through sweeping cluster quantum Monte Carlo (QMC) simulations and field theory analysis that the triangular lattice quantum loop model (QLM) hosts a rich ground state phase diagram with lattice nematic, vison plaquette (VP) crystals, and the $\mathbb{Z}_2$ quantum spin liquid (QSL) close to the Rokhsar-Kivelson point. Here, we focus on the continuous quantum critical point separating the VP and QSL phases and demonstrate via both static and dynamic probes in QMC simulations that this transition is of the (2+1)D cubic* universality. In this transition, the fractionalized visons in QSL condense to give rise to the crystalline VP phase, while leaving their trace in the anomalously large anomalous dimension exponent and pronounced continua in the dimer and vison spectra compared with those at the conventional cubic or O(3) quantum critical points.
△ Less
Submitted 11 June, 2024; v1 submitted 11 September, 2023;
originally announced September 2023.
-
Stable computation of entanglement entropy for 2D interacting fermion systems
Authors:
Gaopei Pan,
Yuan Da Liao,
Weilun Jiang,
Jonathan D'Emidio,
Yang Qi,
Zi Yang Meng
Abstract:
There is no doubt that the information hidden in entanglement entropy (EE), for example, the $n$-th order Rényi EE, i.e., $S^{A}_n=\frac{1}{1-n}\ln \Tr (ρ_A^n)$ where $ρ_A=\mathrm{Tr}_{\overline{A}}ρ$ is the reduced density matrix, can be used to infer the organizing principle of 2D interacting fermion systems, ranging from spontaneous symmetry breaking phases, quantum critical points to topologic…
▽ More
There is no doubt that the information hidden in entanglement entropy (EE), for example, the $n$-th order Rényi EE, i.e., $S^{A}_n=\frac{1}{1-n}\ln \Tr (ρ_A^n)$ where $ρ_A=\mathrm{Tr}_{\overline{A}}ρ$ is the reduced density matrix, can be used to infer the organizing principle of 2D interacting fermion systems, ranging from spontaneous symmetry breaking phases, quantum critical points to topologically ordered states. It is far from clear, however, whether the EE can actually be obtained with the precision required to observe these fundamental features -- usually in the form of universal finite size scaling behavior. Even for the prototypical 2D interacting fermion model -- the Hubbard model, to all existing numerical algorithms, the computation of the EE has not been succeeded with reliable data that the universal scaling regime can be accessed. Here we explain the reason for these unsuccessful attempts in EE computations in quantum Monte Carlo simulations in the past decades and more importantly, show how to overcome the conceptual and computational barrier with the incremental algorithm, such that the stable computation of the EE in 2D interacting fermion systems can be achieved and universal scaling information can be extracted. Relevance towards the experimental 2D interacting fermion systems is discussed.
△ Less
Submitted 28 August, 2023; v1 submitted 24 March, 2023;
originally announced March 2023.
-
The teaching from entanglement: 2D SU(2) antiferromagnet to valence bond solid deconfined quantum critical points are not conformal
Authors:
Yuan Da Liao,
Gaopei Pan,
Weilun Jiang,
Yang Qi,
Zi Yang Meng
Abstract:
The deconfined quantum critical point (DQCP) -- the enigmatic incarnation of the quantum phase transition beyond the Landau-Ginzburg-Wilson paradigm of symmetries and their spontaneous breaking -- has been proposed and actively pursued for more than two decades. Various 2D quantum many-body lattice models, both in spin/boson and fermion representations have been tested with the state-of-the-art nu…
▽ More
The deconfined quantum critical point (DQCP) -- the enigmatic incarnation of the quantum phase transition beyond the Landau-Ginzburg-Wilson paradigm of symmetries and their spontaneous breaking -- has been proposed and actively pursued for more than two decades. Various 2D quantum many-body lattice models, both in spin/boson and fermion representations have been tested with the state-of-the-art numerical techniques and field-theoretical analyses, and yet, the conclusion is still controversial. Experimental realizations of DQCP in the quantum magnet SrCu$_2$(BO$_3$)$_2$ and superconducting quantum criticality in 2D material have either shown first order transition or intermediate phase. The tension between the lattice scale details and the requirement from continuum limit, manifested in the form of the inconsistent critical scaling behavior and violations of generic conformal bootstrap bound, has not been resolved. Here we solve these decades-long controversies from the new and fundamental perspective of the quantum entanglement. We develop the incremental algorithm to compute the entanglement entropy at a fermionic DQCP with unprecedentedly accurate data and reveal the universal coefficient of the logarithmic correction therein is negative and at odds with positivity requirement of the conformal field theory. Together with results in other 2D DQCP lattice models (both in fermion and spin systems), our discoveries clearly demonstrate the 2D SU(2) antiferromagnet to valence bond solid DQCPs are not conformal fixed point and naturally explain the experimental difficulties in finding them. This marks the end of the beginning of unambiguous finding of the quantum phase transitions truely beyond the Landau-Ginzburg-Wilson paradigm, since its suggestion two decades ago.
△ Less
Submitted 1 May, 2023; v1 submitted 22 February, 2023;
originally announced February 2023.
-
Emergent glassy behavior in a kagome Rydberg atom array
Authors:
Zheng Yan,
Yan-Cheng Wang,
Rhine Samajdar,
Subir Sachdev,
Zi Yang Meng
Abstract:
We present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal \textit{emergent} glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glass…
▽ More
We present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal \textit{emergent} glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids -- as well as the crossover towards a trivial paramagnetic phase -- are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.
△ Less
Submitted 18 May, 2023; v1 submitted 17 January, 2023;
originally announced January 2023.
-
Different temperature-dependence for the edge and bulk of entanglement Hamiltonian
Authors:
Menghan Song,
Jiarui Zhao,
Zheng Yan,
Zi Yang Meng
Abstract:
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 p…
▽ More
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 $β=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.
△ Less
Submitted 20 July, 2023; v1 submitted 18 October, 2022;
originally announced October 2022.
-
Detecting Subsystem Symmetry Protected Topological Order Through Strange Correlators
Authors:
Chengkang Zhou,
Meng-Yuan Li,
Zheng Yan,
Peng Ye,
Zi Yang Meng
Abstract:
We employ strange correlators to detect 2D subsystem symmetry-protected topological (SSPT) phases which are nontrivial topological phases protected by subsystem symmetries. Specifically, we analytically construct efficient strange correlators in the 2D cluster model in the presence of a uniform magnetic field and then perform the projector Quantum Monte Carlo simulation within the quantum annealin…
▽ More
We employ strange correlators to detect 2D subsystem symmetry-protected topological (SSPT) phases which are nontrivial topological phases protected by subsystem symmetries. Specifically, we analytically construct efficient strange correlators in the 2D cluster model in the presence of a uniform magnetic field and then perform the projector Quantum Monte Carlo simulation within the quantum annealing scheme. We find that strange correlators show the long-range correlation in the SSPT phase, from which we define strange order parameters to characterize the topological phase transition between the SSPT phase at low fields and the trivial paramagnetic phase at high fields. Thus, the detection of the fully localized zero modes on the 1D physical boundary of SSPT phase has been transformed into the bulk correlation measurement about the local operators with the periodic boundary condition. We also find interesting spatial anisotropy of a strange correlator, which can be intrinsically traced back to the nature of spatial anisotropy of subsystem symmetries that protect SSPT order in the 2D cluster model. By simulating strange correlators, we, therefore, provide the first unbiased large-scale quantum Monte Carlo simulation on the easy and efficient detection in the SSPT phase and open the avenue of the investigation of the subtle yet fundamental nature of the novel interacting topological phases.
△ Less
Submitted 23 December, 2022; v1 submitted 26 September, 2022;
originally announced September 2022.
-
Fully packed quantum loop model on the square lattice: phase diagram and application for Rydberg atoms
Authors:
Xiaoxue Ran,
Zheng Yan,
Yan-Cheng Wang,
Junchen Rong,
Yang Qi,
Zi Yang Meng
Abstract:
The quantum dimer and loop models attract great attentions, partially because the fundamental importance in the phases and phase transitions emerging in these prototypical constrained systems, and partially due to their intimate relevance toward the on-going experiments on Rydberg atom arrays in which the blockade mechanism naturally enforces the local constraint. Here we show, by means of the swe…
▽ More
The quantum dimer and loop models attract great attentions, partially because the fundamental importance in the phases and phase transitions emerging in these prototypical constrained systems, and partially due to their intimate relevance toward the on-going experiments on Rydberg atom arrays in which the blockade mechanism naturally enforces the local constraint. Here we show, by means of the sweeping cluster quantum Monte Carlo method, the complete ground state phase diagram of the fully packed quantum loop model on the square lattice. We find between the lattice nematic (LN) phase with strong dimer attraction and the staggered phase (SP) with strong dimer repulsion, there emerges a resonating plaquette (RP) phase with off-diagonal translational symmetry breaking. Such a quantum phase is separated from the LN via a first order transition and from the SP by the famous Rokhsar-Kivelson point. Our renormalization group analysis reveals the different flow directions, fully consistent with the order parameter histogram in Monte Carlo simulations. The realization and implication of our phase diagram in Rydberg experiments are proposed.
△ Less
Submitted 20 March, 2023; v1 submitted 21 September, 2022;
originally announced September 2022.
-
Dirac fermions with plaquette interactions. III. SU(N) phase diagram with Gross-Neveu criticality and first-order phase transition
Authors:
Yuan Da Liao,
Xiao Yan Xu,
Zi Yang Meng,
Yang Qi
Abstract:
Inspired by our recent works[1, 2] of SU(2) and SU(4) Dirac fermions subjected to plaquette interactions on square lattice, here we extend the large-scale quantum Monte Carlo investigations to the phase digram of correlated Dirac fermions with SU(6) and SU(8) symmetries subjected to the plaquette interaction on the same lattice. From SU(2) to SU(8), the rich phase diagram exhibits a plethora of em…
▽ More
Inspired by our recent works[1, 2] of SU(2) and SU(4) Dirac fermions subjected to plaquette interactions on square lattice, here we extend the large-scale quantum Monte Carlo investigations to the phase digram of correlated Dirac fermions with SU(6) and SU(8) symmetries subjected to the plaquette interaction on the same lattice. From SU(2) to SU(8), the rich phase diagram exhibits a plethora of emerging quantum phases such as the Dirac semimetal, the antiferromagnetic Mott insulator, valence bond solid (VBS) and the Dirac spin liquid and phase transitions including the Gross-Neveu chiral transitions with emergent continuous symmetry, the deconfined quantum criticality and the first order transition between interaction-driven columnar VBS and plaquette VBS. These rich phenomena coming from the simple-looking lattice models, firmly convey the message that the interplay between the $SU(N)$ Dirac fermions -- with enhanced internal symmetries -- and extended plaquette interactions -- beyond the on-site Hubbard type -- is the new playground to synthesise novel highly entangled quantum matter both at the model level and with experimental feasibilities.
△ Less
Submitted 1 November, 2022; v1 submitted 27 July, 2022;
originally announced July 2022.
-
Hidden orders and phase transitions for the fully packed quantum loop model on the triangular lattice
Authors:
Xiaoxue Ran,
Zheng Yan,
Yan-Cheng Wang,
Rhine Samajdar,
Junchen Rong,
Subir Sachdev,
Yang Qi,
Zi Yang Meng
Abstract:
Quantum loop and dimer models are prototypical correlated systems with local constraints, which are not only intimately connected to lattice gauge theories and topological orders but are also widely applicable to the broad research areas of quantum materials and quantum simulation. Employing our sweeping cluster quantum Monte Carlo algorithm, we reveal the complete phase diagram of the triangular-…
▽ More
Quantum loop and dimer models are prototypical correlated systems with local constraints, which are not only intimately connected to lattice gauge theories and topological orders but are also widely applicable to the broad research areas of quantum materials and quantum simulation. Employing our sweeping cluster quantum Monte Carlo algorithm, we reveal the complete phase diagram of the triangular-lattice fully packed quantum loop model. Apart from the known lattice nematic (LN) solid and the even $\mathbb{Z}_2$ quantum spin liquid (QSL) phases, we discover a hidden vison plaquette (VP) phase, which had been overlooked and misinterpreted as a QSL for more than a decade. Moreover, the VP-to-QSL continuous transition belongs to the $(2+1)$D cubic* universality class, which offers a lattice realization of the (fractionalized) cubic fixed point that had long been considered as irrelevant towards the O($3$) symmetry until corrected recently by conformal bootstrap calculations. Our results are therefore of relevance to recent developments in both experiments and theory, and facilitate further investigations of hidden phases and transitions.
△ Less
Submitted 28 June, 2024; v1 submitted 9 May, 2022;
originally announced May 2022.
-
Evolution of Dynamical Signature in the X-cube Fracton Topological Order
Authors:
Chengkang Zhou,
Meng-Yuan Li,
Zheng Yan,
Peng Ye,
Zi Yang Meng
Abstract:
As an unconventional realization of topological orders with an exotic interplay of topology and geometry, fracton (topological) orders feature subextensive topological ground state degeneracy and subdimensional excitations that are movable only within a certain subspace. It has been known in the exactly solvable three-dimensional X-cube model that universally represents the type-I fracton orders,…
▽ More
As an unconventional realization of topological orders with an exotic interplay of topology and geometry, fracton (topological) orders feature subextensive topological ground state degeneracy and subdimensional excitations that are movable only within a certain subspace. It has been known in the exactly solvable three-dimensional X-cube model that universally represents the type-I fracton orders, that mobility constraints on subdimensional excitations originate from the absence of spatially deformable string-like operators. To unveil the interplay of topology and geometry, in this paper, we study the dynamical signature in the X-cube model in the presence of external Zeeman fields via large-scale quantum Monte Carlo simulation and stochastic analytic continuation. We compute both real-space correlation functions and dynamic structure factors of subdimensional excitations (i.e., fractons, lineons, and planons) in the fracton phase and their evolution into the trivial paramagnetic phase by increasing external fields. We find in the fracton phase, that the correlation functions and the spectral functions show clear anisotropy exactly caused by the underlying mobility constraints. On the other hand, the external fields successfully induce quantum fluctuations and offer mobility to excitations along the subspace allowed by mobility constraints. These numerical results provide the evolution of a dynamical signature of subdimensional particles in fracton orders, indicating that the mobility constraints on local dynamical properties of subdimensional excitations are deeply related to the existence of fracton topological order. The results will also be helpful in potential experimental identifications in spectroscopy measurements such as neutron scattering and nuclear magnetic resonance.
△ Less
Submitted 10 August, 2022; v1 submitted 24 March, 2022;
originally announced March 2022.
-
Triangular lattice quantum dimer model with variable dimer density
Authors:
Zheng Yan,
Rhine Samajdar,
Yan-Cheng Wang,
Subir Sachdev,
Zi Yang Meng
Abstract:
Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We…
▽ More
Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We find distinct odd and even $\mathbb{Z}_2$ spin liquids, along with several phases with no topological order: a staggered crystal, a nematic phase, and a trivial symmetric phase with no obvious broken symmetry. We also present dynamic spectra of the phases, and note implications for experiments on Rydberg atoms.
△ Less
Submitted 2 October, 2022; v1 submitted 22 February, 2022;
originally announced February 2022.
-
Unlocking the general relationship between energy and entanglement spectra via the wormhole effect
Authors:
Zheng Yan,
Zi Yang Meng
Abstract:
Based on the path integral formulation of the reduced density matrix, we develop a scheme to overcome the exponential growth of computational complexity in reliably extracting low-lying entanglement spectrum from quantum Monte Carlo simulations. We test the method on the Heisenberg spin ladder with long entangled boundary between two chains and the results support the Li and Haldane's conjecture o…
▽ More
Based on the path integral formulation of the reduced density matrix, we develop a scheme to overcome the exponential growth of computational complexity in reliably extracting low-lying entanglement spectrum from quantum Monte Carlo simulations. We test the method on the Heisenberg spin ladder with long entangled boundary between two chains and the results support the Li and Haldane's conjecture on entanglement spectrum of topological phase. We then explain the conjecture via the wormhole effect in the path integral and show that it can be further generalized for systems beyond gapped topological phases. Our further simulation results on the bilayer antiferromagnetic Heisenberg model with 2D entangled boundary across the $(2+1)$D O(3) quantum phase transition clearly demonstrate the correctness of the wormhole picture. Finally, we state that since the wormhole effect amplifies the bulk energy gap by a factor of $β$, the relative strength of that with respect to the edge energy gap will determine the behavior of low-lying entanglement spectrum of the system.
△ Less
Submitted 26 April, 2023; v1 submitted 10 December, 2021;
originally announced December 2021.
-
Quantum optimization within lattice gauge theory model on a quantum simulator
Authors:
Zheng Yan,
Zheng Zhou,
Yan-Hua Zhou,
Yan-Cheng Wang,
Xingze Qiu,
Zi Yang Meng,
Xue-Feng Zhang
Abstract:
Simulating lattice gauge theory (LGT) Hamiltonian and its nontrivial states by programmable quantum devices has attracted numerous attention in recent years. Rydberg atom arrays constitute one of the most rapidly developing arenas for quantum simulation and quantum computing. The $\mathbb{Z}_2$ LGT and topological order has been realized in experiments while the $U(1)$ LGT is being worked hard on…
▽ More
Simulating lattice gauge theory (LGT) Hamiltonian and its nontrivial states by programmable quantum devices has attracted numerous attention in recent years. Rydberg atom arrays constitute one of the most rapidly developing arenas for quantum simulation and quantum computing. The $\mathbb{Z}_2$ LGT and topological order has been realized in experiments while the $U(1)$ LGT is being worked hard on the way. States of LGT have local constraint and are fragmented into several winding sectors with topological protection. It is therefore difficult to reach the ground state in target sector for experiments, and it is also an important task for quantum topological memory. Here, we propose a protocol of sweeping quantum annealing (SQA) for searching the ground state among topological sectors. With the quantum Monte Carlo method, we show that this SQA has linear time complexity of size with applications to the antiferromagnetic transverse field Ising model, which has emergent $U(1)$ gauge fields. This SQA protocol can be realized easily on quantum simulation platforms such as Rydberg array and D-wave annealer. We expect this approach would provide an efficient recipe for resolving the topological hindrances in quantum optimization and the preparation of quantum topological state.
△ Less
Submitted 18 September, 2023; v1 submitted 15 May, 2021;
originally announced May 2021.
-
Evidence of the Berezinskii-Kosterlitz-Thouless Phase in a Frustrated Magnet
Authors:
Ze Hu,
Zhen Ma,
Yuan-Da Liao,
Han Li,
Chunsheng Ma,
Yi Cui,
Yanyan Shangguan,
Zhentao Huang,
Yang Qi,
Wei Li,
Zi Yang Meng,
Jinsheng Wen,
Weiqiang Yu
Abstract:
The Berezinskii-Kosterlitz-Thouless (BKT) mechanism, building upon proliferation of topological defects in 2D systems, is the first example of phase transition beyond the Landau-Ginzburg paradigm of symmetry breaking. Such a topological phase transition has long been sought yet undiscovered directly in magnetic materials. Here, we pin down two transitions that bound a BKT phase in an ideal 2D frus…
▽ More
The Berezinskii-Kosterlitz-Thouless (BKT) mechanism, building upon proliferation of topological defects in 2D systems, is the first example of phase transition beyond the Landau-Ginzburg paradigm of symmetry breaking. Such a topological phase transition has long been sought yet undiscovered directly in magnetic materials. Here, we pin down two transitions that bound a BKT phase in an ideal 2D frustrated magnet TmMgGaO$_4$, via nuclear magnetic resonance under in-plane magnetic fields, which do not disturb the low-energy electronic states and allow BKT fluctuations to be detected sensitively. Moreover, by applying out-of-plane fields, we find a critical scaling behaviour of the magnetic susceptibility expected for the BKT transition. The experimental findings can be explained by quantum Monte Carlo simulations applied on an accurate triangular-lattice Ising model of the compound which hosts a BKT phase. These results provide a concrete example for the BKT phase and offer an ideal platform for future investigations on the BKT physics in magnetic materials.
△ Less
Submitted 13 October, 2020;
originally announced October 2020.
-
Quantum Simulation of Lattice Gauge Theories on Superconducting Circuits: Quantum Phase Transition and Quench Dynamics
Authors:
Zi-Yong Ge,
Rui-Zhen Huang,
Zi Yang Meng,
Heng Fan
Abstract:
Recently, quantum simulation of low-dimensional lattice gauge theories (LGTs) has attracted many interests, which may improve our understanding of strongly correlated quantum many-body systems. Here, we propose an implementation to approximate $\mathbb{Z}_2$ LGT on superconducting quantum circuits, where the effective theory is a mixture of a LGT and a gauge-broken term. Using matrix product state…
▽ More
Recently, quantum simulation of low-dimensional lattice gauge theories (LGTs) has attracted many interests, which may improve our understanding of strongly correlated quantum many-body systems. Here, we propose an implementation to approximate $\mathbb{Z}_2$ LGT on superconducting quantum circuits, where the effective theory is a mixture of a LGT and a gauge-broken term. Using matrix product state based methods, both the ground state properties and quench dynamics are systematically investigated. With an increase of the transverse (electric) field, the system displays a quantum phase transition from a disordered phase to a translational symmetry breaking phase. In the ordered phase, an approximate Gauss law of the $\mathbb{Z}_2$ LGT emerges in the ground state. Moreover, to shed light on the experiments, we also study the quench dynamics, where there is a dynamical signature of the spontaneous translational symmetry breaking. The spreading of the single particle of matter degree is diffusive under the weak transverse field, while it is ballistic with small velocity for the strong field. Furthermore, due to the emergent Gauss law under the strong transverse field, the matter degree can also exhibit confinement dynamics which leads to a strong suppression of the nearest-neighbor hopping. Our results pave the way for simulating the LGT on superconducting circuits, including the quantum phase transition and quench dynamics.
△ Less
Submitted 9 February, 2022; v1 submitted 28 September, 2020;
originally announced September 2020.
-
Charge-Density-Wave Transitions of Dirac Fermions Coupled to Phonons
Authors:
Chuang Chen,
Xiao Yan Xu,
Zi Yang Meng,
Martin Hohenadler
Abstract:
The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as…
▽ More
The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as well as a line of thermal critical points. Finite-size scaling appears consistent with fermionic Gross-Neveu-Ising universality for the quantum phase transition, and bosonic Ising universality for the thermal phase transition. The critical temperature has a maximum at intermediate couplings. Our findings motivate experimental efforts to identify or engineer Dirac systems with sufficiently strong and tunable electron-phonon coupling.
△ Less
Submitted 20 February, 2019; v1 submitted 20 September, 2018;
originally announced September 2018.
-
Quantum phase transitions in the Kane-Mele-Hubbard model
Authors:
M. Hohenadler,
Z. Y. Meng,
T. C. Lang,
S. Wessel,
A. Muramatsu,
F. F. Assaad
Abstract:
We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adiabatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricte…
▽ More
We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adiabatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricted to the transverse direction. The transition from the topological band insulator to the antiferromagnetic Mott insulator is in the universality class of the three-dimensional XY model. The numerical data suggest that the spin liquid to topological insulator and spin liquid to Mott insulator transitions are both continuous.
△ Less
Submitted 29 March, 2012; v1 submitted 16 November, 2011;
originally announced November 2011.