Superconductivity and the pseudogap in the two-dimensional Hubbard model

Sensitive dependence of pairing symmetry on Ni-eg crystal field splitting in the nickelate superconductor La3Ni2O7

The discovery of high-temperature superconductivity in La3Ni2O7 under pressure has drawn great attention. However, consensus has not been reached on its pairing symmetry in theory. By combining density-functional-theory (DFT), maximally-localized-Wannier-function, and linearized gap equation with random-phase-approximation, we find that the pairing symmetry of La3Ni2O7 is dxy, if its DFT band structure is accurately reproduced by a downfolded bilayer two-orbital model. More importantly, we reveal that the pairing symmetry of La3Ni2O7 sensitively depends on the crystal field splitting between two Ni-eg orbitals. A slight increase in Ni-eg crystal field splitting alters the pairing symmetry from dxy to s±. Such a transition is associated with the change in inverse Fermi velocity and susceptibility, while the shape of Fermi surface remains almost unchanged. Our work highlights the sensitive dependence of pairing symmetry on low-energy electronic structures in multi-orbital superconductors, which calls for care in the downfolding procedure when one calculates their pairing symmetry.

Superconductivity and Mott Physics in Organic Charge Transfer Materials

The phase diagrams of quasi two-dimensional organic superconductors display a plethora of fundamental phenomena associated with strong electron correlations, such as unconventional superconductivity, metal-insulator transitions, frustrated magnetism and spin liquid behavior. We analyze a minimal model for these compounds, the Hubbard model on an anisotropic triangular lattice, using cutting-edge quantum embedding methods respecting the lattice symmetry. We demonstrate the existence of unconventional superconductivity by directly entering the symmetry-broken phase. We show that the crossover from the Fermi liquid metal to the Mott insulator is associated with the formation of a pseudogap. The predicted momentum-selective destruction of the Fermi surface into hot and cold regions provides motivation for further spectroscopic studies. Our theoretical results agree with experimental phase diagrams of κ-BEDT organics.

Origin and fate of the pseudogap in the doped Hubbard model

The relationship between the pseudogap and underlying ground-state phases has not yet been rigorously established. We investigated the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations, allowing for the computation of spectral properties in the infinite-size limit and with arbitrary momentum resolution. We found three distinct regimes as a function of doping and interaction strength: a weakly correlated metal, a correlated metal with strong interaction effects, and a pseudogap regime at low doping. We show that the pseudogap forms both at weak coupling, when the magnetic correlation length is large, and at strong coupling, when it is shorter. As the temperature goes to zero, the pseudogap regime extrapolates precisely to the ordered stripe phase found by ground-state methods.

Unambiguous Fluctuation Decomposition of the Self-Energy: Pseudogap Physics beyond Spin Fluctuations

Correlated electron systems may give rise to multiple effective interactions whose combined impact on quasiparticle properties can be difficult to disentangle. We introduce an unambiguous decomposition of the electronic self-energy which allows us to quantify the contributions of various effective interactions simultaneously. We use this tool to revisit the hole-doped Hubbard model within the dynamical cluster approximation, where commonly spin fluctuations are considered to be the origin of the pseudogap. While our fluctuation decomposition confirms that spin fluctuations indeed suppress antinodal electronic spectral weight, we show that they alone cannot capture the pseudogap self-energy quantitatively. Nonlocal multiboson Feynman diagrams yield substantial contributions and are needed for a quantitative description of the pseudogap.

Selected Nonorthogonal Configuration Interaction with Compressed Single and Double Excitations

Addressing both dynamic and static correlations accurately is a primary goal in electronic structure theory. Nonorthogonal configuration interaction (NOCI) is a versatile tool for treating static correlation, offering chemical insights by combining diverse reference states. Nevertheless, achieving quantitative accuracy requires the inclusion of the missing dynamic correlation. This work introduces a framework for compressing orthogonal single and double excitations into a NOCI of a much smaller dimension. This compression is repeated with each Slater determinant in a reference NOCI, resulting in another NOCI that includes all of its single and double excitations (NOCISD), effectively recovering the missing dynamic correlations from the reference. This compressed NOCISD is further refined through a selection process using metric and energy tests (SNOCISD). We validate the effectiveness of SNOCISD through its application to the dissociation of the nitrogen molecule and the hole-doped two-dimensional Hubbard model at various interaction strengths.

Mechanism of superconductivity in the Hubbard model at intermediate interaction strength

We study the fluctuations responsible for pairing in the d-wave superconducting state of the two-dimensional Hubbard model at intermediate coupling within a cluster dynamical mean-field theory with a numerically exact quantum impurity solver. By analyzing how momentum- and frequency-dependent fluctuations generate the d-wave superconducting state in different representations, we identify antiferromagnetic fluctuations as the pairing glue of superconductivity in both the underdoped and the overdoped regime. Nevertheless, in the intermediate coupling regime, the predominant magnetic fluctuations may differ significantly from those described by conventional spin fluctuation theory.

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

The origin of the pseudogap behavior, found in many high-T_c superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a “gap-filling" rather than “gap-closing" behavior, similar to the one observed in cuprate superconductors. Our results present direct evidence of the pseudogap state, driven by superconducting fluctuations.

The origin of pseudogap in high-T_c superconductors remains a big puzzle. Here, the authors report numerical evidence of pseudogap behavior employing Quantum Monte Carlo algorithm emerging from pairing fluctuations in a quantum-critical non-Fermi liquid, similar to the pseudogap phase observed in cuprate superconductors.

Piecewise Interaction Picture Density Matrix Quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

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

SignificanceMost metals display an electron-scattering rate [Formula: see text] that follows [Formula: see text] at low temperatures, as prescribed by Fermi liquid theory. But there are important exceptions. One of the most prominent examples is the "strange" metal regime in overdoped cuprate supercondcutors, which exhibits a linear T dependence of the scattering rate [Formula: see text] that reaches a putative Planckian limit. Here, using cutting-edge computational approaches, we show that T-linear scattering rate can emerge from the overdoped Hubbard model at low temperatures. Our results agree with cuprate experiments in various aspects but challenge the Planckian limit. Finally, by identifying antiferromagnetic fluctuations as the physical origin of the T-linear scattering rate, we discover the microscopic mechanism of strange metallicity in cuprates.

The Sign Problem in Density Matrix Quantum Monte Carlo

Density matrix quantum Monte Carlo (DMQMC) is a recently developed method for stochastically sampling the N-particle thermal density matrix to obtain exact-on-average energies for model and ab initio systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product basis of Slater determinants. In principle, this means that DMQMC needs to sample a space that scales in the system size, N, as O[(exp(N))²]. In practice, removing the sign problem requires a total walker population that exceeds a system-dependent critical walker population (N_c), imposing limitations on both storage and compute time. We establish that N_c for DMQMC is the square of N_c for FCIQMC. By contrast, the minimum N_c in the interaction picture modification of DMQMC (IP-DMQMC) is only linearly related to the N_c for FCIQMC. We find that this difference originates from the difference in propagation of IP-DMQMC versus canonical DMQMC: the former is asymmetric, whereas the latter is symmetric. When an asymmetric mode of propagation is used in DMQMC, there is a much greater stochastic error and is thus prohibitively expensive for DMQMC without the interaction picture adaptation. Finally, we find that the equivalence between IP-DMQMC and FCIQMC seems to extend to the initiator approximation, which is often required to study larger systems with large basis sets. This suggests that IP-DMQMC offers a way to ameliorate the cost of moving between a Slater determinant space and an outer product basis.

Optical Indistinguishability via Twinning Fields

Here we introduce the concept of the twinning field-a driving electromagnetic pulse that induces an identical optical response from two distinct materials. We show that for a large class of pairs of generic many-body systems, a twinning field which renders the systems optically indistinguishable exists. The conditions under which this field exists are derived, and this analysis is supplemented by numerical calculations of twinning fields for both the 1D Fermi-Hubbard model, and tight-binding models of graphene and hexagonal boron nitride. The existence of twinning fields may lead to new research directions in nonlinear optics, materials science, and quantum technologies.

Valley interference and spin exchange at the atomic scale in silicon

Tunneling is a fundamental quantum process with no classical equivalent, which can compete with Coulomb interactions to give rise to complex phenomena. Phosphorus dopants in silicon can be placed with atomic precision to address the different regimes arising from this competition. However, they exploit wavefunctions relying on crystal band symmetries, which tunneling interactions are inherently sensitive to. Here we directly image lattice-aperiodic valley interference between coupled atoms in silicon using scanning tunneling microscopy. Our atomistic analysis unveils the role of envelope anisotropy, valley interference and dopant placement on the Heisenberg spin exchange interaction. We find that the exchange can become immune to valley interference by engineering in-plane dopant placement along specific crystallographic directions. A vacuum-like behaviour is recovered, where the exchange is maximised to the overlap between the donor orbitals, and pair-to-pair variations limited to a factor of less than 10 considering the accuracy in dopant positioning. This robustness remains over a large range of distances, from the strongly Coulomb interacting regime relevant for high-fidelity quantum computation to strongly coupled donor arrays of interest for quantum simulation in silicon.

Coupled donor wavefunctions in silicon are spatially resolved to evidence valley interference processes. An atomic-scale understanding of the interplay between interference, envelope anisotropy and crystal symmetries unveils a placement strategy compatible with existing technology where the exchange is insensitive to interference.

Detecting quantum critical points in the t-[Formula: see text] Fermi-Hubbard model via complex network theory

A considerable success in phenomenological description of [Formula: see text] superconductors has been achieved within the paradigm of Quantum Critical Point (QCP)-a parental state of a variety of exotic phases that is characterized by dense entanglement and absence of well-defined quasiparticles. However, the microscopic origin of the critical regime in real materials remains an open question. On the other hand, there is a popular view that a single-band t-[Formula: see text] Hubbard model is the minimal model to catch the main relevant physics of superconducting compounds. Here, we suggest that emergence of the QCP is tightly connected with entanglement in real space and identify its location on the phase diagram of the hole-doped t-[Formula: see text] Hubbard model. To detect the QCP we study a weighted graph of inter-site quantum mutual information within a four-by-four plaquette that is solved by exact diagonalization. We demonstrate that some quantitative characteristics of such a graph, viewed as a complex network, exhibit peculiar behavior around a certain submanifold in the parametric space of the model. This method allows us to overcome difficulties caused by finite size effects and to identify precursors of the transition point even on a small lattice, where long-range asymptotics of correlation functions cannot be accessed.

Spin and Charge Correlations across the Metal-to-Insulator Crossover in the Half-Filled 2D Hubbard Model

The 2D Hubbard model with nearest-neighbor hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using the diagrammatic Monte Carlo approach in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and nonlocal fluctuations in both channels are essential.

Extended Crossover from a Fermi Liquid to a Quasiantiferromagnet in the Half-Filled 2D Hubbard Model

The ground state of the Hubbard model with nearest-neighbor hopping on the square lattice at half filling is known to be that of an antiferromagnetic (AFM) band insulator for any on-site repulsion. At finite temperature, the absence of long-range order makes the question of how the interaction-driven insulator is realized nontrivial. We address this problem with controlled accuracy in the thermodynamic limit using self-energy diagrammatic determinant Monte Carlo and dynamical cluster approximation methods and show that development of long-range AFM correlations drives an extended crossover from Fermi liquid to insulating behavior in the parameter regime that precludes a metal-to-insulator transition. The intermediate crossover state is best described as a non-Fermi liquid with a partially gapped Fermi surface.

Anomalies in the pseudogap phase of the cuprates: competing ground states and the role of umklapp scattering

Over the past two decades, advances in computational algorithms have revealed a curious property of the two-dimensional Hubbard model (and related theories) with hole doping: the presence of close-in-energy competing ground states that display very different physical properties. On the one hand, there is a complicated state exhibiting intertwined spin, charge, and pair density wave orders. We call this 'type A'. On the other hand, there is a uniform d-wave superconducting state that we denote as 'type B'. We advocate, with the support of both microscopic theoretical calculations and experimental data, dividing the high-temperature cuprate superconductors into two corresponding families, whose properties reflect either the type A or type B ground states at low temperatures. We review the anomalous properties of the pseudogap phase that led us to this picture, and present a modern perspective on the role that umklapp scattering plays in these phenomena in the type B materials. This reflects a consistent framework that has emerged over the last decade, in which Mott correlations at weak coupling drive the formation of the pseudogap. We discuss this development, recent theory and experiments, and open issues.

Quantum Entanglement in Deep Learning Architectures

Modern deep learning has enabled unprecedented achievements in various domains. Nonetheless, employment of machine learning for wave function representations is focused on more traditional architectures such as restricted Boltzmann machines (RBMs) and fully connected neural networks. In this Letter, we establish that contemporary deep learning architectures, in the form of deep convolutional and recurrent networks, can efficiently represent highly entangled quantum systems. By constructing tensor network equivalents of these architectures, we identify an inherent reuse of information in the network operation as a key trait which distinguishes them from standard tensor network-based representations, and which enhances their entanglement capacity. Our results show that such architectures can support volume-law entanglement scaling, polynomially more efficiently than presently employed RBMs. Thus, beyond a quantification of the entanglement capacity of leading deep learning architectures, our analysis formally motivates a shift of trending neural-network-based wave function representations closer to the state-of-the-art in machine learning.

Correlation-Driven Lifshitz Transition at the Emergence of the Pseudogap Phase in the Two-Dimensional Hubbard Model

We study the relationship between the pseudogap and Fermi-surface topology in the two-dimensional Hubbard model by means of the cellular dynamical mean-field theory. We find two possible mean-field metallic solutions on a broad range of interactions, doping, and frustration: a conventional renormalized metal and an unconventional pseudogap metal. At half filling, the conventional metal is more stable and displays an interaction-driven Mott metal-insulator transition. However, for large interactions and small doping, a region that is relevant for cuprates, the pseudogap phase becomes the ground state. By increasing doping, we show that a first-order transition from the pseudogap to the conventional metal is tied to a change of the Fermi surface from hole- to electronlike, unveiling a correlation-driven mechanism for a Lifshitz transition. This explains the puzzling link between the pseudogap phase and Fermi surface topology that has been pointed out in recent experiments.

Dynamics of correlation-frozen antinodal quasiparticles in superconducting cuprates

Many puzzling properties of high-critical temperature (Tc) superconducting (HTSC) copper oxides have deep roots in the nature of the antinodal quasiparticles, the elementary excitations with wave vector parallel to the Cu-O bonds. These electronic states are most affected by the onset of antiferromagnetic correlations and charge instabilities, and they host the maximum of the anisotropic superconducting gap and pseudogap. We use time-resolved extreme-ultraviolet photoemission with proper photon energy (18 eV) and time resolution (50 fs) to disclose the ultrafast dynamics of the antinodal states in a prototypical HTSC cuprate. After photoinducing a nonthermal charge redistribution within the Cu and O orbitals, we reveal a dramatic momentum-space differentiation of the transient electron dynamics. Whereas the nodal quasiparticle distribution is heated up as in a conventional metal, new quasiparticle states transiently emerge at the antinodes, similarly to what is expected for a photoexcited Mott insulator, where the frozen charges can be released by an impulsive excitation. This transient antinodal metallicity is mapped into the dynamics of the O-2p bands, thus directly demonstrating the intertwining between the low- and high-energy scales that is typical of correlated materials. Our results suggest that the correlation-driven freezing of the electrons moving along the Cu-O bonds, analogous to the Mott localization mechanism, constitutes the starting point for any model of high-Tc superconductivity and other exotic phases of HTSC cuprates.

Coherent organization of electronic correlations as a mechanism to enhance and stabilize high-T C cuprate superconductivity

Strong diffusive or incoherent electronic correlations are the signature of the strange-metal normal state of the cuprate superconductors, with these correlations considered to be undressed or removed in the superconducting state. A critical question is if these correlations are responsible for the high-temperature superconductivity. Here, utilizing a development in the analysis of angle-resolved photoemission data, we show that the strange-metal correlations don't simply disappear in the superconducting state, but are instead converted into a strongly renormalized coherent state, with stronger normal state correlations leading to stronger superconducting state renormalization. This conversion begins well above T C at the onset of superconducting fluctuations and it greatly increases the number of states that can pair. Therefore, there is positive feedback--the superconductive pairing creates the conversion that in turn strengthens the pairing. Although such positive feedback should enhance a conventional pairing mechanism, it could potentially also sustain an electronic pairing mechanism.

Fierz Convergence Criterion: A Controlled Approach to Strongly Interacting Systems with Small Embedded Clusters

We present an embedded-cluster method, based on the triply irreducible local expansion formalism. It turns the Fierz ambiguity, inherent to approaches based on a bosonic decoupling of local fermionic interactions, into a convergence criterion. It is based on the approximation of the three-leg vertex by a coarse-grained vertex computed from a self-consistently determined cluster impurity model. The computed self-energies are, by construction, continuous functions of momentum. We show that, in three interaction and doping regimes of the two-dimensional Hubbard model, self-energies obtained with clusters of size four only are very close to numerically exact benchmark results. We show that the Fierz parameter, which parametrizes the freedom in the Hubbard-Stratonovich decoupling, can be used as a quality control parameter. By contrast, the GW+extended dynamical mean field theory approximation with four cluster sites is shown to yield good results only in the weak-coupling regime and for a particular decoupling. Finally, we show that the vertex has spatially nonlocal components only at low Matsubara frequencies.

Anomalous metallic behaviour in the doped spin liquid candidate κ-(ET)4Hg2.89Br8

Quantum spin liquids are exotic Mott insulators that carry extraordinary spin excitations. Therefore, when doped, they are expected to afford metallic states with unconventional magnetic excitations. Here, we report experimental results which are suggestive of a doped spin liquid with anomalous metallicity in a triangular-lattice organic conductor. The spin susceptibility is nearly perfectly scaled to that of a non-doped spin liquid insulator in spite of the metallic state. Furthermore, the charge transport that is confined in the layer at high temperatures becomes sharply deconfined on cooling, coinciding with the rapid growth of spin correlations or coherence as signified by a steep decrease in spin susceptibility. The present results substantiate the desired doped spin liquid and suggest a strange metal, in which the coherence of the underlying spin liquid promotes the deconfinement of charge from the layers while preserving the non-Fermi-liquid nature.

It is expected that introducing charge carriers into an exotic quantum spin liquid state may lead to an unconventional metal but there are no clear realizations of a metallic spin liquid. Here, the authors present a spin liquid candidate that also shows evidence of strange metal behavior.

Spin-charge transformation of lattice fermion models: duality approach for diagrammatic simulation of strongly correlated systems

I derive a dual description of lattice fermions, specifically focusing on the t-J and Hubbard models, that allow diagrammatic techniques to be employed efficiently in the strongly correlated regime, as well as for systems with a restricted Hilbert space. These constructions are based on spin-charge transformation, where the lattice fermions of the original model are mapped onto spins and spin-less fermions. This mapping can then be combined with Popov-Fedotov fermionisation, where the spins are mapped onto lattice fermions with imaginary chemical potential. The resulting models do not contain any large expansion parameters, even for strongly correlated systems. Also, they exhibit dramatically smaller corrections to the density matrix from nonlinear terms in the Hamiltonian. The combination of these two properties means that they can be addressed with diagrammatic methods, including simulation techniques based on stochastic sampling of diagrammatic expansions.

Holon Wigner Crystal in a Lightly Doped Kagome Quantum Spin Liquid

We address the problem of a lightly doped spin liquid through a large-scale density-matrix renormalization group study of the t-J model on a kagome lattice with a small but nonzero concentration δ of doped holes. It is now widely accepted that the undoped (δ=0) spin-1/2 Heisenberg antiferromagnet has a spin-liquid ground state. Theoretical arguments have been presented that light doping of such a spin liquid could give rise to a high temperature superconductor or an exotic topological Fermi liquid metal. Instead, we infer that the doped holes form an insulating charge-density wave state with one doped hole per unit cell, i.e., a Wigner crystal. Spin correlations remain short ranged, as in the spin-liquid parent state, from which we infer that the state is a crystal of spinless holons, rather than of holes. Our results may be relevant to kagome lattice herbertsmithite upon doping.

Correlation-induced superconductivity dynamically stabilized and enhanced by laser irradiation

Studies on out-of-equilibrium dynamics have paved a way to realize a new state of matter. Superconductor-like properties above room temperatures recently suggested to be in copper oxides achieved by selectively exciting vibrational phonon modes by laser have inspired studies on an alternative and general strategy to be pursued for high-temperature superconductivity. We show that the superconductivity can be enhanced by irradiating laser to correlated electron systems owing to two mechanisms: First, the effective attractive interaction of carriers is enhanced by the dynamical localization mechanism, which drives the system into strong coupling regions. Second, the irradiation allows reaching uniform and enhanced superconductivity dynamically stabilized without deteriorating into equilibrium inhomogeneities that suppress superconductivity. The dynamical superconductivity is subject to the Higgs oscillations during and after the irradiation. Our finding sheds light on a way to enhance superconductivity that is inaccessible in equilibrium in strongly correlated electron systems.

Simulation of the NMR response in the pseudogap regime of the cuprates

The pseudogap in the cuprate high-temperature superconductors was discovered as a suppression of the Knight shift and spin relaxation time measured in nuclear magnetic resonance (NMR) experiments. However, theoretical understanding of this suppression in terms of the magnetic susceptiblility of correlated itinerant fermion systems was so far lacking. Here we study the temperature and doping evolution of these quantities on the two-dimensional Hubbard model using cluster dynamical mean field theory. We recover the suppression of the Knight shift and the linear-in-T spin echo decay that increases with doping. The relaxation rate shows a marked increase as T is lowered but no indication of a pseudogap on the Cu site, and a clear downturn on the O site, consistent with experimental results on single layer materials but different from double layer materials. The consistency of these results with experiment suggests that the pseudogap is well described by strong short-range correlation effects.

Charge orders, magnetism and pairings in the cuprate superconductors

We review the recent developments in the field of cuprate superconductors with special focus on the recently observed charge order in the underdoped compounds. We introduce new theoretical developments following the study of the antiferromagnetic quantum critical point in two dimensions, in which preemptive orders in both charge and superconducting (SC) sectors emerge, that are in turn related by an SU(2) symmetry. We consider the implications of this proliferation of orders in the underdoped region, and provide a study of the type of fluctuations which characterize the SU(2) symmetry. We identify an intermediate energy scale where the SC fluctuations are dominant and argue that they are unstable towards the formation of a resonant excitonic state at the pseudogap temperature T (*). We discuss the implications of this scenario for a few key experiments.

Self-optimized superconductivity attainable by interlayer phase separation at cuprate interfaces

Stabilizing superconductivity at high temperatures and elucidating its mechanism have long been major challenges of materials research in condensed matter physics. Meanwhile, recent progress in nanostructuring offers unprecedented possibilities for designing novel functionalities. Above all, thin films of cuprate and iron-based high-temperature superconductors exhibit remarkably better superconducting characteristics (for example, higher critical temperatures) than in the bulk, but the underlying mechanism is still not understood. Solving microscopic models suitable for cuprates, we demonstrate that, at an interface between a Mott insulator and an overdoped nonsuperconducting metal, the superconducting amplitude is always pinned at the optimum achieved in the bulk, independently of the carrier concentration in the metal. This is in contrast to the dome-like dependence in bulk superconductors but consistent with the astonishing independence of the critical temperature from the carrier density x observed at the interfaces of La2CuO4 and La2-x Sr x CuO4. Furthermore, we identify a self-organization mechanism as responsible for the pinning at the optimum amplitude: An emergent electronic structure induced by interlayer phase separation eludes bulk phase separation and inhomogeneities that would kill superconductivity in the bulk. Thus, interfaces provide an ideal tool to enhance and stabilize superconductivity. This interfacial example opens up further ways of shaping superconductivity by suppressing competing instabilities, with direct perspectives for designing devices.

Pairing in a dry Fermi sea

In the traditional Bardeen–Cooper–Schrieffer theory of superconductivity, the amplitude for the propagation of a pair of electrons with momentum k and −k has a log singularity as the temperature decreases. This so-called Cooper instability arises from the presence of an electron Fermi sea. It means that an attractive interaction, no matter how weak, will eventually lead to a pairing instability. However, in the pseudogap regime of the cuprate superconductors, where parts of the Fermi surface are destroyed, this log singularity is suppressed, raising the question of how pairing occurs in the absence of a Fermi sea. Here we report Hubbard model numerical results and the analysis of angular-resolved photoemission experiments on a cuprate superconductor. In contrast to the traditional theory, we find that in the pseudogap regime the pairing instability arises from an increase in the strength of the spin–fluctuation pairing interaction as the temperature decreases rather than the Cooper log instability.

Pairing interaction appears at room temperature in traditional superconductors with a Cooper instability in the Fermi sea. Here, Maier et al. report that in the pseudogap phase of cuprate, where this instability is absent, superconductivity arises from an increase in the strength of the spin fluctuation pairing interaction as the temperature decreases.

Unconventional High-Energy-State Contribution to the Cooper Pairing in the Underdoped Copper-Oxide Superconductor HgBa_{2}Ca_{2}Cu_{3}O_{8+δ}

We study the temperature-dependent electronic B_{1g} Raman response of a slightly underdoped single crystal HgBa_{2}Ca_{2}Cu_{3}O_{8+δ} with a superconducting critical temperature T_{c}=122  K. Our main finding is that the superconducting pair-breaking peak is associated with a dip on its higher-energy side, disappearing together at T_{c}. This result reveals a key aspect of the unconventional pairing mechanism: spectral weight lost in the dip is transferred to the pair-breaking peak at lower energies. This conclusion is supported by cellular dynamical mean-field theory on the Hubbard model, which is able to reproduce all the main features of the B_{1g} Raman response and explain the peak-dip behavior in terms of a nontrivial relationship between the superconducting gap and the pseudogap.

Quantum simulation of the Hubbard model with dopant atoms in silicon

In quantum simulation, many-body phenomena are probed in controllable quantum systems. Recently, simulation of Bose-Hubbard Hamiltonians using cold atoms revealed previously hidden local correlations. However, fermionic many-body Hubbard phenomena such as unconventional superconductivity and spin liquids are more difficult to simulate using cold atoms. To date the required single-site measurements and cooling remain problematic, while only ensemble measurements have been achieved. Here we simulate a two-site Hubbard Hamiltonian at low effective temperatures with single-site resolution using subsurface dopants in silicon. We measure quasi-particle tunnelling maps of spin-resolved states with atomic resolution, finding interference processes from which the entanglement entropy and Hubbard interactions are quantified. Entanglement, determined by spin and orbital degrees of freedom, increases with increasing valence bond length. We find separation-tunable Hubbard interaction strengths that are suitable for simulating strongly correlated phenomena in larger arrays of dopants, establishing dopants as a platform for quantum simulation of the Hubbard model.

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 , the maximum of the condensation energy as a function of doping, the correlation between maximum 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.

Coexistence of Incommensurate Magnetism and Superconductivity in the Two-Dimensional Hubbard Model

We analyze the competition of magnetism and superconductivity in the two-dimensional Hubbard model with a moderate interaction strength, including the possibility of incommensurate spiral magnetic order. Using an unbiased renormalization group approach, we compute magnetic and superconducting order parameters in the ground state. In addition to previously established regions of Néel order coexisting with d-wave superconductivity, the calculations reveal further coexistence regions where superconductivity is accompanied by incommensurate magnetic order.

Hidden Fermionic Excitation Boosting High-Temperature Superconductivity in Cuprates

The dynamics of a microscopic cuprate model, namely, the two-dimensional Hubbard model, is studied with a cluster extension of the dynamical mean-field theory. We find a nontrivial structure of the frequency-dependent self-energies, which describes an unprecedented interplay between the pseudogap and superconductivity. We show that these properties are well described by quasiparticles hybridizing with (hidden) fermionic excitations, emergent from the strong electronic correlations. The hidden fermion enhances superconductivity via a mechanism distinct from a conventional boson-mediated pairing, and originates the normal-state pseudogap. Though the hidden fermion is elusive in experiments, it can solve many experimental puzzles.

Self-energy behavior away from the Fermi surface in doped Mott insulators

We analyze self-energies of electrons away from the Fermi surface in doped Mott insulators using the dynamical cluster approximation to the Hubbard model. For large onsite repulsion, U, and hole doping, the magnitude of the self-energy for imaginary frequencies at the top of the band ([Formula: see text]) is enhanced with respect to the self-energy magnitude at the bottom of the band ([Formula: see text]). The self-energy behavior at these two [Formula: see text]-points is switched for electron doping. Although the hybridization is much larger for (0, 0) than for [Formula: see text], we demonstrate that this is not the origin of this difference. Isolated clusters under a downward shift of the chemical potential, [Formula: see text], at half-filling reproduce the overall self-energy behavior at (0, 0) and [Formula: see text] found in low hole doped embedded clusters. This happens although there is no change in the electronic structure of the isolated clusters. Our analysis shows that a downward shift of the chemical potential which weakly hole dopes the Mott insulator can lead to a large enhancement of the [Formula: see text] self-energy for imaginary frequencies which is not associated with electronic correlation effects, even in embedded clusters. Interpretations of the strength of electronic correlations based on self-energies for imaginary frequencies are, in general, misleading for states away from the Fermi surface.

Superconducting fluctuations in the normal state of the two-dimensional Hubbard model

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the d_{x^{2}-y^{2}} superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half filling. We find a change in sign of the vertex contributions to d_{xy} superconductivity from repulsive near half filling to attractive at large doping. p-wave superconductivity is not found at the parameters we study, and s-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity, which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase T_{c}.

Pressure-induced Mott transition in an organic superconductor with a finite doping level

We report the pressure study of a doped organic superconductor with a Hall coefficient and conductivity measurements. We find that maximally enhanced superconductivity and a marginal-Fermi liquid appear around a certain pressure where mobile carriers increase critically, suggesting a possible quantum phase transition between strongly and weakly correlated regimes. This observation points to the presence of a criticality in Mottness for a doped Mott insulator with tunable correlation.

Pseudogap in cuprates in the loop-current ordered state

Scanning tunneling microscopy (STM) has revealed that the magnitude of the pseudo-gap in under-doped cuprates varies spatially and is correlated with disorder. The loop-current order, characterized by the anapole vector Ω, discovered in under-doped cuprates occurs in the same region of the temperature and doping as the pseudo gap observed in STM and ARPES experiments. Since translational symmetry remains unchanged in the pure limit, no gap occurs at the chemical potential. On the other hand for disorder coupling linearly to the different possible orientations of Ω, there can only be a finite temperature dependent static correlation length for the loop-current state at any temperature. This leads to formation of domains of the ordered state with different orientation and magnitude of Ω in each. For the characteristic size of the domains much larger than the Fermi-vectors [Formula: see text], the boundary of the domains leads to forward scattering of the Fermions. Such forward scattering is shown to push states near the chemical potential to energies both above and below it leading to a pseudo-gap with an angular dependence which is maximum in the [Formula: see text] directions because the single-particle energies are degenerate in these directions for all domains. The magnitude of the average gap systematically increases with the square of the average loop order parameter measured by polarized neutron scattering. This result is tested. A unique result of the gap due to forward scattering is the lack of a bump in the density of states at the 'edge' of the pseudo-gap so that the depletion of states near the chemical potential is recovered only in integration up to the edge of the band. This is also in agreement with a variety of experiments. Some predictions for further experiments are provided. Due to the finite correlation length, low frequency excitations are expected at long wavelength at all temperatures in the 'ordered' phase. Such fluctuations motionally average over the shifts in frequencies of local probes such as NMR and muon resonance expected for a truly static order.

Size, dimensionality, and strong electron correlation in nanoscience

In electronic structure theory, electron-electron repulsion is normally considered only in an average (or mean field) sense, for example, in a single Hartree-Fock determinant. This is the simple molecular orbital model, which is often a good approximation for molecules. In infinite systems, this averaging treatment leads to delocalized electronic bands, an excellent description of bulk 3D sp(3) semiconductors. However, in reality electrons try to instantaneously avoid each other; their relative motion is correlated. Strong electron-electron repulsion and correlation create new collective states and cause new femtosecond kinetic processes. This is especially true in 1D and 2D systems. The quantum size effect, a single electron property, is widely known: the band gap increases with decreasing size. This Account focuses on the experimental consequences of strong correlation. We first describe π-π* excited states in carbon nanotubes (CNTs). To obtain the spectra of individual CNTs, we developed a white-light, right-angle resonant Rayleigh scattering method. Discrete exciton transitions dominate the optical absorption spectra of both semiconducting and metallic tubes. Excitons are bound neutral excited states in which the electron and hole tightly orbit each other due to their mutual Coulomb attraction. We then describe more generally the independent roles of size and dimensionality in nanoelectronic structure, using additional examples from graphene, trans-polyacetylene chains, transition metal dichalcogenides, organic/inorganic Pb iodide perovskites, quantum dots, and pentacene van der Waals crystals. In 1D and 2D chemical systems, the electronic band structure diagram can be a poor predictor of properties if explicit correlation is not considered. One- and two-dimensional systems show quantum confinement and especially strong correlation as compared with their 3D parent systems. The Coulomb interaction is enhanced because the electrons are on the surface. One- and two-dimensional systems can exhibit essentially molecular properties even though they are infinite in size. Zero-dimensional Qdots show quantum confinement and modest electron correlation. Correlation is weak in 3D bulk semiconductors. Strongly correlated electronic states can behave as if they have fractional charge and effectively separate the spin and charge of the electron. This is apparent in the "soliton" state of polyacetylene, the fractional charge quantum Hall state of graphene, and the Luttinger electrical conductivity of metallic CNTs.

Normal-state nodal electronic structure in underdoped high-Tc copper oxides

An outstanding problem in the field of high-transition-temperature (high-Tc) superconductivity is the identification of the normal state out of which superconductivity emerges in the mysterious underdoped regime. The normal state uncomplicated by thermal fluctuations can be studied using applied magnetic fields that are sufficiently strong to suppress long-range superconductivity at low temperatures. Proposals in which the normal ground state is characterized by small Fermi surface pockets that exist in the absence of symmetry breaking have been superseded by models based on the existence of a superlattice that breaks the translational symmetry of the underlying lattice. Recently, a charge superlattice model that positions a small electron-like Fermi pocket in the vicinity of the nodes (where the superconducting gap is minimum) has been proposed as a replacement for the prevalent superlattice models that position the Fermi pocket in the vicinity of the pseudogap at the antinodes (where the superconducting gap is maximum). Although some ingredients of symmetry breaking have been recently revealed by crystallographic studies, their relevance to the electronic structure remains unresolved. Here we report angle-resolved quantum oscillation measurements in the underdoped copper oxide YBa2Cu3O6 + x. These measurements reveal a normal ground state comprising electron-like Fermi surface pockets located in the vicinity of the nodes, and also point to an underlying superlattice structure of low frequency and long wavelength with features in common with the charge order identified recently by complementary spectroscopic techniques.

Thermodynamics and magnetic properties of the anisotropic 3D Hubbard model

We study the anisotropic 3D Hubbard model with increased nearest-neighbor tunneling amplitudes along one direction using the dynamical cluster approximation and compare the results to a quantum simulation experiment of ultracold fermions in an optical lattice. We find that the short-range spin correlations are significantly enhanced in the direction with stronger tunneling amplitudes. Our results agree with the experimental observations and show that the experimental temperature is lower than the strong tunneling amplitude. We characterize the system by examining the spin correlations beyond neighboring sites and determine the distribution of density, entropy, and spin correlation in the trapped system. We furthermore investigate the dependence of the critical entropy at the Néel transition on anisotropy.

Dielectric versus magnetic pairing mechanisms in high-temperature cuprate superconductors investigated using Raman scattering

We suggest, and demonstrate, a systematic approach to the study of cuprate superconductors, namely, progressive change of ion size in order to systematically alter the interaction strength and other key parameters. R(Ba,Sr)2Cu3Oy (R={La,…,Lu,Y}) is such a system where potentially obscuring structural changes are minimal. We thereby systematically alter both dielectric and magnetic properties. Dielectric fluctuation is characterized by ionic polarizability while magnetic fluctuation is characterized by exchange interactions measurable by Raman scattering. The range of transition temperatures is 70-107 K, and we find that these correlate only with the dielectric properties, a behavior which persists with external pressure. The ultimate significance may remain to be proven, but it highlights the role of dielectric screening in the cuprates and adds support to a previously proposed novel pairing mechanism involving exchange of quantized waves of electronic polarization.

Raman-scattering measurements and theory of the energy-momentum spectrum for underdoped Bi2Sr2CaCuO(8+δ) superconductors: evidence of an s-wave structure for the pseudogap

We reveal the full energy-momentum structure of the pseudogap of underdoped high-Tc cuprate superconductors. Our combined theoretical and experimental analysis explains the spectral-weight suppression observed in the B2g Raman response at finite energies in terms of a pseudogap appearing in the single-electron excitation spectra above the Fermi level in the nodal direction of momentum space. This result suggests an s-wave pseudogap (which never closes in the energy-momentum space), distinct from the d-wave superconducting gap. Recent tunneling and photoemission experiments on underdoped cuprates also find a natural explanation within the s-wave pseudogap scenario.