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

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

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

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

Spin nematic phases in models of correlated electron systems: a numerical study

Strongly interacting systems are known to often spontaneously develop exotic ground states under certain conditions. For instance, spin nematic phases have been discovered in various magnetic models. Such phases, which break spin symmetry but have no net local magnetization, have also been proposed by Nersesyan et al. (J. Phys.: Cond. Matt. 3, 3353 (1991)) in the context of electronic models. We introduce a N-flavor microscopic model that interpolates from the large-N limit, where mean-field is valid and such a nematic phase occurs, to the more realistic N=1 case. By using a sign-free quantum Monte-Carlo, we show the existence of a spin nematic phase (analogous to a spin flux phase) for finite N; when N decreases, quantum fluctuations increase and this phase ultimately disappears in favor of an s-wave superconducting state. We also show that this nematic phase extends up to a finite critical charge doping. Dynamical studies allow us to clarify the Fermi surface property: in the nematic phase at half-filling, it consists of 4 points and the low-energy structure has a Dirac cone-like shape. Under doping, we observe clear signatures of Fermi pockets around these points. This is one of the few examples where numerical simulations show how quantum fluctuations can destroy a large-N phase.Comment: 9 pages, 19 figures. Problem with figures has been fixe

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