Strong Coupling Superconductivity, Pseudogap, and Mott Transition

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.

With high-Tc cuprates in mind, properties of correlated dx2-y2-wave superconducting (SC) and antiferromagnetic (AF) states are studied for the Hubbard (t-t'-U) model on square lattices, using a variational Monte Carlo method. We employ simple trial wave functions including only crucial parameters, such as a doublon-holon binding factor indispensable to describe correlated SC and normal states as doped Mott insulators. U/t, t'/t and \delta (doping rate) dependence of relevant quantities are systematically calculated. As U/t increases, a sharp crossover of SC properties occurs at U_co/t \sim 10 from a conventional BCS type to a kinetic-energy-driven type for any t'/t. As \delta decreases, U_co/t is smoothly connected to the Mott transition point at half filling. For U/t\lsim 5, steady superconductivity corresponding to the cuprates is not found, whereas the d-wave SC correlation function Pd^\infty rapidly increases for U/t\gsim 6 and becomes maximum at U=U_co. Comparing the \delta dependence of Pd^\infty with experimentally observed dome-shaped Tc and condensation energy, we find that the effective value of $U$ for the cuprates should be larger than the band width, for which the t-J model is valid. Analyzing the kinetic energy, we reveal that for U>U_co only doped holes (electrons) become charge carriers, which will make a small Fermi surface (hole pocket), but for U<U_co all the electrons (holes) contribute to conduction and will make an ordinary large Fermi surface, which is contradictory to the feature of cuprates. By introducing a proper negative (positive) t'/t, the SC (AF) state is stabilized. In the underdoped regime, the strength of SC for U>U_co is determined by two factors, i.e., the AF spin correlation, which creates singlet pairs (pseudogap), and the charge mobility dominated by Mott physics.

A cellular dynamical mean field theory approach to Mottness

A half-filled-band Hubbard model on an anisotropic triangular lattice ( t in two bond directions and t ' in the other) is studied using an optimization variational Monte Carlo method, to consider a Mott transition and superconductivity arising in κ-(BEDT-TTF) 2 X. Adopting wave functions with doublon–holon binding factors, we reveal that a first-order Mott (conductor-to-nonmagnetic insulator) transition takes place at U = U c , approximately of the bandwidth, for a wide t '/ t range. This transition is not directly connected to magnetism. Robust d -wave superconductivity appears in a restricted parameter range: immediately below U c and the moderate strength of frustration (\(0.4 {\lesssim} t'/t {\lesssim} 0.7\)), where short-range antiferromagnetic correlation sufficiently develops but does not come to a long-range order. The relevance to experiments is also discussed.

The newly discovered 122-iron-chalcogenide superconductors show clear signatures of Mott-insulating and strange, incoherent metal normal states. Motivated thereby, we use local-density approximation plus dynamical mean-field theory (LDA + DMFT) to study the electronic structure and transport properties of K${}_{x}$Fe${}_{2\ensuremath{-}y}$Se${}_{2}$. We find that the undoped KFe${}_{1.6}$Se${}_{2}$ system is a new kind of Mott-Kondo insulator (MKI). Electron doping this MKI drives a Mott transition to an orbital-selective non-Fermi-liquid metal. Good agreement with spectral and transport responses supports our view, implying that superconductivity arises from a doped Mott insulator, as in the high-${T}_{c}$ cuprates.

The single-particle spectral properties near the Mott transition in the two-dimensional Hubbard model with next-nearest-neighbor hopping are investigated by using cluster perturbation theory. Complicated spectral features of this model are simply interpreted, by considering how the next-nearest-neighbor hopping shifts the spectral weights of the two-dimensional Hubbard model. Various anomalous features observed in hole-doped and electron-doped cuprate high-temperature superconductors are explained in a unified manner as properties near the Mott transition in a two-dimensional system whose spectral weights are shifted by next-nearest-neighbor hopping.

A hole injected into a Mott insulator will gain an internal structure as recently identified by exact numerics, which is characterized by a nontrivial quantum number whose nature is of central importance in understanding the Mott physics. In this work, we show that a spin texture associated with such an internal degree of freedom can explicitly manifest after the spin degeneracy is lifted by a \emph{weak} Rashba spin-orbit coupling (SOC). It is described by an emergent angular momentum $J_{z}=\pm3/2$ as shown by both exact diagonalization (ED) and variational Monte Carlo (VMC) calculations, which are in good agreement with each other at a finite size. In particular, as the internal structure such a spin texture is generally present in the hole composite even at high excited energies, such that a corresponding texture in momentum space, extending deep inside the Brillouin zone, can be directly probed by the spin-polarized angle-resolved photoemission spectroscopy (ARPES). This is in contrast to a Landau quasiparticle under the SOC, in which the spin texture induced by SOC will not be protected once the excited energy is larger than the weak SOC coupling strength, away from the Fermi energy. We point out that the spin texture due to the SOC should be monotonically enhanced 1with reducing spin-spin correlation length in the superconducting/pseudogap phase at finite doping. A brief discussion of a recent experiment of the spin-polarized ARPES will be made.

A correlated material in the vicinity of an insulator-metal transition (IMT) exhibits rich phenomenology and a variety of interesting phases. A common avenue to induce IMTs in Mott insulators is doping, which inevitably leads to disorder. While disorder is well known to create electronic inhomogeneity, recent theoretical studies have indicated that it may play an unexpected and much more profound role in controlling the properties of Mott systems. Theory predicts that disorder might play a role in driving a Mott insulator across an IMT, with the emergent metallic state hosting a power-law suppression of the density of states (with exponent close to 1; V-shaped gap) centered at the Fermi energy. Such V-shaped gaps have been observed in Mott systems, but their origins are as-yet unknown. To investigate this, we use scanning tunneling microscopy and spectroscopy to study isovalent Ru substitutions in Sr3(Ir1-xRux)2O7 (0 ≤ x ≤ 0.5) which drive the system into an antiferromagnetic, metallic state. Our experiments reveal that many core features of the IMT, such as power-law density of states, pinning of the Fermi energy with increasing disorder, and persistence of antiferromagnetism, can be understood as universal features of a disordered Mott system near an IMT and suggest that V-shaped gaps may be an inevitable consequence of disorder in doped Mott insulators.

Antiferromagnetism and $d$-wave superconductivity are the most important competing ground-state phases of cuprate superconductors. Using cellular dynamical mean-field theory (CDMFT) for the Hubbard model, we revisit the question of the coexistence and competition of these phases in the one-band Hubbard model with realistic band parameters and interaction strengths. With an exact diagonalization solver, we improve on previous works with a more complete bath parametrization which is carefully chosen to grant the maximal possible freedom to the hybridization function for a given number of bath orbitals. Compared with previous incomplete parametrizations, this general bath parametrization shows that the range of microscopic coexistence of superconductivity and antiferromagnetism is reduced for band parameters for NCCO, and confined to electron-doping with parameters relevant for YBCO.

We study the non-equilibrium dynamics of a simple model for V2O3 that consists of a quarter-filled Hubbard model for two orbitals that are split by a weak crystal field. Peculiarities of this model are: (1) a Mott insulator whose gap corresponds to transferring an electron from the occupied lower orbital to the empty upper one, rather than from the lower to the upper Hubbard sub-bands; (2) a Mott transition generically of first order even at zero temperature. We simulate by means of time-dependent Gutzwiller approximation the evolution within the insulating phase of an initial state endowed by a non-equilibrium population of electrons in the upper orbital and holes in the lower one. We find that the excess population may lead, above a threshold, to a gap-collapse and drive the insulator into the metastable metallic phase within the coexistence region around the Mott transition. This result foresees a non-thermal pathway to revert a Mott insulator into a metal. Even though this physical scenario is uncovered in a very specific toy-model, we argue it might apply to other Mott insulating materials that share similar features.

We analyze cellular dynamical mean-field theory (CDMFT) and the dynamical cluster approximation (DCA). We derive exact sum-rules for the hybridization functions and give examples for DMFT, CDMFT, and DCA. For impurity solvers based on a Hamiltonian, these sum-rules can be used to monitor convergence of the bath-parametrization. We further discuss how the symmetry of the cluster naturally leads to a decomposition of the bath Green matrix into irreducible components, which can be parametrized independently, and give an explicit recipe for finding the optimal bath-parametrization. As a benchmark we revisit the one-dimensional Hubbard model. We carefully analyze the evolution of the density as a function of chemical potential and find that, close to the Mott transition, convergence with cluster size is unexpectedly slow. In two dimensions we find, that we need so many bath-sites to obtain a reliable parametrization that Lanczos calculations are hardly feasible with current computers. For such large baths our symmetry-adapted approach should prove crucial for finding a reliable bath-parametrization.

We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest coupled cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.

Spin response and weak pseudogap behavior in electron-doped cobaltates

Pyqcm is a Python/C++ library that implements a few quantum cluster methods with an exact diagonalization impurity solver. Quantum cluster methods are used in the study of strongly correlated electrons to provide an approximate solution to Hubbard-like models. The methods covered by this library are Cluster Perturbation Theory (CPT), the Variational Cluster Approach (VCA) and Cellular (or Cluster) Dynamical Mean Field Theory (CDMFT). The impurity solver (the technique used to compute the cluster’s interacting Green function) is exact diagonalization from sparse matrices, using the Lanczos algorithm and variants thereof. The core library is written in C++ for performance, but the interface is in Python, for ease of use and inter-operability with the numerical Python ecosystem. The library is distributed under the GPL license.

Pyqcm is a Python/C++ library that implements a few quantum cluster methods with an exact diagonalization impurity solver. Quantum cluster methods are used in the study of strongly correlated electrons to provide an approximate solution to Hubbard-like models. The methods covered by this library are Cluster Perturbation Theory (CPT), the Variational Cluster Approach (VCA) and Cellular (or Cluster) Dynamical Mean Field Theory (CDMFT). The impurity solver (the technique used to compute the cluster’s interacting Green function) is exact diagonalization from sparse matrices, using the Lanczos algorithm and variants thereof. The core library is written in C++ for performance, but the interface is in Python, for ease of use and inter-operability with the numerical Python ecosystem. The library is distributed under the GPL license.