BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zi Yang Meng

Zi Yang Meng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.stat-mechcond-mat.supr-concond-mat.mes-hallhep-latquant-phcond-mat.quant-gashep-thcond-mat.dis-nnphysics.comp-ph

Catalog footprint

What is connected

59works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2026arXiv

Chiral Graviton Modes in Fermionic Fractional Chern Insulators

Chiral graviton modes are hallmark collective excitations of Fractional Quantum Hall (FQH) liquids. However, their existence on the lattice, where continuum symmetries that protect them from decay are lost, is still an open and urgent question, especially considering the recent advances in the realization of Fractional Chern Insulators (FCI) in transition metal dichalcogenides and rhombohedral pentalayer graphene. Here we present a comprehensive theoretical and numerical study of graviton-modes in fermionic FCI, and thoroughly demonstrate their existence. We first derive a lattice stress tensor operator in the context of the fermionic Harper-Hofstadter(HH) model which captures the graviton in the flat band limit. Importantly, we discover that such lattice stress-tensor operators are deeply connected to lattice quadrupolar density correlators, readily generalizable to generic Chern bands. We then explicitly show the adiabatic connection between FQH and FCI chiral graviton modes by interpolating from a low flux HH model to a Checkerboard lattice model that hosts a topological flat band. In particular, using state-of-the-art matrix product state and exact diagonalization simulations, we provide strong evidence that chiral graviton modes are long-lived excitations in FCIs despite the lack of continuous symmetries and the scattering with a two-magnetoroton continuum. By means of a careful finite-size analysis, we show that the lattice generates a finite but small intrinsic decay rate for the graviton mode. We discuss the relevance of our results for the exploration of graviton modes in FCI phases realized in solid state settings, as well as cold atom experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Discrete Shift and Polarization from Response to Symmetry Defects in Interacting Topological Phases

We extend the previous study of extracting crystalline symmetry-protected topological invariants to the correlated regime. We construct the interacting Hofstadter model defined on square lattice with the rotation and translation symmetry defects: disclination and dislocation. The model realizes Chern insulator and the charge density wave state as one tunes interactions. Employing the density matrix renormalization group (DMRG) method, we calculate the excess charge around the defects and find that the topological invariants remain quantized in both phases, with the topological quantity extracted to great precision. This study paves the way for utilizing matrix product state, and potentially other quantum many-body computation methods, to efficiently study crystalline symmetry defects on 2D interacting lattice systems.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Precise computation of universal corner entanglement entropy at 2+1 dimension: From Ising to Gaussian quantum critical points

Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Spectra of Magnetoroton and Chiral Graviton Modes of Fractional Chern Insulator

Employing the state-of-the-art time-dependent variational principle (TDVP) algorithm, we compute the spectra of charge-neutral excitations in the $ν=1/2$ (bosonic) \updated{ and $1/3$ (fermionic) fractional Chern insulator (FCI)} on the Haldane honeycomb lattice model. The magnetoroton visualized from the dynamic density structure factor acquires a minimum gap at finite momentum that can go soft with increasing interaction and give rise to a charge density wave (CDW) at the same wavevector. As the system approaches the FCI-to-CDW transition point, we observe a pronounced sharpening of the roton mode, suggesting that the magnetoroton behaves more like a quasiparticle as it softens. Notably, this occurs while the single-particle gap remains finite. Besides the magnetoroton at finite momentum, we also construct quadrupolar chiral operators in a discrete lattice and resolve the chiral graviton mode around the $Γ$ point of the Brillouin zone. Furthermore, we show the different chiralities of the gravitons of FCIs with opposite-sign Hall conductance for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamical properties of collective excitations in twisted bilayer Graphene

Employing the recently developed momentum-space quantum Monte Carlo scheme, we study the dynamic response of single-particle and collective excitations in realistic continuum models of twisted bilayer graphene. At charge neutrality, this unbiased numerical method reveals strong competition between different symmetry breaking channels with a leading instability towards the intervalley coherent state. Single-particle spectra indicate that repulsive interactions push the fermion spectral weight away from the Fermi energy and open up an insulating gap. The spectra of collective excitations suggest an approximate valley $SU(2)$ symmetry. At low-energy, long-lived valley waves are observed, which resemble spin waves of Heisenberg ferromagnetism. At high-energy, these sharp modes quickly become over-damped, when their energy reaches the fermion particle-hole continuum.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evolution of Dynamical Signature in the X-cube Fracton Topological Order

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exciton Proliferation and Fate of the Topological Mott Insulator in a Twisted Bilayer Graphene Lattice Model

Topological Mott insulator (TMI) with spontaneous time-reversal symmetry breaking and nonzero Chern number has been discovered in a real-space effective model for twisted bilayer graphene (TBG) at 3/4 filling in the strong coupling limit. However, the finite temperature properties of such a TMI state remain illusive. In this work, employing the state-of-the-art thermal tensor network and the perturbative field-theoretical approaches, we obtain the finite-$T$ phase diagram and the dynamical properties of the TBG model. The phase diagram includes the quantum anomalous Hall and charge density wave phases at low $T$, and a Ising transition separating them from the high-$T$ symmetric phases. Due to the proliferation of excitons -- particle-hole bound states -- the transitions take place at a significantly reduced temperature than the mean-field estimation. The exciton phase is accompanied with distinctive experimental signatures in such as charge compressibilities and optical conductivities close to the transition. Our work explains the smearing of the many-electron state topology by proliferating excitons and opens the avenue for controlled many-body investigations on finite-temperature states in the TBG and other quantum moiré systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fermion sign bounds theory in quantum Monte Carlo simulation

Sign problem in fermion quantum Monte Carlo (QMC) simulation appears to be an extremely hard problem. Traditional lore passing around for years tells people that when there is a sign problem, the average sign in QMC simulation approaches zero exponentially fast with the space-time volume of the configurational space. We, however, analytically show this is not always the case and manage to find physical bounds for the average sign. Our understanding is based on a direct connection between the sign bounds and a well-defined partition function of reference system and could distinguish when the bounds have the usual exponential scaling, and when they are bestowed on an algebraic scaling at low temperature limit. We analytically explain such algebraic sign problems found in flat band moiré lattice models at low temperature limit. At finite temperature, a domain size argument based on sign bounds also explains the connection between sign behavior and finite temperature phase transition. Sign bounds, as a well-defined observable, may have ability to ease or even make use of the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Height-conserving quantum dimer models

We propose a height-conserving quantum dimer model (hQDM) such that the lattice sum of its associated height field is conserved, and that it admits a Rokhsar-Kivelson (RK) point. The hQDM with minimal interaction range on the square lattice exhibits Hilbert space fragmentation and maps exactly to the XXZ spin model on the square lattice in certain Krylov subspaces. We obtain the ground-state phase diagram of hQDM via quantum Monte Carlo simulations, and demonstrate that a large portion of it is within the Krylov subspaces which admit the exact mapping to the XXZ model, with dimer ordered phases corresponding to easy-axis and easy-plane spin orders. At the RK point, the apparent dynamical exponents obtained from the single mode approximation and the height correlation function show drastically different behavior across the Krylov subspaces, exemplifying Hilbert space fragmentation and emergent glassy phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult $(2+1)d$ quantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code $\Z_2$ topological ordered state and the Kagome $\Z_2$ quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars and minimal computational costs. Our method therefore establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

The origin of the pseudogap behavior, found in many high-$T_c$ superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a "gap-filling" rather than "gap-closing" behavior, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nematic Quantum Criticality in Dirac Systems

We investigate nematic quantum phase transitions in two different Dirac fermion models. The models feature twofold and fourfold, respectively, lattice rotational symmetries that are spontaneously broken in the ordered phase. Using negative-sign-free quantum Monte Carlo simulations and an $ε$-expansion renormalization group analysis, we show that both models exhibit continuous phase transitions. In contrast to generic Gross-Neveu dynamical mass generation, the quantum critical regime is characterized by large velocity anisotropies, with fixed-point values being approached very slowly. Hence both experimental and numerical investigations will not be representative of the infrared fixed point, but of a crossover regime characterized by drifting exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Network-Initialized Monte Carlo Based on Generative Neural Networks

We design generative neural networks that generate Monte Carlo configurations with complete absence of autocorrelation from which only short Markov chains are needed before making measurements for physical observables, irrespective of the system locating at the classical critical point, fermionic Mott insulator, Dirac semimetal, or quantum critical point. We further propose a network-initialized Monte Carlo scheme based on such neural networks, which provides independent samplings and can accelerate the Monte Carlo simulations by significantly reducing the thermalization process. We demonstrate the performance of our approach on the two-dimensional Ising and fermion Hubbard models, and expect it can systematically speed up the Monte Carlo simulations especially for the very challenging many-electron problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Possible Dirac quantum spin liquid in a kagome quantum antiferromagnet YCu$_3$(OH)$_6$Br$_2$[Br$_{x}$(OH)$_{1-x}$]

We studied the magnetic properties of YCu$_3$(OH)$_6$Br$_2$[Br$_{1-x}$(OH)$_{x}$] ($x$ = 0.33), where Cu$^{2+}$ ions form two-dimensional kagome layers. There is no magnetic order down to 50 mK while the Curie-Weiss temperature is on the order of -100 K. At zero magnetic field, the low-temperature specific heat shows a $T^2$ dependence. Above 2 T, a linear temperature dependence term in specific heat emerges, and the value of $γ= C/T$ increases linearly with the field. Furthermore, the magnetic susceptibility tends to a constant value at $T = 0$. Our results suggest that the magnetic ground state of YCu$_3$(OH)$_6$Br$_2$[Br$_{1-x}$(OH)$_{x}$] is consistent with a Dirac quantum-spin-liquid state with a linearly dispersing spinon strongly coupled to an emergent gauge field, which has long been theoretically proposed as a candidate ground state in the two-dimensional kagome Heisenberg antiferromagnetic system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Simulation of Lattice Gauge Theories on Superconducting Circuits: Quantum Phase Transition and Quench Dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scaling of entanglement entropy at deconfined quantum criticality

We develop a nonequilibrium increment method to compute the Rényi entanglement entropy and investigate its scaling behavior at the deconfined critical (DQC) point via large-scale quantum Monte Carlo simulations. To benchmark the method, we first show that at an conformally-invariant critical point of O(3) transition, the entanglement entropy exhibits universal scaling behavior of area law with logarithmic corner corrections and the obtained correction exponent represents the current central charge of the critical theory. Then we move on to the deconfined quantum critical point, where although we still observe similar scaling behavior but with a very different exponent. Namely, the corner correction exponent is found to be negative. Such a negative exponent is in sharp contrast with positivity condition of the Rényi entanglement entropy, which holds for unitary conformal field theories(CFT). Our results unambiguously reveal fundamental differences between DQC and quantum critical points(QCPs) described by unitary CFTs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solving quantum rotor model with different Monte Carlo techniques

We systematically test the performance of several Monte Carlo update schemes for the $(2+1)$d XY phase transition of quantum rotor model. By comparing the local Metropolis (LM), LM plus over-relaxation (OR), Wolff-cluster (WC), hybrid Monte Carlo (HM), hybrid Monte Carlo with Fourier acceleration (FA) scheme, it is clear that among the five different update schemes, at the quantum critical point, the WC and FA schemes acquire the smallest autocorrelation time and cost the least amount of CPU hours in achieving the same level of relative error, and FA enjoys a further advantage of easily implementable for more complicated interactions such as the long-range ones. These results bestow one with the necessary knowledge of extending the quantum rotor model, which plays the role of ferromagnetic/antiferromagnetic critical bosons or Z$_2$ topological order, to more realistic and yet challenging models such as Fermi surface Yukawa-coupled to quantum rotor models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The dynamical exponent of a quantum critical itinerant ferromagnet: a Monte Carlo study

We consider the effect of the coupling between 2D quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions.A common belief is that near a $q=0$ ferromagnetic transition, fermions induce an $Ω/q$ Landau damping of rotors (i.e., the dynamical critical exponent is $z=3$) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy $Σ\propto ω^{2/3}$. This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies.Here we show that for the XY case the behavior is different.We report the results of large scale quantum Monte Carlo simulations,which show that at small frequencies $z=2$ and $Σ\propto ω^{1/2}$. We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces $1/q$ in the Landau damping by a constant. We discuss the implication of these results to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermodynamic characteristic for correlated flat-band system with quantum anomalous Hall ground state

While the ground state phase diagram of the correlated flat-band systems have been intensively investigated, the dynamic and thermodynamic properties of such lattice models are less explored, but it is the latter which is most relevant to the experimental probes (transport, quantum capacitance and spectroscopy) of the quantum moiré materials such as twisted bilayer graphene and transition metal dichalcogenides. Here we show, by means of momentum-space quantum Monte Carlo and exact diagonalization, there exists a unique thermodynamic characteristic for the correlated flat-band models with interaction-driven quantum anomalous Hall (QAH) ground state, namely, the transition from the QAH insulator to the metallic state takes place at a much lower temperature compared with the zero-temperature single-particle gap generated by the long-range Coulomb interaction. Such low transition temperature comes from the proliferation of excitonic particle-hole excitations, which "quantum teleport" the electrons across the gap between different topological bands to restore the broken time-reversal symmetry and give rise to a pronounced enhancement in the charge compressibility. Future experiments, to verify such generic thermodynamic characteristics, are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Topological Disorder Parameter

We introduce a many-body topological invariant, called the topological disorder parameter (TDP), to characterize gapped quantum phases with global internal symmetry in (2+1)d. TDP is defined as the constant correction that appears in the ground state expectation value of a partial symmetry transformation applied to a connected spatial region $M$, the absolute value of which scales generically as $\exp(-αl+γ)$ where $l$ is the perimeter of $M$ and $γ$ is the TDP. Motivated by a topological quantum field theory interpretation of the operator, we show that $e^γ$ can be related to the quantum dimension of the symmetry defect, and provide a general formula for $γ$ when the entanglement Hamiltonian of the topological phase can be described by a (1+1)d conformal field theory (CFT). A special case of TDP is equivalent to the topological Rényi entanglement entropy when the symmetry is the cyclic permutation of the replica of the gapped phase. We then investigate several examples of lattice models of topological phases, both analytically and numerically, in particular when the assumption of having a CFT edge theory is not satisfied. We also consider an example of partial translation symmetry in Wen's plaquette model and show that the result can be understood using the edge CFT. Our results establish a new tool to detect quantum topological order.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlated Insulating Phases in the Twisted Bilayer Graphene

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fermion enhanced first-order phase transition and chiral Gross-Neveu tricritical point

The fluctuations of massless Dirac fermion can not only turn a first-order bosonic phase transition (in the Landau sense) to a quantum critical point, but also work reversely to enhance the first-order transition itself, depending on the implementation of finite size effects in the coupling corrections. Here, we report a case study of the latter by employing quantum Monte Carlo simulation upon a lattice model in which the bosonic part featuring the Landau-Devonshire first-order phase transition and Yukawa coupled to the Dirac fermions. We find that the parameter range for the first-order phase transition becomes larger as the Yukawa coupling increases and the microscopic mechanism of this phenomena is revealed, at a quantitative level, as the interplay between the critical fluctuations and the finite-size effects. Moreover, the scaling behavior at the separation point between the first-order and the continuous phase transitions is found to belong to the chiral tricritical Gross-Neveu universality. Our result demonstrates that the interplay of massless Dirac fermions, critical fluctuations and the finite size effects could trigger a plethora of interesting phenomena and therefore great care is called for when making generalizations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Vestigial anyon condensation in kagome quantum spin liquids

We construct a lattice model of topological order (kagome quantum spin liquids) and solve it with unbiased quantum Monte Carlo simulations. A three-stage anyon condensation with two transitions from a $\mathbb Z_2\boxtimes\mathbb Z_2$ topological order to a $\mathbb Z_2$ topological order and eventually to a trivial symmetric phase is revealed. These results provide concrete examples of phase transitions between topological orders in quantum magnets. The designed quantum spin liquid model and its numerical solution offer a playground for further investigations on vestigial anyon condensation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Confinement transition in the QED$_3$-Gross-Neveu-XY universality class

The coupling between fermionic matter and gauge fields plays a fundamental role in our understanding of nature, while at the same time posing a challenging problem for theoretical modeling. In this situation, controlled information can be gained by combining different complementary approaches. Here, we study a confinement transition in a system of $N_f$ flavors of interacting Dirac fermions charged under a U(1) gauge field in 2+1 dimensions. Using Quantum Monte Carlo simulations, we investigate a lattice model that exhibits a continuous transition at zero temperature between a gapless deconfined phase, described by three-dimensional quantum electrodynamics, and a gapped confined phase, in which the system develops valence-bond-solid order. We argue that the quantum critical point is in the universality class of the QED$_3$-Gross-Neveu-XY model. We study this field theory within a $1/N_f$ expansion in fixed dimension as well as a renormalization group analysis in $4-ε$ space-time dimensions. The consistency between numerical and analytical results is revealed from large to intermediate flavor number.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Correlated states in twisted double bilayer graphene

Electron-electron interactions play an important role in graphene and related systems and can induce exotic quantum states, especially in a stacked bilayer with a small twist angle. For bilayer graphene where the two layers are twisted by a "magic angle", flat band and strong many-body effects lead to correlated insulating states and superconductivity. In contrast to monolayer graphene, the band structure of untwisted bilayer graphene can be further tuned by a displacement field, providing an extra degree of freedom to control the flat band that should appear when two bilayers are stacked on top of each other. Here, we report the discovery and characterization of such displacement-field tunable electronic phases in twisted double bilayer graphene. We observe insulating states at a half-filled conduction band in an intermediate range of displacement fields. Furthermore, the resistance gap in the correlated insulator increases with respect to the in-plane magnetic fields and we find that the g factor according to spin Zeeman effect is ~2, indicating spin polarization at half filling. These results establish the twisted double bilayer graphene as an easily tunable platform for exploring quantum many-body states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Designer Monte Carlo Simulation for Gross-Neveu Transition

In this manuscript, we study quantum criticality of Dirac fermions via large-scale numerical simulations, focusing on the Gross-Neveu-Yukawa(GNY) chiral-Ising quantum critical point with critical bosonic modes coupled with Dirac fermions. We show that finite-size effects at this quantum critical point can be efficiently minimized via model design, which maximizes the ultraviolet cutoff and at the same time places the bare control parameters closer to the nontrivial fixed point to better expose the critical region. Combined with the efficient self-learning quantum Monte Carlo algorithm, which enables non-local update of the bosonic field, we find that moderately-large system size (up to $16\times 16$) is already sufficient to produce robust scaling behavior and critical exponents.The conductance of the Dirac fermions is also calculated and its frequency dependence is found to be consistent with the scaling behavior predicted by the conformal field theory. The methods and model-design principles developed for this study can be generalized to other fermionic QCPs, and thus provide a promising direction for controlled studies of strongly-correlated itinerant systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data

Quantum Monte Carlo (QMC) simulations of correlated electron systems provide unbiased information about system behavior at a quantum critical point (QCP) and can verify or disprove the existing theories of non-Fermi liquid (NFL) behavior at a QCP. However, simulations are carried out at a finite temperature, where quantum-critical features are masked by finite temperature effects. Here we present a theoretical framework within which it is possible to separate thermal and quantum effects and extract the information about NFL physics at $T=0$. We demonstrate our method for a specific example of 2D fermions near a Ising-ferromagnetic QCP. We show that one can extract from QMC data the zero-temperature form of fermionic self-energy $Σ(ω)$ even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of $Σ(ω)$ agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys $ω^{2/3}$ scaling at low frequencies. Our results open up an avenue for QMC studies of quantum-critical metals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetism of Finite Graphene Samples: Mean-Field Theory compared with Exact Diagonalization and Quantum Monte Carlo Simulation

The magnetic properties of graphene on finite geometries are studied using a self-consistent mean-field theory of the Hubbard model. This approach is known to predict ferromagnetic edge states close to the zig-zag edges in single-layer graphene quantum dots and nanoribbons. In order to assess the accuracy of this method, we perform complementary exact diagonalization and quantum Monte Carlo simulations. We observe good quantitative agreement for all quantities investigated provided that the Coulomb interaction is not too strong.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Metal to Orthogonal Metal Transition

Orthogonal metal is a new quantum metallic state that conducts electricity but acquires no Fermi surface (FS) or quasiparticles, and hence orthogonal to the established paradigm of Landau's Fermi-liquid (FL). Such a state may hold the key of understanding the perplexing experimental observations of quantum metals that are beyond FL, i.e., dubbed non-Fermi-liquid (nFL), ranging from the Cu- and Fe-based oxides, heavy fermion compounds to the recently discovered twisted graphene heterostructures. However, to fully understand such an exotic state of matter, at least theoretically, one would like to construct a lattice model and to solve it with unbiased quantum many-body machinery. Here we achieve this goal by designing a 2D lattice model comprised of fermionic and bosonic matter fields coupled with dynamic $\mathbb Z_2$ gauge fields, and obtain its exact properties with sign-free quantum Monte Carlo simulations. We find that as the bosonic matter fields become disordered, with the help of deconfinement of the $\mathbb Z_2$ gauge fields, the system reacts with changing its nature from the conventional normal metal with an FS to an orthogonal metal of nFL without FS and quasiparticles and yet still responds to magnetic probe like an FL. Such a quantum phase transition from a normal metal to an orthogonal metal, with its electronic and magnetic spectral properties revealed, is calling for the establishment of new paradigm of quantum metals and their transition with conventional ones.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum Many-Body Simulations of the 2D Fermi-Hubbard Model in Ultracold Optical Lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, has been recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping, to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Phys. Rev. X 8, 031082 (2018)], complemented with independent determinant quantum Monte Carlo (DQMC) offer a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to unprecedented low temperature of few percents of the fermion tunneling energy scale. We observe excellent agreement with ultracold fermion measurements at both half-filling and finite-doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behavior of the antimoments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum phases of SrCu2(BO3)2 from high-pressure thermodynamics

We report heat capacity measurements of SrCu$_2$(BO$_3$)$_2$ under high pressure along with simulations of relevant quantum spin models and map out the $(P,T)$ phase diagram of the material. We find a first-order quantum phase transition between the low-pressure quantum dimer paramagnet and a phase with signatures of a plaquette-singlet state below T = $2$ K. At higher pressures, we observe a transition into a previously unknown antiferromagnetic state below $4$ K. Our findings can be explained within the two-dimensional Shastry-Sutherland quantum spin model supplemented by weak inter-layer couplings. The possibility to tune SrCu$_2$(BO$_3$)$_2$ between the plaquette-singlet and antiferromagnetic states opens opportunities for experimental tests of quantum field theories and lattice models involving fractionalized excitations, emergent symmetries, and gauge fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bona fide interaction-driven topological phase transition in correlated SPT states

It is expected that the interplay between non-trivial band topology and strong electron correlation will lead to very rich physics. Thus a controlled study of the competition between topology and correlation is of great interest. Here, employing large-scale quantum Monte Carlo (QMC) simulations, we provide a concrete example of the Kane-Mele-Hubbard (KMH) model on an AA stacking bilayer honeycomb lattice with inter-layer antiferromagnetic interaction. Our simulation identified several different phases: a quantum spin-Hall insulator (QSH), a $xy$-plane antiferromagnetic Mott insulator ($xy$-AFM) and an inter-layer dimer-singlet insulator (dimer-singlet). Most importantly, a bona fide topological phase transition between the QSH and the dimer-singlet insulators, purely driven by the inter-layer antiferromagnetic interaction is found. At the transition, the spin and charge gap of the system close while the single-particle excitations remain gapped, which means that this transition has no mean field analogue and it can be viewed as a transition between bosonic SPT states. At one special point, this transition is described by a $(2+1)d$ $O(4)$ nonlinear sigma model (NLSM) with {\it exact} $SO(4)$ symmetry, and a topological term at {\it exactly} $Θ= π$. Relevance of this work towards more general interacting SPT states is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Competing pairing channels in the doped honeycomb lattice Hubbard model

Proposals for superconductivity emerging from correlated electrons in the doped Hubbard model on the honeycomb lattice range from chiral $d+id$ singlet to $p+ip$ triplet pairing, depending on the considered range of doping and interaction strength, as well as the approach used to analyze the pairing instabilities. Here, we consider these scenarios using large-scale dynamic cluster approximation (DCA) calculations to examine the evolution in the leading pairing symmetry from weak to intermediate coupling strength. These calculations focus on doping levels around the van Hove singularity (VHS) and are performed using DCA simulations with an interaction-expansion continuous-time quantum Monte Carlo cluster solver. We calculated explicitly the temperature dependence of different uniform superconducting pairing susceptibilities and found a consistent picture emerging upon gradually increasing the cluster size: while at weak coupling the $d+id$ singlet pairing dominates close to the VHS filling, an enhanced tendency towards $p$-wave triplet pairing upon further increasing the interaction strength is observed. The relevance of these systematic results for existing proposals and ongoing pursuits of odd-parity topological superconductivity are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Diagnosis of interaction-driven topological phase via exact diagonalization

We propose a general scheme for diagnosing interaction-driven topological phases in the weak interaction regime using exact diagonalization (ED). The scheme comprises the analysis of eigenvalues of the point-group operators for the many-body eigenstates and the correlation functions for physical observables to extract the symmetries of the order parameters and the topological numbers of the underlying ground states at the thermodynamic limit from a relatively small size system afforded by ED. As a concrete example, we investigate the interaction effects on the half-filled spinless fermions on the checkerboard lattice with a quadratic band crossing point. Numerical results support the existence of a spontaneous quantum anomalous Hall phase purely driven by a nearest-neighbor weak repulsive interaction, separated from a nematic Mott insulator phase at strong repulsive interaction by a first-order phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effect of nematic order on the low-energy spin fluctuations in detwinned BaFe$_{1.935}$Ni$_{0.065}$As$_2$

The origin of nematic order remains one of the major debates in iron-based superconductors. In theories based on spin nematicity, one major prediction is that the spin-spin correlation length at (0,$π$) should decrease with decreasing temperature below the structural transition temperature $T_s$. Here we report inelastic neutron scattering studies on the low-energy spin fluctuations in BaFe$_{1.935}$Ni$_{0.065}$As$_2$ under uniaxial pressure. Both intensity and spin-spin correlation start to show anisotropic behavior at high temperature, while the reduction of the spin-spin correlation length at (0,$π$) happens just below $T_s$, suggesting strong effect of nematic order on low-energy spin fluctuations. Our results favor the idea that treats the spin degree of freedom as the driving force of the electronic nematic order.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum critical point of Dirac fermion mass generation without spontaneous symmetry breaking

We study a lattice model of interacting Dirac fermions in $(2+1)$ dimension space-time with an SU(4) symmetry. While increasing interaction strength, this model undergoes a {\it continuous} quantum phase transition from the weakly interacting Dirac semimetal to a fully gapped and nondegenerate phase without condensing any Dirac fermion bilinear mass operator. This unusual mechanism for mass generation is consistent with recent studies of interacting topological insulators/superconductors, and also consistent with recent progresses in lattice QCD community.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Visualizing a Bosonic Symmetry Protected Topological Phase in an Interacting Fermion Model

Symmetry protected topological (SPT) phases in free fermion and interacting bosonic systems have been classified, but the physical phenomena of interacting fermionic SPT phases have not been fully explored. Here, employing large-scale quantum Monte Carlo simulation, we investigate the edge physics of a bilayer Kane-Mele-Hubbard model with zigzag ribbon geometry. Our unbiased numerical results show that the fermion edge modes are gapped out by interaction, while the bosonic edge modes remain gapless at the $(1+1)d$ boundary, before the bulk quantum phase transition to a topologically trivial phase. Therefore, finite fermion gaps both in the bulk and on the edge, together with the robust gapless bosonic edge modes, prove that our system becomes an emergent bosonic SPT phase at low energy, which is, for the first time, directly observed in an interacting fermion lattice model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

$i$QIST: An open source continuous-time quantum Monte Carlo impurity solver toolkit

Quantum impurity solvers have a broad range of applications in theoretical studies of strongly correlated electron systems. Especially, they play a key role in dynamical mean-field theory calculations of correlated lattice models and realistic materials. Therefore, the development and implementation of efficient quantum impurity solvers is an important task. In this paper, we present an open source interacting quantum impurity solver toolkit (dubbed $i$QIST). This package contains several highly optimized quantum impurity solvers which are based on the hybridization expansion continuous-time quantum Monte Carlo algorithm, as well as some essential pre- and post-processing tools. We first introduce the basic principle of continuous-time quantum Monte Carlo algorithm and then discuss the implementation details and optimization strategies. The software framework, major features, and installation procedure for $i$QIST are also explained. Finally, several simple tutorials are presented in order to demonstrate the usage and power of $i$QIST.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coulomb liquid phases of bosonic cluster Mott insulators on a pyrochlore lattice

Employing large-scale quantum Monte Carlo simulations, we reveal the full phase diagram of the extended Hubbard model of hard-core bosons on the pyrochlore lattice with partial fillings. When the intersite repulsion is dominant, the system is in a cluster Mott insulator phase with an integer number of bosons localized inside the tetrahedral units of the pyrochlore lattice. We show that the full phase diagram contains three cluster Mott insulator phases with 1/4, 1/2, and 3/4 boson fillings, respectively. We further demonstrate that all three cluster Mott insulators are Coulomb liquid phases and its low-energy property is described by the emergent compact U(1) quantum electrodynamics. In addition to measuring the specific heat and entropy of the cluster Mott insulators, we investigate the correlation function of the emergent electric field and verify it is consistent with the compact U(1) quantum electrodynamics description. Our result sheds light on the magnetic properties of various pyrochlore systems, as well as the charge physics of the cluster magnets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for spin-triplet odd-parity superconductivity close to type-II van Hove singularities

Searching for unconventional Cooper pairing states has been at the heart of superconductivity research since the discovery of BCS superconductors. In particular, spin-triplet odd-parity pairing states were recently revisited due to the possibility of tuning towards topological superconductors. In this context, it is interesting to note a recent proposal that such a spin-triplet pairing instability occurs when the band filling is near van Hove singularities (vHS) associated with momenta away from time-reversal invariant momenta named type-II vHS. However, this result was obtained within a weak coupling renormalization group with Fermi surface patch approximation. To explore superconducting instabilities beyond this weak coupling Fermi surface patch approximation, we perform systematic study on Hubbard model in a two-dimensional square lattice using three different methods: random phase approximation, large-scale dynamical mean field theory simulations with continuous time quantum Monte Carlo (CTQMC) impurity solver, and large-scale dynamical cluster simulations with CTQMC cluster solver. We find, in a wide doping range centered around the type-II van Hove filling, a two-fold degenerate, spin-triplet, odd-parity $p$-wave pairing state emerges due to repulsive interaction, when the Fermi surface is not sufficiently nested. Possible relevance of our findings to the recently discovered superconductors LaO$_{1-x}$F$_{x}$BiS$_{2}$, Ir$_{1-x}$Pt$_{x}$Te$_{2}$ and proposed doped BC$_{3}$ are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mott insulating states and quantum phase transitions of correlated SU(2N) Dirac fermions

The interplay between charge and spin degrees of freedom in strongly correlated fermionic systems, in particular of Dirac fermions, is a long-standing problem in condensed matter physics. We investigate the competing orders in the half-filled SU(2N) Hubbard model on a honeycomb lattice, which can be accurately realized in optical lattices with large-spin ultra-cold alkaline-earth fermions. Employing large-scale projector determinant quantum Monte Carlo simulations, we have explored quantum phase transitions from the gapless Dirac semi-metals to the gapped Mott-insulating phases in the SU(4) and SU(6) cases. Both of these Mott-insulating states are found to be columnar valence bond solid (cVBS) and to be absent of the antiferromagnetic Neel ordering and the loop current ordering. Inside the cVBS phases, the dimer ordering is enhanced by increasing fermion components and behaves non-monotonically as the interaction strength increases. Although the transitions generally should be of first order due to a cubic invariance possessed by the cVBS order, the coupling to gapless Dirac fermions can soften the transitions to second order through a non-analytic term in the free energy. Our simulations provide important guidance for the experimental exploration of novel states of matter with ultra-cold alkaline earth fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mott transition in the triangular lattice Hubbard model: a dynamical cluster approximation study

Based on dynamical cluster approximation (DCA) quantum Monte Carlo simulations, we study the interaction-driven Mott metal-insulator transition (MIT) in the half-filled Hubbard model on the anisotropic two-dimensional triangular lattice, where the degree of frustration is varied between the unfrustrated case and the fully frustrated, isotropic triangular lattice. Upon increasing the DCA cluster size, we analyze the evolution of the MIT phase boundary as a function of frustration in the phase diagram spanned by the interaction strength and temperature, and provide a quantitative description of the MIT phase boundary in the triangular lattice Hubbard model. Qualitative differences in the phase boundary between the unfrustrated and fully frustrated cases are exhibited. In particular, a change in the sign of the phase boundary slope is observed, which via an impurity cluster eigenstate analysis, may be related to a change in the nature of the insulating state. We discuss our findings within the scenario that the triangular lattice electron system might exhibit a quantum critical Mott MIT with a possible quantum spin liquid insulating state, such as considered for the organic charge transfer salts $κ$-(BEDT-TTF)$_{2}$Cu$_{2}$(CN)$_{3}$ and EtMe$_{3}$Sb[Pd(dmit)$_{2}$]$_{2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiplicative logarithmic corrections to quantum criticality in three-dimensional dimerized antiferromagnets

We investigate the quantum phase transition in an $S = 1/2$ dimerized Heisenberg antiferromagnet in three spatial dimensions. By performing large-scale quantum Monte Carlo simulations and detailed finite-size scaling analyses, we obtain high-precision results for the quantum critical properties at the transition from the magnetically disordered dimer-singlet phase to the antiferromagnetically ordered Néel phase. This transition breaks O($N$) symmetry with $N = 3$ in $D = 3 + 1$ dimensions. This is the upper critical dimension, where multiplicative logarithmic corrections to the leading mean-field critical properties are expected; we extract these corrections, establishing their precise forms for both the zero-temperature staggered magnetization, $m_s$, and the Néel temperature, $T_N$. We present a scaling Ansatz for $T_N$, including logarithmic corrections, which agrees with our data and indicates exact linearity with $m_s$, implying a complete decoupling of quantum and thermal fluctuation effects even arbitrarily close to the quantum critical point. We also demonstrate the predicted $N$-independent leading and subleading logarithmic corrections in the size-dependence of the staggered magnetic susceptibility. These logarithmic scaling forms have not previously been identified or verified by unbiased numerical methods and we discuss their relevance to experimental studies of dimerized quantum antiferromagnets such as TlCuCl$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum Monte Carlo study of strange correlator in interacting topological insulators

Distinguishing the nontrivial symmetry-protected topological (SPT) phase from the trivial insulator phase in the presence of electron-electron interaction is an urgent question to the study of topological insulators, due to the fact that most of the topological indices defined for free electron systems are very likely unsuitable for interacting cases. In this work, we demonstrate that the strange correlator is a sensitive diagnosis to detect SPT states in interacting systems. Employing large-scale quantum Monte Carlo (QMC) simulations, we investigate the interaction-driven quantum phase transition in the Kane-Mele-Hubbard model. The transition from the quantum spin Hall insulator at weak interaction to an antiferromagnetic Mott insulator at strong interaction can be readily detected by the momentum space behavior of the strange correlator in single-particle, spin, and pairing sectors. The interaction effects on the symmetry-protected edge states in various sectors, i.e., the helical Luttinger liquid behavior, are well captured in the QMC measurements of strange correlators. Moreover, we demonstrate that the strange correlator is technically easier to implement in QMC and more robust in performance than other proposed numerical diagnoses for interacting topological states, as only static correlations are needed. The attempt in this work paves the way for using the strange correlator to study interaction-driven topological phase transitions in fermionic as well as bosonic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological invariants for interacting topological insulators: I. Efficient numerical evaluation scheme and implementations

The aim of this series of two papers is to discuss topological invariants for interacting topological insulators (TIs). In the first paper (I), we provide a paradigm of efficient numerical evaluation scheme for topological invariants, in which we demystify the procedures and techniques employed in calculating Z2 invariant and spin Chern number via zero-frequency single-particle Green's function in quantum Monte Carlo (QMC) simulations. Here we introduce a periodization process to overcome the ubiquitous finite-size effect, so that the calculated spin Chern number shows ideally quantized values. We also show that making use of symmetry properties of the underlying systems can greatly reduce the computational effort. To demonstrate the effectiveness of our numerical evaluation scheme, especially the periodization process, of topological invariants, we apply it on two independent two-dimensional models of interacting topological insulators. In the subsequent paper (II), we apply the scheme developed here to wider classes of models of interacting topological insulators, for which certain limitation of constructing topological invariant via single-particle Green's functions will be presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological invariants for interacting topological insulators: II. Breakdown of the Green's function formalism

Topological phase transitions in free fermion systems can be characterized by closing of single-particle gap and change in topological invariants. However, in the presence of electronic interactions, topological phase transitions are more complicated. In paper I of this series (arXiv:1510.07816), we have developed an efficient scheme to evaluate the topological invariants based on Green's function formalism. Here, in paper II, we demonstrate four interaction-drive topological phase transitions (TPTs) in two-dimensional (2D) interacting topological insulators (TIs) via large-scale quantum Monte Carlo (QMC) simulations, based on the scheme of evaluating topological invariants presented in paper I. Across these transitions, the defining symmetries of the TIs have been neither explicitly nor spontaneously broken. In the first two models, the topological invariants calculated from Green's function formalism succeed in characterizing interaction-driven TPTs. However, in the second two models, we find single-particle gap does not close and the topological invariants constructed from single-particle Green's function acquire no change across the TPTs. Unexpected breakdown of the Green's function formalism in constructing topological invariants is thus discovered. We thence classify the TPTs in interacting TIs into two categories: those have noninteracting correspondence can be characterized successfully by the topological invariants constructed from Green's functions, while for the others that do not have noninteracting correspondence, the Green's function formalism experiences a breakdown but more interesting and exciting phenomena, such as emergent collective critical modes at the transition, arise. Discussion on the success and breakdown of topological invariants constructed from the Green's function formalism in the context of symmetry protected topological (SPT) states is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Direct observation of fragile Mott insulators on plaquette Hubbard lattices

Employing extensive cellular dynamical mean-field theory (CDMFT) calculations with exact diagonalization impurity solver, we investigate the ground state phase diagrams and non-magnetic metal-insulator transitions of the half-filled Hubbard model on two plaquette -- the 1/5 depleted and checkerboard -- square lattices. We identify three different insulators in the phase diagrams: dimer insulator, antiferromagnetic insulator, and plaquette insulator. And we demonstrate that the plaquette insulator is a novel fragile Mott insulator (FMI) which features a nontrivial one-dimensional irreducible representation of the $C_{4v}$ crystalline point-group and cannot be adiabatically connected to any band insulator with time-reversal symmetry. Furthermore, we study the non-magnetic quantum phase transitions from the metal to the FMI and find that this Mott metal-insulator transition is characterized by the splitting of the non-interacting bands due to interaction effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effective Cluster Typical Medium Theory for Diagonal Anderson Disorder Model in One- and Two-Dimensions

We develop a cluster typical medium theory to study localization in disordered electronic systems. Our formalism is able to incorporate non-local correlations beyond the local typical medium theory in a systematic way. The cluster typical medium theory utilizes the momentum resolved typical density of states and hybridization function to characterize the localization transition. We apply the formalism to the Anderson model of localization in one- and two-dimensions. In one dimension, we find that the critical disorder strength scales inversely with the linear cluster size with a power-law, $W_c \sim (1/L_c)^{1/ν}$; whereas in two dimensions, the critical disorder strength decreases logarithmically with the linear cluster size. Our results are consistent with previous numerical work and in agreement with the one-parameter scaling theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Odd-parity triplet superconductivity in multi-orbital materials with strong spin-orbit coupling: applications to doped Sr2IrO4

We explore possible superconducting states in $t_{2g}$ multi-orbital correlated electron systems with strong spin-orbit coupling (SOC). In order to study such systems in a controlled manner, we employ large-scale dynamical mean-field theory (DMFT) simulations with the hybridization expansion continuous-time Quantum Monte Carlo (CTQMC) impurity solver. To determine the pairing symmetry, we go beyond the local DMFT formalism using parquet equations to introduce the momentum dependence in the two-particle vertex and correlation functions. In the strong SOC limit, a singlet, $d$-wave pairing state in the electron-doped side of the phase diagram is observed at weak Hund's coupling, which is triggered by antiferromagnetic fluctuations. When the Hund's coupling is comparable to SOC, a two-fold degenerate, triplet $p$-wave pairing state with relatively high $T_c$ emerges in the hole-doped side of the phase diagram, which is associated with enhanced charge fluctuations. Experimental implications to doped Sr$_{2}$IrO$_{4}$ are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dimerized Solids and Resonating Plaquette Order in SU(N)-Dirac Fermions

We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolution of the Superconductivity Dome in the two dimensional Hubbard Model

In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two dimensional Hubbard model. Here we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the d-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (t'), chemical potential and temperature as control parameters. To obtain the superconducting transition temperature Tc, we employ two-particle measurements of the pairing susceptibilities. As t' goes from positive to negative values, we find the d-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling (n close to 0.85) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two dimensional Hubbard model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unconventional superconductivity on the triangular lattice Hubbard model

Using large-scale dynamical cluster quantum Monte Carlo simulations, we explore the unconventional superconductivity in the hole-doped Hubbard model on the triangular lattice. Due to the interplay of electronic correlations, geometric frustration, and Fermi surface topology, we find a doubly degenerate singlet pairing state at an interaction strength close to the bare bandwidth. Such an unconventional superconducting state is mediated by antiferromagnetic spin fluctuations along the $Γ$-$K$ direction, where the Fermi surface is nested. An exact decomposition of the irreducible particle-particle vertex further confirms the dominant component of the effective pairing interaction comes from the spin channel. Our findings provide support for chiral $d +i d$ superconductivity in water-intercalated sodium cobaltates Na$_{x}$CoO$_{2} \cdot y$H$_{2}$O, as well as insight into the superconducting phases of the organic compounds $κ$-(ET)$_{2}$X and Pd(dmit)$_{2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Antiferromagnetism in the Hubbard Model on the Bernal-stacked Honeycomb Bilayer

Using a combination of quantum Monte Carlo simulations, functional renormalization group calculations and mean-field theory, we study the Hubbard model on the Bernal-stacked honeycomb bilayer at half-filling as a model system for bilayer graphene. The free bands consisting of two Fermi points with quadratic dispersions lead to a finite density of states at the Fermi level, which triggers an antiferromagnetic instability that spontaneously breaks sublattice and spin rotational symmetry once local Coulomb repulsions are introduced. Our results reveal an inhomogeneous participation of the spin moments in the ordered ground state, with enhanced moments at the three-fold coordinated sites. Furthermore, we find the antiferromagnetic ground state to be robust with respect to enhanced interlayer couplings and extended Coulomb interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lifshitz Transition in the Two Dimensional Hubbard Model

Using large-scale dynamical cluster quantum Monte Carlo simulations, we study the Lifshitz transition of the two dimensional Hubbard model with next-nearest-neighbor hopping ($t'$), chemical potential and temperature as control parameters. At $t'\le0$, we identify a line of Lifshitz transition points associated with a change of the Fermi surface topology at zero temperature. In the overdoped region, the Fermi surface is complete and electron-like; across the Lifshitz transition, the Fermi surface becomes hole-like and develops a pseudogap. At (or very close to) the Lifshitz transition points, a van Hove singularity in the density of states crosses the Fermi level. The van Hove singularity occurs at finite doping due to correlation effects, and becomes more singular when $t'$ becomes more negative. The resulting temperature dependence on the bare d-wave pairing susceptibility close to the Lifshitz points is significantly different from that found in the traditional van Hove scenarios. Such unambiguous numerical observation of the Lifshitz transition at $t'\le0$ extends our understanding of the quantum critical region in the phase diagram, and shines lights on future investigations of the nature of the quantum critical point in the two dimensional Hubbard model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Signatures of Edge-State Magnetism on Graphene Nanoribbons

We investigate the edge-state magnetism of graphene nanoribbons using projective quantum Monte Carlo simulations and a self-consistent mean-field approximation of the Hubbard model. The static magnetic correlations are found to be short ranged. Nevertheless, the correlation length increases with the width of the ribbon such that already for ribbons of moderate widths we observe a strong trend towards mean-field-type ferromagnetic correlations at a zigzag edge. These correlations are accompanied by a dominant low-energy peak in the local spectral function and we propose that this can be used to detect edge-state magnetism by scanning tunneling microscopy. The dynamic spin structure factor at the edge of a ribbon exhibits an approximately linearly dispersing collective magnonlike mode at low energies that decays into Stoner modes beyond the energy scale where it merges into the particle-hole continuum.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Phases and magnetization process of an anisotropic Shastry-Sutherland model

We examine ground state properties of the spin-1/2 easy-axis Heisenberg model on the Shastry-Sutherland lattice with ferromagnetic transverse spin exchange using quantum Monte Carlo and degenerate perturbation theory. For vanishing transverse exchange, the model reduces to an antiferromagnetic Ising model that besides Neel order harbors regions of extensive ground state degeneracy. In the quantum regime, we find a dimerized phase of triplet states, separated from the Neel ordered phase by a superfluid. The quantum phase transitions between these phases are characterized. The magnetization process shows a magnetization plateau at 1/3 of the saturation value, that persists down to the Ising limit, and a further plateau at 1/2 only in the quantum regime. For both plateaus, we determine the crystalline patterns of the localized triplet excitations. No further plateaus nor supersolid phases are found in this model.

0 citations0 reviewsNo code linkedOpen access