Source author record

Han-Qing Wu

Han-Qing Wu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

14 published item(s)

preprint2026arXiv

Ground-state phase diagram of S = 1/2 Heisenberg model on 2D square-hexagon-octagon lattice

Using stochastic series expansion quantum Monte Carlo and density matrix renormalization group methods, we investigate the ground-state phase diagram of the $S=1/2$ Heisenberg model on the two-dimensional square-hexagon-octagon (SHO) lattice. The model incorporates nearest-neighbor interactions $J_1$ (intrahexagon interaction) and $J_2$ (interhexagon), as well as a selected third-neighbor interaction $J_3$ along the $x$ direction. We identify five distinct phases in the parameter regime $0<λ_1=J_2/J_1<4, 0<λ_2=J_3/J_1<4$: a Néel antiferromagentic phase, two dimer phases (orthogonal and ladder staggered dimers), a hexagon singlet phase, and notably a Haldane-like symmetry-protected topological (SPT) phase. The topological nature of the Haldane-like phase is confirmed by the degeneracy of the ground-state energy under open boundary conditions and the twofold degeneracy of the entanglement spectrum. Phase boundaries are accurately determined using finite-size scaling of the spin stiffness and Binder cumulant. Data collapse analysis reveals that all transitions from nonmagnetic phases to the antiferromagnetic phase belong to the three-dimensional $O(3)$ Heisenberg universality class. In addition, we investigate the robustness of the SPT phase to other interactions, such as those that could arise in experimental materials. Our work establishes a comprehensive theoretical framework for understanding magnetic and topological phases on the SHO lattice.

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.

preprint2023arXiv

Characterizing Random-singlet State in Two-dimensional Frustrated Quantum Magnets and Implications for the Double Perovskite Sr$_2$CuTe$_{1-x}$W$_{x}$O$_6$

Motivated by experimental observation of the non-magnetic phase in the compounds with frustration and disorder, we study the ground state of the spin-$1/2$ square-lattice Heisenberg model with randomly distributed nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ couplings. By using the density matrix renormalization group (DMRG) calculation on cylinder system with circumference up to $10$ lattice sites, we identify a disordered phase between the Néel and stripe magnetic phase with growing $J_2 / J_1$ in the presence of strong bond randomness. The vanished spin-freezing parameter indicates the absence of spin glass order. The large-scale DMRG results unveil the size-scaling behaviors of the spin-freezing parameter, the power-law decay of the average spin correlation, and the exponential decay of the typical spin correlation, which all agree with the corresponding behavior in the one-dimensional random singlet (RS) state and characterize the RS nature of this disordered phase. The DMRG simulation also provides insights and opportunities for characterizing a class of non-magnetic states in two-dimensional frustrated magnets with disorder. We also compare with existing experiments and suggest more measurements for understanding the spin-liquid-like behaviors in the double perovskite Sr$_2$CuTe$_{1-x}$W$_{x}$O$_6$.

preprint2022arXiv

Fractional and composite excitations of antiferromagnetic quantum spin trimer chains

Using quantum Monte Carlo, exact diagonalization and perturbation theory, we study the spectrum of the $S=1/2$ antiferromagnetic Heisenberg trimer chain by varying the ratio $g=J_2/J_1$ of the intertrimer and intratrimer coupling strengths. The doublet ground states of trimers form effective interacting $S=1/2$ degrees of freedom described by a Heisenberg chain. Therefore, the conventional two-spinon continuum of width $\propto J_1$ when $g=1$ evolves into to a similar continuum of width $\propto J_2$ when $g\to 0$. The intermediate-energy and high-energy modes are termed \emph{doublons} and \emph{quartons} which fractionalize with increasing $g$ to form the conventional spinon continuum. In particular, at $g \approx 0.716$, the gap between the low-energy spinon branch and the high-energy band with mixed doublons, quartons, and spinons closes. These features should be observable in inelastic neutron scattering experiments if a quasi-one-dimensional quantum magnet with the linear trimer structure and $J_2<J_1$ can be identified. Our results may open a window for exploring the high-energy fractional excitations.

preprint2022arXiv

Ground States and Dynamical Properties of $S>1/2$ Quantum Heisenberg Model on the 1/5-Depleted Square Lattice

We study the $S>1/2$ antiferromagnetic Heisenberg model on the 1/5-depleted square lattice as a function of the ratio of the intra-plaquette coupling to the inter-plaquette coupling. Using stochastic series expansion quantum Monte Carlo simulations, we numerically identify three quantum phases, including the dimer phase, Néel phase and plaquette valence bond solid phase. We also obtain the accurate quantum critical points that belong to the O(3) universality class using the large-scale finite-size scaling. Most importantly, we study the dynamic spin structure factors of different phases, which can be measured by inelastic neutron scattering experiments. The low-energy excitations can be explained as triplons in the dimer phase and plaquette valence bond solid phase. While in the Néel phase, the more prominent magnon mode can be found as the spin magnitude increases. Furthermore, we find a broader continuum at smaller $S$, which may be the dynamical signature of nearly deconfined spinon excitations.

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.

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.

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.

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.

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.

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.

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.

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.

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.