Unconventional Discontinuous Transitions in Isospin Systems
Not Available
Raines, Zachary M. (ORCID:0000000307309089)↗
Engineering topics
Chubukov, Andrey V.
Publications and source records attributed to Chubukov, Andrey V..
Location and thermal evolution of the pseudogap due to spin fluctuations
We study pseudogap behavior in a metal near a spin density wave (SDW) instability due to thermal magnetic fluctuations. We consider the t–t' Hubbard model on a square lattice at a finite doping, at intermediate coupling strength, and analyze the thermal evolution of the electron spectral function between a SDW ordered state at low temperatures and a normal Fermi liquid at high temperatures. We argue that for proper description of the pseudogap one needs to sum up infinite series of diagrams for both the fermionic self-energy and the SDW order parameter in the SDW state or the magnetic correlation length in the paramagnetic state. We use the eikonal approach to sum up an infinite series of diagrammatic contributions from thermal fluctuations. Earlier studies found that in the SDW state the spectral function A k (ω) of a hot fermion at a finite T is exponentially small below the energy scale Δ(T), which scales with SDW order and vanishes at the ordering temperature T N , and has a hump at a larger frequency Δ PG , comparable to the zero-temperature SDW gap Δ(T = 0). We argue that the hump, which we associate with the pseudogap, survives in some T range above T N . We show that this range is split by regions of strong and weak pseudogap behavior. In the first region, Δ PG is weakly temperature dependent, despite the fact that it comes from thermal fluctuations. Such a behavior has been seen in numerical studies of the Hubbard model. We show that to obtain it one needs to go beyond the one-loop approximation and sum up the infinite series of diagrams. In the second regime, Δ PG decreases with increasing T and eventually vanishes. Here, we further argue that a magnetic pseudogap at a finite T emerges only if the ground state is magnetically ordered. We present the phase diagram and apply the results to high-T c cuprates.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Crucial role of thermal fluctuations and vertex corrections for the magnetic pseudogap
It is generally believed that in a two-dimensional metal, whose ground state is antiferromagnetically ordered with Q = (π,π), thermal (static) magnetic fluctuations give rise to precursor behavior above TN in which the spectral function of a hot fermion (the one for which k and k + Q are Fermi surface points) contains two peaks, separated by roughly the same energy as in the antiferromagnetically ordered state. The two peaks persist in some range of T > T N and eventually merge into a single peak at zero frequency. This behavior is obtained theoretically by departing from free fermions in a paramagnet and evaluating the dressed fermionic Green's function by summing up infinite series of diagrams with contributions from thermal magnetic fluctuations. We show, following that keeping vertex renormalization diagrams in these series is crucial as other terms only broaden the spectral function of a hot fermion but do not shift its maximum away from zero frequency. As the consequence, the magnetic pseudogap should be treated as an input for theories that neglect vertex corrections, such as, e.g., Eliashberg theory for magnetically mediated superconductivity. We also analyze the potential pseudogap behavior at T = 0. Here, we argue that it may exist, but only at a finite correlation length, and not as a precursor to antiferromagnetism.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Density of States and Spectral Function of a Superconductor out of a Quantum-Critical Metal
We analyze the validity of a quasiparticle description of a superconducting state above a metallic quantum-critical point (QCP). A normal state at a QCP is a non-Fermi liquid with no coherent quasiparticles. A superconducting order gaps out low-energy excitations, except for a sliver of states for non-s-wave gap symmetry, and at a first glance, restores coherent quasiparticle behavior. We argue that this does not necessarily hold as the fermionic self-energy may remain singular above the gap edge. This singularity gives rise to markedly non-BCS behavior of the density of states and to the appearance of a nondispersing mode at the gap edge in the spectral function. We analyze the set of quantum-critical models with an effective dynamical four-fermion interaction V(Ω) ∝ 1/Ω γ , where Ω is a frequency of a boson, which mediates the interaction. We show that coherent quasiparticle behavior in a superconducting state holds for γ < 1/2, but breaks down for larger γ. Here, we discuss signatures of quasiparticle breakdown and compare our results with the photoemission data for Bi2201 and Bi2212.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Transformer spin-triplet superconductivity at the onset of isospin order in bilayer graphene
We consider the origin of superconductivity found recently in Bernal bilayer graphene at the onset of isospin-polarized order, trying to infer the pairing mechanism and superconducting order from the measurements available to date. The superconductivity is induced by a parallel magnetic field and persists well above the Pauli limit, indicating an unconventional scenario of quantum-critical pairing, where soft fluctuations of isospin give rise to spin-triplet superconductivity. We consider the scenario in which the pairing interaction is entirely repulsive, which stands in contrast to the typical quantum-critical pairing mechanisms. Superconductivity emerges through a “transformer” mechanism where, in the presence of an in-plane magnetic field, the incipient valley polarization converts a frequency-independent repulsion into one with a strong nonmonotonic frequency dependence. Such an interaction enables a nonzero solution for the pairing gap function that changes sign as a function of frequency. Furthermore, the same mechanism holds at zero field in the presence of spin-orbit coupling, providing a likely explanation for the recently observed superconductivity in bilayer graphene on the WSe 2 monolayer.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Following the Higgs mode across the BCS-BEC crossover in two dimensions
Although substantial effort has been dedicated to analyzing the Higgs (amplitude) mode in superconducting systems, there are relatively few studies of the Higgs peak's dispersion and width, quantities which are observable in spectroscopic measurements. These properties can be obtained from the location of the pole of the amplitude (Higgs) susceptibility in the lower half plane of complex frequency. We analyze the behavior of the Higgs mode in a two-dimensional neutral fermionic superfluid at T = 0 throughout the crossover from Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensation (BEC) behavior. Here, this occurs when the dressed chemical potential μ changes sign from positive to negative. For μ > 0, we find a pole in the Higgs susceptibility in the lower half plane of frequency and demonstrate that it leads to a well-defined peak in spectroscopic probes. For μ < 0, the pole still exists but is “hidden,” not giving rise to a peak in spectroscopic probes. Extending this analysis to a charged superconductor, we find that the Higgs mode is unaffected by the long-range Coulomb interaction.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Free energy and specific heat near a quantum critical point of a metal
We analyze the free energy and the specific heat for fermions interacting with a gapless boson at a quantum-critical point (QCP) in a metal. We use the Luttinger-Ward-Eliashberg formula for the free energy in the normal state, which includes contributions from bosons, fermions, and their interaction, all expressed via fully dressed fermionic and bosonic propagators. The sum of the last two contributions is the free energy F γ of an effective low-energy model of fermions with boson-mediated dynamical 4-fermion interaction V(Ω m ) ∝ 1/|Ω m | γ (the γ model). This purely electronic model has been used to analyze the interplay between non-Fermi liquid behavior and pairing near a QCP, which are both independent of the upper energy cutoff Λ. However, the specific heat C γ (T), obtained from F γ , does depend on Λ. We argue that this dependence is spurious and cancels out, once we include the contribution from bosons. We further argue that the full C(T) is the sum of the contribution from free fermions and the one from a critical boson, with the fully dressed propagator, other terms cancel out. We compare the full C(T) with the C γ (T), obtained using recently proposed regularization of F γ . We argue that for γ < 1, the full C(T) and the regularized C γ (T) differ by a γ-dependent prefactor, while for γ > 1, the full C(T) and C γ (T) differ by the positive contribution from free massless fermions (a positive constant for the electron-phonon case γ = 2). For these γ, C γ (T) is negative, but the full C(T) is positive. We argue that only the full C(T) matters as the positive and the negative contributions originate from the term in C(T) which contains the fully dressed bosonic propagator. Here, we then argue that the normal state remains stable until the pairing instability develops.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Superconductivity out of a non-Fermi liquid: Free energy analysis
In this paper, we present an in-depth analysis of the condensation energy E c for a superconductor in a situation when superconductivity emerges out of a non-Fermi liquid due to pairing mediated by a massless boson. This is the case for electronic-mediated pairing near a quantum-critical point in metal, for pairing in SYK-type models, and for phonon-mediated pairing in the properly defined limit, when the dressed Debye frequency vanishes. We consider a subset of these quantum-critical models, in which the pairing in a channel with a proper spatial symmetry is described by an effective 0 + 1 dimensional model with the effective dynamical interaction V(Ω m ) = g¯ γ /|Ω m | γ , where γ is model-specific (the γ model). In previous papers, we argued that the pairing in the γ model is qualitatively different from that in a Fermi liquid, and the gap equation at T = 0 has an infinite number of topologically distinct solutions, Δ n (ω m ), where an integer n, running between 0 and infinity, is the number of zeros of Δ n (ω m ) on the positive Matsubara axis. This gives rise to the set of extrema of E c at E c,n , of which E c,0 is the global minimum. Here, the spectrum E c,n is discrete for a generic γ < 2 but becomes continuous at γ = 2–0. Here, we discuss in more detail the profile of the condensation energy near each E c,n and the transformation from a discrete to a continuous spectrum at γ → 2. We also discuss the free energy and the specific heat of the γ model in the normal state.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
SU(4) Symmetry in Twisted Bilayer Graphene: An Itinerant Perspective
We study symmetry-broken phases in twisted bilayer graphene at small filling above charge neutrality and at van Hove filling. We argue that the Landau functionals for the particle-hole order parameters at these fillings both have an approximate SU(4) symmetry, but differ in the sign of quartic terms. We determine the order parameter manifold of the ground state and analyze its excitations. For small fillings, we find a strong first-order transition to an SU(3) Ⓧ U(1) manifold of orders that break spin-valley symmetry and induce a 3-1 splitting of fermionic excitations. For van Hove filling, we find a weak first-order transition to an SU(2) Ⓧ SU(2) Ⓧ U(1) manifold of orders that preserves the twofold band degeneracy. Here, we discuss the effect of particle-hole orders on superconductivity and compare with strong-coupling approaches.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Effect of repulsion on superconductivity at low density
Here we examine the effect of repulsion on superconductivity in a three-dimensional system with a Bardeen-Pines-like interaction in the low-density limit, where the chemical potential μ is much smaller than the phonon frequency ω L . We parametrize the strength of the repulsion by a dimensionless parameter f, and find that the superconducting transition temperature T c approaches a nonzero value in the μ = 0 limit as long as f is below a certain threshold f * . In this limit, we find that T c goes to zero as a power of f * - f, in contrast to the high density limit, where T c goes to zero exponentially quickly as f approaches f * . For all nonzero f, the gap function Δ(ω m ) changes sign along the Matsubara axis, which allows the system to partially overcome the repulsion at high frequencies. We trace the position of the gap node with f and show that it approaches zero frequency as f approaches f * . To investigate the robustness of our conclusions, we then go beyond the Bardeen-Pines model and include full dynamical screening of the interaction, finding that T c still saturates to a nonzero value at μ = 0 when f < f * .
36 MATERIALS SCIENCE↗
Dynamical exponent of a quantum critical itinerant ferromagnet: A Monte Carlo study
Here we consider the effect of the coupling between two-dimensional (2D) quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions. A common belief is that near a q = 0 ferromagnetic transition, fermions induce an Ω/q Landau damping of rotors (i.e., the dynamical critical exponent is z = 3) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy Σ∝ω 2/3 . This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies. Here we show that for the XY case the behavior is different. We report the results of large-scale quantum Monte Carlo simulations, which show that at small frequencies z = 2 and Σ∝ω 1/2 . We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces 1/q in the Landau damping by a constant. We discuss the implication of these results to experiments.
36 MATERIALS SCIENCE↗
Quasiparticle Scattering in a Superconductor near a Nematic Critical Point: Resonance Mode and Multiple Attractive Channels
We analyze the scattering rate for 2D fermions interacting via soft nematic fluctuations. The ground state is an s-wave superconductor, but other pairing channels are almost equally attractive. This strongly alters the scattering rate: At energies beyond the pairing gap Δ, it is renormalized by contributions from all pairing channels. At energies of order Δ, it is determined by the competition between scattering into a gapped continuum and dispersing nematic resonance. The outcome is a “peak-peak-dip-hump” spectrum, similar, but not identical, to the “peak-dip-hump” structure in the cuprates.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Specific heat and gap structure of a nematic superconductor: Application to FeSe
We report the results of our in-depth analysis of spectroscopic and thermodynamic properties of a multiorbital metal, like FeSe, which first develops a nematic order and then undergoes a transition into a superconducting state, which coexists with nematicity. We analyze the angular dependence of the gap function and specific heat C v (T) of such a nematic superconductor. We specifically address three issues: (i) the angular dependence of the gap in light of the competition between the nematicity-induced s - d mixture and the orbital transmutation of low-energy excitations in the nematic state, (ii) the effect of nematicity on the magnitude of the jump of the specific heat C v (T) at T c and the temperature dependence of C v (T) below T c , and (iii) a potential transition at T c1 < T c from an s + d state to an s + e iη d state that breaks time-reversal symmetry. We consider two scenarios for a nematic order: scenario A, in which this order develops between d xz and d yz orbitals on hole and electron pockets, and scenario B, in which there is an additional component of the nematic order for d xy fermions on the two electron pockets.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
One-dimensional scattering of two-dimensional fermions near quantum criticality
Forward scattering and backscattering play an exceptional role in the physics of two-dimensional interacting fermions. In a Fermi liquid, both give rise to a nonanalytic ω 2 ln (ω) form of the fermionic scattering rate at second order in the interaction. Here we argue that higher powers of ln (ω) show up in the backscattering contribution at higher orders. We show that these terms come from “planar” processes, which are effectively one-dimensional. This is explicitly demonstrated by extending a Fermi liquid to the limit of N >> 1 fermionic flavors, when only planar processes survive. We sum the leading logarithms for the case of a two-dimensional Fermi liquid near a nematic transition, and we obtain an expression for the scattering rate at T = 0 to all orders in the interaction. For a repulsive interaction, the resulting rate is logarithmically suppressed, and the result is valid down to ω = 0. For an attractive interaction, the ground state is an s-wave superconductor with a gap Δ 0 . We show that in this case the scattering rate increases as ω is reduced toward Δ 0 . At ω ≥ Δ 0 , the behavior of the scattering rate is rather unconventional as many pairing channels compete near a nematic critical point, and the s-wave wins only by a narrow margin. We take superconductivity into consideration and obtain the scattering rate also at smaller ω ≃ Δ 0 .
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Competing orders at higher-order Van Hove points
Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. In this study we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization-group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive interaction, we find an unconventional spin Pomeranchuk order, in wich the spin oder parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case.
36 MATERIALS SCIENCE↗
Valley magnetism, nematicity, and density wave orders in twisted bilayer graphene
We analyze density wave and Pomeranchuk orders in twisted bilayer graphene. This complements our earlier analysis of the pairing instabilities. We assume that near half filling of either conduction or valence band, the Fermi level is close to Van Hove points, where the density of states diverges, and study potential instabilities in the particle-hole channel within a patch model with two valley degrees of freedom. The hexagonal symmetry of twisted bilayer graphene allows for either six or twelve Van Hove points. We consider both cases and find the same two leading candidates for particle-hole order. One is an SU(2)-breaking spin state with ferromagnetism within a valley. A subleading intervalley hopping induces antiferromagnetism between the valleys. The same state has also been obtained in strong-coupling approaches, indicating that this order is robust. Here, the other is a mixed state with 120° complex spin order and orthogonal complex charge order. In addition, we find a weaker but still attractive interaction in nematic channels, and discuss the type of a nematic order.
36 MATERIALS SCIENCE↗
Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data
Quantum Monte Carlo (QMC) simulations of correlated electron systems provide unbiased information about system behavior at a quantum critical point (QCP) and can verify or disprove the existing theories of non-Fermi liquid (NFL) behavior at a QCP. However, simulations are carried out at a finite temperature, where quantum critical features are masked by finite-temperature effects. Here, we present a theoretical framework within which it is possible to separate thermal and quantum effects and extract the information about NFL physics at T = 0. We demonstrate our method for a specific example of 2D fermions near an Ising ferromagnetic QCP. We show that one can extract from QMC data the zero-temperature form of fermionic self-energy Σ(ω) even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of Σ(ω) agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys ω 2/3 scaling at low frequencies. Our results open up an avenue for QMC studies of quantum critical metals.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome
Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the T dependence of the specific heat C(T) deep inside a fully gapped s +– superconducting state as the AFM-QCP is approached. We find that at the AFM-QCP, the specific heat C(T) vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully gapped BCS superconductors. This robust result is due to a nonanalytic contribution to the free energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, C(T) shows a crossover from a T2 behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. As a result, we argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.
36 MATERIALS SCIENCE↗