---
title: Bardasis–Schrieffer Mode in Superconductors
url: https://www.emergentmind.com/topics/bardasis-schrieffer-mode
type: topic
---

# Bardasis–Schrieffer Mode in Superconductors

Searching arXiv for recent Bardasis-Schrieffer mode papers to support the article.
The Bardasis–Schrieffer mode is a collective, exciton-like bound state associated with fluctuations in a subdominant attractive pairing channel of a superconductor or, more generally, of a condensate with competing internal channels. In the canonical superconducting setting, a dominant condensate forms in one symmetry channel while a residual attraction survives in another; the latter does not condense, but supports a sharp in-gap excitation below the pair-breaking threshold \(2\Delta\) [2011.01042]. In an \(s\)-wave ground state, a \(d\)-wave subdominant channel yields the standard Bardasis–Schrieffer mode; in a \(d\)-wave ground state, an \(s\)-wave residual channel can generate an analogous excitation, although nodal quasiparticles can strongly damp it unless additional structure, such as nematicity, suppresses that damping [2011.01042]. The mode was originally proposed as an internal fluctuation of the superconducting order parameter and has since become a diagnostic of pairing competition, collective-mode hybridization, and unconventional superconducting symmetry [2109.11275].

## 1. Definition and microscopic origin

The defining condition for a Bardasis–Schrieffer mode is the coexistence of a dominant pairing instability and a distinct subdominant attractive interaction. In the formulation used for unconventional superconductors with multiple attractive channels, the condensate selects a ground-state symmetry \(l_0\), while fluctuations in a different channel \(l \neq l_0\) form a bound-state collective mode below the pair-breaking edge \(2\Delta\) [1907.10464]. In this sense the mode is a particle–particle exciton internal to the paired state rather than a quasiparticle continuum feature.

A standard example is an \(s\)-wave ground state with residual attraction in a \(d\)-wave channel. The \(d\)-wave channel is not strong enough to condense, but it binds a subgap pair fluctuation. In the notation of the non-equilibrium square-lattice model, this is the mode in the subdominant \(d_{x^2-y^2}\) sector when the equilibrium state is extended \(s\)-wave [1907.10464]. In the nematicity-coupled model, the same logic applies to \(s\)- and \(d\)-wave competition on a two-dimensional single-band square lattice, where the mode sits at \(\omega_{\mathrm{BS}}<2\Delta\) and is explicitly associated with the subdominant pairing symmetry [2011.01042].

The mode softens as the subdominant channel approaches degeneracy with the dominant one. In the pumped unconventional-superconductor study, the Bardasis–Schrieffer mode approaches the Higgs continuum far from the phase boundary, while near the transition between competing ground states it softens toward zero frequency [1907.10464]. The same softening appears in the \(s\)- versus \(d\)-competition phase diagram of superconductors with competing \(s\)- and \(d\)-wave interactions, where the one-band \(s\)-wave BS mode collapses at the \(s+id\) boundary and evolves into a mixed-symmetry collective mode inside the \(s+id\) phase [1507.07501].

The basic pole condition is a susceptibility equation in the subdominant channel. In the nematicity-coupled formulation, the single-channel limit recovers the familiar condition \(1-V_d \Pi_d(\omega)=0\) for an \(s\)-wave ground state [2011.01042]. In the current-enabled optical formalism, the collective frequencies are poles of the effective interaction kernel,
\[
\det V_\text{eff}^{-1} \overset{!}{=}0,
\]
with the single-band \(s\)-phase Bardasis–Schrieffer mode satisfying
\[
\frac{1}{V_d}-L_d+\frac{(i\Omega)^2}{2}I_2=0
\]
[2504.06642]. In Raman-based multiband formulations, the same physics appears as zeros of denominators such as \([2/\tilde U-\Pi_{22}(\Omega)]\) in the subleading pairing sector [1611.04541].

## 2. Field-theoretic and pseudospin formulations

Several complementary microscopic descriptions are used in the literature. A common approach starts from a BCS or mean-field Hamiltonian with separable pairing interactions in multiple angular-momentum channels, then derives collective modes from Gaussian fluctuations around the saddle point [2107.02834]. In the time-domain studies of unconventional superconductors, the problem is recast in terms of Anderson pseudospins,
\[
\mathbf{s}_{\mathbf{k}}=\frac{1}{2}(c^\dagger_{\mathbf{k}\uparrow},c_{-\mathbf{k}\downarrow})\,\boldsymbol{\sigma}\,\begin{pmatrix}c_{\mathbf{k}\uparrow}\\ c^\dagger_{-\mathbf{k}\downarrow}\end{pmatrix},
\]
whose dynamics obey
\[
\frac{d}{dt}\langle \mathbf{s}_{\mathbf{k}}\rangle=\mathbf{B}_{\mathbf{k}}(t)\times \langle \mathbf{s}_{\mathbf{k}}\rangle
\]
[1907.10464].

In the nematicity-coupled model, the mean-field Hamiltonian is
\[
H = \sum_{k,\sigma} \big(\xi_k + \Delta_n \gamma_k^n\big)\, c_{k\sigma}^\dagger c_{k\sigma}
    + \sum_k \Big[\Delta_k\, c_{k\uparrow}^\dagger c_{-k\downarrow}^\dagger + \text{h.c.}\Big],
\]
with
\[
\Delta_k = \Delta_s \gamma_k^s + e^{i\alpha}\, \Delta_d \gamma_k^d,
\]
and the choice \(\gamma_k^s=1\), \(\gamma_k^d=\gamma_k^n=\sqrt{2}\cos(2\phi)\) [2011.01042]. Linearizing the pseudospin equations around equilibrium gives
\[
\delta \boldsymbol{\Delta}(\omega)
=
\big(\mathbf{1}-\boldsymbol{\chi}(\omega)\big)^{-1}
\big(2\alpha\,\mathcal{F}[A^2](\omega)\big)
\sum_k V G_k^T M_k \hat{e}_z,
\]
and collective modes satisfy
\[
\det[\mathbf{1}-\boldsymbol{\chi}(\omega)]=0
\]
[2011.01042].

In the Raman-response framework for multicomponent superconductors, one starts from a general BdG Hamiltonian with separable pairing interaction,
\[
H_{\mathrm{BdG}}(\mathbf{k})=
\begin{bmatrix}
H_0(\mathbf{k}) & \Delta(\mathbf{k})\\
\Delta^\dagger(\mathbf{k}) & -H_0^{\mathrm{T}}(-\mathbf{k})
\end{bmatrix},
\quad
\Delta(\mathbf{k})=\sum_{\alpha=1}^{N_{\Delta}}\Delta^{\mathrm{MF}}_\alpha \hat{\Delta}_\alpha(\mathbf{k}),
\]
and derives a gauge-invariant Raman susceptibility in which collective modes appear as poles of an effective fluctuation kernel \(U_{\mathrm{eff},\phi}\) [2602.05607]. In that treatment the Bardasis–Schrieffer mode is one member of a broader class of relative pairing-channel modes, distinct from Leggett and clapping modes by its intraband, subdominant-channel character [2602.05607].

## 3. Energetics, softening, and relation to competing phases

A central property of the Bardasis–Schrieffer mode is that it lies below the pair-breaking edge of the dominant condensate. In the pumped square-lattice model this is expressed directly as \(\omega_{\mathrm{BS}}<\omega_H\), where \(\omega_H=2\Delta_{l_0,\max}\) is the Higgs energy for the dominant gap [1907.10464]. In rotationally symmetric \(s\)-wave states, the mode is undamped because it sits below the continuum. In nodal \(d\)-wave states, however, the analogous subdominant-channel oscillation is strongly damped by nodal quasiparticles [1907.10464].

The frequency tracks the competition between pairing channels. In the one-band \(s\)-wave phase of the \(s\)- versus \(d\)-competition model, the phase-sector equation for the BS mode is
\[
-\left(\frac{\Omega}{2\Delta_0^s}\right)^2 I_2(\Omega)=\frac{2}{u_d}-\frac{2}{u_s},
\]
and the solution softens as \(u_d\to u_s\), marking the approach to the \(s+id\) phase boundary [1507.07501]. In the mixed \(s+id\) regime the mode no longer remains a pure \(d\)-channel phase oscillation, but becomes a “mixed-symmetry Bardasis–Schrieffer mode” determined by a coupled equation involving \(I_0\), \(I_2\), and \(I_4\) [1507.07501].

The multiband case adds further structure. In the three-pocket model for hole-doped Fe-based superconductors, an \(s\)-wave ground state supports one well-defined Bardasis–Schrieffer mode below the lowest gap edge and a second damped BS-like mode between the two gap energies [1507.07501]. These coexist with a damped Leggett mode between the two gap scales [1507.07501]. This separation of sharp subgap and damped intergap resonances is important experimentally because a subgap peak is not automatically unique to one collective degree of freedom.

A closely related issue is the distinction from other relative modes. In the Raman theory of multicomponent superconductors, the Bardasis–Schrieffer mode is explicitly contrasted with the Leggett mode, which is interband and corresponds to relative-phase oscillation between condensates on different bands, and with clapping modes, which arise in chiral states as fluctuations of Cooper-pair angular momentum [2602.05607]. This distinction becomes blurred in systems with strong symmetry mixing. In FeSe with electronic nematicity, the observed low-energy mode is attributed to a fluctuation between an \(s+d\)-wave-like ground state and a subleading pairing channel; it is described as corresponding to the Bardasis–Schrieffer mode but also resembling an intraband Leggett mode [2507.14466]. This suggests that in lowered-symmetry environments the modal taxonomy can become basis dependent even when the underlying subdominant-channel origin remains intact.

## 4. Hybridization with other collective modes and with symmetry breaking

The Bardasis–Schrieffer mode need not remain isolated. In the nematicity-coupled model, the most important result is that in a rotationally symmetric \(s\)-wave ground state the \(d\)-wave Bardasis–Schrieffer fluctuation hybridizes with the nematic collective mode, producing a single in-gap excitation of mixed nematic–BS character [2011.01042]. The relevant \(2\times2\) kernel couples the imaginary \(d\)-wave component \(\Delta_d''\) and the nematic amplitude \(\Delta_n\),
\[
\chi_{\mathrm{N\text{-}BS}}(\omega)=
\begin{pmatrix}
-2V_d(I_1+I_2) & iV_d\,\omega\,I_2/\Delta \\
-2iV_n\,\omega\,I_2/\Delta & -4V_n I_2
\end{pmatrix},
\]
and only one in-gap solution exists for \(0<x<1\), with \(x=\omega/(2\Delta)\) [2011.01042].

This mixing is symmetry controlled. For \(B_{1g}\) nematicity, where \(\gamma_k^n=\sqrt{2}\cos 2\phi\), the nematic and \(d\)-wave pairing sectors hybridize because they transform identically [2011.01042]. If the nematicity instead has \(B_{2g}\) form factor, \(\gamma_k^n\propto \sin 2\phi\), the off-diagonal couplings vanish and BS and nematic modes remain distinct [2011.01042]. The mixed mode softens both at the nematic critical point \(V_n\to V_{n,c}\) and near \(V_d\to V_s\), where \(s+id\) degeneracy is approached [2011.01042].

In bilayer and locally non-centrosymmetric settings, the subdominant-channel fluctuation can take a phase-like form and couple strongly to Coulomb or electromagnetic degrees of freedom. In a superconducting bilayer with strong spin-orbit coupling, the competing odd-parity pairing instability gives rise to a Bardasis–Schrieffer-like phase mode inside the excitation gap when Coulomb interaction is neglected [2306.16611]. There the mode satisfies
\[
G^{-1}_{\mathrm{BS}}(\omega)=\frac{1}{g}+\chi_{\mathrm{BS}}(\omega),
\]
with the pole condition \((1/g)+\chi_{\mathrm{BS}}(\omega_{\mathrm{BS}})=0\) [2306.16611]. Once long-range Coulomb interaction is included, however, the mode is converted into an antisymmetric plasmon and is generally pushed into the quasiparticle continuum [2306.16611]. This is unlike conventional single-layer BS modes, which are not primarily longitudinal phase oscillations.

Hybridization also appears in optical and cavity contexts. In cavity proposals, a supercurrent produces a linear light–BS coupling, allowing one photon polarization to hybridize with the BS mode and form upper and lower BS polaritons [1807.07501]. In graphene–superconductor heterostructures, the graphene plasmon hybridizes with the Bardasis–Schrieffer mode of a two-dimensional superconductor, producing clear anticrossings in the near-field reflection coefficient \( \mathrm{Im}\,r_p(q,\omega)\) [2109.11275]. These examples show that the BS mode can enter composite bosonic sectors well beyond isolated superconducting fluctuation spectra.

## 5. Spectroscopic signatures: Raman, pump–probe, THG, and linear optics

Raman spectroscopy remains the canonical probe because the Bardasis–Schrieffer mode is symmetry selective and typically appears below the pair-breaking threshold in the channel corresponding to the subleading pairing symmetry. In the multiband Raman theory, the weights of multiple BS modes are set by overlaps between the Raman vertex and the subleading gap harmonics, allowing several orthogonal subleading eigenfunctions within the same irreducible representation to produce several distinct Bardasis–Schrieffer resonances [1611.04541]. In centrosymmetric crystals, Raman activity requires
\[
\Gamma_{\hat{\Delta}^{\dagger}_\alpha}\otimes\Gamma_{\gamma}\otimes\Gamma_{\hat{\Delta}_\beta}\ni A_{\mathrm{tot}},
\]
which gives the group-theoretical selection rule for Raman-active collective modes, including BS excitations [2602.05607].

For iron-based superconductors, the \(B_{1g}\) Raman channel is especially important. In hole-doped BaFe\(_2\)As\(_2\), a general multiband Raman treatment identifies multiple Bardasis–Schrieffer-type modes and relates their spectral weights to subleading gap structures [1611.04541]. In systems with nearby \(d\)-wave and nematic fluctuations, a single \(B_{1g}\) in-gap feature can instead be a mixed nematic–BS mode rather than a pure BS excitation, which makes assigning a Raman peak solely to subdominant pairing unreliable except very close to phase boundaries where mixing vanishes [2011.01042].

Ultrafast pump–probe and THz time-domain methods reveal the mode in non-equilibrium dynamics. In the pumped unconventional-superconductor model, a pump with polarization along the crystal axis induces a finite \(d_{x^2-y^2}\) component of the pseudomagnetic field and excites the subdominant channel, while diagonal polarization suppresses the BS excitation because \(\xi_{\mathbf{k},A}^{d}=0\) [1907.10464]. In an \(s\)-wave ground state, \(\Delta_d(t)\) oscillates at \(\omega_{\mathrm{BS}}\) whereas \(|\Delta_s(t)|\) exhibits Higgs oscillations at \(\omega_H\) [1907.10464]. The quasiparticle distribution function \(s_{\mathbf{k}}^z\) can carry especially strong BS signatures, suggesting visibility in time-resolved ARPES [1907.10464].

Third-harmonic generation provides another nonlinear route. In the THG analysis of unconventional superconductors, the Bardasis–Schrieffer mode remains unaffected by long-range Coulomb interaction, while the phase mode is pushed to the plasmon frequency [2107.02834]. The BS mode contributes a sizable signal to third-harmonic currents and has a characteristic polarization dependence. Unlike the Higgs mode, it can generate a current perpendicular to the applied field, with angular structure
\[
j_{3,\perp}(\phi)\propto \sin(4\phi),
\]
and a resonance when the effective drive frequency satisfies \(2\Omega=\Omega_{\mathrm{BS}}\) [2107.02834]. This provides a clean way to distinguish BS contributions from charge-density-fluctuation backgrounds near \(2\Delta\).

Linear optical access is usually forbidden in centrosymmetric single-band systems, but several mechanisms can lift that restriction. A supercurrent can enable linear optical conductivity from BS modes in superconductors with competing pairing channels [2504.06642]. Rashba spin-orbit coupling can render BS modes optically active even without a supercurrent in interband \(p\)-wave systems [2504.06642]. Trigonal warping in single-valley superconductors can also make BS and clapping modes visible in longitudinal and Hall optical response [2508.09387]. These results suggest that optical darkness is not intrinsic to the mode itself, but rather to a restrictive set of symmetry conditions.

## 6. Extensions beyond conventional superconductors and current directions

The Bardasis–Schrieffer concept extends beyond conventional superconductors. In excitonic insulators, the same basic mechanism produces subgap bound states in noncondensed internal channels. For an \(s\)-wave excitonic insulator, a \(p\)-wave Bardasis–Schrieffer mode exists below the gap energy and is optically active at \(q=0\), because a uniform electric field acts directly on the internal structure of the electron–hole pair [2003.02997]. The mode satisfies a bound-state condition
\[
1+g_p\,\Pi_p(\omega,\mathbf{q})=0,
\]
and can hybridize strongly with photons to form Bardasis–Schrieffer polaritons observable in far-field and near-field optics [2003.02997].

A first-principles extension has also been developed for excitonic insulators using a quasiparticle Bethe–Salpeter equation. In biased WSe\(_2\)–MoSe\(_2\) bilayers, this approach reveals a gapless phase mode, subgap Bardasis–Schrieffer modes, and above-gap scattering states, with a particularly strong optically bright \(2p\) BS mode dominating the subgap optical conductivity [2512.20969]. This suggests that the BS concept is not limited to separable weak-coupling models, but can be carried into ab initio descriptions of symmetry-broken electronic phases.

Recent superconducting experiments continue to sharpen the mode’s materials relevance. In FeSe, terahertz nonlinear spectroscopy detects a collective mode substantially below the superconducting gap energy, distinct from the amplitude Higgs mode, and theory attributes it to a fluctuation between an \(s+d\)-wave-like ground state and a subleading pairing channel activated by electronic nematicity [2507.14466]. In kagome superconductors, scanning tunneling spectroscopy reveals a distinct collective mode below \(2\Delta\) whose interpretation remains between a Leggett mode and a Bardasis–Schrieffer mode due to a subleading superconducting component [2403.09395]. These observations indicate that the BS framework is increasingly relevant in multicomponent, symmetry-lowered, and topologically nontrivial superconductors.

A recurring theme in current work is that the Bardasis–Schrieffer mode is best understood not as a single universal spectral line, but as a class of subdominant-channel collective excitations whose detailed manifestation depends on nodal structure, multiband composition, symmetry lowering, and coupling to additional bosonic sectors. This suggests that identifying a BS mode experimentally requires combined information from symmetry selection rules, location relative to \(2\Delta\), polarization dependence, and evolution with tuning parameters such as doping, pressure, nematicity, or supercurrent [2011.01042].

Source: https://www.emergentmind.com/topics/bardasis-schrieffer-mode