Showing 1–17 of 17 results for author: Da Liao, Y

Extract Universal Corner Entanglement Entropy in a Single 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 of stable computation of EE in (2+1)D quantum many body systems, we develop a new method for directly measuring the corner contribution in EE with less computational… ▽ 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 of stable computation of EE in (2+1)D quantum many body systems, we develop 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 corner entanglement boundaries in \textit{a single} 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. 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 corner opening angles and obtain the accurate values of the corresponding universal corner log-coefficients. These findings will encourage further theoretical investigations in this field. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 6+5 pages, 4+3 figures

Universal collective Larmor-Silin mode emerging in magnetized correlated Dirac fermions

Authors: Chuang Chen, Yuan Da Liao, Chengkang Zhou, Gaopei Pan, Zi Yang Meng, Yang Qi

Abstract: Employing large-scale quantum Monte Carlo simulations, we find in magnetized interacting Dirac fermion model, there emerges a new and universal collective Larmor-Silin spin wave mode in the transverse dynamical spin susceptibility. Such mode purely originates from the interaction among Dirac fermions and distinguishes itself from the usual particle-hole continuum with finite lifetime and clear dis… ▽ More Employing large-scale quantum Monte Carlo simulations, we find in magnetized interacting Dirac fermion model, there emerges a new and universal collective Larmor-Silin spin wave mode in the transverse dynamical spin susceptibility. Such mode purely originates from the interaction among Dirac fermions and distinguishes itself from the usual particle-hole continuum with finite lifetime and clear dispersion, both at small and large momenta in a large portion of the Brillouin zone. Our unbiased numerical results offer the dynamic signature of this new collective excitations in interacting Dirac fermion systems, and provide experimental guidance for inelastic neutron scattering, electron spin resonance and other spectroscopic approaches in the investigation of such universal collective modes in quantum Moire materials, topological insulators and quantum spin liquid materials under magnetic field, with quintessential interaction nature beyond the commonly assumed noninteracting Dirac fermion or spinon approximations. △ Less

Submitted 29 January, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 11 pages, 5 figures

Incremental SWAP Operator for Entanglement Entropy: Application for Exponential Observables in Quantum Monte Carlo Simulation

Authors: Xuan Zhou, Zi Yang Meng, Yang Qi, Yuan Da Liao

Abstract: We propose a new method to efficiently compute the entanglement entropy (EE) of quantum many-body systems. Our approach, called the incremental SWAP operator method , combines the simplicity of the SWAP operator used in projector quantum Monte Carlo simulations with recent advances in precisely computing exponential observables using incremental algorithms. We apply this technique to obtain accura… ▽ More We propose a new method to efficiently compute the entanglement entropy (EE) of quantum many-body systems. Our approach, called the incremental SWAP operator method , combines the simplicity of the SWAP operator used in projector quantum Monte Carlo simulations with recent advances in precisely computing exponential observables using incremental algorithms. We apply this technique to obtain accurate EE data at reduced computational cost for 1d and 2d antiferromagnetic Heisenberg models with different bipartition schemes. Using the computed EE data, we extract the area law coefficient, universal logarithmic corrections from Goldstone modes, and the geometric constant, finding quantitative agreement with analytical predictions. Moreover, for the first time in an unbiased numerical simulation of 2d antiferromagnetic Heisenberg model, we successfully obtain reliable universal logarithmic corrections from sharp corners that match expected theoretical values. The consistency between our numerical results and theoretical calculations demonstrates the power of our approach for accessing challenging universal entanglement properties. The extensions of our method to other quantum spin/boson models and the interacting fermion models, are outlined. △ Less

Submitted 2 April, 2024; v1 submitted 14 January, 2024; originally announced January 2024.

Journal ref: Phys. Rev. B 109, 165106 (2024)

Disorder Operator and Rényi Entanglement Entropy of Symmetric Mass Generation

Authors: Zi Hong Liu, Yuan Da Liao, Gaopei Pan, Menghan Song, Jiarui Zhao, Weilun Jiang, Chao-Ming Jian, Yi-Zhuang You, Fakher F. Assaad, Zi Yang Meng, Cenke Xu

Abstract: 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 supposedl… ▽ More 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ényi 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. △ Less

Submitted 17 March, 2024; v1 submitted 14 August, 2023; originally announced August 2023.

Comments: 16 pages, 12 figures

Journal ref: PhysRevLett.132.156503(2024)

Controllable Incremental Algorithm for Entanglement Entropy in Quantum Monte Carlo Simulations

Authors: Yuan Da Liao

Abstract: Researchers in physical science persistently seek universal features in strongly interacting many-body systems. However, these traits often remain concealed within highly complex observables like entanglement entropy (EE). The non-local essence of such observables renders experimental measurement and numerical evaluation challenging. Thus, developing an accurate and convenient algorithm to obtain… ▽ More Researchers in physical science persistently seek universal features in strongly interacting many-body systems. However, these traits often remain concealed within highly complex observables like entanglement entropy (EE). The non-local essence of such observables renders experimental measurement and numerical evaluation challenging. Thus, developing an accurate and convenient algorithm to obtain these observables is imperative. In this work, aided by quantum Monte Carlo, we unveil an exponentially exploding coefficient of variation with system size in the direct calculation of EE, making the evaluation of EE inaccurate. We thus introduce a controllable incremental algorithm to solve this conundrum. We demonstrate the efficacy and convenience of our algorithm by accurately evaluating the EE of a 2D Hubbard model.Moreover, our broadly adaptable technique can calculate numerically unstable observables with exploding coefficient of variation, including partition function, entanglement spectra, and negativity in correlated bosonic/fermionic systems. Precisely measuring these complex observables may motivate theoretical and experimental △ Less

Submitted 27 December, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: 5 pages, 4 figures

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.

Comments: 6+4 pages, 4+4 figures

Journal ref: Phys. Rev. B 108, L081123 (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.

Comments: 9+5 pages, 3+1 figures

Caution on Gross-Neveu criticality with a single Dirac cone: Violation of locality and its consequence of unexpected finite-temperature transition

Authors: Yuan Da Liao, Xiao Yan Xu, Zi Yang Meng, Yang Qi

Abstract: Lately there are many SLAC fermion investigations on the (2+1)D Gross-Neveu criticality of a single Dirac cone [1,2]. While the SLAC fermion construction indeed gives rise to the linear energy-momentum relation for all lattice momenta at the non-interacting limit, the long-range hopping and its consequent violation of locality on the Gross-Neveu quantum critical point (GN-QCP) -- which a priori re… ▽ More Lately there are many SLAC fermion investigations on the (2+1)D Gross-Neveu criticality of a single Dirac cone [1,2]. While the SLAC fermion construction indeed gives rise to the linear energy-momentum relation for all lattice momenta at the non-interacting limit, the long-range hopping and its consequent violation of locality on the Gross-Neveu quantum critical point (GN-QCP) -- which a priori requires short-range interaction -- has not been verified. Here we show, by means of large-scale quantum Monte Carlo simulations, that the interaction-driven antiferromagnetic insulator in this case is fundamentally different from that on a purely local $πぱい$-flux Hubbard model on the square lattice. In particular, we find the antiferromagnetic long-range order in the SLAC fermion model has a finite temperature continuous phase transition, which violates the Mermin-Wagner theorem, and smoothly connects to the previously determined GN-QCP. The magnetic excitations inside the antiferromagnetic insulator are gapped without Goldstone mode, even though the state spontaneously breaks continuous $SU(2)$ symmetry. These unusual results proclaim caution on the interpretation of the quantum phase transition in SLAC fermion model as that of GN-QCP with short-range interaction. △ Less

Submitted 7 November, 2023; v1 submitted 9 October, 2022; originally announced October 2022.

Journal ref: Phys.Rev.B 108,195112(2023)

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.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. B 106, 155159 (2022)

Dirac fermions with plaquette interactions. II. SU(4) phase diagram with Gross-Neveu criticality and quantum spin liquid

Authors: Yuan Da Liao, Xiao Yan Xu, Zi Yang Meng, Yang Qi

Abstract: At sufficiently low temperatures, interacting electron systems tend to develop orders. Exceptions are quantum critical point (QCP) and quantum spin liquid (QSL), where fluctuations prevent the highly entangled quantum matter to an ordered state down to the lowest temperature. While the ramification of these states may have appeared in high-temperature superconductors, ultra-cold atoms, frustrated… ▽ More At sufficiently low temperatures, interacting electron systems tend to develop orders. Exceptions are quantum critical point (QCP) and quantum spin liquid (QSL), where fluctuations prevent the highly entangled quantum matter to an ordered state down to the lowest temperature. While the ramification of these states may have appeared in high-temperature superconductors, ultra-cold atoms, frustrated magnets and quantum Moiré materials, their unbiased presence remain elusive in microscopic two-dimensional lattice models. Here, we show by means of large-scale quantum Monte Carlo simulations of correlated electrons on the $πぱい$-flux square lattice subjected to extended Hubbard interaction, that a Gross-Neveu QCP separating massless Dirac fermions and an columnar valence bond solid at finite interaction, and a possible Dirac QSL at the infinite yet tractable interaction limit emerge in a coherent sequence. These unexpected novel quantum states resides in this simple-looking model, unifying ingredients including emergent symmetry, deconfined fractionalization and the dynamic coupling between emergent matter and gauge fields, will have profound implications both in quantum many-body theory and understanding of the aforementioned experimental systems. △ Less

Submitted 29 September, 2022; v1 submitted 14 May, 2022; originally announced May 2022.

Comments: 12 pages, 7 figures

Journal ref: PhysRevB.106.115149(2022)

Dirac fermions with plaquette interactions. I. SU(2) phase diagram with Gross-Neveu and deconfined quantum criticalities

Authors: Yuan Da Liao, Xiao Yan Xu, Zi Yang Meng, Yang Qi

Abstract: We investigate the ground state phase diagram of an extended Hubbard model with $πぱい$-flux hopping term at half-filling on a square lattice, with unbiased large-scale auxiliary-field quantum Monte Carlo simulations. As a function of interaction strength, there emerges an intermediate phase which realizes two interaction-driven quantum critical points, with the first between the Dirac semimetal and a… ▽ More We investigate the ground state phase diagram of an extended Hubbard model with $πぱい$-flux hopping term at half-filling on a square lattice, with unbiased large-scale auxiliary-field quantum Monte Carlo simulations. As a function of interaction strength, there emerges an intermediate phase which realizes two interaction-driven quantum critical points, with the first between the Dirac semimetal and an insulating phase of weak valence bond solid (VBS) order, and the second separating the VBS order and an antiferromagnetic insulating phase. These intriguing quantum critical points are respectively bestowed with Gross-Neveu and deconfined quantum criticalities, and the critical exponents $ηいーた_\text{VBS}=0.6(1)$ and $ηいーた_\text{AF}=0.58(3)$ at deconfined quantum critical point satisfy the CFT Bootstrap bound. We also investigate the dynamical properties of the spin excitation and find the spin gap open near the first transition and close at the second. The relevance of our findings in realizing deconfined quantum criticality in fermion systems and the implication to lattice models with further extended interactions such as those in quantum Moiré systems, are discussed. △ Less

Submitted 4 August, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 6+2 pages, 5+2 figures

Journal ref: Phys. Rev. B 106,075111(2022)

Phase Diagram of Triangular Lattice Quantum Ising Model under External Field

Authors: Yuan Da Liao, Han Li, Zheng Yan, Hao-Tian Wei, Wei Li, Yang Qi, Zi Yang Meng

Abstract: Quantum Ising model on a triangular lattice hosts a finite temperature Berezinskii-Kosterlitz-Thouless (BKT) phase with emergent U(1) symmetry, and it will transit into an up-up-down (UUD) phase with $C_3$ symmetry breaking upon an infinitesimal external field along the longitudinal direction, but the overall phase diagram spanned by the axes of external field and temperature remains opaque due to… ▽ More Quantum Ising model on a triangular lattice hosts a finite temperature Berezinskii-Kosterlitz-Thouless (BKT) phase with emergent U(1) symmetry, and it will transit into an up-up-down (UUD) phase with $C_3$ symmetry breaking upon an infinitesimal external field along the longitudinal direction, but the overall phase diagram spanned by the axes of external field and temperature remains opaque due to the lack of systematic invesitgations with controlled methodologies. By means of quantum Monte Carlo at finite temperature and ground state density matrix renormalization group simulations, we map out the phase diagram of triangular quantum Ising model. Stemming from the upper BKT temperature at zero field, we obtain the phase boundary between the UUD and paramagnetic phases with its 2D $q=3$ Potts universality at weak field and weakly first order transition at strong field. Originated from the lower BKT temperature at zero field, we analyze the low temperature phase boundary between the clock phase and the UUD phase with Ising symmetry breaking at weak fields and the quantum phase transition between the UUD and fully polarized phases at strong fields. The accurate many-body numerical results are consistent with our field theoretical analysis. The experimental relevance towards the BKT magnet TmMgGaO$_4$ and programmable quantum simulators are also discussed. △ Less

Submitted 16 March, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 10 pages, 8 figures

Journal ref: Phys. Rev. B 103, 104416 (2021)

Realization of Topological Mott Insulator in a Twisted Bilayer Graphene Lattice Model

Authors: Bin-Bin Chen, Yuan Da Liao, Ziyu Chen, Oskar Vafek, Jian Kang, Wei Li, Zi Yang Meng

Abstract: Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG s… ▽ More Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of $\pm 1$. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit. △ Less

Submitted 16 September, 2021; v1 submitted 15 November, 2020; originally announced November 2020.

Comments: 6+7 pages, 3+7 figures

Journal ref: Nature Communications 12, 5480 (2021)

Correlated Insulating Phases in the Twisted Bilayer Graphene

Authors: Yuan Da Liao, Xiao Yan Xu, Zi Yang Meng, Jian Kang

Abstract: We review analytical and numerical studies of correlated insulating states in twisted bilayer graphene, focusing on real-space lattice models constructions and their unbiased quantum many-body solutions. We show that by constructing localized Wannier states for the narrow bands, the projected Coulomb interactions can be approximated by interactions of cluster charges with assisted nearest neighbor… ▽ More We review analytical and numerical studies of correlated insulating states in twisted bilayer graphene, focusing on real-space lattice models constructions and their unbiased quantum many-body solutions. We show that by constructing localized Wannier states for the narrow bands, the projected Coulomb interactions can be approximated by interactions of cluster charges with assisted nearest neighbor hopping terms. With the interaction part only, the Hamiltonian is $SU(4)$ symmetric considering both spin and valley degrees of freedom. In the strong coupling limit where the kinetic terms are neglected, the ground states are found to be in the $SU(4)$ manifold with degeneracy. The kinetic terms, treated as perturbation, break this large $SU(4)$ symmetry and propel the appearance of intervalley coherent state, quantum topological insulators and other symmetry-breaking insulating states. We first present the theoretical analysis of moiré lattice model construction and then show how to solve the model with large-scale quantum Monte Carlo simulations in an unbiased manner. We further provide potential directions such that from the real-space model construction and its quantum many-body solutions how the perplexing yet exciting experimental discoveries in the correlation physics of twisted bilayer graphene can be gradually understood. This review will be helpful for the readers to grasp the fast growing field of the model study of twisted bilayer graphene. △ Less

Submitted 9 January, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 15 pages, 12 figures, invited review article

Journal ref: Chin. Phys. B, 2021, Vol. 30(1): 017305

Correlation-induced insulating topological phases at charge neutrality in twisted bilayer graphene

Authors: Yuan Da Liao, Jian Kang, Clara N. Breiø, Xiao Yan Xu, Han-Qing Wu, Brian M. Andersen, Rafael M. Fernandes, Zi Yang Meng

Abstract: Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlation… ▽ More Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlations is fundamental to shed light on the nature of the superconducting and correlated insulating states observed in the TBG experiments. Here, we use unbiased, sign-problem-free quantum Monte Carlo simulations to solve an effective interacting lattice model for TBG at charge neutrality. Besides the usual cluster Hubbard-like repulsion, this model also contains an assisted hopping interaction that emerges due to the non-trivial topological properties of TBG. Such a non-local interaction fundamentally alters the phase diagram at charge neutrality, gapping the Dirac cones even for infinitesimally small interaction. As the interaction strength increases, a sequence of different correlated insulating phases emerge, including a quantum valley Hall state with topological edge states, an intervalley-coherent insulator, and a valence bond solid. The charge-neutrality correlated insulating phases discovered here provide the sought-after reference states needed for a comprehensive understanding of the insulating states at integer fillings and the proximate superconducting states of TBG. △ Less

Submitted 23 January, 2021; v1 submitted 26 April, 2020; originally announced April 2020.

Comments: 15 pages, 9 figures, 2 tables

Journal ref: Phys. Rev. X 11, 011014 (2021)

Kosterlitz-Thouless Melting of Magnetic Order in the Triangular Quantum Ising Material TmMgGaO$_4$

Authors: Han Li, Yuan Da Liao, Bin-Bin Chen, Xu-Tao Zeng, Xian-Lei Sheng, Yang Qi, Zi Yang Meng, Wei Li

Abstract: Frustrated magnets host the promises of material realizations of new paradigm of quantum matter, while direct comparison of unbiased model calculations with experimental measurements is still very challenging. Here, we design and implement a protocol of employing many-body computation methodologies for accurate model calculation -- both equilibrium and dynamical properties -- of a frustrated rare-… ▽ More Frustrated magnets host the promises of material realizations of new paradigm of quantum matter, while direct comparison of unbiased model calculations with experimental measurements is still very challenging. Here, we design and implement a protocol of employing many-body computation methodologies for accurate model calculation -- both equilibrium and dynamical properties -- of a frustrated rare-earth magnet TmMgGaO$_4$ (TMGO), which perfectly explains the corresponding experimental findings. Our results confirm TMGO is an ideal realization of triangular-lattice Ising model with an intrinsic transverse field. The magnetic order of TMGO is predicted to melt through two successive Kosterlitz-Thouless (KT) phase transitions, with a floating KT phase in between. The dynamical spectra calculated suggest remnant images of a vanishing magnetic stripe order that represent vortex-antivortex pairs, resembling rotons in a superfluid helium film. TMGO therefore constitutes a rare quantum magnet for realizing KT physics and we further propose experimental detections of its intriguing properties. △ Less

Submitted 1 March, 2020; v1 submitted 18 July, 2019; originally announced July 2019.

Comments: 18 pages, 5 figures, supplementary information

Journal ref: Nature Communications 11, 1111 (2020)

Is magic-angle twisted bilayer graphene near a quantum critical point?

Authors: Yuan Da Liao, Zi Yang Meng, Xiao Yan Xu

Abstract: An extended Hubbard model on a honeycomb lattice with two orbitals per site at charge neutrality is investigated with unbiased large-scale quantum Monte Carlo simulations. The Fermi velocity of the Dirac fermions is renormalized as the cluster charge interaction increases, until a mass term emerges and a quantum phase transition from Dirac semi-metal to valence bond solid (VBS) insulator is establ… ▽ More An extended Hubbard model on a honeycomb lattice with two orbitals per site at charge neutrality is investigated with unbiased large-scale quantum Monte Carlo simulations. The Fermi velocity of the Dirac fermions is renormalized as the cluster charge interaction increases, until a mass term emerges and a quantum phase transition from Dirac semi-metal to valence bond solid (VBS) insulator is established. The quantum critical point is discovered to belong to 3D $N=4$ Gross-Neveu chiral XY universality with the critical exponents obtained at high precision. Further enhancement of the interaction drives the system into two different VBS phases, the properties and transition between them are also revealed. Since the model is related to magic-angle twisted bilayer graphene, our results may have relevance towards the symmetry breaking order at the charge neutrality point of the material, and associate the wide range of universal strange metal behavior around it with quantum critical fluctuations. △ Less

Submitted 9 October, 2019; v1 submitted 31 January, 2019; originally announced January 2019.

Comments: 5+5 pages, 5+3 figures, 0+1 table

Journal ref: Phys. Rev. Lett. 123, 157601 (2019)