Source author record

A. -M. S. Tremblay

A. -M. S. Tremblay 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

50works
10topics
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

50 published item(s)

preprint2022arXiv

Disorder effects on hot spots in electron-doped cuprates

Antiferromagnetic fluctuations in two dimensions cause a decrease in spectral weight at so-called hot spots associated with the pseudogap in electron-doped cuprates. In the 2D Hubbard model, these hot spots occur when the Vilk criterion is satisfied, namely when the spin correlation length exceeds the thermal de Broglie wavelength. Using the two-particle self-consistent approach, here we show that this criterion may be violated near the antiferromagnetic quantum critical point when sufficient disorder is added to the model. Static disorder decreases inelastic scattering, in contradiction with Matthiessen's rule, leading to a shift in the position of the quantum critical point and a modification of the conditions for the appearance of hot spots. This opens the road to a study of the interplay between disorder, antiferromagnetic fluctuations and superconductivity in electron-doped cuprates.

preprint2021arXiv

Fermi arcs vs hole pockets: periodization of a cellular two-band model

It is still debated whether the low-doping Fermi surface of cuprates is composed of hole pockets or of disconnected Fermi arcs. Results from cellular dynamical mean field theory (c-DMFT) support the Fermi arcs hypothesis by predicting corresponding Fermi arcs for the Hubbard model. Here, we introduce a simple parametrization of the self-energy, in the spirit of Yang-Rice-Zhang theory, and show that state of the art c-CDMFT calculations cannot give a definitive answer to the question of Fermi arcs vs holes pockets, and this, independently of the periodization (cumulant or Green's function) used to display spectral weights of the infinite lattice. Indeed, when our model is restricted to a cluster and periodized like in c-DMFT, only two adjustable parameters suffice to reproduce the qualitative details of the frequency and momentum dependence of the low energy c-DMFT spectral weight for both periodizations. In other words, even though our starting model has a Fermi surface composed of hole and electron pockets, it leads to Fermi arcs when restricted to a cluster and periodized like in c-DMFT. We provide a new "compact tiling" scheme to recover the hole and electron pockets of our starting non-interacting lattice model, suggesting that better periodization schemes might exist.

preprint2021arXiv

Information-theoretic measures of superconductivity in a two-dimensional doped Mott insulator

A key open issue in condensed matter physics is how quantum and classical correlations emerge in an unconventional superconductor from the underlying normal state. We study this problem in a doped Mott insulator with information theory tools on the two-dimensional Hubbard model at finite temperature with cluster dynamical mean-field theory. We find that the local entropy detects the superconducting state and that the difference in the local entropy between the superconducting and normal states follows the same difference in the potential energy. We find that the thermodynamic entropy is suppressed in the superconducting state and monotonically decreases with decreasing doping. The maximum in entropy found in the normal state above the overdoped region of the superconducting dome is obliterated by superconductivity. The total mutual information, which quantifies quantum and classical correlations, is amplified in the superconducting state of the doped Mott insulator for all doping levels, and shows a broad peak versus doping, as a result of competing quantum and classical effects.

preprint2021arXiv

Non-Fermi liquid phase and linear-in-temperature scattering rate in overdoped two dimensional Hubbard model

Understanding electronic properties that violate the Landau Fermi liquid paradigm in cuprate superconductors remains a major challenge in condensed matter physics. The strange metal state in overdoped cuprates that exhibits linear-in-temperature scattering rate and dc resistivity is a particularly puzzling example. Here, we compute the electronic scattering rate in the two-dimensional Hubbard model using cluster generalization of dynamical mean-field theory. We present a global phase diagram documenting an apparent non-Fermi liquid phase, in between the pseudogap and Fermi liquid phase in the doped Mott insulator regime. We discover that in this non-Fermi liquid phase, the electronic scattering rate $γ_k(T)$ can display linear temperature dependence as temperature $T$ goes to zero. In the temperature range that we can access, the $T-$ dependent scattering rate is isotropic on the Fermi surface, in agreement with recent experiments. Using fluctuation diagnostic techniques, we identify antiferromagnetic fluctuations as the physical origin of the $T-$ linear electronic scattering rate.

preprint2021arXiv

Oxygen hole content, charge-transfer gap, covalency, and cuprate superconductivity

Experiments have shown that the families of cuprate superconductors that have the largest transition temperature at optimal doping also have the largest oxygen hole content at that doping. They have also shown that a large charge-transfer gap, a quantity accessible in the normal state, is detrimental to superconductivity. We solve the three-band Hubbard model with cellular dynamical mean-field theory and show that both of these observations follow from the model. Cuprates play a special role amongst doped charge-transfer insulators of transition metal oxides because copper has the largest covalent bonding with oxygen.

preprint2021arXiv

Resilient Fermi liquid and strength of correlations near an antiferromagnetic quantum critical point

Near the antiferromagnetic quantum critical point (QCP) of electron-doped cuprate superconductors, angle-resolved photoemission experiments detect hot spots where the Fermi surface disappears. Here we demonstrate, using the two-particle self-consistent theory, that in the antinodal region the Fermi liquid remains stable for a broad range of angles on the Fermi surface and for all dopings near the QCP. We show how the quasiparticle weight Z and effective mass m* change and then abruptly become meaningless as the hot spots are approached. We propose a dimensionless number, easily accessible in ARPES experiments, that can be used to gauge the strength of correlations.

preprint2020arXiv

Analysis of the magnetic and magnetocaloric properties of ALaFeMnO6 (A= Sr, Ba and Ca) double perovskites

In previous studies, we have reported that double perovskite La2NiMnO6 presents a non-negligible potential for room temperature magnetocaloric tasks. With the aim of improving even further the cooling performances and the working temperature range of double perovskites, we report the magnetic and magnetocaloric properties of La2MnFeO6 and ALaMnFeO6 (A = Sr, Ba, Ca) compounds. X-ray diffraction (XRD) and Rietveld refinement show that La2MnFeO6 (LMFO) and CaLaMnFeO6 (CLMFO) samples crystallize in an orthorhombic structure with the Pnma space group. However, a rhombohedral structure with the R3C space group is obtained for BaLaMnFeO6 (BLMFO) and SrLaMnFeO6 (SLMFO) samples. Substituting La by Ba or Sr in LMFO leads to a clear increase of the Curie temperature (Tc) compared to LMFO from 150 K for BLMFO up to 350 K for SLMFO. Moreover, CLMFO shows the smallest Tc down to 70 K. Ferromagnetic-like behavior is observed for SLMFO and BLMFO while CLMFO's magnetism resembles that of LMFO. A clear connection between the structural parameters and the magnetic properties of these doped LMFO samples is unveiled as the highest Tc and the largest magnetization are observed for SLMFO which shows also bond angles closest to 180° and the smallest bond lengths, thus optimizing the superexchange interaction. The partial substitution of Sr for La leads in fact to a significant magnetocaloric effect over a wide operating temperature range extending beyond 300 K. For some optimal growth conditions, its entropy change varies slowly over an unusually large temperature range, which is of clear interest from a practical point of view.

preprint2020arXiv

Interaction and temperature effects on the magneto-optical conductivity of Weyl liquids

Negative magnetoresistance is one of the manifestations of the chiral anomaly in Weyl semimetals. The magneto-optical conductivity also shows transitions between Landau levels that are not spaced as in an ordinary electron gas. How are such topological properties modified by interactions and temperature? We answer this question by studying a lattice model of Weyl semimetals with an on-site Hubbard interaction. Such an interacting Weyl semimetal, dubbed as Weyl liquid, may be realized in Mn$_3$Sn. We solve that model with single-site dynamical mean-field theory. We find that in a Weyl liquid, quasiparticles can be characterized by a quasiparticle spectral weight $Z$, although their lifetime increases much more rapidly as frequency approaches zero than in an ordinary Fermi liquid. The negative magnetoresistance still exists, even though the slope of the linear dependence of the DC conductivity with respect to the magnetic filed is decreased by the interaction. At elevated temperatures, a Weyl liquid crossesover to bad metallic behavior where the Drude peak becomes flat and featureless.

preprint2020arXiv

Mott transition and high-temperature crossovers at half-filling

The interaction-driven Mott transition in the half-filled Hubbard model is a first-order phase transition that terminates at a critical point $(T_\mathrm{c},U_\mathrm{c})$ in the temperature-interaction plane $T-U$. A number of crossovers occur along lines that extend for some range above $(T_\mathrm{c},U_\mathrm{c})$. Asymptotically close to $(T_\mathrm{c},U_\mathrm{c})$, these lines coalesce into the so-called Widom line. The existence of $(T_\mathrm{c},U_\mathrm{c})$ and of the associated crossovers becomes unclear when long-wavelength fluctuations or long-range order occur above $(T_\mathrm{c},U_\mathrm{c})$. We study this problem using continuous-time quantum Monte Carlo methods as impurity solvers for both Dynamical Mean-Field Theory (DMFT) and Cellular Dynamical Mean-Field Theory (CDMFT). We contrast the cases of the square lattice, where antiferromagnetic fluctuations dominate in the vicinity of the Mott transition, and the triangular lattice where they do not. The inflexion points and maxima found near the Widom line for the square lattice can serve as proxy for the triangular lattice case. But the only crossover observable in all cases at sufficiently high temperature is that associated with the opening of the Mott gap. The same physics also controls an analog crossover in the resistivity called the "Quantum Widom line".

preprint2019arXiv

Intrinsic Cluster-Shaped Density Waves in Cellular Dynamical Mean-Field Theory

It is well known that cellular dynamical mean-field theory (CDMFT) leads to the artificial breaking of translation invariance. In spite of this, it is one of the most successful methods to treat strongly correlated electrons systems. Here, we investigate in more detail how this broken translation invariance manifests itself. This allows to disentangle artificial broken translation invariance effects from the genuine strongly correlated effects captured by CDMFT. We report artificial density waves taking the shape of the cluster---cluster density waves---in all our zero temperature CDMFT solutions, including pair density waves in the superconducting state. We discuss the limitations of periodization regarding this phenomenon, and we present mean-field density-wave models that reproduce CDMFT results at low energy in the superconducting state. We then discuss how these artificial density waves help the agreement of CDMFT with high temperature superconducting cuprates regarding the low-energy spectrum, in particular for subgap structures observed in tunnelling microscopy. We relate these subgap structures to nodal and anti-nodal gaps in our results, similar to those observed in photoemission experiments. This fortuitous agreement suggests that spatial inhomogeneity may be a key ingredient to explain some features of the low-energy underdoped spectrum of cuprates with strongly correlated methods. This work deepens our understanding of CDMFT and clearly identifies signatures of broken translation invariance in the presence of strong correlations.

preprint2016arXiv

Algorithms for optimized maximum entropy and diagnostic tools for analytic continuation

Analytic continuation of numerical data obtained in imaginary time or frequency has become an essential part of many branches of quantum computational physics. It is, however, an ill-conditioned procedure and thus a hard numerical problem. The maximum-entropy approach, based on bayesian inference, is the most widely used method to tackle that problem. Although the approach is well established and among the most reliable and efficient ones, useful developments of the method and of its implementation are still possible. In addition, while a few free software implementations are available, a well-documented, optimized, general purpose and user-friendly software dedicated to that specific task is still lacking. Here we analyze all aspects of the implementation that are critical for accuracy and speed, and present a highly optimized approach to maximum-entropy. Original algorithmic and conceptual contributions include (1) numerical approximations that yield a computational complexity that is almost independent of temperature and spectrum shape (including sharp Drude peaks in broad background for example) while ensuring quantitative accuracy of the result whenever precision of the data is sufficient, (2) a robust method of choosing the entropy weight $α$ that follows from a simple consistency condition of the approach and the observation that information- and noise-fitting regimes can be identified clearly from the behavior of $χ^2$ with respect to $α$, and (3) several diagnostics to assess the reliability of the result. Benchmarks with test spectral functions of different complexity and an example with an actual physical simulation are presented. Our implementation, which covers most typical cases for fermions, bosons and response functions, is available as an open source, user friendly software.

preprint2016arXiv

An organizing principle for two-dimensional strongly correlated superconductivity

Superconductivity in the cuprates exhibits many unusual features. We study the two-dimensional Hubbard model with plaquette dynamical mean-field theory to address these unusual features and relate them to other normal-state phenomena, such as the pseudogap. Previous studies with this method found that upon doping the Mott insulator at low temperature a pseudogap phase appears. The low-temperature transition between that phase and the correlated metal at higher doping is first-order. A series of crossovers emerge along the Widom line extension of that first-order transition in the supercritical region. Here we show that the highly asymmetric dome of the dynamical mean-field superconducting transition temperature $T_c^d$, the maximum of the condensation energy as a function of doping, the correlation between maximum $T_c^d$ and normal-state scattering rate, the change from potential-energy driven to kinetic-energy driven pairing mechanisms can all be understood as remnants of the normal state first-order transition and its associated crossovers that also act as an organizing principle for the superconducting state.

preprint2016arXiv

Antagonistic effects of nearest-neighbor repulsion on the superconducting pairing dynamics in the doped Mott insulator regime

The nearest-neighbor superexchange-mediated mechanism for d_{x^2-y^2}-wave superconductivity in the one-band Hubbard model faces the challenge that nearest-neighbor Coulomb repulsion can be larger than superexchange. To answer this question, we use cellular dynamical mean-field theory (CDMFT) with a continuous-time quantum Monte Carlo solver to determine the superconducting phase diagram as a function of temperature and doping for on-site repulsion $U=9t$ and nearest-neighbor repulsion $V=0,2t,4t$. In the underdoped regime, $V$ increases the CDMFT superconducting transition temperature $T_c^d$ even though it decreases the superconducting order parameter at low temperature for all dopings. However, $V$ decreases $T_c^d$ in the overdoped regime. We gain insight into these paradoxical results through a detailed study of the frequency dependence of the anomalous spectral function, extracted at finite temperature via the MaxEntAux method for analytic continuation. A systematic study of dynamical positive and negative contributions to pairing reveals that even though $V$ has a high-frequency depairing contribution, it also has a low frequency pairing contribution since it can reinforce superexchange through $J=4t^2/(U-V)$. Retardation is thus crucial to understand pairing in doped Mott insulators, as suggested by previous zero-temperature studies. We also comment on the tendency to charge order for large $V$ and on the persistence of d-wave superconductivity over extended-$s$ or s+d-wave.

preprint2016arXiv

Correlation-enhanced odd-parity inter-orbital singlet pairing in the iron-pnictide superconductor LiFeAs

The rich variety of iron-based superconductors and their complex electronic structure lead to a wide range of possibilities for gap symmetry and pairing components. Here we solve in the 2-Fe Brillouin zone the full frequency-dependent linearized Eliashberg equations for LiFeAs with spin-fluctuation mediated pairing interactions. The magnetic excitations are calculated with the random phase approximation on a correlated electronic structure obtained with density functional theory and dynamical mean field theory. Correlations induce long-lived local moments with orbital-dependent dynamics. The interaction between electrons through Hund's coupling promotes inter-orbital magnetic susceptibility. As a consequence, the leading pairing channel, conventional $s^{+-}$, acquires sizeable inter-orbital $d_{xy}-d_{xz(yz)}$ singlet pairing with odd parity under glide-plane symmetry. These components reduce the superconducting gap magnitude induced by the intra-orbital components of the gap function at the electron pockets intersection where the Fe-d $t_{2g}$ orbitals strongly mix. This in turn makes the results consistent with available experiments on the angular dependence of the gaps observed on the different Fermi surfaces.

preprint2016arXiv

Mott transition and magnetism on the anisotropic triangular lattice

Spin-liquid behavior was recently suggested experimentally in the moderately one-dimensional organic compound $κ$-H$_3$(Cat-EDT-TTF)$_2$. This compound can be modeled by the one-band Hubbard model on the anisotropic triangular lattice with $t^\prime/t \simeq 1.5$, where $t'$ is the minority hopping. It thus becomes important to extend previous studies, that were performed in the range $0 \leq t'/t \leq 1.2$, to find out whether there is a regime where Mott insulating behavior can be found without long-range magnetic order. To this end, we study the above model in the range $1.2 \leq t'/t \leq 2$ using cluster dynamical mean-field theory (CDMFT). We argue that it is important to choose a symmetry-preserving cluster rather than a quasi one-dimensional cluster. We find that, upon increasing $t'/t$ beyond $t^\prime/t \approx 1.3$, the Mott transition at zero-temperature is replaced by a first-order transition separating a metallic state from a collinear magnetic insulating state. Nevertheless, at the physically relevant value $t^\prime/t \simeq 1.5$, the transitions toward the magnetic and the Mott insulating phases are very close. The phase diagram obtained in this study can provide a working basis for moderately one-dimensional compounds on the anisotropic triangular lattice.

preprint2016arXiv

Pseudogap and superconductivity in two-dimensional doped charge-transfer insulators

High-temperature superconductivity emerges in the CuO$_2$ plane upon doping a Mott insulator. To ascertain the influence of Mott physics plus short-range correlations, we solve a three-band copper-oxide model in the charge-transfer regime using cellular dynamical mean-field theory with continuous-time quantum Monte Carlo as an impurity solver. We report the normal and superconducting phase diagram of this model as a function of doping, interaction strength and temperature. Upon hole doping of the charge-transfer insulator, the phase boundary between pseudogap and correlated metal consists of a first-order transition line at finite doping ending at a critical point, as in the one-band model. Beyond the endpoint, the phase boundary continues as a Widom crossover line, across which thermodynamic quantities peak. This phase boundary determines changes in the pairing mechanism and is an emergent phenomenon characteristic of doped Mott insulators, independent of many microscopic details. Broader implications are discussed.

preprint2015arXiv

Antiferromagnetism in the Hubbard model on the honeycomb lattice: a Two-Particle Self-Consistent study

The semimetal to antiferromagnet quantum phase transition of the Hubbard model on the honeycomb lattice has come to the forefront in the context of the proposal that a semimetal to spin liquid transition can occur before the transition to the antiferromagnetic phase. To study the semimetal to antiferromagnet transition, we generalize the two-particle self-consistent (TPSC) approach to the honeycomb lattice (a structure that can be realized in graphene for example). We show that the critical interaction strength where the transition occurs is $U_c/t=3.79\pm 0.01$ quite close to the value $U_c/t=3.869 \pm 0.013$ reported using large-scale quantum Monte Carlo simulations. This reinforces the conclusion that the semimetal to spin liquid transition is pre-empted by the transition to the antiferromagnet. Since TPSC satisfies the Mermin-Wagner theorem, we find temperature-dependent results for the antiferromagnetic and ferromagnetic correlation lengths as well as the dependence of double occupancy and of the renormalized spin and charge interactions on the bare interaction strength. We also estimate the value of the crossover temperature to the renormalized classical regime as a function of interaction strength.

preprint2015arXiv

Maximum Entropy Analytic Continuation for Spectral Functions with Non-Positive Spectral Weight

Information about the pairing mechanism for superconductivity is contained in the spectral weight for the anomalous (Gorkov) Green function. In the most general case, this spectral weight can change sign on the positive real axis or even be complex in the presence of broken time-reversal symmetry. This renders impossible the direct analytic continuation with Maximum Entropy methods of numerical results obtained for the anomalous Green function either in Matsubara frequency or in imaginary time. Here we show that one can define auxiliary Green functions that allow one to extract the anomalous spectral weight in the most general multi-orbital spin-singlet or triplet cases with no particular symmetry. As a simple example, we treat the case of lead in Eliashberg theory: Maximum Entropy analytic continuation with auxiliary Green functions (MaxEntAux) allows us to recover spectral functions that change sign on the positive real axis. The approach can be extended to transport quantities.

preprint2015arXiv

Superconducting dome in doped quasi-2d organic Mott insulators: a paradigm for strongly-correlated superconductivity

Layered organic superconductors of the BEDT family are model systems for the interplay of the Mott transition with superconductivity, magnetic order and frustration. Recent experimental studies on a hole-doped version of BEDT compounds reveal an enhancement of superconductivity and a rapid crossover between two different conducting phases above the superconducting dome. One of these phases is a Fermi liquid, the other not. Using plaquette cellular dynamical mean field theory with state of the art continuous-time quantum Monte Carlo calculations, we study this problem with the two-dimensional Hubbard model on the anisotropic triangular lattice. Phase diagrams as a function of temperature $T$ and interaction strength $U/t$ are obtained for anisotropy parameters $t'=0.4t$, $t'=0.8t$ and various fillings. As for cuprates, we find, at finite doping, a first-order transition between two normal-state phases. For $T$ above the critical point of the first-order transition, there is a Widom line where crossovers occur. The results are in broad agreement with experiment. This suggests that for compounds with intermediate to high frustration, very light-doping should reveal the first-order transition and associated crossovers. These crossovers could leave traces in the superconducting phase. We also predict that destroying the superconducting phase by a magnetic field should reveal the first-order transition. Finally, we predict that electron-doping should also lead to an increased range of $U/t$ for superconductivity but with a reduced maximum $T_c$. This work also clearly shows that the superconducting dome here is tied to the Mott transition and its continuation as a transition separating pseudogap phase from correlated metal in doped compounds, as in the cuprates. Contrary to heavy fermions for example, the maximum $T_c$ is definitely not attached to an antiferromagnetic quantum critical point.

preprint2014arXiv

Electric polarization of Sr$_{0.5}$Ba$_{0.5}$MnO$_{3}$: a multiferroic Mott insulator

Ab initio calculations of the electric polarization of correlation-driven insulating materials, namely Mott insulators, have not been possible so far. Using a combination of density functional theory and dynamical-mean-field theory we study the electric polarization of the Mott insulator Sr$_{0.5}$Ba$_{0.5}$MnO$_{3}$. We predict a ferroelectric polarization of $\simeq 16.5 μC/cm^2$ in the high temperature paramagnetic phase and recover the measured value of $\simeq 13.3 μC/cm^2$ in the low temperature antiferromagnetic phase. Our calculations reveal that the driving force for the ferroelectricity, the hybridization between Mn e$_g$ and O p orbitals, is suppressed by correlations, in particular by the Hund coupling and by the onset of magnetic order. They also confirm that the half-filled Mn $t_{2g}$ orbitals give rise to the antiferromagnetic Mott phase. This magnetic ordering leads to changes in the ionic polar displacement and in turn to the electronic polarization. In addition, for a fixed ionic displacement, we find that there is a reduction in the electronic contribution due to partial magnetic polarization of the Mn e$_g$ orbitals. The reduction of the polarization due to ionic displacement dominates over the additional electronic part, hence the net magneto-electric coupling is negative.

preprint2014arXiv

Ergodicity of the Hybridization-Expansion Monte Carlo Algorithm for Broken-Symmetry States

With the success of dynamical mean field theories, solvers for quantum-impurity problems have become an important tool for the numerical study of strongly correlated systems. Continuous-time Quantum Monte Carlo sampling of the expansion in powers of the hybridization between the ``impurity'' and the bath provides a powerful solver when interactions are strong. Here we show that the usual updates that add or remove a pair of creation-annihilation operators are rigorously not ergodic for several classes of broken-symmetries that involve spatial components. We show that updates with larger numbers of simultaneous updates of pairs of creation-annihilation operators remedy this problem. As an example, we apply the four operator updates that are necessary for ergodicity to the case of d-wave superconductivity in plaquette dynamical mean-field theory for the one-band Hubbard model. While the results are qualitatively similar to those previously published, they are quantitatively better that previous ones, being closer to those obtained by other approaches.

preprint2014arXiv

Impurity-induced magnetic moments on the graphene-lattice Hubbard model: an inhomogeneous cluster DMFT study

Defect-induced magnetic moments are at the center of the research effort on spintronic applications of graphene. Here we study the problem of a nonmagnetic impurity in graphene with a new theoretical method, inhomogeneous cluster dynamical mean field theory (I-CDMFT), which takes into account interaction-induced short-range correlations while allowing long-range inhomogeneities. The system is described by a Hubbard model on the honeycomb lattice. The impurity is modeled by a local potential. For a large enough potential, interactions induce local antiferromagnetic correlations around the impurity and a net total spin $\frac12$ appears, in agreement with Lieb's theorem. Bound states caused by the impurity are visible in the local density of states (LDOS) and have their energies shifted by interactions in a spin-dependent way, leading to the antiferromagnetic correlations. Our results take into account dynamical correlations; nevertheless they qualitatively agree with previous mean-field and density functional theory (DFT) studies. Moreover, they provide a relation between impurity potential and on-site repulsion $U$ that could in principle be used to determine experimentally the value of $U$.

preprint2014arXiv

Lazy skip lists, a new algorithm for fast hybridization-expansion quantum Monte Carlo

The solution of a generalized impurity model lies at the heart of electronic structure calculations with dynamical mean-field theory (DMFT). In the strongly-correlated regime, the method of choice for solving the impurity model is the hybridization expansion continuous time quantum Monte Carlo(CT-HYB). Enhancements to the CT-HYB algorithm are critical for bringing new physical regimes within reach of current computational power. Taking advantage of the fact that the bottleneck in the algorithm is a product of hundreds of matrices, we present optimizations based on the introduction and combination of two concepts of more general applicability: a) skip lists and b) fast rejection of proposed configurations based on matrix bounds. Considering two very different test cases with $d$ electrons, we find speedups of $\sim 25$ up to $\sim 500$ compared to the direct evaluation of the matrix product. Even larger speedups are likely with $f$ electron systems and with clusters of correlated atoms.

preprint2014arXiv

Orbital magnetization of correlated electrons with arbitrary band topology

Spin-orbit coupling introduces chirality into electronic structure. This can have profound effects on the magnetization induced by orbital motion of electrons. Here we derive a formula for the orbital magnetization of interacting electrons in terms of the full Green's function and vertex functions. The formula is applied within dynamical mean-field theory to the Kane-Mele-Hubbard model that allows both topological and trivial insulating phases. We study the insulating and metallic phases in the presence of an exchange magnetic field. In the presence of interactions, the orbital magnetization of the quantum spin Hall insulating phase with inversion symmetry is renormalized by the bulk quasi-particle weight. The orbital magnetization vanishes for the in-plane antiferromagnetic phase with trivial topology. In the metallic phase, the enhanced effective spin-orbit coupling due to the interaction sometimes leads to an enhancement of the orbital magnetization. However, at low doping, magnetization is suppressed at large interaction strengths.

preprint2013arXiv

Strongly correlated superconductivity

Band theory and BCS theory are arguably the most successful theories of condensed matter physics. Yet, in a number of materials, in particular the high-temperature superconductors and the layered organic superconductors, they fail. In these lecture notes for an international school, I emphasize that even though the low energy properties of a phase of matter are generally emergent and entirely determined by the broken symmetry, there are many differences between a strongly correlated and a weakly correlated state of matter. For example, spin waves are an emergent property for antiferromagnets, but for weak correlations the normal phase of a (Slater) antiferromagnet is metallic, whereas it is insulating (Heisenberg) for strong correlations. As a function of interaction strength, above the antiferromagnetic phase, the crossover between a metal and a local moment paramagnetic insulator is described by the Mott transition, whose mean-field theory down to T=0 is best formulated with dynamical mean-field theory. A similar situation occurs for superconductors: despite similar emergent properties, there are many differences between both kinds of superconductors. Experimental evidence suggests that hole-doped cuprates are strongly correlated whereas the correlations are weaker in the electron-doped case, especially near optimal doping. Evidence for a pseudogap arising from antiferromagnetic fluctuations is strong in the latter case, whereas in the hole-doped case the pseudogap temperature appears as a consequence of Mott physics. In cluster dynamical mean-field theory, the pseudogap line T* is a Widom line arising from a T=0 first-order transition that terminates at a finite T critical point. Strongly correlated superconductors are much more resilient to near-neighbor repulsion than their weakly correlated counterpart. Many different methods to attack these problems theoretically are described.

preprint2013arXiv

Transport functions for hypercubic and Bethe lattices

In calculations of transport quantities, such as the electrical conductivity, thermal conductivity, Seebeck, Peltier, Nernst, Ettingshausen, Righi-Leduc, or Hall coefficients, sums over the Brillouin zone of wave-vector derivatives of the dispersion relation commonly appear. When the self-energy depends only on frequency, as in single-site dynamical mean-field theory, it is advantageous to perform these sums once and for all. We show here that in the case of a hypercubic lattice in d dimensions, the sums needed for any of the transport coefficients can be expressed as integrals over powers of the energy weighted by the energy-dependent non-interacting density of states. It is also shown that our exact expressions for the transport functions can be obtained from differential equations that follow from sum rules. By substituting the Bethe lattice density of states, one can obtain the previously unknown transport function for the electrical or thermal Hall coefficients and for the Nernst coefficient of the Bethe lattice.

preprint2012arXiv

Benchmark of a modified Iterated Perturbation Theory approach on the 3d FCC lattice at strong coupling

The Dynamical Mean-Field theory (DMFT) approach to the Hubbard model requires a method to solve the problem of a quantum impurity in a bath of non-interacting electrons. Iterated Perturbation Theory (IPT) has proven its effectiveness as a solver in many cases of interest. Based on general principles and on comparisons with an essentially exact Continuous-Time Quantum Monte Carlo (CTQMC) solver, here we show that the standard implementation of IPT fails away from half-filling when the interaction strength is much larger than the bandwidth. We propose a slight modification to the IPT algorithm that replaces one of the equations by the requirement that double occupancy calculated with IPT gives the correct value. We call this method IPT-$D$. We recover the Fermi liquid ground state away from half-filling. The Fermi liquid parameters, density of states, chemical potential, energy and specific heat on the FCC lattice are calculated with both IPT-$D$ and CTQMC as benchmark examples. We also calculated the resistivity and the optical conductivity within IPT-$D$. Particle-hole asymmetry persists even at coupling twice the bandwidth. Several algorithms that speed up the calculations are described in appendices.

preprint2012arXiv

Breakdown of Fermi liquid behavior at the (π,π)=2k_F spin-density wave quantum-critical point: the case of electron-doped cuprates

Many correlated materials display a quantum critical point between a paramagnetic and a SDW state. The SDW wave vector connects points (hot spots) on opposite sides of the Fermi surface. The Fermi velocities at these pairs of points are in general not parallel. Here we consider the case where pairs of hot spots coalesce, and the wave vector (π,π) of the SDW connects hot spots with parallel Fermi velocities. Using the specific example of electron-doped cuprates, we first show that Kanamori screening and generic features of the Lindhard function make this case experimentally relevant. The temperature dependence of the correlation length, the spin susceptibility and the self-energy at the hot spots are found using the Two-Particle-Self-Consistent theory and specific numerical examples worked out for parameters characteristic of the electron-doped cuprates. While the curvature of the Fermi surface at the hot spots leads to deviations from perfect nesting, the pseudo-nesting conditions lead to drastic modifications to the temperature dependence of these physical observables: Neglecting logarithmic corrections, the correlation length ξscales like 1/T, i.e. z=1 instead of the naive z=2, the (π,π) static spin susceptibility χlike $1/\sqrt T$, and the imaginary part of the self-energy at the hot spots like $T^{3/2}$. The correction $T_1^{-1}\sim T^{3/2}$ to the Korringa NMR relaxation rate is subdominant. We also consider this problem at zero temperature, or for frequencies larger than temperature, using a field-theoretical model of gapless SDW fluctuations interacting with fermions. The imaginary part of the fermionic self-energy close to the hot spots scales as $-ω^{3/2}\logω$. This is less singular than earlier predictions of the form $-ω\logω$. The difference arises from the effects of umklapp terms that were not included in previous studies.

preprint2012arXiv

c-axis resistivity, pseudogap, superconductivity and Widom line in doped Mott insulators

Layered doped Mott insulators, such as the cuprates, show unusual temperature dependence of the resistivity. Intriguingly, the resistivity perpendicular to the CuO$_2$ planes, $ρ_c(T)$, shows both metallic ($dρ_c/dT > 0$) and semi-conducting ($dρ_c/dT<0$) behavior. We shed light on this puzzle by calculating $ρ_c$ for the two-dimensional Hubbard model within plaquette cellular dynamical mean-field theory and strong-coupling continuous-time quantum Monte Carlo as the impurity solver. The temperature, $T$, and doping, $δ$, dependencies of $ρ_c$ are controlled by the first-order transition between pseudogap and correlated metal phases from which superconductivity can emerge. On the large doping side of the transition $ρ_c(T)$ is metallic, while on the low-doping side $ρ_c(T)$ changes from metallic to semi-conducting behavior with decreasing $T$. As a function of doping, the jump in $ρ_c$ across the first-order transition evolves into a sharp crossover at higher temperatures. This crossover coincides with the pseudogap temperature $T^*$ in the single-particle density of states, the spin susceptibility and other observables. Such coincidence in crossovers is expected along the continuation of the first-order transition into the super-critical regime, called the Widom line. This implies that not only the dynamic and the thermodynamic properties but also the DC transport in the normal state are governed by the hidden first-order transition. $ρ_c(T)$ has a high-temperature quasi-linear regime where it can exceed the Mott-Ioffe-Regel limit and when it has a minimum it is nearly parallel to the Widom line.

preprint2012arXiv

Entropy, frustration and large thermopower of doped Mott insulators on the fcc lattice

Electronic frustration and strong correlations may lead to large Seebeck coefficients. To understand this physics on general grounds, we compute the thermopower of the one-band Hubbard model on the 3-dimensional fcc lattice over the whole range of fillings for intermediate and large interaction strength. Dynamical mean-field theory shows that when the density approaches half-filling, the fcc lattice at strong coupling exhibits a large low temperature Seebeck coefficient $S$. The largest effect occurs as one approaches $n=1$ from dopings where electronic frustration is maximized. The high-frequency limit of the thermopower and the Kelvin limit are both used to provide physical insight as well as practical tools to estimate the thermopower. The high-frequency limit gives a reliable estimate of the DC limit at low temperature when the metal becomes coherent. By contrast, the Kelvin approach is useful in the strongly interacting case at high temperature when transport is incoherent. The latter result shows that in doped Mott insulators at high temperature and strong coupling the thermopower can be understood on entropic grounds.

preprint2012arXiv

Mott p-n Junctions in layered materials

The p-n junction has provided the basis for the semiconductor-device industry. Investigations of p-n junctions based on Mott insulators is still in its infancy. Layered Mott insulators, such as the cuprates or other transition metal-oxides, present a special challenge since strong in-plane correlations are important. Here we model the planes carefully using plaquette Cellular Dynamical Mean Field Theory with an exact diagonalization solver. The energy associated with inter-plane hopping is neglected compared with the long-range Coulomb interaction that we treat in the Hartree-Fock approximation. Within this new approach, "Dynamical Layer Theory", the charge redistribution is obtained at the final step from minimization of a function of the layer fillings. A simple analytical description of the solution, in the spirit of Thomas-Fermi theory, reproduces quite accurately the numerical results. Various interesting charge reconstructions can be obtained by varying the Fermi energy differences between both sides of the junction. One can even obtain quasi-two dimensional charge carriers at the interface, in the middle of a Mott insulating layer. The density of states as a function of position does not follow the simple band bending picture of semiconductors.

preprint2012arXiv

Optimal doping and entropic origin of giant thermopower in doped Mott insulators

This paper has been withdrawn by the authors because the results were incorrect due to the use of an impurity solver that fails at large interaction strength and above half-filling. This led, first, to the development of a new impurity solver (arXiv:1202.5814 or Phys.Rev. B 86, 085133 (2012)) and, second, to a complete revision of the thermopower results with different conclusions. Therefore a new paper on thermopower was submitted (arXiv:1209.4349) and the present one withdrawn.

preprint2012arXiv

Pseudogap temperature as a Widom line in doped Mott insulators

The pseudogap refers to an enigmatic state of matter with unusual physical properties found below a characteristic temperature $T^*$ in hole-doped high-temperature superconductors. Determining $T^*$ is critical for understanding this state. Here we study the simplest model of correlated electron systems, the Hubbard model, with cluster dynamical mean-field theory to find out whether the pseudogap can occur solely because of strong coupling physics and short nonlocal correlations. We find that the pseudogap characteristic temperature $T^*$ is a sharp crossover between different dynamical regimes along a line of thermodynamic anomalies that appears above a first-order phase transition, the Widom line. The Widom line emanating from the critical endpoint of a first-order transition is thus the organizing principle for the pseudogap phase diagram of the cuprates. No additional broken symmetry is necessary to explain the phenomenon. Broken symmetry states appear in the pseudogap and not the other way around.

preprint2012arXiv

Resilience of d-wave superconductivity to nearest-neighbor repulsion

Many theoretical approaches find d-wave superconductivity in the prototypical one-band Hubbard model for high-temperature superconductors. At strong-coupling (U > W, where U is the on-site repulsion and W=8t the bandwidth) pairing is controlled by the exchange energy J=4t^2/U. One may then surmise, ignoring retardation effects, that near-neighbor Coulomb repulsion V will destroy superconductivity when it becomes larger than J, a condition that is easily satisfied in cuprates for example. Using Cellular Dynamical Mean-Field theory with an exact diagonalization solver for the extended Hubbard model, we show that pairing, at strong coupling, is preserved even when V>>J, as long as V<U/2. While at weak coupling V always reduces the spin fluctuations and hence d-wave pairing, at strong coupling, in the underdoped regime, the increase of J=4t^2/(U-V) caused by V increases binding at low frequency while the pair-breaking effect of V is pushed to high frequency. These two effects compensate in the underdoped regime, in the presence of a pseudogap. While the pseudogap competes with superconductivity, the proximity to the Mott transition that leads to the pseudogap, and retardation effects, protect d-wave superconductivity from V.

preprint2012arXiv

Strong coupling superconductivity, pseudogap and Mott transition

An intricate interplay between superconductivity, pseudogap and Mott transition, either bandwidth driven or doping driven, occurs in materials. Layered organic conductors and cuprates offer two prime examples. We provide a unified perspective of this interplay in the two-dimensional Hubbard model within cellular dynamical mean-field theory on a $2\times 2$ plaquette and using the continuous-time quantum Monte Carlo method as impurity solver. Both at half filling and at finite doping, the metallic normal state close to the Mott insulator is unstable to d-wave superconductivity. Superconductivity can destroy the first-order transition that separates the pseudogap phase from the overdoped metal, yet that normal state transition leaves its marks on the dynamic properties of the superconducting phase. For example, as a function of doping one finds a rapid change in the particle-hole asymmetry of the superconducting density of states. In the doped Mott insulator, the dynamical mean-field superconducting transition temperature $T_c^d$ does not scale with the order parameter when there is a normal-state pseudogap. $T_c^d$ corresponds to the local pair formation temperature observed in tunneling experiments and is distinct from the pseudogap temperature.

preprint2011arXiv

d-Wave superconductivity on the checkerboard Hubbard model at weak and strong coupling

It has been argued that inhomogeneity generally can enhance superconductivity in the cuprate high-Tc materials. To check the validity of this claim, we study d-wave superconductivity on the checkerboard Hubbard model on a square lattice using the Cellular Dynamical Mean Field theory method with an exact diagonalization solver at zero temperature. The d-wave order parameter is computed for various inhomogeneity levels over the entire doping range of interest in both strong and weak coupling regimes. At a given doping, the size of the d-wave order parameter manifests itself directly in the height of the coherence peaks and hence is an appropriate measure of the strength of superconductivity. The weak coupling results reveal a suppression of the order parameter in the presence of inhomogeneity for small to intermediate hole dopings, while it is enhanced for large dopings. In contrast, for strong coupling there is a monotonic decrease in the maximum amplitude of the superconducting order parameter with inhomogeneity over the entire doping range of interest. Furthermore, at moderately high inhomogeneity, the system undergoes a first-order transition from the superconducting to the normal state in the underdoped regime. In the overdoped regime, the change in the value of the superconducting order parameter correlates with the height of the lowest energy peak in the spectral weight of antiferromagnetic spin fluctuations, confirming the connection between antiferromagnetic fluctuations and d-wave superconductivity found in earlier studies on the homogeneous case. Our results are benchmarked by comparisons with numerically exact results on the checkerboard Hubbard ladder.

preprint2011arXiv

Importance of subleading corrections for the Mott critical point

The interaction-induced metal-insulator transition should be in the Ising universality class. Experiments on layered organic superconductors suggest that the observed critical endpoint of the first-order Mott transition belongs instead to a different universality class. To address this question, we use dynamical mean-field theory and a cluster generalization that is necessary to account for short-range spatial correlations in two dimensions. Such calculations can give information on crossover effects, in particular quantum ones, that are not included in the simplest mean-field. In the cluster calculation, a canonical transformation that minimizes the sign problem in continuous-time quantum Monte Carlo calculations allows us to obtain very accurate results for double occupancy. These results show that there are important subleading corrections that can lead to apparent exponents that are different from mean-field. Experiments on optical lattices could verify our predictions.

preprint2011arXiv

Mott physics and first-order transition between two metals in the normal state phase diagram of the two-dimensional Hubbard model

For doped two-dimensional Mott insulators in their normal state, the challenge is to understand the evolution from a conventional metal at high doping to a strongly correlated metal near the Mott insulator at zero doping. To this end, we solve the cellular dynamical mean-field equations for the two-dimensional Hubbard model using a plaquette as the reference quantum impurity model and continuous-time quantum Monte Carlo method as impurity solver. The normal-state phase diagram as a function of interaction strength $U$, temperature $T$, and filling $n$ shows that, upon increasing $n$ towards the Mott insulator, there is a surface of first-order transition between two metals at nonzero doping. That surface ends at a finite temperature critical line originating at the half-filled Mott critical point. Associated with this transition, there is a maximum in scattering rate as well as thermodynamic signatures. These findings suggest a new scenario for the normal-state phase diagram of the high temperature superconductors. The criticality surmised in these systems can originate not from a T=0 quantum critical point, nor from the proximity of a long-range ordered phase, but from a low temperature transition between two types of metals at finite doping. The influence of Mott physics therefore extends well beyond half-filling.

preprint2011arXiv

Optical and DC conductivity of the two-dimensional Hubbard model in the pseudogap regime and across the antiferromagnetic quantum critical point, including vertex corrections

The conductivity of the two-dimensional Hubbard model is particularly relevant for high-temperature superconductors. Vertex corrections are expected to be important because of strongly momentum dependent self-energies. We use the Two-Particle Self-Consistent approach that satisfies crucial constraints such as the Mermin-Wagner theorem, the Pauli principle and sum rules in order to reach non-perturbative regimes. This approach is reliable from weak to intermediate coupling. A functional derivative approach ensures that vertex corrections are included in a way that satisfies the f sum-rule. The two types of vertex corrections that we find are the antiferromagnetic analogs of the Maki-Thompson and Aslamasov-Larkin contributions of superconducting fluctuations to the conductivity but, contrary to the latter, they include non-perturbative effects. The resulting analytical expressions must be evaluated numerically. The calculations are impossible unless a number of advanced numerical algorithms are used. A maximum entropy approach is specially developed for analytical continuation of our results. The numerical results are for nearest neighbor hoppings. In the pseudogap regime induced by two-dimensional antiferromagnetic fluctuations, the effect of vertex corrections is dramatic. Without vertex corrections the resistivity increases as we enter the pseudogap regime. Adding vertex corrections leads to a drop in resistivity, as observed in some high temperature superconductors. At high temperature, the resistivity saturates at the Ioffe-Regel limit. At the quantum critical point and beyond, the resistivity displays both linear and quadratic temperature dependence and there is a correlation between the linear term and the superconducting transition temperature. A hump is observed in the mid-infrared range of the optical conductivity in the presence of antiferromagnetic fluctuations.

preprint2011arXiv

Two-Particle-Self-Consistent Approach for the Hubbard Model

Even at weak to intermediate coupling, the Hubbard model poses a formidable challenge. In two dimensions in particular, standard methods such as the Random Phase Approximation are no longer valid since they predict a finite temperature antiferromagnetic phase transition prohibited by the Mermin-Wagner theorem. The Two-Particle-Self-Consistent (TPSC) approach satisfies that theorem as well as particle conservation, the Pauli principle, the local moment and local charge sum rules. The self-energy formula does not assume a Migdal theorem. There is consistency between one- and two-particle quantities. Internal accuracy checks allow one to test the limits of validity of TPSC. Here I present a pedagogical review of TPSC along with a short summary of existing results and two case studies: a) the opening of a pseudogap in two dimensions when the correlation length is larger than the thermal de Broglie wavelength, and b) the conditions for the appearance of d-wave superconductivity in the two-dimensional Hubbard model.

preprint2010arXiv

Dynamical electronic nematicity from Mott physics

Very large anisotropies in transport quantities have been observed in the presence of very small in-plane structural anisotropy in many strongly correlated electron materials. By studying the two-dimensional Hubbard model with dynamical-mean-field theory for clusters, we show that such large anisotropies can be induced without static stripe order if the interaction is large enough to yield a Mott transition. Anisotropy decreases at large frequency. The maximum effect on conductivity anisotropy occurs in the underdoped regime, as observed in high temperature superconductors.

preprint2010arXiv

Finite doping signatures of the Mott transition in the two-dimensional Hubbard model

Experiments on layered materials call for a study of the influence of short-range spin correlations on the Mott transition. To this end, we solve the cluster dynamical mean-field equations for the Hubbard model on a plaquette with continuous-time quantum Monte Carlo. The normal state phase diagram as a function of temperature $T$, interaction strength $U$ and filling $n$ reveals that upon increasing $n$ towards the insulator, there is a surface of first-order transition between two metals at non-zero doping. For $T$ above the critical end line there is a maximum in scattering rate.

preprint2010arXiv

Generalized DC and AC Josephson effects in antiferromagnets and in antiferromagnetic d-wave superconductors

The Josephson effect is generally described as Cooper pair tunneling, but it can also be understood in a more general context. The DC Josephson effect is the pseudo-Goldstone boson of two coupled systems with a broken continuous abelian U(1) symmetry. Hence, an analog should exist for systems with broken continuous non-Abelian symmetries. To exhibit the generality of the phenomenon and make predictions from a realistic model, we study tunneling between antiferromagnets and also between antiferromagnetic d-wave superconductors. Performing a calculation analogous to that of Ambegaokar and Baratoff for the Josephson junction, we find an equilibrium current of the staggered magnetization through the junction that, in antiferromagnets, is proportional to s_L X s_R where s_L and s_R are the Neel vectors on either sides of the junction. Microscopically, this effect exists because of the coherent tunneling of spin-one particle-hole pairs. In the presence of a magnetic field which is different on either sides of the junction, we find an analog of the AC Josephson effect where the angle between Neel vectors depends on time. In the case of antiferromagnetic d-wave superconductors we predict that there is a contribution to the critical current that depends on the antiferromagnetic order and a contribution to the spin-critical current that depends on superconducting order. The latter contributions come from tunneling of the triplet Cooper pair that is necessarily present in the ground state of an antiferromagnetic d-wave superconductor. All these effects appear to leading order in the square of the tunneling matrix elements.

preprint2009arXiv

Mixed pairing symmetry in κ-(BEDT-TTF)_2 X organic superconductors from ultrasonic velocity measurements

Discontinuities in elastic constants are detected at the superconducting transition of layered organic conductors κ-(BEDT-TTF)_{2}X by longitudinal and transverse ultrasonic velocity measurements. Symmetry arguments show that discontinuities in shear elastic constants can be explained in the orthorhombic compound only if the superconducting order parameter has a mixed character that can be of two types, either A_{1g}+B_{1g} or B_{2g}+B_{3g} in the classification of irreducible representations of the orthorhombic point group D_{2h}. Consistency with other measurements suggests that the A_{1g}+B_{1g} (d_{xy}+d_{z(x+y)}) possibility is realized. Such clear symmetry-imposed signatures of mixed order parameters have not been observed in other superconducting compounds.

preprint2009arXiv

Pairing dynamics in strongly correlated superconductivity

Confirmation of the phononic origin of Cooper pair formation in superconductors came with the demonstration that the interaction was retarded and that the corresponding energy scales were associated with phonons. Using cellular dynamical mean-field theory for the two-dimensional Hubbard model, we identify such retardation effects in d-wave pairing and associate the corresponding energy scales with short-range spin fluctuations. We find which frequencies are relevant for pairing as a function of interaction strength and doping and show that the disappearance of superconductivity on the overdoped side coincides with the disappearance of the low energy feature in the antiferromagnetic fluctuations, as observed in neutron scattering experiments.

preprint2007arXiv

Supersolidity, entropy and frustration

We study the properties of t-t'-V model of hard-core bosons on the triangular lattice that can be realized in optical lattices. By mapping to the spin-1/2 XXZ model in a field, we determine the phase diagram of the t-V model where the supersolid characterized by the ordering pattern (x,x,-2x') ("ferrimagnetic" or SS A) is a ground state for chemical potential μ>3V. By turning on either temperature or t' at half-filling μ=3V, we find a first order transition from SS A to the elusive supersolid characterized by the (x,-x,0) ordering pattern ("antiferromagnetic" or SS C). In addition, we find a large region where a superfluid phase becomes a solid upon raising temperature at fixed chemical potential. This is an analog of the Pomeranchuk effect driven by the large entropic effects associated with geometric frustration on the triangular lattice.

preprint2006arXiv

Pseudogap induced by short-range spin correlations in a doped Mott insulator

We study the evolution of a Mott-Hubbard insulator into a correlated metal upon doping in the two-dimensional Hubbard model using the Cellular Dynamical Mean Field Theory. Short-range spin correlations create two additional bands apart from the familiar Hubbard bands in the spectral function. Even a tiny doping into this insulator causes a jump of the Fermi energy to one of these additional bands and an immediate momentum dependent suppression of the spectral weight at this Fermi energy. The pseudogap is closely tied to the existence of these bands. This suggests a strong-coupling mechanism that arises from short-range spin correlations and large scattering rates for the pseudogap phenomenon seen in several cuprates.

preprint1996arXiv

Crossover from two- to three-dimensional critical behavior for nearly antiferromagnetic itinerant electrons

The crossover from two- to three-dimensional critical behavior of nearly antiferromagnetic itinerant electrons is studied in a regime where the inter-plane single-particle motion of electrons is quantum-mechanically incoherent because of thermal fluctuations. This is a relevant regime for very anisotropic materials like the cuprates. The problem is studied within the Two-Particle Self-Consistent approach (TPSC), that has been previously shown to give a quantitative description of Monte Carlo data for the Hubbard model. It is shown that TPSC belongs to the $n\rightarrow \infty $ limit of the $O\left( n\right) $ universality class. However, contrary to the usual approaches, cutoffs appear naturally in the microscopic TPSC theory so that parameter-free calculations can be done for Hubbard models with arbitrary band structure. A general discussion of universality in the renormalized-classical crossover from $d=2$ to $d=3$ is also given.

preprint1995arXiv

Magnetic and pair correlations of the Hubbard model with next-nearest-neighbor hopping

A combination of analytical approaches and quantum Monte Carlo simulations is used to study both magnetic and pairing correlations for a version of the Hubbard model that includes second-neighbor hopping $t^{\prime }=-0.35t$ as a model for high-temperature superconductors. Magnetic properties are analyzed using the Two-Particle Self-Consistent approach. The maximum in magnetic susceptibility as a function of doping appears both at finite $% t^{\prime }$ and at $t^{\prime }=0$ but for two totally different physical reasons. When $t^{\prime }=0$, it is induced by antiferromagnetic correlations while at $t^{\prime }=-0.35t$ it is a band structure effect amplified by interactions. Finally, pairing fluctuations are compared with $% T $-matrix results to disentangle the effects of van Hove singularity and of nesting on superconducting correlations. The addition of antiferromagnetic fluctuations increases slightly the $d$-wave superconducting correlations despite the presence of a van Hove singularity which tends to decrease them in the repulsive model. Some aspects of the phase diagram and some subtleties of finite-size scaling in Monte Carlo simulations, such as inverted finite-size dependence, are also discussed.

preprint1995arXiv

Theory of single-particle properties of the Hubbard model

It is shown that it is possible to quantitatively explain quantum Monte Carlo results for the Green's function of the two-dimensional Hubbard model in the weak to intermediate coupling regime. The analytic approach includes vertex corrections in a paramagnon-like self-energy. All parameters are determined self-consistently. This approach clearly shows that in two dimensions Fermi-liquid quasiparticles disappear in the paramagnetic state when the antiferromagnetic correlation length becomes larger than the electronic thermal de Broglie wavelength.