To reduce the rapidly growing computational cost of the dual fermion lattice calculation with increasing system size, we introduce two embedding schemes. One is the real fermion embedding, and the other is the dual fermion embedding. Our numerical te
sts show that the real fermion and dual fermion embedding approaches converge to essentially the same result. The application on the Anderson disorder and Hubbard models shows that these embedding algorithms converge more quickly with system size as compared to the conventional dual fermion method, for the calculation of both single-particle and two-particle quantities.
While the coherent potential approximation (CPA) is the prevalent method for the study of disordered electronic systems, it fails to capture non-local correlations and Anderson localization. To incorporate such effects, we extend the dual fermion app
roach to disordered non-interacting systems using the replica method. Results for single- and two- particle quantities show good agreement with cluster extensions of the CPA; moreover, weak localization is captured. As a natural extension of the CPA, our method presents an alternative to the existing cluster theories. It can be used in various applications, including the study of disordered interacting systems, or for the description of non-local effects in electronic structure calculations.
We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adi
abatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricted to the transverse direction. The transition from the topological band insulator to the antiferromagnetic Mott insulator is in the universality class of the three-dimensional XY model. The numerical data suggest that the spin liquid to topological insulator and spin liquid to Mott insulator transitions are both continuous.
At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. While such states are possibly reali
zed in two-dimensional organic compounds, they have remained elusive in experimentally relevant microscopic two-dimensional models. Here, we show by means of large-scale quantum Monte Carlo simulations of correlated fermions on the honeycomb lattice, a structure realized in graphene, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Therefore, the possibility of unconventional superconductivity through doping arises. We foresee its realization with ultra-cold atoms or with honeycomb lattices made with group IV elements.
