Pair Susceptibility Research Articles

Density response and collective modes of semiholographic non-Fermi liquids

Semi-holographic models of non-Fermi liquids have been shown to have generically stable generalised quasi-particles on the Fermi surface. Although these excitations are broad and exhibit particle-hole asymmetry, they were argued to be stable from interactions at the Fermi surface. In this work, we use this observation to compute the density response and collective behaviour in these systems. Compared to the Fermi liquid case, we find that the boundaries of the particle-hole continuum are blurred by incoherent contributions. However, there is a region inside this continuum, that we call inner core, within which salient features of the Fermi liquid case are preserved. A particularly striking prediction of our work is that these systems support a plasmonic collective excitation which is well-defined at large momenta, has an approximately linear dispersion relation and is located in the low-energy tail of the particle-hole continuum. Furthermore, the dynamic screening potential shows deep attractive regions as a function of the distance at higher frequencies which might lead to long-lived pair formation depending on the behaviour of the pair susceptibility. We also find that Friedel oscillations are present in these systems but are highly suppressed.

Pairing susceptibility of iron-based superconductors within a two-layer Hubbard model

By using the determinant quantum Monte Carlo method, we studied the dominant pairing susceptibility of iron-based superconductors within an extended Hubbard model, which describes the underlying electronic structure of both iron pnictides and iron chalcogenides. The extended Hubbard model is constructed by two iron layers, each of which forms two sublattices on a square structure. Although the coupling between the two layers has different effects on the behavior of pairings in iron pnictides and iron chalcogenides, our non-biased numerical simulations reveal that the pairing with Sxy symmetry dominates over the studied parameter for both materials.

Variational cluster approach to superconductivity and magnetism in the Kondo lattice model

We investigate in detail antiferromagnetic (AF) and superconducting (SC) phases as well as their coexistence in the two-dimensional Kondo lattice model on a square lattice, which is a paradigmatic model for heavy fermion materials. The results presented are mainly obtained using the variational cluster approximation (VCA) and are complemented by analytical findings for the equations of motion of pairing susceptibilities. A particularly interesting aspect is the possibility to have s-wave SC near half filling as reported by Bodensiek \textit{et al.} [Phys. Rev. Lett. \textbf{110}, 146406 (2013)]. When doping the system, we identify three regions which correspond to an AF metallic phase with small Fermi surface at weak coupling, an AF metal with a different Fermi surface topology at intermediate coupling, and a paramagnetic metal with a large Fermi surface at strong coupling. The transition between these two AF phases is found to be discontinuous at lower fillings, but turns to a continuous one when approaching half-filling. In the quest for s-wave superconductivity, only solutions are found which possess mean-field character. No true superconducting solutions caused by correlation effects are found in the s-wave channel. In contrast, we clearly identify robust d-wave pairing away from half-filling. However, we show that only by treating antiferromagnetism and superconductivity on equal footing artificial superconducting solutions at half-filling can be avoided. Our VCA findings support scenarios previously identified by variational Monte Carlo approaches and are a starting point for future investigations with VCA and further approaches such as cluster-embedding methods.

The ab initio study of unconventional superconductivity in CeCoIn5 and FeSe

The electronic structure and shape of the Fermi surface are known to be of fundamental importance for the superconducting instability in real materials. We demonstrate that such an instability may be explored by static Cooper pair susceptibility renormalized by pairing interaction and present an efficient method of its evaluation using Wannier orbitals derived from ab initio calculation. As an example, this approach is used to search for an unconventional superconducting phase of the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) type in a heavy-fermion compound CeCoIn5 and an iron-based superconductor FeSe. The results suggest that the FFLO superconducting phase occurs at finite magnetic field in both materials.

Superconductivity mediated by quantum critical antiferromagnetic fluctuations: The rise and fall of hot spots

In several unconventional superconductors, the highest superconducting transition temperature $T_{c}$ is found in a region of the phase diagram where the antiferromagnetic transition temperature extrapolates to zero, signaling a putative quantum critical point. The elucidation of the interplay between these two phenomena - high-$T_{c}$ superconductivity and magnetic quantum criticality - remains an important piece of the complex puzzle of unconventional superconductivity. In this paper, we combine sign-problem-free Quantum Monte Carlo simulations and field-theoretical analytical calculations to unveil the microscopic mechanism responsible for the superconducting instability of a general low-energy model, called spin-fermion model. In this approach, low-energy electronic states interact with each other via the exchange of quantum critical magnetic fluctuations. We find that even in the regime of moderately strong interactions, both the superconducting transition temperature and the pairing susceptibility are governed not by the properties of the entire Fermi surface, but instead by the properties of small portions of the Fermi surface called hot spots. Moreover, $T_{c}$ increases with increasing interaction strength, until it starts to saturate at the crossover from hot-spots dominated to Fermi-surface dominated pairing. Our work provides not only invaluable insights into the system parameters that most strongly affect $T_{c}$, but also important benchmarks to assess the origin of superconductivity in both microscopic models and actual materials.

Competing pairing channels in the doped honeycomb lattice Hubbard model

Proposals for superconductivity emerging from correlated electrons in the doped Hubbard model on the honeycomb lattice range from chiral $d+id$ singlet to $p+ip$ triplet pairing, depending on the considered range of doping and interaction strength, as well as the approach used to analyze the pairing instabilities. Here, we consider these scenarios using large-scale dynamic cluster approximation (DCA) calculations to examine the evolution in the leading pairing symmetry from weak to intermediate coupling strength. These calculations focus on doping levels around the van Hove singularity (VHS) and are performed using DCA simulations with an interaction-expansion continuous-time quantum Monte Carlo cluster solver. We calculated explicitly the temperature dependence of different uniform superconducting pairing susceptibilities and found a consistent picture emerging upon gradually increasing the cluster size: while at weak coupling the $d+id$ singlet pairing dominates close to the VHS filling, an enhanced tendency towards $p$-wave triplet pairing upon further increasing the interaction strength is observed. The relevance of these systematic results for existing proposals and ongoing pursuits of odd-parity topological superconductivity are also discussed.

Pairing correlations near a Kondo-destruction quantum critical point

Motivated by the unconventional superconductivity observed in heavy-fermion metals, we investigate pairing susceptibilities near a continuous quantum phase transition of the Kondo-destruction type. We solve two-impurity Bose-Fermi Anderson models with Ising and Heisenberg forms of the interimpurity exchange interaction using continuous-time quantum Monte-Carlo and numerical renormalization-group methods. Each model exhibits a Kondo-destruction quantum critical point separating Kondo-screened and local-moment phases. For antiferromagnetic interimpurity exchange interactions, singlet pairing is found to be enhanced in the vicinity of the transition. Implications of this result for heavy-fermion superconductivity are discussed.

Odd-frequency pairing and Ising spin susceptibility in time-reversal-invariant superfluids and superconductors

We here examine the relation between odd-frequency spin-triple even-parity (OTE) Cooper pairs and anomalous surface magnetic response in time-reversal invariant (TRI) spin-triplet superfluids and superconductors. The spin susceptibility generally consists of two contributions: Even-frequency odd-parity pair amplitudes and odd-frequency even-parity pair amplitudes. The OTE pair amplitudes are absent in the bulk region, but ubiquitously exist in the surface and interface region as Andreev bound states. We here clarify that additional discrete symmetries, originating from the internal symmetry and point group symmetry, impose strong constraint on the OTE pair amplitudes emergent in the surface of TRI superfluids and superconductors. As a result of the symmetry constraint, the magnetic response of the OTE pairs yields Ising-like anisotropy. For the topological phase of the $^3$He-B in a restricted geometry, the coupling of the OTE pair amplitudes to an applied field is prohibited by an additional discrete symmetry. Once the discrete symmetry is broken, however, the OTE pairs start to couple to the applied field, which anomalously enhances surface spin susceptibility. Furthermore, we extend this theory to TRI superconductors, where the corresponding discrete symmetry is the mirror reflection symmetry.

Observing the inhomogeneity in the holographic models of superconductors

We study the gravity duals of striped holographic superconductors in the AdS black hole and AdS soliton backgrounds. We show the dependences of the condensation and the critical temperature/critical chemical potential on the inhomogeneity in these two different spacetimes. By exploring the dynamics of the normal phase through the scalar field perturbation, we argue that the pair susceptibility and the conductivity can be possible phenomenological indications to disclose the property of inhomogeneity.

Superconductivity of composite particles in a two-channel Kondo lattice.

Emergence of odd-frequency s-wave superconductivity is demonstrated in the two-channel Kondo lattice by means of the dynamical mean-field theory combined with the continuous-time quantum MonteCarlo method. Around half filling of the conduction bands, divergence of an odd-frequency pairing susceptibility is found, which signals instability toward the superconductivity. The corresponding order parameter is equivalent to a staggered composite-pair amplitude with even frequencies, which involves both localized spins and conduction electrons. A model wave function is constructed for the composite order with the use of symmetry operations such as charge conjugation and channel rotations. Given a certain asymmetry of the conduction bands, another s-wave superconductivity is found that has a uniform order parameter. The Kondo effect in the presence of two channels is essential for both types of unconventional superconductivity.

Influence of s± symmetry on unconventional superconductivity in pnictides above the Pauli limit – two-band model study

The theoretical analysis of the Cooper pair susceptibility shows the two-band Fe-based superconductors (FeSC) to support the existence of the phase with nonzero Cooper pair momentum (called the Fulde--Ferrel--Larkin--Ovchinnikov phase or shortly FFLO), regardless of the order parameter symmetry. Moreover this phase for the FeSC model with $s_{\pm}$ symmetry is the ground state of the system near the Pauli limit. This article discusses the phase diagram $h-T$ for FeSC in the two-band model and its physical consequences. We compare the results for the superconducting order parameter with s-wave and $s_{\pm}$-wave symmetry -- in first case the FFLO phase can occur in both bands, while in second case only in one band. We analyze the resulting order parameter in real space -- showing that the FeSC with $s_{\pm}$-wave symmetry in the Pauli limit have typical properties of one-band systems, such as oscillations of the order parameter in real space with constant amplitude, whereas with s-wave symmetry the oscillations have an amplitude modulation. Discussing the free energy in the superconducting state we show that in absence of orbital effects, the phase transition from the BCS to the FFLO state is always first order, whereas from the FFLO phase to normal state is second order.

Evolution of the superconductivity dome in the two-dimensional Hubbard model

In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the $d$-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature ${T}_{c}$, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the $d$-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.

Numerical Assessment of Radiated Susceptibility of Twisted-Wire Pairs With Random Nonuniform Twisting

The influence of twist-pitch nonuniformity on the radiated susceptibility of a twisted-wire pair (TWP) running above ground and illuminated by a plane-wave field is investigated via a transmission-line model. The geometrical representation of the TWP involves a bifilar helix with a variable pitch function, which readily allows for the description of any unwanted/unknown deviations from the ideal twist geometry. Basic examples are used to show that such a deformation plays a fundamental role in the pickup of differential-mode (DM) noise, thus substantiating the need for the proposed prediction model. By describing the pitch random variations as a spatial stochastic process, a realistic representation of the twisting features is obtained (both for TWPs with imperfect uniform twists and for TWPs intentionally manufactured with random nonuniform twists). Repeated runs are used to characterize the sensitivity of the induced DM noise to random nonuniform twisting. Especially, it is shown that the twist-pitch nonuniformity, sometimes intentionally used to reduce pair-to-pair crosstalk in TWP cables, may however lead to a severe decrease of the TWP immunity to external fields.

Exchange and spin states in quantum dots under strong spatial correlations. Computer simulation by the Feynman path integral method

The fundamental laws in the behavior of electrons in model quantum dots that are caused by exchange and strong Coulomb correlations are studied. The ab initio path integral method is used to numerically simulate systems of two, three, four, and six interacting identical electrons confined in a three-dimensional spherical potential well with a parabolic confining potential against the background of thermal fluctuations. The temperature dependences of spin and collective spin magnetic susceptibility are calculated for model quantum dots of various spatial sizes. A basically exact procedure is proposed for taking into account the permutation symmetry and the spin state of electrons, which makes it possible to perform numerical calculations using modern computer facilities. The conditions of applicability of a virial energy estimator and its optimum form in exchange systems are determined. A correlation estimator of kinetic energy, which is an alternative to a basic estimator, is suggested. A fundamental relation between the kinetic energy of a quantum particle and the character of its virtual diffusion in imaginary time is demonstrated. The process of natural “pairing” of electron spins during the compression of a quantum dot and cooling of a system is numerically reproduced in terms of path integrals. The temperature dependences of the spin magnetic susceptibility of electron pairs with a characteristic maximum caused by spin pairing are obtained.

Numerical study of magnetic and pairing correlation in a bilayer triangular lattice

By using the determinant Quantum Monte Carlo method, the magnetic and pairing correlation of the NaxCoO2⋅yH2O system are studied within the Hubbard model on a bilayer triangular lattice. The temperature dependence of the spin correlation function and pairing susceptibility with several kinds of symmetries at different electron fillings and interlayer coupling terms are investigated. It is found that the system shows antiferromagnetic correlation around half filling, and the fn-wave pairing correlation dominates over other kinds of pairing symmetry in the low doping region. As the electron filling decreases from half filling, both ferromagnetic correlation and the f-wave pairing susceptibility are enhanced and tend to dominate. It is also shown that both the magnetic susceptibility and pairing susceptibility decrease as the interlayer coupling increases.

BCS-BEC crossover at finite temperature in spin-orbit-coupled Fermi gases

By adopting a $T$-matrix-based method within the ${G}_{0}G$ approximation for the pair susceptibility, we study the effects of the pairing fluctuation on the three-dimensional spin-orbit-coupled Fermi gases at finite temperature. The critical temperatures of the superfluid to normal phase transition are determined for three different types of spin-orbit coupling (SOC): (1) the extreme oblate (EO) or Rashba SOC, (2) the extreme prolate or equal Rashba-Dresselhaus SOC, and (3) the spherical (S) SOC. For EO- and S-type SOC, the SOC dependence of the critical temperature signals a crossover from BCS to BEC state; at strong SOC limit, the critical temperature recovers those of ideal BEC of rashbons. The pairing fluctuation induces a pseudogap in the fermionic excitation spectrum in both superfluid and normal phases. We find that, for EO- and S-type SOC, even at weak coupling, sufficiently strong SOC can induce sizable pseudogap. Our research suggests that the spin-orbit-coupled Fermi gases may open new means to the study of the pseudogap formation in fermionic systems.

Competing electronic orders on kagome lattices at van Hove filling

The electronic orders in Hubbard models on a Kagome lattice at van Hove filling are of intense current interest and debate. We study this issue using the singular-mode functional renormalization group theory. We discover a rich variety of electronic instabilities under short range interactions. With increasing on-site repulsion $U$, the system develops successively ferromagnetism, intra unit-cell antiferromagnetism, and charge bond order. With nearest-neighbor Coulomb interaction $V$ alone (U=0), the system develops intra-unit-cell charge density wave order for small $V$, s-wave superconductivity for moderate $V$, and the charge density wave order appears again for even larger $V$. With both $U$ and $V$, we also find spin bond order and chiral $d_{x^2 - y^2} + i d_{xy}$ superconductivity in some particular regimes of the phase diagram. We find that the s-wave superconductivity is a result of charge density wave fluctuations and the squared logarithmic divergence in the pairing susceptibility. On the other hand, the d-wave superconductivity follows from bond order fluctuations that avoid the matrix element effect. The phase diagram is vastly different from that in honeycomb lattices because of the geometrical frustration in the Kagome lattice.

Quantum Monte Carlo Study of a Dominants-Wave Pairing Symmetry in Iron-Based Superconductors

We perform a systematic quantum Monte Carlo study of the pairing correlation in the S(4) symmetric microscopic model for iron-based superconductors. It is found that the pairing with an extensive s-wave symmetry robustly dominates over other pairings at low temperature in a reasonable parameter region regardless of the change of Fermi surface topologies. The pairing susceptibility, the effective pairing interaction, and the (π, 0) antiferromagnetic correlation strongly increase as the on-site Coulomb interaction increases, indicating the importance of the effect of electron-electron correlation. Our nonbiased numerical results provide a unified understanding of the superconducting mechanism in iron pnictides and iron chalcogenides and demonstrate that the superconductivity is driven by strong electron-electron correlation effects.

The Fulde–Ferrell–Larkin–Ovchinnikov State in Pnictides

Fe-based superconductors (FeSC) exhibit all the properties of systems that allow the formation of a superconducting phase with oscillating order parameter, called the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase. By the analysis of the Cooper pair susceptibility in two-band FeSC, such systems are shown to support the existence of a FFLO phase, regardless of the exhibited order parameter symmetry. We also show the state with nonzero Cooper pair momentum, in superconducting FeSC with ∼cos(k x )⋅cos(k y ) symmetry, to be the ground state of the system in a certain parameter range.

Lifshitz transition in the two-dimensional Hubbard model

Using large-scale dynamical cluster quantum Monte Carlo simulations, we study the Lifshitz transition of the two dimensional Hubbard model with next-nearest-neighbor hopping ($t'$), chemical potential and temperature as control parameters. At $t'\le0$, we identify a line of Lifshitz transition points associated with a change of the Fermi surface topology at zero temperature. In the overdoped region, the Fermi surface is complete and electron-like; across the Lifshitz transition, the Fermi surface becomes hole-like and develops a pseudogap. At (or very close to) the Lifshitz transition points, a van Hove singularity in the density of states crosses the Fermi level. The van Hove singularity occurs at finite doping due to correlation effects, and becomes more singular when $t'$ becomes more negative. The resulting temperature dependence on the bare d-wave pairing susceptibility close to the Lifshitz points is significantly different from that found in the traditional van Hove scenarios. Such unambiguous numerical observation of the Lifshitz transition at $t'\le0$ extends our understanding of the quantum critical region in the phase diagram, and shines lights on future investigations of the nature of the quantum critical point in the two dimensional Hubbard model.