We have made significant progress in applying the sign-problem free QMC simulations to explore various fundamental problems of strong correlation problems, including interacting topological insulators, doped Mott insulators, itinerant ferromagnetism and Curie-Weiss metals, novel quantum magnetism with SU(N) symmetries.
2D topological insulators
We performed one of the two earliest QMC simulations for topological
insulators and the helical edge states. It was based on the Kane-Mele
Hubbard model and we proved that it is sign-problem free.
See
the helical edge liquid and [Ref 1.]
Curie-Weiss metal and magnetic criticality
We performed, to our knowledge, the first sign-problem free QMC
simulations on itinerant ferromagnetism and the Curie-Weiss metal
state. Magnetic criticality is explored.
See
Cuire-Weiss metal and magnetic criticality and
[Ref 2.]
Doping an Antiferromagnetic Mott Insulator
The emergence of superconductivity in the vicinity of antiferromagnetism
is a universal phenomenon frequently observed in strongly correlated
materials.
We study the transition from an antiferromagnetic insulator to a
superconductor driven by hole-doping via the projector QMC simulations
[Ref. 3] , which are sign-problem-free both at
and away from the half-filling.
It is based on a bilayer lattice model with the onsite Hubbard
interaction and intra-rung charge and spin exchange interactions.
An Ising anisotropic antiferromagnetic Mott insulating phase occurs
at half-filling, which is weakened by hole doping.
Below a critical doping value xc=0.11, the antiferromagnetism coexists
with the singlet superconductivity, which is an intra-rung pairing
with an extended s-wave symmetry.
As further increasing doping, the antiferromagnetic order vanishes,
leaving only a pure superconducting phase.
These results provide important information how superconductivity
appears upon doping the parent Mott insulating state.
Inter-layer current phase
The interlayer staggered current phase Strongly interacting
systems have been conjectured to spontaneously develop current
carrying ground states under certain conditions. However, most results in
2D are based on mean field approximations whose validity is hard to justify.
In contrast, we (with S. Capponi and S. C. Zhang) performed the QMC
simulation to an extended bilayer Hubbard model without the sign problem,
finding a commensurate staggered interlayer current phase
[Ref. 4] .
To our knowledge, it is the first work to conclusively demonstrate
the existence of such a phase in 2D.
Novel quantum magnetism of SU(2N) fermions -- Slater v.s. Mott, Dirac
fermions, new type of quantum phase transition
There exist two basic pictures for interacting insulators --
the weak coupling Slater mechanism based on Fermi surface nesting, and
the strong coupling Mott one based on local moments and super-exchange.
In the SU(2) case, the evolution between these two limits is smooth
without a quantum phase transition.
How about in the SU(2N) case? We have found fundamentally different features by the sign-problem free QMC based on SU(2N) Hubbard models on the square lattice [Ref. 5] . The AFM order appears in weak and intermediate coupling regions agreeing with the Slater insulator picture. However, as further increasing interaction, it becomes suppressed. In the SU(6) case, the AFM is completely suppressed in the strong interacting region, i.e., the Mott insulator, where the quantum paramagnetic columnar dimer state appears. This shows a clear quantum phase transition between the Slater and Mott insulating regions.
We investigated the Mott-insulating state of the SU(2N) Dirac fermions via the sign-problem free QMC simulations [Ref. 6 ] . It is based on the SU(2N) Hubbard model on a honeycomb lattice. Unlike the SU(2) case, the SU(4) and SU(6) Mott-insulating phases are identified with the columnar valence-bond-solid (cVBS) order in the absence of the Neel ordering. Although the Dirac semimetal-to-cVBS transitions are typically first order due to a cubic invariance possessed by the cVBS order, the coupling to gapless Dirac fermions can soften these transitions to the second order at zero temperature.
Detect edge degeneracy in interacting topological insulators using entanglement entropy
The existance of degenerate or gapless edge states is a characteristic
feature of topological insulators, but is difficult to detect in the
presence of interactions.
We propose a new method to obtain the degeneracy of the edge states
from the perspective of entanglement entropy
[Ref. 7] , which is very useful
to identify interacting topological states.
Employing the determinant quantum Monte Carlo technique, we investigate
the interaction effect on two representative models of fermionic
topological insulators in one and two dimensions, respectively.
In the two topologically nontrivial phases, the edge degeneracies
are reduced by interactions but remain to be nontrivial.
References and Talks
Phys. Rev. Lett. 112, 156403 (2014), See pdf file
Phys. Rev. B 93, 245157 (2016) . See pdf file
Phys. Rev. B 91, 115118 (2015) . See pdf file

Last modified: July 15, 2007.