Source author record

Andrey V. Chubukov

Andrey V. Chubukov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2026arXiv

Spatially resolved collective modes in d-wave superconductors

We analyze the dispersion of collective modes in a superconductor with $d-$wave symmetry of the order parameter in the presence of long-range Coulomb interaction. We use diagrammatic technique and quasiclassical theory in Keldysh-Nambu formalism to compute longitudinal and transverse pair susceptibilities and extract from them the dispersion of the longitudinal and transverse collective mode. We show that at T=0, the dispersion of the transverse (plasma) mode is the same as in an s-wave superconductor, but at a finite temperature it is softer and has a much larger decay rate due to the partial screening of the Coulomb potential by nodal quasiparticles. We show that the dispersion of the longitudinal mode depends on the direction of momentum with respect to the positions of the nodes of the d-wave gap, while the decay rate of this mode does not depend on momentum. We discuss experimental implications of our results.

preprint2022arXiv

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, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations.

preprint2022arXiv

Quantum phase transition in a clean superconductor with repulsive dynamical interaction

We consider a model of electrons at zero temperature, with a repulsive interaction which is a function of the energy transfer. Such an interaction can arise from the combination of electron-electron repulsion at high energies and the weaker electron-phonon attraction at low energies. As shown in previous works, superconductivity can develop despite the overall repulsion due to the energy dependence of the interaction, but the gap $Δ(ω)$ must change sign at some (imaginary) frequency $ω_0$ to counteract the repulsion. However, when the constant repulsive part of the interaction is increased, a quantum phase transition towards the normal state occurs. We show that, as the phase transition is approached, $Δ$ and $ω_0$ must vanish in a correlated way such that $1/|\log[Δ(0)]| \sim ω_0^2$. We discuss the behavior of phase fluctuations near this transition and show that the correlation between $Δ(0)$ and $ω_0$ locks the phase stiffness to a non-zero value.

preprint2022arXiv

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.

preprint2022arXiv

The dynamical exponent of a quantum critical itinerant ferromagnet: a Monte Carlo study

We consider the effect of the coupling between 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 $Σ\propto ω^{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 $Σ\propto ω^{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.

preprint2022arXiv

Twists and turns of superconductivity from a repulsive dynamical interaction

We review recent theoretical progress in understanding spatially uniform $s$-wave superconductivity which arises from a fermion-fermion interaction, which is repulsive on the Matsubara axis, where it is real, but does depend on the transferred frequency. Such a situation holds, e.g., for systems with a screened Coulomb and retarded electron-phonon interaction. We show that despite repulsion, superconductivity is possible in a certain range of system parameters. However, at $T=0$ the gap function on the Matsubara axis, $Δ(ω_m)$, must pass through zero and change sign at least once. These zeros of $Δ(ω_m)$ have a topological interpretation in terms of dynamical vortices, and their presence imposes a constraint on the variation of the phase of the gap function along the real frequency axis, which can potentially be extracted from ARPES and other measurements. We discuss how superconductivity vanishes when the repulsion becomes too strong, and obtain a critical line which terminates at $T=0$ at a quantum-critical point for superconductivity. We show that the behavior of the gap function near this point is highly non-trivial. In particular, an infinitesimally small $Δ(ω_m)$ contains a singular $δ-$function piece $ω_m δ(ω_m)$. We argue that near the critical point superconductivity may be a mixed state with even-frequency and odd-frequency gap components.

preprint2021arXiv

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 1st-order transition to an SU(3)$\otimes$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 1st-order transition to an SO(4)$\otimes$U(1) manifold of orders that preserves the two-fold band degeneracy. We discuss the effect of particle-hole orders on superconductivity and compare with strong-coupling approaches.

preprint2021arXiv

The Effect of Repulsion on Superconductivity at Low Density

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 parameterize 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 non-zero value at $μ= 0$ when $f < f^*$.

preprint2020arXiv

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 a 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.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum critical point in a metal. II. The $γ$-model at a finite $T$ for $0<γ<1$

In this paper we continue the analysis of the interplay between non-Fermi liquid and superconductivity for quantum-critical systems, the low-energy physics of which is described by an effective model with dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$ model). In paper I [A. Abanov and A. V. Chubukov, Phys Rev B. 102, 024524 (2020)] two of us analyzed the $γ$ model at $T=0$ for $0<γ<1$ and argued that there exist a discrete, infinite set of topologically distinct solutions for the superconducting gap, all with the same spatial symmetry. The gap function $Δ_n (ω_m)$ for the $n$th solution changes sign $n$ times as the function of Matsubara frequency. In this paper we analyze the linearized gap equation at a finite $T$. We show that there exist an infinite set of pairing instability temperatures, $T_{p,n}$, and the eigenfunction $Δ_n (ω_{m})$ changes sign $n$ times as a function of a Matsubara number $m$. We argue that $Δ_n (ω_{m})$ retains its functional form below $T_{p,n}$ and at $T=0$ coincides with the $n$th solution of the nonlinear gap equation. Like in paper I, we extend the model to the case when the interaction in the pairing channel has an additional factor $1/N$ compared to that in the particle-hole channel. We show that $T_{p,0}$ remains finite at large $N$ due to special properties of fermions with Matsubara frequencies $\pm πT$, but all other $T_{p,n}$ terminate at $N_{cr} = O(1)$. The gap function vanishes at $T \to 0$ for $N > N_{cr}$ and remains finite for $N < N_{cr}$. This is consistent with the $T =0$ analysis.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal. I: The $γ$-model and its phase diagram at $T=0$. The case $0 < γ<1$

We analyze a class of quantum-critical models, in which momentum integration and the selection of a particular pairing symmetry can be done explicitly, and the competition between non-Fermi liquid and pairing can be analyzed within an effective model with dynamical electron-electron interaction $V(Ω_m)\sim 1/|Ω_m|^γ$ (the $γ$-model). In this paper, the first in the series, we consider the case $T=0$ and $0<γ<1$. We argue that tendency to pairing is stronger, and the ground state is a superconductor. We argue, however, that superconducting state is highly non-trivial as there exists a discrete set of topologically distinct solutions for the pairing gap $Δ_n (ω_m)$ ($n = 0, 1, 2..., \infty$). All solutions have the same spatial pairing symmetry, but differ in the time domain: $Δ_n (ω_m)$ changes sign $n$ times as a function of Matsubara frequency $ω_m$. The $n =0$ solution $Δ_0 (ω_m)$ is sign-preserving and tends to a finite value at $ω_m =0$, like in BCS theory. The $n = \infty$ solution corresponds to an infinitesimally small $Δ(ω_m)$. As a proof, we obtain the exact solution of the linearized gap equation at $T=0$ on the entire frequency axis for all $0<γ<1$, and an approximate solution of the non-linear gap equation.We argue that the presence of an infinite set of solutions opens up a new channel of gap fluctuations. We extend the analysis to the case where the pairing component of the interaction has additional factor $1/N$ and show that there exists a critical $N_{cr} >1$, above which superconductivity disappears, and the ground state becomes a non-Fermi liquid.We show that all solutions develop simultaneously once $N$ gets smaller than $N_{cr}$.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal. IV: The $γ$ model and its phase diagram at $1<γ<2$

In this paper we continue our analysis of the interplay between the pairing and the non-Fermi liquid behavior in a metal for a set of quantum-critical (QC) systems with an effective dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$-model). In previous papers we studied the cases $0<γ<1$ and $γ\approx 1$. We argued that the pairing by a gapless boson is fundamentally different from BCS/Eliashberg pairing by a massive boson as for the former there exists an infinite number of topologically distinct solutions for the gap function $Δ_n (ω_m)$ at $T=0$ ($n=0,1,2...$), each with its own condensation energy $E_{c,n}$. Here we extend the analysis to larger $1< γ<2$. We argue that the discrete set of solutions survives, and the spectrum of $E_{c,n}$ gets progressively denser as $γ$ approaches $2$ and eventually becomes continuous at $γ\to 2$. This increases the strength of "longitudinal" gap fluctuations, which tend to reduce the actual superconducting $T_c$ and give rise to a pseudogap region of preformed pairs. We also detect two features on the real axis for $γ>1$ which become critical at $γ\to 2$. First, the density of states evolves towards a set of discrete $δ-$functions. Second, an array of dynamical vortices emerges in the upper frequency half-plane. These two features come about because on a real axis, the real part of the interaction, $V' (Ω) \propto \cos(πγ/2)/|Ω|^γ$, becomes repulsive for $γ>1$, and the imaginary $V^{''} (Ω) \propto \sin(πγ/2)/|Ω|^γ$, gets progressively smaller at $γ\to 2$. The features on the real axis are consistent with the development of a continuum spectrum of $E_{c,n}$ obtained using $Δ_n (ω_m)$ on the Matsubara axis. We consider the case $γ=2$ separately in the next paper.

preprint2020arXiv

Nematic superconductivity in twisted bilayer graphene

Twisted bilayer graphene displays insulating and superconducting phases caused by exceptional flattening of its lowest energy bands. Superconductivity with highest $T_c$ appears at hole and electron dopings, near half-filling for valence or conduction bands. In the hole-doped case, the data show that three-fold lattice rotation symmetry is broken in the superconducting phase, i.e., a superconductor is also a nematic. We present the mechanism for nematic superconductivity. We take as an input the fact that at relevant dopings the Fermi energy lies in the vicinity of twist-induced Van Hove singularities in the density of states and argue that the low-energy physics can be properly described by patch models with six Van Hove points for electron doping and twelve Van Hove points for hole doping. We obtain pairing interactions for the patch models in terms of parameters of the microscopic model for the flat bands, which contains both local and twist-induced non-local interactions and show that the latter gives rise to attraction in different superconducting channels. For hole-doping, we find two attractive channels, $g$ and $i$-waves, with almost equal coupling constants. We show that in the co-existence state, where both order parameters are non-zero, the three-fold lattice rotation symmetry is broken, i.e., a superconductor is also a nematic. We find two possible nematic states, one is time-reversal symmetric, the other additionally breaks time-reversal symmetry. Our scenario for nematic superconductivity is based on generic symmetry considerations, and we expect it to be applicable also to other systems with two (or more) attractive channels with similar couplings.

preprint2020arXiv

Normal state properties of quantum critical metals at finite temperature

We study the effects of finite temperature on normal state properties of a metal near a quantum critical point to an antiferromagnetic or Ising-nematic state. At $T = 0$ bosonic and fermionic self-energies are traditionally computed within Eliashberg theory and obey scaling relations with characteristic power-laws. Quantum Monte Carlo (QMC) simulations have shown strong systematic deviations from these predictions, casting doubt on the validity of the theoretical analysis. We extend Eliashberg theory to finite $T$ and argue that for the $T$ range accessible in the QMC simulations, the scaling forms for both fermionic and bosonic self energies are quite different from those at $T = 0$. We compare finite $T$ results with QMC data and find good agreement for both systems. This, we argue, resolves the key apparent contradiction between the theory and the QMC simulations.

preprint2020arXiv

Orbital transmutation and the electronic spectrum of FeSe in the nematic phase

We consider the electronic spectrum near $M=(π,π)$ in the nematic phase of FeSe ($T<T_{\rm nem}$) and make a detailed comparison with recent ARPES and STM experiments. Our main focus is the unexpected temperature dependence of the excitations at the $M$ point. These have been identified as having $xz$ and $yz$ orbital character well below $T_{\rm nem}$, but remain split at $T>T_{\rm nem}$, in apparent contradiction to the fact that in the tetragonal phase the $xz$ and $yz$ orbitals are degenerate. Here we present two scenarios which can describe the data. In both scenarios, hybridization terms present in the tetragonal phase leads to an orbital transmutation, a change in the dominant orbital character of some of the bands, between $T > T_{\rm nem}$ and $T \ll T_{\rm nem}$. The first scenario relies on the spin-orbit coupling at the $M$ point. We show that a finite spin-orbit coupling gives rise to orbital transmutation, in which one of the modes, identified as $xz$ ($yz)$ at $T \ll T_{\rm nem}$, becomes predominantly $xy$ at $T > T_{\rm nem}$ and hence does not merge with the predominantly $yz$ ($xz$) mode. The second scenario, complementary to the first, takes into consideration the fact that both ARPES and STM are surface probes. In the bulk, a direct hybridization between the $xz$ and $yz$ orbitals is not allowed at the $M$ point, however, it is permitted on the surface. In the presence of a direct $xz/yz$ hybridization, the orbital character of the $xz/yz$ modes changes from pure $xz$ and pure $yz$ at $T \ll T_{\rm nem}$ to $xz \pm yz$ at $T > T_{\rm nem}$, i.e., the two modes again have mono-orbital character at low $T$, but do not merge at $T_{\rm nem}$. We discuss how these scenarios can be distinguished in polarized ARPES experiments.

preprint2020arXiv

Quantum Phase Transition in the Yukawa-SYK Model

We study the quantum phase transition upon variation of the fermionic density $ν$ in a solvable model with random Yukawa interactions between $N$ bosons and $M$ fermions, dubbed the Yukawa-SYK model. We show that there are two distinct phases in the model: an incompressible state with gapped excitations and an exotic quantum-critical, non-Fermi liquid state with exponents varying with $ν$. We show analytically and numerically that the quantum phase transition between these two states is first-order, as for some range of $ν$ the NFL state has a negative compressibility. In the limit $N/M\to \infty$ the first-order transition gets weaker and asymptotically becomes second-order, with an exotic quantum-critical behavior. We show that fermions and bosons display highly unconventional spectral behavior in the transition region.

preprint2020arXiv

The pairing glue in cuprate superconductors from the self-energy revealed via machine learning

Recently, machine learning was applied to extract both the normal and the anomalous components of the self-energy from photoemission data at the antinodal points in Bi-based cuprate high-temperature superconductors [Y. Yamaji {\it et al.}, arXiv:1903.08060]. It was argued that both components do show prominent peaks near 50 meV, which hold information about the pairing glue, but the peaks are hidden in the actual data, which measure only the total self-energy. We analyze the self-energy within an effective fermion-boson theory. We show that soft thermal fluctuations give rise to peaks in both components of the self-energy at a frequency comparable to superconducting gap, while they cancel in the total self-energy; all irrespective of the nature of the pairing boson. However, in the quantum limit $T \to 0$ prominent peaks survive only for a very restricted subclass of pairing interactions. We argue that the way to potentially nail down the pairing boson is to determine the thermal evolution of the peaks.

preprint2020arXiv

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 non-analytic 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 $T^2$ behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. 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.

preprint2020arXiv

Valley magnetism, nematicity, and density wave orders in twisted bilayer graphene

We analyze density-wave and Pomeranchuk orders in twisted bilayer graphene. This compliments 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 inter-valley hopping induces antiferromagnetism between the valleys. The same state has also been obtained in strong coupling approaches, indicating that this order is robust. The other is a mixed state with $120^\circ$ 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.

preprint2019arXiv

Raman response in the nematic phase of FeSe

Raman experiments on bulk FeSe revealed that the low-frequency part of $B_{1g}$ Raman response $R_{B_{1g}}$, which probes nematic fluctuations, rapidly decreases below the nematic transition at $T_n \sim 85$K. Such behavior is usually associated with the gap opening and at a first glance is inconsistent with the fact that FeSe remains a metal below $T_n$, with sizable hole and electron pockets. We argue that the drop of $R_{B_{1g}}$ in a nematic metal comes about because the nematic order drastically changes the orbital content of the pockets and makes them nearly mono-orbital. In this situation $B_{1g}$ Raman response gets reduced by the same vertex corrections that enforce charge conservation. The reduction holds at low frequencies and gives rise to gap-like behavior of $R_{B_{1g}}$, in full agreement with the experimental data.

preprint2019arXiv

The interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal

Near a quantum-critical point, a metal reveals two competing tendencies: destruction of fermionic coherence and attraction in one or more pairing channels. We analyze the competition within Eliashberg theory for a class of quantum-critical models with an effective dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$-model) for $0 < γ<1$. We argue that the two tendencies are comparable in strength, yet the one towards pairing is stronger, and the ground state is a superconductor. We show, however, that there exist two distinct regimes of system behavior below the onset temperature of the pairing $T_p$. In the range $T_{cr} < T < T_p$ fermions remain incoherent and the density of states $N(ω)$ displays "gap filling" behavior in which the position of the maximum in $N(ω)$ is set by temperature rather than the pairing gap. At lower $T < T_{cr}$, fermions acquire coherence, and $N(ω)$ displays a conventional "gap closing" behavior. We argue that the existence of the two regimes comes about because of special behavior of fermions with frequencies $ω= \pm πT$ along the Matsubara axis. Specifically, for these fermions, the component of the self-energy, which competes with the pairing, vanishes in the normal state. We further argue that the crossover at $T \sim T_{cr}$ comes about because Eliashberg equations allow an infinite number of topologically distinct solutions for the onset temperature of the pairing within the same gap symmetry. Only one solution, with the highest $T_p$, actually emerges, but other solutions are generated and modify the form of the gap function at $T \leq T_{cr}$. Finally, we argue that the actual $T_c$ is comparable to $T_{cr}$, while at $T_{cr} < T < T_{p}$ phase fluctuations destroy superconducting long-range order, and the system displays pseudogap behavior.

preprint2016arXiv

Displacement and annihilation of Dirac gap-nodes in d-wave iron-based superconductors

Several experimental and theoretical arguments have been made in favor of a $d-$wave symmetry for the superconducting state in some Fe-based materials. It is a common belief that a $d-$wave gap in the Fe-based superconductors must have nodes on the Fermi surfaces centered at the $Γ$ point of the Brillouin zone. Here we show that, while this is the case for a single Fermi surface made out of a single orbital, the situation is more complex if there is an even number of Fermi surfaces made out of different orbitals. In particular, we show that for the two $Γ$-centered hole Fermi surfaces made out of $d_{xz}$ and $d_{yz}$ orbitals, the nodal points still exist near $T_{c}$ along the symmetry-imposed directions, but are are displaced to momenta between the two Fermi surfaces. If the two hole pockets are close enough, pairs of nodal points can merge and annihilate at some $T<T_{c}$, making the $d-$wave state completely nodeless. These results imply that photoemission evidence for a nodeless gap on the $d_{xz}/d_{yz}$ Fermi surfaces of KFe$_{2}$As$_{2}$ does not rule out $d-$wave gap symmetry in this material, while a nodeless gap observed on the $d_{xy}$ pocket in K$_{x}$Fe$_{2-y}$Se$_{2}$ is truly inconsistent with the $d-$wave gap symmetry.

preprint2016arXiv

Interplay between uni-directional and bi-directional charge-density-wave orders in underdoped cuprates

We analyze the interplay between charge-density-wave (CDW) orders with axial momenta $(Q, 0)$ and $(0,Q)$ ($Δ_x$ and $Δ_y$ respectively), detected in the underdoped cuprates. The CDW order in real space can be uni-directional (either $Δ_x$ or $Δ_y$ is non-zero) or bi-directional (both $Δ_x$ and $Δ_y$ are non-zero). To understand which of the two orders develop, we adopt the magnetic scenario, in which the CDW order appears due to spin-fluctuation exchange, and derive the Ginzburg-Landau action to the sixth order in $Δ_x$ and $Δ_y$. We argue that, at the mean-field level, the CDW order is bi-directional at the onset, with equal amplitudes of $Δ_x$ and $Δ_y$, but changes to uni-directional inside the CDW phase. This implies that, at a given temperature, CDW order is uni-directional at smaller dopings, but becomes bi-directional at larger dopings. This is consistent with recent x-ray data on YBCO, which detected tendency towards bi-directional order at larger dopings. We discuss the role of discrete symmetry breaking at a higher temperature for the interplay between bi-directional and uni-directional CDW orders and also discuss the role of pair-density-wave (PDW) order, which may appear along with CDW. We argue that PDW with the same momentum as CDW changes the structure of the bi-directional charge order by completely replacing either $Δ_x$ or $Δ_y$ CDW components by PDW. However, if an "Amperean" PDW order, which pairs fermions with approximately the same momenta, is also present, both $Δ_x$ and $Δ_y$ remain non-zero in the bi-directional phase, albeit with non-equal amplitudes. This is again consistent with x-ray experiments, which at larger doping found non-equal $Δ_x$ and $Δ_y$ in every domain.

preprint2016arXiv

Itinerant scenario for Fe-pnictides: comparison with quantum Monte Carlo

Recent applications of Quantum Monte Carlo (QMC) technique to Fe-based superconductors opened a way to directly verify the applicability of the itinerant scenario for these systems. Fe-based superconductors undergo various instabilities upon lowering temperature (magnetism, superconductivity, nematicity/orbital order), and one can check whether the hierarchy of instabilities obtained within the itinerant approach is the same as in unbiased QMC simulations. In a recent paper [arXiv:1512:08523] the authors considered the simplest two-band model with interaction tailored to favor orbital order. The type of the orbital order found in QMC is different from the one found in earlier itinerant analysis. We report the results of our calculations within the itinerant scenario and argue that they are in perfect agreement with QMC.

preprint2016arXiv

Magnetism, superconductivity, and spontaneous orbital order in iron-based superconductors: who comes first and why?

Magnetism and nematic order are the two non-superconducting orders observed in iron-based superconductors. To elucidate the interplay between them and ultimately unveil the pairing mechanism, several models have been investigated. In models with quenched orbital degrees of freedom, magnetic fluctuations promote stripe magnetism which induces orbital order. In models with quenched spin degrees of freedom, charge fluctuations promote spontaneous orbital order which induces stripe magnetism. Here we develop an unbiased approach, in which we treat magnetic and orbital fluctuations on equal footing. Key to our approach is the inclusion of the orbital character of the low-energy electronic states into renormalization group analysis. Our results show that in systems with large Fermi energies, such as BaFe2As2, LaFeAsO, and NaFeAs, orbital order is induced by stripe magnetism. However, in systems with small Fermi energies, such as FeSe, the system develops a spontaneous orbital order, while magnetic order does not develop. Our results provide a unifying description of different iron-based materials.

preprint2016arXiv

Raman resonance in iron-based superconductors: The magnetic scenario

We perform theoretical analysis of polarization-sensitive Raman spectroscopy on NaFe$_{1-x}$Co$_x$As, EuFe$2$As$_2$, SrFe$_2$As$_2$, and Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$, focusing on two features seen in the $B_{1g}$ symmetry channel (in one Fe unit cell notation): the strong temperature dependence of the static, uniform Raman response in the normal state and the existence of a collective mode in the superconducting state. We show that both features can be explained by the coupling of fermions to pairs of magnetic fluctuations via the Aslamazov-Larkin process. We first analyze magnetically-mediated Raman intensity at the leading two-loop order and then include interactions between pairs of magnetic fluctuations. We show that the full Raman intensity in the $B_{1g}$ channel can be viewed as the result of the coupling of light to Ising-nematic susceptibility via Aslamazov-Larkin process. We argue that the singular temperature dependence in the normal state is the combination of the temperature dependencies of the Aslamazov-Larkin vertex and of Ising-nematic susceptibility. We discuss two scenarios for the resonance below $T_c$. In one, the resonance is due to the development of a pole in the fully renormalized Ising-nematic susceptibility. The other is the orbital excitonic scenario, in which spin fluctuations generate an attractive interaction between low-energy fermions.

preprint2016arXiv

Superconductivity near a quantum-critical point: The special role of the first Matsubara frequency

Near a quantum-critical point in a metal strong fermion-fermion interaction mediated by a soft collective boson gives rise to incoherent, non-Fermi liquid behavior. It also often gives rise to superconductivity which masks the non-Fermi liquid behavior. We analyze the interplay between superconductivity and fermionic incoherence for a set of quantum-critical models with effective dynamical interaction between low-energy fermions. We argue that superconductivity always develops above a quantum-critical point, no matter how strong the fermionic self-energy is. We argue that superconductivity should not be viewed as a pairing of incoherent fermions, as previously thought. Rather, $T_c$ is non-zero due to the fact that the self-energy is suppressed at the two lowest fermionic Matsubara frequencies $ω_m = \pm πT$. We obtain the analytic formula for $T_c$ which reproduces well numerical results for the electron-phonon model at vanishing Debye frequency.

preprint2015arXiv

Enhancement of superconductivity at the onset of charge-density-wave order in a metal

We analyze superconductivity in the cuprates near the onset of an incommensurate charge density wave (CDW) order with momentum ${\bf Q} = (Q,0)/(0,Q)$, as observed in the experiments. We first consider a semi-phenomenological charge-fermion model in which hot fermions, separated by ${\bf Q}$, attract each other by exchanging soft CDW fluctuations. We find that in a quantum-critical region near CDW transition, $T_c = A {\bar g}_c$, where ${\bar g}_c$ is charge-fermion coupling and $A$ is the prefactor which we explicitly compute. We then consider the particular microscopic scenario in which CDW order parameter emerges as a composite field made out of primary spin-density-wave fields. We show that charge-fermion coupling ${\bar g}_c$ is of order of spin-fermion coupling ${\bar g}_s$. As the consequence, superconducting $T_c$ is substantially enhanced near the onset of CDW order. Finally we analyze the effect of an external magnetic field $H$. We show that, as $H$ increases, optimal $T_c$ decreases and the superconducting dome becomes progressively more confined to the CDW quantum-critical point. These results are consistent with the experiments.

preprint2015arXiv

Fluctuating charge order in the cuprates: spatial anisotropy and feedback from superconductivity

We analyze the form of static charge susceptibility $χ(q)$ in underdoped cuprates near axial momenta $(Q,0)$ and $(0,Q)$ at which short-range static charge order has been observed. We show that the momentum dependence of $χ(q)$ is anisotropic, and the correlation length in the longitudinal direction is larger than in the transverse direction. We show that correlation lengths in both directions decrease once the system evolves into a superconductor, as a result of the competition between superconductivity and charge order. These results are in agreement with resonant x-ray scattering data [R. Comin et al., Science 347, 1335 (2015)]. We also argue that density and current components of the charge order parameter are affected differently by superconductivity --- the charge-density component is reduced less than the current component and hence extends deeper into the superconducting state. This gives rise to two distinct charge order transitions at zero temperature.

preprint2015arXiv

Gap structure in Fe-based superconductors with accidental nodes: the role of hybridization

We study the effects of hybridization between the two electron pockets in Fe-based superconductors with s-wave gap with accidental nodes. We argue that hybridization reconstructs the Fermi surfaces and also induces an additional inter-pocket pairing component. We analyze how these two effects modify the gap structure by tracing the position of the nodal points of the energy dispersions in the superconducting state. We find three possible outcomes. In the first, the nodes simply shift their positions in the Brillouin zone; in the second, the nodes merge and disappear, in which case the gap function has either equal or opposite signs on the electron pockets; in the third, a new set of nodal points emerges, doubling the original number of nodes.

preprint2015arXiv

Itinerant electron scenario for Fe-based superconductors

I review works on Fe-based superconductors which depart from a metal with well defined Fermi surfaces and Fermi liquid-type quasiparticles. I consider normal state instabilities -- SDW magnetism and nematic order, and superconductivity, all three as the consequences of the instability of a Fermi surface due to interactions between low-energy fermionic quasiparticles. This approach assumes that renormalizations coming from fermions from high energies, of order bandwidth, modify but do not destroy Fermi liquid behavior in the normal state and can be absorbed into the effective low-energy model of interacting fermions located near hole and electron-type Fermi surfaces. I argue that the interactions between these fermions are responsible for (i) a stripe-type SDW magnetic order (and, in some cases, a checkerboard order), (ii) a preliminary nematic-type instability, in which magnetic fluctuations break $C_4$ lattice rotational symmetry down to $C_2$, but magnetic order does not yet develop, and (iii) a superconductivity, which competes with these two orders. The experimental data on superconductivity show rich behavior with potentially different symmetry of a superconducting state even for different compositions of the same material. I argue that, despite all this, the physics of superconductivity in the itinerant scenario for Fe-based materials is governed by a single underlying pairing mechanism.

preprint2015arXiv

Magnetic fluctuations and specific heat in NaxCoO2 near a Lifshitz topological transition

We analyze the temperature and doping dependence of the specific heat $C(T)$ in Na$_x$CoO$_2$. This material was conjectured to undergo a Lifshitz -type topological transition at $x =x_c =0.62$, in which a new electron Fermi pocket emerges at the $Γ$ point, in addition to the existing hole pocket with large $k_F$. The data show that near $x =x_c$, the temperature dependence of $C(T)/T$ at low $T$ gets stronger as $x$ approaches $x_c$ from below and then reverses the trend and changes sign at $x \geq x_c$. We argue that this behavior can be quantitatively explained within the spin-fluctuation theory. We show that magnetic fluctuations are enhanced near $x_c$ at momenta around $k_F$ and their dynamics changes between $x \leq x_c$ and $x >x_c$, when the new pocket forms. We demonstrate that this explains the temperature dependence of $C(T)/T$. We show that at larger $x$ ($x > 0.65$) the system enters a magnetic quantum critical regime where $C(T)/T$ roughly scales as $\log T$. This behavior extends to progressively lower $T$ as $x$ increases towards a magnetic instability at $x \approx 0.75$.

preprint2015arXiv

The origin of nematic order in FeSe

The origin of the 90 K nematic transition in the chalcogenide FeSe, which displays no magnetic order down to T=0, remains a major puzzle for a unifying theory for the iron-based superconductors. We analyze this problem in light of recent experimental data which reveal very small Fermi pockets in this material. We show that the smallness of the Fermi energy leads to a near-degeneracy between magnetic fluctuations and fluctuations in the charge-current density-wave channel. While the two fluctuation modes cooperate to promote the same preemptive Ising-nematic order, they compete for primary order. We argue that this explains why in FeSe the nematic order emerges when the magnetic correlation length is smaller than in other Fe-based materials, and why no magnetism is observed. We discuss how pressure lifts this near-degeneracy, resulting in a non-monotonic dependence of the nematic transition with pressure, in agreement with experiments.

preprint2014arXiv

Antiferromagnetism in Iron-Based Superconductors: Selection of Magnetic Order and Quasiparticle Interference

The recent discovery of superconductivity in the iron-based layered pnictides with T_c ranging between 26 and 56K generated enormous interest in the physics of these materials. Here, we review some of the peculiarities of the antiferromagnetic order in the iron pnictides, including the selection of the stripe magnetic order and the formation of the Ising-nematic state in the unfolded BZ within an itinerant description. In addition we analyze the properties of the quasiparticle interference spectrum in the parent antiferromagnetic phase.

preprint2014arXiv

Charge-density-wave order with momentum $(2Q, 0)$ and $(0, 2Q)$ within the spin-fermion model: continuous and discrete symmetry breaking, preemptive composite order, and relation to pseudogap in hole-doped cuprates

We analyze charge order within the the spin-fermion model. We show that magnetically-mediated interaction gives rise to charge order $Δ_k^Q = \langle c^\dagger_{{\bf k}+{\bf Q}} c_{{\bf k}-{\bf Q}}\rangle$ with momenta ${\bf Q}=Q_x =(2Q,0)$ and ${\bf Q}=Q_y =(0,2Q)$, if the magnetic correlation length $ξ$ exceeds some critical value. We argue that $Δ_k^Q$ and $Δ_{-k}^Q$ are not equivalent, and their symmetric and antisymmetric combinations describe density modulations and bond current. We derive GL functional for four-component $U(1)$ order parameters $Δ^Q_{\pm k}$ with ${\bf Q} = Q_x$ or $Q_y$. Within mean-field we find two types of CDW states, I and II, depending on system parameters. In state I density and current modulations emerge with the same ${\bf Q} = Q_x$ or $Q_y$, breaking $Z_2$ lattice rotational symmetry, and differ in phase by $\pmπ/2$. The selection of $π/2$ or $-π/2$ additionally breaks $Z_2$ time-reversal symmetry, such that the total order parameter manifold is $U(1) \times Z_2 \times Z_2$. In state II density and current modulations emerge with different $\bf Q$ and the order parameter manifold is $U(1) \times U (1) \times Z_2$. We go beyond mean-field and show that discrete symmetries get broken before long-range charge order sets in. For state I, the system first breaks $Z_2$ lattice rotational symmetry ($C_4 \to C_2$) at $T= T_n$ and develops a nematic order, then breaks $Z_2$ time-reversal symmetry at $T_t < T_n$, and finally breaks $U(1)$ symmetry of a common phase of even and odd components of $Δ^Q_{k}$ at $T= T_{\rm cdw} < T_t < T_n$ and develops a true charge order. We argue that the preemptive orders lift $T_{\rm cdw}$ and reduces $T_{\rm sc}$ such that at large $ξ$ charge order may develop prior to superconductivity. We obtain the phase diagram and present quantitative comparison with ARPES data for hole-doped cuprates.

preprint2014arXiv

Effect of pairing fluctuations on the spin resonance in Fe-based superconductors

The spin resonance observed in the inelastic neutron scattering data on Fe-based superconductors has played a prominent role in the quest for determining the symmetry of the order parameter in these compounds. Most theoretical studies of the resonance employ an RPA-type approach in the particle-hole channel and associate the resonance in the spin susceptibility $χ_S (\mathbf{q}, ω)$ at momentum ${\bf Q} = (π,π)$ with the spin-exciton of an $s^{+-}$ superconductor, pulled below $2Δ$ by residual attraction associated with the sign change of the gap between Fermi points connected by ${\bf Q}$. Here we explore the effect of fluctuations in the particle-particle channel on the spin resonance. Particle-particle and particle-hole channels are coupled in a superconductor and to what extent the spin resonance can be viewed as a particle-hole exciton needs to be addressed. In the case of purely repulsive interactions we find that the particle-particle channel at total momentum ${\bf Q}$ (the $π$channel) contributes little to the resonance. However, if the interband density-density interaction is attractive and the $π-$resonance is possible on its own, along with the spin-exciton, we find a much stronger shift of the resonance frequency from the position of the would-be spin-exciton resonance. We also show that the expected double-peak structure in this situation does not appear because of the strong coupling between particle-hole and particle-particle channels, and $\mathrm{Im} χ_S ({\bf Q}, ω)$ displays only a single peak.

preprint2014arXiv

Non-Landau damping of magnetic excitations in systems with localized and itinerant electrons

We discuss the form of the damping of magnetic excitations in a metal near a ferromagnetic instability. The paramagnon theory predicts that the damping term should have the form $Ω/Γ(q)$ with $Γ(q) \propto q$ (the Landau damping). However, the experiments on uranium metallic compounds UGe$_2$ and UCoGe showed that $Γ(q)$ tends to a constant value at vanishing $q$. A non-zero $Γ(0)$ is impossible in systems with one type of carriers (either localized or itinerant) because it would violate the spin conservation. It has been conjectured recently that a non-zero $Γ(q)$ in UGe$_2$ and UCoGe may be due to the presence of both localized and itinerant electrons in these materials, with ferromagnetism involving predominantly localized spins. We present microscopic analysis of the damping of near-critical localized excitations due to interaction with itinerant carriers. We show explicitly how the presence of two types of electrons breaks the cancellation between the contributions to $Γ(0)$ from self-energy and vertex correction insertions into the spin polarization bubble and discuss the special role of the Aslamazov-Larkin processes. We show that $Γ(0)$ increases with $T$ both in the paramagnetic and ferromagnetic regions, but in-between it has a peak at $T_c$. We compare our theory with the available experimental data.

preprint2014arXiv

Optical conductivity of a two-dimensional metal at the onset of spin-density-wave order

We consider the optical conductivity of a clean two-dimensional metal near a quantum spin-density-wave transition. Critical magnetic fluctuations are known to destroy fermionic coherence at "hot spots" of the Fermi surface but coherent quasiparticles survive in the rest of the Fermi surface. A large part of the Fermi surface is not really "cold" but rather "lukewarm" in a sense that coherent quasiparticles in that part survive but are strongly renormalized compared to the non-interacting case. We discuss the self-energy of lukewarm fermions and their contribution to the optical conductivity, $σ(Ω)$, focusing specifically on scattering off composite bosons made of two critical magnetic fluctuations. Recent study [S.A. Hartnoll et al., Phys. Rev. B {\bf 84}, 125115 (2011)] found that composite scattering gives the strongest contribution to the self-energy of lukewarm fermions and suggested that this may give rise to non-Fermi liquid behavior of the optical conductivity at the lowest frequencies. We show that the most singular term in the conductivity coming from self-energy insertions into the conductivity bubble, $σ'(Ω)\propto \ln^3Ω/Ω^{1/3}$, is canceled out by the vertex-correction and Aslamazov-Larkin diagrams. However, the cancelation does not hold beyond logarithmic accuracy, and the remaining conductivity still diverges as $1/Ω^{1/3}$. We further argue that the $1/Ω^{1/3}$ behavior holds only at asymptotically low frequencies, well inside the frequency range affected by superconductivity. At larger $Ω$, up to frequencies above the Fermi energy, $σ'(Ω)$ scales as $1/Ω$, which is reminiscent of the behavior observed in the superconducting cuprates.

preprint2014arXiv

Phases of triangular lattice antiferromagnet near saturation

We consider 2D Heisenberg antiferromagnets on a triangular lattice with spatially anisotropic interactions in a high magnetic field close to the saturation. We show that this system possess rich phase diagram in field/anisotropy plane due to competition between classical and quantum orders: an incommensurate non-coplanar spiral state, which is favored classically, and a commensurate co-planar state, which is stabilized by quantum fluctuations. We show that the transformation between these two states is highly non-trivial and involves two intermediate phases -- the phase with co-planar incommensurate spin order and the one with non-coplanar double-${\bf Q}$ spiral order. The transition between the two co-planar states is of commensurate- incommensurate type, not accompanied by softening of spin-wave excitations. We show that a different sequence of transitions holds in triangular antiferromagnets with exchange anisotropy, such as Ba3CoSb2O9.

preprint2014arXiv

Polar Kerr effect from chiral-nematic charge order

We analyze the polar Kerr effect in an itinerant electron system on a square lattice in the presence of a composite charge order proposed for the pseudogap state in underdoped cuprates. This composite charge order preserves discrete translational symmetries, and is "chiral-nematic" in the sense that it breaks time-reversal symmetry, mirror symmetries in $x$ and $y$ directions, and $C_4$ lattice rotation symmetry. The Kerr angle $θ_K$ in $C_4$-symmetric system is proportional to the antisymmetric component of the anomalous Hall conductivity $σ_{xy}-σ_{yx}$. We show that this result holds when $C_4$ symmetry is broken. We show that in order for $σ_{xy}$ and $σ_{yx}$ to be non-zero the mirror symmetries in $x$ and $y$ directions have to be broken, and that for $σ_{xy}-σ_{yx}$ to be non-zero time-reversal symmetry has to be broken. The chiral-nematic charge order satisfies all these conditions, such that a non-zero signal in a polar Kerr effect experiment is symmetry allowed. We further show that to get a non-zero $θ_K$ in a one-band spin-fluctuation scenario, in the absence of disorder, one has to extend the spin-mediated interaction to momenta away from $(π,π)$ and has to include particle-hole asymmetry. Alternatively, in the presence of disorder one can get a non-zero $θ_K$ from impurity scattering: either due to skew scattering (with non-Gaussian disorder) or due to particle-hole asymmetry in case of Gaussian disorder. The impurity analysis in our case is similar to that in earlier works on Kerr effect in $p_x+ip_y$ superconductor, however in our case the magnitude of $θ_K$ is enhanced by the flattening of the Fermi surface in the "hot" regions which mostly contribute to charge order.

preprint2014arXiv

Superconductivity from weak repulsion in hexagonal lattice systems

We analyze the pairing instabilities for fermions on hexagonal lattices (both honeycomb and triangular ones) in a wide range of fermionic densities. We argue that for a generic doping in this range, superconductivity at weak coupling is of Kohn-Luttinger type, and, due to the presence of electronic interactions beyond on-site repulsion, is a threshold phenomenon, with superconductivity emerging only if the attraction generated by the Kohn-Luttinger mechanism exceeds the bare repulsion in some channel. For disconnected Fermi pockets, we predict that Kohn-Luttinger superconductivity, if it occurs, is likely to be $f$-wave. While the Kohn-Luttinger analysis is adequate over most of the doping range, a more sophisticated analysis is needed near Van Hove doping. We treat Van Hove doping using a parquet renormalization group, the equations for which we derive and analyze. Near this doping level, superconductivity is a universal phenomenon, arising from any choice of bare repulsive interactions. The strongest pairing instability is into a chiral $d-$wave state ($d+id$). At a truly weak coupling, the strongest competitor is a spin-density-wave instability, however, $d-$wave superconductivity still wins. Moreover, the feedback of the spin density fluctuations into the Cooper channel significantly enhances the critical temperature over the estimates of the Kohn Luttinger theory. We analyze renormalization group equations at stronger couplings and find that the main competitor to $d-$wave supoerconductivity away from weak coupling is actually ferromagnetism. We also discuss the effect of the edge fermions and show that they are unimportant in the asymptotic weak coupling limit, but may give rise to, e.g., a charge-density-wave order at moderate coupling strengths.

preprint2014arXiv

Time-reversal symmetry breaking superconductivity in the coexistence phase with magnetism in Fe-pnictides

We argue that superconductivity in the coexistence region with spin-density-wave (SDW) order in weakly doped Fe-pnictides differs qualitatively from the ordinary $s^{+-}$ state outside the coexistence region, as it develops an additional gap component which is a mixture of intra-pocket singlet ($s^{++}$) and inter-pocket spin-triplet pairings (the $t-$state). The coupling constant for the $t-$channel is proportional to the SDW order and involves interactions that do not contribute to superconductivity outside of the SDW region. We argue that the $s^{+-}$ and $t-$type superconducting orders coexist at low temperatures, and the relative phase between the two is in general different than $0$ or $π$, manifesting explicitly the breaking of the time-reversal symmetry promoted by long-range SDW order. We show that this exotic state emerges already in the simplest model of Fe-pnictides, with one hole pocketand two symmetry-related electron pockets. We argue that in some parameter range time-reversal gets broken even before long-range superconducting order develops.

preprint2013arXiv

How many quantum phase transitions exist inside the superconducting dome of the iron pnictides?

Recent experiments on two iron-pnictide families suggest the existence of a single quantum phase transition (QPT) inside the superconducting dome despite the fact that two separate transition lines - magnetic and nematic - cross the superconducting dome at $T_{c}$. Here we argue that these two observations are actually consistent. We show, using a microscopic model, that each order coexists with superconductivity for a wide range of parameters, and both transition lines continue into the superconducting dome below $T_{c}$. However, at some $T_{\mathrm{merge}}<T_{c}$, the two transitions merge and continue down to T=0 as a single simultaneous first-order nematic/magnetic transition. We show that superconductivity has a profound effect on the character of this first-order transition, rendering it weakly first-order and allowing strong fluctuations to exist near the QPT.

preprint2013arXiv

s+is State with Broken Time Reversal Symmetry in Fe-Based Superconductors

We analyze the evolution of the superconducting gap structure in strongly hole doped Ba_{1-x}K_xFe_2As_2 between x=1 and x ~ 0.4 (optimal doping). In the latter case, the pairing state is most likely s+-, with different gap signs on hole and electron pockets, but with the same signs of the gap on the two Gamma-centered hole pockets (a ++ state on hole pockets). In a pure KFe_2As_2 (x=1), which has only hole pockets, laser ARPES data suggested another s+- state, in which the gap changes sign between hole pockets (a +- state). We analyze how ++ gap transforms into a +- gap as x -> 1. We found that this transformation occurs via an intermediate s+is, state in which the gaps on the two hole pockets differ in phase by `phi', which gradually involves from `phi' = pi (the +- state) to phi =0 (the ++ state). This state breaks time-reversal symmetry and has huge potential for applications. We compute the dispersion of collective excitations and show that two different Leggett-type phase modes soften at the two end points of TRSB state.

preprint2013arXiv

Spin-current order in anisotropic triangular antiferromagnets

We analyze instabilities of the collinear up-up-down state of a two-dimensional quantum spin-S spatially anisotropic triangular lattice antiferromagnet in a magnetic field. We find, within large-S approximation, that near the end point of the plateau, the collinear state becomes unstable due to condensation of two-magnon bound pairs rather than single magnons. The two-magnon instability leads to a novel 2D vector chiral phase with alternating spin currents but no magnetic order in the direction transverse to the field. This phase breaks a discrete Z_2 symmetry but preserves a continuous U(1) one of rotations about the field axis. It possesses orbital antiferromagnetism and displays a magnetoelectric effect.

preprint2013arXiv

Superconductivity at the onset of spin-density-wave order in a metal

We revisit the issue of superconductivity at the quantum-critical point (QCP) between a 2D paramagnet and a spin-density-wave metal with ordering momentum (π,π). This problem is highly non-trivial because the system at criticality displays a non-Fermi liquid behavior and because the effective coupling constant λfor the pairing is generally of order one, even when the actual interaction is smaller than fermionic bandwidth. Previous study [M. A. Metlitski, S. Sachdev, Phys.Rev.B 82, 075128 (2010)] has found that the renormalizations of the pairing vertex are stronger than in BCS theory and hold in powers of \log^2 (1/T), like in color superconductivity. We analyze the full gap equation and argue that, for QCP problem, summing up of the leading logarithms does not lead to a pairing instability. Yet, we show that superconductivity has no threshold and appears even if the coupling is set to be small, because subleading logarithmical renormalizations diverge and give rise to BCS-like \log(1/T_c) \propto 1/λ. We argue that the analogy with BCS is not accidental as at small coupling superconductivity at a QCP predominantly comes from fermions which retain Fermi liquid behavior at criticality. We computed T_c for the actual λ\sim 1, and found that both Fermi-liquid and non-Fermi liquid fermions contribute to the pairing. The value of T_c agrees well with the numerical results.

preprint2013arXiv

Superconductivity from repulsive interaction

The BCS theory of superconductivity named electron-phonon interaction as a glue that overcomes Coulomb repulsion and binds fermions into pairs which then condense and superconduct. We review recent and not so recent works aiming to understand whether a nominally repulsive Coulomb interaction can by itself give rise to a superconductivity. We first discuss a generic scenario of the pairing by electron-electron interaction, put forward by Kohn and Luttinger back in 1965, and then turn to modern studies of the electronic mechanism of superconductivity in the lattice models for the cuprates, the Fe-pnictides, and the doped graphene. We show that the pairing in all three classes of materials can be viewed as lattice version of Kohn-Luttinger physics, despite that the pairing symmetries are different. We discuss under what conditions the pairing occurs and rationalize the need to do parquet renormalization-group analysis. We also discuss the interplay between superconductivity and density-instabilities.

preprint2012arXiv

Broken translational symmetry in an emergent paramagnetic phase of graphene

We show that the spin-density wave state on the partially filled honeycomb and triangular lattices is preempted by a paramagnetic phase that breaks an emergent Z_4 symmetry of the system, associated with the four inequivalent arrangements of spins in the quadrupled unit cell. Unlike other emergent paramagnetic phases in itinerant and localized-spin systems, this state preserves the rotational symmetry of the lattice but breaks its translational symmetry, giving rise to a super-lattice structure that can be detected by scanning tunneling microscopy. This emergent phase also has distinctive signatures in the magnetic spectrum that can be probed experimentally.

preprint2012arXiv

First-Matsubara-frequency rule in a Fermi liquid. Part I: Fermionic self-energy

We analyze in detail the fermionic self-energy Σ(ω, T) in a Fermi liquid (FL) at finite temperature T and frequency ω. We consider both canonical FLs -- systems in spatial dimension D >2, where the leading term in the fermionic self-energy is analytic [the retarded ImΣ^R(ω,T) = C(ω^2 +π^2 T^2)], and non-canonical FLs in 1<D <2, where the leading term in ImΣ^R(ω,T) scales as T^D or ω^D. We relate the ω^2 + π^2 T^2 form to a special property of the self-energy -"the first-Matsubara-frequency rule", which stipulates that Σ^R(iπT,T) in a canonical FL contains an O(T) but no T^2 term. We show that in any D >1 the next term after O(T) in Σ^R(iπT,T) is of order T^D (T^3\ln T in D=3). This T^D term comes from only forward- and backward scattering, and is expressed in terms of fully renormalized amplitudes for these processes. The overall prefactor of the T^D term vanishes in the "local approximation", when the interaction can be approximated by its value for the initial and final fermionic states right on the Fermi surface. The local approximation is justified near a Pomeranchuk instability, even if the vertex corrections are non-negligible. We show that the strength of the first-Matsubara-frequency rule is amplified in the local approximation, where it states that not only the T^D term vanishes but also that Σ^R(iπT,T) does not contain any terms beyond O(T). This rule imposes two constraints on the scaling form of the self-energy: upon replacing ωby iπT, ImΣ^R(ω,T) must vanish and ReΣ^R (ω, T) must reduce to O(T). These two constraints should be taken into consideration in extracting scaling forms of Σ^R(ω,T) from experimental and numerical data.

preprint2012arXiv

Interplay of superconductivity and spin density wave order in doped graphene

We study the interplay between superconductivity and spin density wave order in graphene doped to 3/8 or 5/8 filling (a Van Hove doping). At this doping level, the system is known to exhibit weak coupling instabilities to both chiral d + id superconductivity and to a uniaxial spin density wave. Right at van Hove doping, the superconducting instability is strongest and emerges at the highest Tc, but slightly away from van-Hove doping a spin-density-wave likely emerges first. We investigate whether at some lower temperature superconductivity and spin-density-waves co-exist. We derive the Landau-Ginzburg functional describing interplay of the two order parameters. Our calculations show that superconductivity and spin density wave order do not co-exist and are separated by first-order transitions, either as a function of doping or as a function of T.

preprint2012arXiv

Itinerant half-metal spin-density-wave state on the hexagonal lattice

We consider electrons on a honeycomb or triangular lattice doped to the saddle point of the bandstructure. We assume system parameters are such that spin density wave (SDW) order emerges below a temperature $T_N$ and investigate the nature of the SDW phase. We argue that at $T < T_N$ the system develops a uniaxial SDW phase whose ordering pattern breaks $O(3) \times Z_4$ symmetry and corresponds to an eight site unit cell with non-uniform spin moments on different sites. This state is a half-metal -- it preserves full original Fermi surface, but has gapless charged excitations in one spin branch only. It allows for electrical control of spin currents and is desirable for nano-science.

preprint2011arXiv

Relation between 2Delta/Tc and nodes in Fe-based superconductors

We analyze the interplay between the absence or presence of the nodes in the superconducting gap along electron Fermi surfaces (FSs) in Fe-pnictides, and 2Delta/T_c along hole FSs, measured by ARPES. We solve the set of coupled gap equations for 4-band and 5-band models of Fe-pnictides and relate the presence of the nodes to 2Δ/T_c being below a certain threshold. Using ARPES data for 2Δ/T_c, we find that optimally doped (BaK)Fe_2As_2 and Ba(FeCo)_2As_2 and likely nodeless, but iso-valent BaFe_2(AsP)_2 likely has nodes.

preprint2011arXiv

Renormalization group flow, competing phases, and the structure of superconducting gap in multi-band models of Iron based superconductors

We perform an analytical renormalization group (RG) study to address the role of Coulomb repulsion, the competition between extended s-wave superconducting order (s+-) and the spin-density wave (SDW) order and the angular dependence of the superconducting gap in multi-pocket models of Iron based superconductors. Previous analytic RG studies considered a toy model of one hole and one electron pocket. We consider more realistic models of two electron pockets and either two or three hole pockets, and also incorporate the angular dependence of interaction. We neglect for simplicity one of the two hole pockets centered at k=0, which is less nested with electron pockets, i.e., consider 3-pocket and 4-pocket models. In a toy 2-pocket model, SDW order always wins over s+- order at perfect nesting, and s+- order only appears at a finite doping, and only if RG flow extends long enough to overcome intra-pocket Coulomb repulsion. For multi-pocket models, we find two new effects. First, there always exists an attractive component of the interaction in s+- channel, such that the system necessary becomes a superconductor once it overcomes the competition from the SDW state. Second, in 3-pocket case (but not in 4-pocket case), there are situations when s+- order wins over SDW order even for perfect nesting, suggesting that SDW order is not a necessary pre-condition for the s+- order. Our results are in good agreement with recent numerical functional RG studies by Thomale et al. [arXiv:1002.3599]

preprint2011arXiv

Resistivity of a non-Galilean--invariant Fermi Liquid near Pomeranchuk Quantum Criticality

We analyze the effect of the electron-electron interaction on the resistivity of a metal near a Pomeranchuk quantum phase transition (QPT). We show that Umklapp processes are not effective near a QPT, and one must consider both interactions and disorder to obtain finite and T dependent resistivity. By power counting, the correction to the residual resistivity at low T scales as AT^{(D+2)/3} near a Z=3 QPT. We show, however, that A=0 for a simply connected, convex Fermi surface in 2D, due to hidden integrability of the electron motion. We argue that A >0 in a two-band (s-d) model and propose this model as an explanation for the observed T^{(D+2)/3} behavior.

preprint2010arXiv

Electron self-energy near a nematic quantum critical point

We consider an isotropic Fermi liquid in two dimensions near the n=2 Pomeranchuk instability in the charge channel. The order parameter is a quadrupolar stress tensor with two polarizations, longitudinal and transverse to the quadrupolar momentum tensor. Longitudinal and transverse bosonic modes are characterized by dynamical exponents z_parallel=3 and z_perp=2, respectively. Previous studies have found that such a system exhibits multiscale quantum criticality with two different energy scales omega ~ xi^{-z_{parallel,perp}}, where xi is the correlation length. We study the impact of the multiple energy scales on the electron Green function. The interaction with the critical z_parallel =3 mode is known to give rise to a local self-energy that develops a non-Fermi liquid form, Sigma(omega) ~ omega^{2/3} for frequencies larger than the energy scale omega ~ xi^{-3}. We find that the exchange of transverse z_perp=2 fluctuations leads to a logarithmically singular renormalizations of the quasiparticle residue Z and the vertex Gamma. We derive and solve renormalization group equations for the flow of Z and Gamma and show that the system develops an anomalous dimension at the nematic quantum-critical point (QCP). As a result, the spectral function at a fixed omega and varying k has a non-Lorentzian form. Away from the QCP, we find that the flow of Z is cut at the energy scale omega_{FL} ~ xi^{-1}, associated with the z=1 dynamics of electrons. The z_perp=2 energy scale, omega ~ xi^{-2}, affects the flow of Z only if one includes into the theory self-interaction of transverse fluctuations.

preprint2010arXiv

Optical Integral and Sum Rule Violation

The purpose of this work is to investigate the role of the lattice in the optical Kubo sum rule in the cuprates. We compute conductivities, optical integrals W, and ΔW between superconducting and normal states for 2-D systems with lattice dispersion typical of the cuprates for four different models -- a dirty BCS model, a single Einstein boson model, a marginal Fermi liquid model, and a collective boson model with a feedback from super-conductivity on a collective boson. The goal of the paper is two-fold. First, we analyze the dependence of W on the upper cut-off w_c placed on the optical integral because in experiments W is measured up to frequencies of order bandwidth. For a BCS model, the Kubo sum rule is almost fully reproduced at w_c equal to the bandwidth. But for other models only 70%-80% of Kubo sum rule is obtained up to this scale and even less so for ΔW, implying that the Kubo sum rule has to be applied with caution. Second, we analyze the sign of ΔW. In all models we studied ΔW is positive at small w_c, then crosses zero and approaches a negative value at large w_c, i.e. the optical integral in a superconductor is smaller than in a normal state. The point of zero crossing, however, increases with the interaction strength and in a collective boson model becomes comparable to the bandwidth at strong coupling. We argue that this model exhibits the behavior consistent with that in the cuprates.

preprint2010arXiv

Quantum-critical pairing with varying exponents

We analyse the onset temperature T_p for the pairing in cuprate superconductors at small doping, when tendency towards antiferromagnetism is strong. We consider the model of Moon and Sachdev (MS), which assumes that electron and hole pockets survive in a paramagnetic phase. Within this model, the pairing between fermions is mediated by a gauge boson, whose propagator remains massless in a paramagnet. We relate the MS model to a generic γ-model of quantum-critical pairing with the pairing kernel λ(Ω) \propto 1/Ω^γ. We show that, over some range of parameters, the MS model is equivalent to the γ-model with γ=1/3 (λ(Ω) \propto Ω^{-1/3}). We find, however, that the parameter range where this analogy works is bounded on both ends. At larger deviations from a magnetic phase, the MS model becomes equivalent to the γ-model with varying γ>1/3, whose value depends on the distance to a magnetic transition and approaches γ=1 deep in a paramagnetic phase. Very near the transition, the MS model becomes equivalent to the γ-model with varying γ<1/3. Right at the magnetic QCP, the MS model is equivalent to the γ-model with γ=0+ (λ(Ω) \propto \log Ω), which is the model for color superconductivity. Using this analogy, we verified the formula for T_c derived for color superconductivity.

preprint2010arXiv

Self-consistent Eliashberg theory, Tc, and the gap function in electron-doped cuprates

We consider normal state properties, the pairing instability temperature, and the structure of the pairing gap in electron-doped cuprates. We assume that the pairing is mediated by collective spin excitations, with antiferromagnetism emerging with the appearance of hot spots. We use a low-energy spin-fermion model and Eliashberg theory up to two-loop order. We justify ignoring vertex corrections by extending the model to N >>1 fermionic flavors, with 1/N playing the role of a small Eliashberg parameter. We argue, however, that it is still necessary to solve coupled integral equations for the frequency dependent fermionic and bosonic self-energies, both in the normal and superconducting state. Using the solution of the coupled equations, we find an onset of d-wave pairing at Tc ~ 30 K, roughly three times larger than the one obtained previously [P. Krotkov and A. Chubukov, Phys. Rev. B 74, 014509 (2006)], where it was assumed that the equations for fermionic and bosonic self-energies decouple in the normal state. To obtain the momentum and frequency dependent d-wave superconducting gap D(k,w), we derive and solve the non-linear gap equation together with the modified equation for the bosonic self energy which below Tc depends on the pairing gap. We find that the gap is a non-monotonic function of momentum along the Fermi surface, with its node along the zone diagonal and its maximum some distance away from it. We obtain 2D/Tc ~ 4. We argue that the value of Tc, the non-monotonicity of the gap, and 2D/Tc ratio are all in good agreement with the experimental data on electron-doped cuprates.

preprint2010arXiv

Universal and non-universal renormalizations in Fermi liquids

We discuss an interplay between the Fermi-liquid (FL) theory and diagrammatic perturbative approach to interacting Fermi systems. In the FL theory for Galilean-invariant systems, mass renormalization $m^*/m$ comes exclusively from fermions at the Fermi surface. We show that in a diagrammatic perturbation theory the same result for $m^*/m$ comes from fermions both at and away from the Fermi surface. The equivalence of the FL and pertubative approaches is based on a particular relation between self-energy contributions from high- and low-energy fermions. We argue that care has to be exercised in the renormalization group approach to a FL in order not to miss the high-energy contribution to $m^*/m$. As particular examples, we discuss $m^*/m$ and the quasiparticle residue $Z$ for 2D and 3D systems with both SU(2) and SU(N) symmetries, and with a short-range interaction. We derive an expression for the anisotropic part of the Fermi-liquid vertex in the large-$N$ limit of the SU(N) case.

preprint2009arXiv

Fermi liquid near Pomeranchuk quantum criticality

We analyze the behavior of an itinerant Fermi system near a charge nematic(n=2) Pomeranchuk instability in terms of the Landau Fermi liquid (FL) theory. The main object of our study is the fully renormalized vertex function $ΓΩ$, related to the Landau interaction function. We derive $Γ^Ω$ for a model case of the long-range interaction in the nematic channel. Already within the Random Phase Approximation (RPA), the vertex is singular near the instability. The full vertex, obtained by resumming the ladder series composed of the RPA vertices, differs from the RPA result by a multiplicative renormalization factor $Z_Γ$, related to the single-particle residue $Z$ and effective mass renormalization $m^*/m$. We employ the Pitaevski-Landau identities, which express the derivatives of the self-energy in terms of $Γ^Ω$, to obtain and solve a set of coupled non-linear equations for $Z_Γ$, $Z$, and $m^*/m$. We show that near the transition the system enters a critical FL regime, where $Z_Γ\sim Z \propto (1 + g_{c,2})^{1/2}$ and $m^*/m \approx 1/Z$, where $g_{c,2}$ is the $n=2$ charge Landau component which approaches -1 at the instability. We construct the Landau function of the critical FL and show that all but $g_{c,2}$ Landau components diverge at the critical point. We also show that in the critical regime the one-loop result for the self-energy $Σ(K) \propto \int dP G(P) D (K-P)$ is asymptotically exact if one identifies the effective interaction $D$ with the RPA form of $Γ^Ω$.

preprint2009arXiv

Magnetic degeneracy and hidden metallicity of the spin density wave state in ferropnictides

We analyze spin density wave (SDW) order in iron-based superconductors and electronic structure in the SDW phase. We consider an itinerant model for Fe-pnictides with two hole bands centered at $(0,0)$ and two electron bands centered at $(0,π)$ and $(π,0)$ in the unfolded BZ. A SDW order in such a model is generally a combination of two components with momenta $(0,π)$ and $(π,0)$, both yield $(π,π)$ order in the folded zone. Neutron experiments, however, indicate that only one component is present. We show that $(0,π)$ or $(π,0)$ order is selected if we assume that only one hole band is involved in the SDW mixing with electron bands. A SDW order in such 3-band model is highly degenerate for a perfect nesting and hole-electron interaction only, but we show that ellipticity of electron pockets and interactions between electron bands break the degeneracy and favor the desired $(0,π)$ or $(π,0)$ order. We further show that stripe-ordered system remains a metal for arbitrary coupling. We analyze electronic structure for parameters relevant to the pnictides and argue that the resulting electronic structure is in good agreement with ARPES experiments. We discuss the differences between our model and $J_1-J_2$ model of localized spins.

preprint2009arXiv

Spin conservation and Fermi liquid near a ferromagnetic quantum critical point

We propose a new low-energy theory for itinerant fermions near a ferromagnetic quantum critical point. We show that the full low-energy model includes, in addition to conventional interaction via spin fluctuations, another type of interaction, whose presence is crucial for the theory to satisfy SU(2) spin conservation. We demonstrate the consistency between a loop-wise expansion and a Fermi liquid description for the full model. We further show that, prior to the ferromagnetic instability, the system develops a Pomeranchuk-type instability into a state with zero magnetization but with p-wave deformations of the Fermi surfaces of spin-up and -down electrons (a spin nematic).

preprint2006arXiv

Nonanalytic Magnetic Response of Fermi- and non-Fermi Liquids

We study the nonanalytic behavior of the static spin susceptibility of 2D fermions as a function of temperature and magnetic field. For a generic Fermi liquid, χ_s (T, H)= const+c_1 max (T,μ_B|H|), where c_1 is shown to be expressed via complicated combinations of the Landau parameters, rather than via the backscattering amplitude, contrary to the case of the specific heat. Near a ferromagnetic quantum critical point, the field dependence acquires a universal form χ^{-1}_s(H)= const- c_2|H|^{3/2}$, with $c_2>0$. This behavior implies a first-order transition into a ferromagnetic state. We establish a criterion for such a transition to win over the transition into an incommensurate phase.

preprint2005arXiv

Condensation energy in Eliashberg theory -- from weak to strong coupling

We consider two issues related to the condensation energy in superconductors described by the Eliashberg theory for various forms of the pairing interaction, associated either with phonon or electronic mechanisms of superconductivity. First, we derive a leading correction to the BCS formula for the condensation energy to first order in the coupling $λ$. Second, we show that at a given $λ$, the value of the condensation energy strongly depends on the functional form of the effective pairing interaction $Γ(ω)$.

preprint2005arXiv

Self-generated locality near a ferromagnetic quantum-critical point

We analyze the behavior of interacting fermions near a ferromagnetic Stoner instability. We show that the Landau damping of the spin susceptibility is a relevant perturbation near a ferromagnetic quantum-critical point (FQCP). We argue that, as the system approaches a FQCP, the fermionic self-energy crosses over from predominantly momentum dependent away from the transition to predominantly frequency dependent in the immediate vicinity of the transition. We argue that due to this self-generated locality, the quasiparticle effective mass does not diverge before a FQCP is reached.

preprint2005arXiv

Strong coupling superconductivity due to massless boson exchange

We solve the problem of fermionic pairing mediated by a massless boson in the limit of large coupling constant. At weak coupling, the transition temperature is exponentially small and superconductivity is robust against phase fluctuation. In the strong coupling limit, the pair formation occurs at a temperature of the order of the Fermi energy, however, the actual transition temperature is much smaller due to phase and amplitude fluctuations of the pairing gap. Our model calculations describe superconductivity due to color magnetic interactions in quark matter and in systems close to a ferromagnetic quantum critical point with Ising symmetry. Our strong coupling results are, however, more general and can be applied to other systems as well, including the antiferromagnetic exchange in 2D used for description of the cuprates.

preprint2004arXiv

Ward identities for strongly coupled Eliashberg theories

We discuss Ward identities for strongly interacting fermion systems described by Eliashberg-type theories. We show that Ward identities are not in conflict with Migdal theorem. We derive diagrammatically Ward identity for a charge vertex in a Fermi liquid, and when a Fermi liquid is destroyed at quantum criticality. We argue that Ward identity for a spin vertex cannot be obtained within Eliashberg theory.

preprint2003arXiv

Quantum-critical behavior of itinerant ferromagnets

We study the stability of the Quantum Critical Point (QCP) for itinerant ferromagnets commonly described by the Hertz-Millis-Moriya (HMM) theory. We argue that in $D \leq 3$, long-range spatial correlations associated with the Landau damping of the order parameter field generate a universal {\it negative}, non-analytic $|q|^{(D+1)/2}$ contribution to the static magnetic susceptibility $χ_s (q, 0)$, which makes HMM theory unstable. We argue that the actual transition is either towards incommensurate ordering, or first order. We also show that singular corrections are specific to the spin problem, while charge susceptibility remains analytic at criticality.

preprint2001arXiv

Coherent vs incoherent pairing in 2D systems near magnetic instability

We study the superconductivity in 2D fermionic systems near antiferromagnetic instability, assuming that the pairing is mediated by spin fluctuations. This pairing involves fully incoherent fermions and diffusive spin excitations. We show that the competition between fermionic incoherence and strong pairing interaction yields the pairing instability temperature $T_{ins}$ which increases and saturates as the magnetic correlation length $ξ\to \infty$. We argue that in this quantum-critical regime the pairing problem is qualitatively different from the BCS one.

preprint2001arXiv

The differential sum rule for the relaxation rate in the cuprates

Motivated by recent experiments by Basov {\it et al}, we study the differential sum rule for the effective scattering rate $1/τ(ω)$. We show that in a dirty BCS superconductor, the area under $1/τ(ω)$ does not change between the normal and the superconducting states. For magnetically mediated pairing, similar result holds between $T<T_c$ and $T \geq T_c$, while in the pseudogap phase, $1/τ(ω)$ is just suppressed compared to $1/τ(ω)$ in the normal state. We argue that this violation of the differential sum rule is due to the fact that low-energy fermions in the pseudogap phase are still incoherent.

preprint2000arXiv

Fingerprints of spin-fermion pairing in cuprates

We demonstrate that the feedback effect from bosonic excitations on fermions, which in the past allowed one to verify the phononic mechanism of a conventional, $s-$wave superconductivity, may also allow one to experimentally detect the ``fingerprints'' of the pairing mechanism in cuprates. We argue that for spin-mediated $d-$wave superconductivity, the fermionic spectral function, the density of states, the tunneling conductance through an insulating junction, and the optical conductivity are affected by the interaction with collective spin excitations, which below $T_c$ are propagating, magnon-like quasiparticles with gap $Δ_s$. We show that the interaction with a propagating spin excitation gives rise to singularities at frequencies $Δ+ Δ_s$ for the spectral function and the density of states, and at $2Δ+ Δ_s$ for tunneling and optical conductivities, where $Δ$ is the maximum value of the $d-$wave gap. We further argue that recent optical measurements also allow one to detect subleading singularities at $4Δ$ and $2Δ+ 2Δ_s$. We consider the experimental detection of these singularities as a strong evidence in favor of the magnetic scenario for superconductivity in cuprates.

preprint2000arXiv

Spectral function and conductivity in the normal state of the cuprates: a spin fluctuation study

We study the spectral function $A_k (ω)$ and the optical conductivity $σ_1 (ω)$ for a system of fermions interacting with their collective spin fluctuations (the spin-fermion model), and apply the results to optimally doped cuprates in the normal state. We show that there is no qualitative distinction between hot and cold regions in the Brillouin zone - in both cases, there exists a wide range of frequencies where the width of the peak in $A_k (ω)$ scales linearly with $ω$. We demonstrate that $σ_1 (ω)$ is inversely linear in $ω$ up to very high frequencies. We argue that these results agree {\it quantitatively} with the photoemission and optical data for $Bi2212$.

preprint1997arXiv

Temperature variation of the pseudogap in underdoped cuprates

We consider thermal evolution of the spectral function A(Omega, T) in underdoped cuprates at T >Tc. We find that in the strong coupling limit, the fermionic Green's function near (0,π) behaves as G^{-1} (Omega) ~ \sqrt{Omega} in the whole frequency range relevant to photoemission experiments. The analysis of the pairing problem with this form of G yields a leading edge gap in A(Ω) and a broad maximum at larger frequencies. We find that classical fluctuations predominantly add a frequency independent term to G^{-1} (Omega). This term destroys the leading edge gap at T =T* (k)$which has a maximum value of about 200K near (0,π).

preprint1996arXiv

Crossover from $O(3)$ to $O(4)$ behavior in weakly frustrated antiferromagnets

We consider an anisotropic version of the $CP^{1}$ model which describes frustrated quantum antiferromagnets with incommensurate spin correlations. We extend the two-component spinon field, describing lattice spins, to the $M$-component complex vector, and show, in the $1/M$ expansion, that for arbitrary small incommensurability longitudinal and transverse stiffnesses tend to the same value as the system approaches the quantum critical point. For physical spins ($M=2$), this yields $O(4)$ critical behavior. However, if the spin structure is commensurate, the longitudinal stiffness is identically zero. In this case, the critical behavior is the same as in $O(3)$ sigma model. We show how the critical exponents interpolate between $O(3)$ and $O(4)$ values near the transition. We also show that the competition between these two fixed points leads to a confinement-deconfinement transition at a finite temperature.

preprint1995arXiv

Magnetic Phases of the 2D Hubbard Model at Low Doping

We study the equilibrium spin configuration of the 2D Hubbard model at low doping, when a long-range magnetic order is still present. We use the spin-density-wave formalism and examine the two suggested candidates for a ground state: the planar spiral $(π,Q)$ and $(Q,Q)$ states. In a mean-field theory, each of these two states is a ground state in a certain range of parameters. We however show that both states have negative bosonic modes in a transverse channel which implies that the spiral states are in fact unstable (the $(Q,Q)$ state, in addition, has a negative longitudinal stiffness). We searched for a ground state in a parameter range which on a mean-field level favored the $(π,Q)$ phase, and found that the equilibrium spin configuration is noncoplanar. These findings lead to a novel scenario of spin reorientation upon doping in Hubbard antiferromagnets.

preprint1995arXiv

On the confinement of spinons in the $CP^{M-1}$ model

We use the $1/M$ expansion for the $CP^{M-1}$ model to study the long-distance behaviour of the staggered spin susceptibility in the commensurate, two-dimensional quantum antiferromagnet at finite temperature. At $M=\infty$ this model possesses deconfined spin-1/2 bosonic spinons (Schwinger bosons), and the susceptibility has a branch cut along the imaginary $k$ axis. We show that in all three scaling regimes at finite $T$, the interaction between spinons and gauge field fluctuations leads to divergent $1/M$ corrections near the branch cut. We identify the most divergent corrections to the susceptibility at each order in $1/M$ and explicitly show that the full static staggered susceptibility has a number of simple poles rather than a branch cut. We compare our results with the $1/N$ expansion for the $O(N)$ sigma-model.

preprint1995arXiv

Phase Transition, Longitudinal Spin Fluctuations and Scaling in a Two-Layer Antiferromagnet

We consider a two-layer Heisenberg antiferromagnet which can be either in the Néel-ordered or in the disordered phase at $T=0$, depending on the ratio of the intra- and interlayer exchange constants. We reduce the problem to an interacting Bose-gas and study the sublattice magnetization and the transverse susceptibility in the ordered phase, and the spectrum of quasiparticle excitations in both phases. We compare the results with the spin-wave theory and argue that the longitudinal spin fluctuations, which are not included in the spin-wave description, are small at vanishing coupling between the layers, but increase as the system approaches the transition point. We also compute the uniform susceptibility at the critical point to order $O(T^2)$, and show that the corrections to scaling are numerically small, and the linear behavior of $χ_u$ extends to high temperatures. This is consistent with the results of the recent Monte-Carlo simulations by Sandvik and Scalapino.

preprint1995arXiv

Quantum-Critical Behavior in a Two-Layer Antiferromagnet

We analyze quantum Monte Carlo data in the vicinity of the quantum transition between a Neel state and a quantum paramagnet in a two-layer, square lattice spin 1/2 Heisenberg antiferromagnet. The real-space correlation function and the universal amplitude ratio of the structure factor and the dynamic susceptibility show clear evidence of quantum critical behavior at low temperatures. The numerical results are in good quantitative agreement with $1/N$ calculations for the $O(N)$ non-linear sigma model. A discrepancy, reported earlier, between the critical properties of the antiferromagnet and the sigma model is resolved. We also discuss the values of prefactors of the dynamic susceptibility and the structure factor in a single layer antiferromagnet at low $T$.

preprint1995arXiv

Resonant two-magnon Raman scattering in parent compounds of high-T$_c$ superconductors.

We propose a theory of two-magnon Raman scattering from the insulating parent compounds of high-T$_c$ superconductors, which contains information not only on magnetism, but also on the electronic properties in these materials. We use spin density wave formalism for the Hubbard model, and study diagrammatically the profile of the two-magnon scattering and its intensity dependence on the incoming photon frequency $ω_i$ both for $ω_i \ll U$ and in the resonant regime, in which the energy of the incident photon is close to the gap between conduction and valence bands. In the nonresonant case, we identify the diagrams which contribute to the conventional Loudon-Fleury Hamiltonian. In the resonant regime, where most of the experiments have been done, we find that the dominant contribution to Raman intensity comes from a different diagram, one which allows for a simultaneous vanishing of all three of its denominators (i.e., a triple resonance). We study this diagram in detail and show that the triple resonance, combined with the spin-density-wave dispersion relation for the carriers, explains the unusual features found in the two-magnon profile and in the two-magnon peak intensity dependence on the incoming photon frequency. In particular, our theory predicts a maximum of the two-magnon peak intensity right at the upper edge of the features in the optical data, which has been one of the key experimental puzzles.