Secondary Collective Excitations in Intermediate to Strong-Coupling Superconductors

This paper demonstrates that systematically derived energy-transfer-dependent effective electron-electron interactions in intermediate-to-strong-coupling superconductors give rise to secondary, long-lived collective phase and amplitude excitations below the quasiparticle continuum, which exhibit lattice-independent behavior and eigenoperator structures analogous to Hydrogen wave functions.

The Gist

The Big Picture: A Symphony of Electrons

Imagine a superconductor not as a wire, but as a massive, perfectly synchronized dance floor filled with electrons. In a normal metal, these electrons are like a chaotic crowd bumping into each other. But in a superconductor, they pair up and move in perfect unison, creating a "superfluid" that flows without friction.

This paper investigates what happens when you try to "jiggle" this synchronized dance floor. Specifically, the authors are looking for the specific "notes" or vibrations (called collective excitations) that the electron pairs make when disturbed.

The Two Main Dancers: The Higgs and the Phase

In the world of superconductors, there are two fundamental ways the dance can be disturbed:

1. The Higgs Mode (The Amplitude): Imagine the dancers holding hands. The "Higgs mode" is when they squeeze their grip tighter or loosen it. They are changing the strength of their connection.
2. The Phase Mode (The Rhythm): Imagine the dancers are all stepping in time. The "Phase mode" is when they all shift their steps slightly earlier or later. They aren't changing how hard they hold hands, but they are changing the timing of the dance.

In simple, weak interactions, scientists already knew about these two main dancers. The Higgs mode usually vibrates at a specific high energy (twice the energy gap), and the Phase mode vibrates at zero energy (like a perfect, silent rhythm).

The Discovery: The "Secondary" Dancers

The main discovery of this paper is that when the electrons interact strongly (like a very crowded, energetic dance floor), new, hidden dancers appear.

The authors found that if you crank up the interaction strength, secondary modes emerge. These are like backup dancers that were hiding in the crowd.

• They appear below the main energy limit where the electrons usually break apart.
• They are very long-lived (they don't fade away quickly).
• They appear in a very regular pattern. As the interaction gets stronger, these new modes pop up one by one, like bubbles rising in a pot of boiling water.

The paper shows that this happens regardless of the specific shape of the "dance floor" (whether it's a simple cubic, body-centered cubic, or face-centered cubic lattice). It seems to be a universal rule of strong superconductivity.

The "Hydrogen Atom" of Superconductivity

One of the most fascinating parts of the paper is how the authors figured out what these secondary dancers look like. They calculated the "wave functions" of these modes—the mathematical description of how the electrons move to create these vibrations.

They found a surprising pattern:

• The primary (first) mode looks like a smooth hill with no bumps.
• The second mode has two "nodes" (places where the vibration cancels out to zero, like a wave crossing the water line).
• The third mode has four nodes.
• The fourth mode has six nodes.

The Analogy: This is exactly like the Hydrogen atom in physics. In a hydrogen atom, electrons orbit the nucleus in specific shells. The first shell is a smooth sphere; the second has a node; the third has more. The authors found that these superconducting vibrations follow the exact same mathematical rules as electrons in a Hydrogen atom, but instead of orbiting a nucleus, they are "orbiting" in energy space. It's as if the superconductor has its own internal "quantum number" system for these vibrations.

Why Does This Happen?

The paper explains that this happens because the interaction between electrons isn't a simple, constant rule. It depends on how much energy the electrons exchange (a concept called "retardation").

Think of it like a conversation:

• Weak coupling: You shout a constant message to everyone. The reaction is simple.
• Strong coupling: You only talk to people who are within a certain distance and time window. This complex, time-delayed conversation creates a much richer set of possible reactions (the secondary modes).

The "W-Shape" Surprise

The authors also noticed something odd about the energy of the electrons themselves. Usually, the lowest energy point is right in the middle of the band. But in strong coupling, the energy landscape can twist into a "W" shape.

Imagine a valley that usually has one bottom. In these strong superconductors, the valley splits, creating two side valleys and a small hill in the middle. This means the electrons have multiple "favorite" places to sit, which is a direct result of the complex interactions described above.

Summary

In short, this paper reveals that superconductors are more complex than we thought. When the electrons interact strongly:

1. New vibrations appear: Hidden "secondary" modes emerge below the main energy limit.
2. They are universal: This happens on different types of crystal structures.
3. They have a pattern: These modes look mathematically identical to the energy levels of a Hydrogen atom, with increasing numbers of "nodes" or zeros.
4. They are stable: These new modes don't decay quickly; they are robust features of strong superconductivity.

The authors didn't propose a new device or a medical application. Instead, they provided a deeper theoretical map of how these quantum dances work, showing that even in a superconductor, there is a hidden, structured "universe" of vibrations waiting to be discovered.

Technical Summary

Technical Summary: Secondary Collective Excitations in Intermediate to Strong-Coupling Superconductors

Problem Statement
The paper investigates collective excitations in isotropic superconductors, specifically focusing on the behavior of amplitude (Higgs) and phase (Anderson-Bogoliubov) modes in the intermediate-to-strong-coupling regime. While the weak-coupling limit is well-understood within the Bardeen-Cooper-Schrieffer (BCS) framework—where the Higgs mode sits at the lower edge of the quasiparticle continuum ($\omega = 2\Delta$) and the phase mode is a Goldstone boson at $\omega=0$—the behavior of these modes under more realistic, momentum-dependent interactions derived from electron-phonon coupling remains less explored. Previous work on rotationally invariant continuum systems suggested the existence of "secondary" collective modes at higher energies, but it was unclear how these phenomena depend on the underlying lattice geometry and specific Fermi level configurations.

Methodology
The authors employ a pairing Hamiltonian derived from an electron-phonon interaction via a continuous unitary transformation (CUT). This results in an effective interaction potential $g(k, k')$ that inherently includes retardation effects, restricting interactions to electrons with energy differences smaller than the phonon energy scale. The study focuses on three-dimensional Bravais lattices: simple cubic (sc), body-centered cubic (bcc), and face-centered cubic (fcc).

To analyze collective excitations beyond the static mean-field approximation, the authors utilize the iterated equations of motion (iEoM) approach. This method involves:

1. Defining a suitable operator basis $\mathcal{B}$ containing pair creation/annihilation operators ($\hat{R}, \hat{I}$) and number operators ($\hat{N}$) discretized over single-particle energies.
2. Constructing the dynamical matrix $M$ and the norm matrix $N$ using a pseudo-scalar product defined by the expectation value of commutators with the full Hamiltonian.
3. Computing Fourier-transformed Green's functions to identify spectral functions $A(\omega)$, where collective modes appear as peaks (or $\delta$-functions represented as Gaussians).
4. Developing an algorithm to compute eigenoperators: specific linear combinations of the basis operators that excite a single collective mode. The coefficients of these combinations are termed "operator amplitudes."

Calculations are performed at zero temperature ($T=0$) for two Fermi level scenarios: particle-hole symmetric ($E_F = 0$) and asymmetric ($E_F = -0.5W$). The Coulomb interaction is explicitly excluded, though the authors note its primary effect would be a renormalization of the coupling strength.

Key Results

• Primary Modes and Detachment: In the weak-coupling limit, the primary Higgs mode is found at $\omega = 2\Delta_{true}$. As the coupling strength $g$ increases, the Higgs mode detaches from the lower edge of the quasiparticle continuum, moving to lower energies ($\omega < 2\Delta_{true}$). This detachment renders the mode infinitely long-lived (dissipationless) at zero temperature, a phenomenon driven by the non-trivial momentum dependence of the interaction.
• Emergence of Secondary Modes: Upon further increasing the coupling strength, additional secondary collective excitations emerge from the quasiparticle continuum. These appear in both the amplitude and phase channels.
  • On symmetric lattices (sc, bcc), secondary modes appear in an alternating fashion between the phase and amplitude spectral functions.
  • On the asymmetric fcc lattice, the lack of particle-hole symmetry causes the Higgs and phase modes to couple, resulting in signatures of both modes appearing in both spectral functions.
• Universality of Mode Emergence: The energy gap $\Delta_{true}$ required for a new mode to emerge from the continuum follows a linear relationship with the mode number. The authors find that an additional mode appears approximately every $\Delta_{true} \approx 0.64 \omega_D$ (where $\omega_D$ is the Debye frequency). This slope is found to be largely independent of the lattice type, Fermi level, and Coulomb interaction, suggesting a universal property of the secondary mode emergence.
• Quasiparticle Dispersion Anomalies: The study identifies that the minimum of the quasiparticle dispersion $E(\epsilon)$ does not necessarily lie at the chemical potential $\mu$. In cases where the order parameter $\Delta(\epsilon)$ decays faster than $|\epsilon - \mu|$ rises, side minima emerge, leading to a "W-shape" in the dispersion. This is particularly pronounced in the bcc lattice near specific coupling strengths.
• Structure of Eigenoperators: The computed eigenoperators reveal distinct nodal structures in their coefficients (operator amplitudes).
  • The primary modes exhibit no roots (excluding the node at $\epsilon = \mu$).
  • The $n$-th secondary mode exhibits $2(n-1)$ roots.
  • This behavior is described as reminiscent of the wave functions of bound states in one-dimensional quantum systems (e.g., the Hydrogen atom or harmonic oscillator), suggesting a "radial" quantum number character for these modes.
  • The amplitudes for the Higgs mode are shown to satisfy an analytical relation between pair creation and number operators, derived from the conservation of specific bilinear operators.

Significance and Claims
The paper claims to demonstrate that the existence of secondary collective excitations is a generic feature of superconductors with intermediate-to-strong coupling and momentum-dependent interactions, independent of specific lattice geometry or Fermi level details.

The authors emphasize that these secondary modes are distinct from Bardasis-Schrieffer modes, which correspond to different angular momentum channels. Instead, the identified modes possess the same symmetry as the underlying lattice (extended s-wave) but differ in their "radial" structure in energy space, characterized by the number of nodes in the operator amplitudes.

The study provides a theoretical framework for identifying these modes via eigenoperators and suggests that their detection could be possible through conventional spectroscopic methods (e.g., Raman or Terahertz spectroscopy) capable of probing the primary Higgs mode. The work establishes a foundation for future investigations into competing ordering phenomena and the dispersion of collective modes at finite momentum.