Efficient calculation of imaginary time displaced correlation functions in the projector auxiliary field quantum Monte-Carlo algorithm

The calculation of imaginary time displaced correlation functions with the auxiliary field projector quantum Monte-Carlo algorithm provides valuable insight (such as spin and charge gaps) in the model under consideration. One of the authors and M. Imada [F.F. Assaad and M. Imada, J. Phys. Soc. Jpn. 65 189 (1996).] have proposed a numerically stable method to compute those quantities. Although precise this method is expensive in CPU time. Here, we present an alternative approach which is an order of magnitude quicker, just as precise, and very simple to implement. The method is based on the observation that for a given auxiliary field the equal time Green function matrix, $G$, is a projector: $G^2 = G$.Comment: 4 papes, 1 figure in eps forma

The Fractional Quantum Hall Effect on a Lattice

Starting from the Hofstadter butterfly, we define lattice versions of Landau levels as well as a continuum limit which ensures that they scale to continuum Landau levels. By including a next-neighbor repulsive interaction and projecting onto the lowest lattice Landau level, we show that incompressible ground states exist at filling fractions, $\nu = 1/3, 2/5$ and $3/7$. Already for values of $l_0/a \sim 2$ where $l_0$ ($a$) is the magnetic length (lattice constant), the lattice version of the $\nu = 1/3$ state reproduces with nearly perfect accuracy the the continuum Laughlin state. The numerical data strongly suggests that at odd filling fractions of the lowest lattice Landau level, the lattice constant is an irrelevant length scale. We find a new relation between the hierarchy of incompressible states and the self-similar structure of the Hofstadter butterfly.Comment: 11 pages, Latex, 3 compressed and uuencoded postscript figures appended, FFA-94-0

Edge instabilities of topological superconductors

Nodal topological superconductors display zero-energy Majorana flat bands at generic edges. The flatness of these edge bands, which is protected by time-reversal and translation symmetry, gives rise to an extensive ground-state degeneracy. Therefore, even arbitrarily weak interactions lead to an instability of the flat-band edge states towards time-reversal and translation-symmetry-broken phases, which lift the ground-state degeneracy. We examine the instabilities of the flat-band edge states of d_{xy}-wave superconductors by performing a mean-field analysis in the Majorana basis of the edge states. The leading instabilities are Majorana mass terms, which correspond to coherent superpositions of particle-particle and particle-hole channels in the fermionic language. We find that attractive interactions induce three different mass terms. One is a coherent superposition of imaginary s-wave pairing and current order, and another combines a charge-density-wave and finite-momentum singlet pairing. Repulsive interactions, on the other hand, lead to ferromagnetism together with spin-triplet pairing at the edge. Our quantum Monte Carlo simulations confirm these findings and demonstrate that these instabilities occur even in the presence of strong quantum fluctuations. We discuss the implications of our results for experiments on cuprate high-temperature superconductors.Comment: 4 pages, 3 figure

Metamagnetism and Lifshitz Transitions in Models for Heavy Fermions

We investigate metamagnetic transitions in models for heavy fermions by considering the doped Kondo lattice model in two dimensions. Results are obtained within the framework of dynamical mean field and dynamical cluster approximations. Universal magnetization curves for different temperatures and Kondo couplings develop upon scaling with the lattice coherence temperature. Furthermore, the coupling of the local moments to the magnetic field is varied to take into account the different Land\'e factors of localized and itinerant electrons. The competition between the lattice coherence scale and the Zeeman energy scale allows for two interpretations of the metamagnetism in heavy fermions: Kondo breakdown or Lifshitz transitions. By tracking the single-particle residue through the transition, we can uniquely conclude in favor of the Lifshitz transition scenario. In this scenario, a quasiparticle band drops below the Fermi energy which leads to a change in topology of the Fermi surface.Comment: 8 pages, 7 figure

Dynamical dimer correlations at bipartite and non-bipartite Rokhsar-Kivelson points

We determine the dynamical dimer correlation functions of quantum dimer models at the Rokhsar-Kivelson point on the bipartite square and cubic lattices and the non-bipartite triangular lattice. Based on an algorithmic idea by Henley, we simulate a stochastic process of classical dimer configurations in continuous time and perform a stochastic analytical continuation to obtain the dynamical correlations in momentum space and the frequency domain. This approach allows us to observe directly the dispersion relations and the evolution of the spectral intensity within the Brillouin zone beyond the single-mode approximation. On the square lattice, we confirm analytical predictions related to soft modes close to the wavevectors (pi,pi) and (pi,0) and further reveal the existence of shadow bands close to the wavevector (0,0). On the cubic lattice the spectrum is also gapless but here only a single soft mode at (pi,pi,pi) is found, as predicted by the single mode approximation. The soft mode has a quadratic dispersion at very long wavelength, but crosses over to a linear behavior very rapidly. We believe this to be the remnant of the linearly dispersing "photon" of the Coulomb phase. Finally the triangular lattice is in a fully gapped liquid phase where the bottom of the dimer spectrum exhibits a rich structure. At the M point the gap is minimal and the spectral response is dominated by a sharp quasiparticle peak. On the other hand, at the X point the spectral function is much broader. We sketch a possible explanation based on the crossing of the coherent dimer excitations into the two-vison continuum.Comment: 16 pages, 7 figures, published versio

Spontaneous particle-hole symmetry breaking of correlated fermions on the Lieb lattice

We study spinless fermions with nearest-neighbor repulsive interactions ($t$-$V$ model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However, due to the three-orbital unit cell, commensurate charge order implies an imbalance of electron and hole densities and therefore doping away from half-filling. Our numerical results show that below a finite-temperature Ising transition a charge density wave with one electron and two holes per unit cell and its partner under particle-hole transformation are spontaneously generated. Our calculations are based on recent advances in auxiliary-field and continuous-time quantum Monte Carlo simulations that allow sign-free simulations of spinless fermions at half-filling. It is argued that particle-hole symmetry breaking provides a route to access levels of finite doping, without introducing a sign problem.Comment: 9 pages, 6 figures, added data for strong Coulomb repulsion and classical Ising-limi

Comment on "Quantum Monte Carlo Evidence for Superconductivity in the Three-Band Hubbard Model in Two Dimensions"

In a recent Letter, Kuroki and Aoki [Phys. Rev. Lett. 76, 440 (1996)] presented quantum Monte-Carlo (QMC) results for pairing correlations in the three-band Hubbard model, which describes the Cu-d_{x^2-y^2} and O-p_{x,y} orbitals present in the CuO_2 planes of high-T_c materials. In this comment we argue that (i) the used parameter set is not appropriate for the description of high-T_c materials since it does not satisfy the minimal requirement of a charge-transfer gap at half-filling, and (ii) the observed increase in the d_{x^2-y^2} channel is dominantly produced by the pair-field correlations without the vertex part. Hence, the claim of evidence of ODLRO is not justified.Comment: 1 page latex and 2 eps-figures, uses epsfig, submitted to PR