Measuring intrinsic relaxation rates in superconductors using nonlinear response

This paper proposes a method to experimentally extract intrinsic relaxation rates (including Higgs mode decay, quasiparticle redistribution, and dephasing) in both $s$- and $d$-wave superconductors by analyzing their nonlinear terahertz optical response, such as third harmonic generation and time-dependent gap dynamics, under various polarization conditions and damping scenarios.

The Gist

The Big Picture: Listening to the "Heartbeat" of Superconductors

Imagine a superconductor as a massive, synchronized dance floor. In a normal metal, electrons are like a chaotic crowd bumping into each other. But in a superconductor, they pair up (like dance partners) and move in perfect unison. This collective dance creates a "gap" in their energy, which is what allows electricity to flow without resistance.

The authors of this paper want to know: How long does this perfect dance last before the partners get tired, lose their rhythm, or bump into the floor?

In physics terms, they are trying to measure the intrinsic relaxation rates—how fast the system relaxes back to calm after being disturbed. To do this, they use a technique called nonlinear optical spectroscopy, which is like hitting the dance floor with a specific type of light pulse to see how the dancers react.

---

The Tools of the Trade

1. The Anderson Pseudospin (The "Compass" Analogy)

To understand the math, the authors use a concept called "Anderson pseudospins."

• The Analogy: Imagine every electron pair has a tiny compass needle attached to it.
  • If the needle points down, the pair is "asleep" (no electron).
  • If it points up, the pair is "awake" (two electrons).
  • If it points sideways, the pair is in a superposition (a mix of both).
• The Goal: When you shine light on the superconductor, you twist these compass needles. The way they wobble back to their resting position tells you everything about the material's health.

2. The Light Pulse (The "Whack-a-Mole" Game)

The researchers use a laser pulse (specifically in the Terahertz range) to "whack" the compass needles.

• Linear vs. Nonlinear: If you hit a drum lightly, it makes a simple sound (linear). If you hit it hard or in a specific way, it makes complex harmonics (nonlinear).
• The Trick: The "Higgs mode" (a specific vibration of the superconductor) is invisible to normal light. It's like a ghost that doesn't reflect a flashlight. But if you hit it with a strong, nonlinear pulse, the ghost leaves a shadow. The authors look for this shadow in the form of a Third Harmonic Generation (a signal that vibrates three times faster than the light you sent in).

---

The Two Main Characters: s-wave and d-wave

The paper studies two types of superconductors, which behave like different types of dance floors.

1. s-wave Superconductors (The Round Dance Floor)

• The Shape: Imagine a perfectly circular dance floor. The "dance partners" are evenly spaced all around.
• The Result: When you hit this floor, the compass needles wobble and slowly lose energy.
• The Decay: In a perfect, clean world, the wobble doesn't stop abruptly; it fades away slowly like a bell ringing in a canyon (a "power-law" decay). It takes a long time to die out.
• The Measurement: By measuring how fast the wobble fades, the scientists can calculate $T_2$ (how fast the rhythm is lost) and $T_1$ (how fast the energy is lost).

2. d-wave Superconductors (The Square Dance Floor with Corners)

• The Shape: Imagine a square dance floor. The dancers are strong in the middle of the sides (the "antinodes") but completely stop dancing at the corners (the "nodes").
• The Result: This is trickier. Because the dancers stop at the corners, the "wobble" fades away much faster than in the round floor.
• The Polarization Knob: This is where the paper gets clever. The authors realized that by rotating the polarization of the light (changing the angle of the "whack"), they can target specific parts of the square floor.
  • Shine light one way? You excite the sides.
  • Shine light at a 45-degree angle? You excite the corners.
  • This allows them to isolate different "modes" of the dance and measure their relaxation rates separately.

---

The "Intrinsic" Mystery: Why do they slow down?

In a perfect vacuum, these compass needles would wobble forever. But in real life, they slow down. The paper asks: Why?

1. Uniform Damping: Imagine the whole dance floor is covered in sticky mud. Everyone slows down at the same rate.
2. Energy-Dependent Damping: Imagine the dancers near the edge of the floor are on ice (fast), but those in the middle are on sand (slow). The paper explores how the "slowness" changes depending on how much energy the electron has.

The Key Finding:
If the damping is zero right at the "Fermi surface" (the edge of the dance floor), the signal doesn't fade exponentially (like a battery dying). Instead, it fades as a power law (like a slow leak). By measuring the shape of this leak, the scientists can figure out exactly how the damping depends on energy.

---

The "Aha!" Moment: How to Read the Data

The paper provides a recipe for experimentalists:

1. Hit the superconductor with a strong light pulse.
2. Watch the "Third Harmonic" signal (the echo).
3. Rotate the light polarization to isolate different parts of the material (especially for the square-shaped d-wave superconductors).
4. Analyze the decay:
   • If it decays exponentially, you've found the "intrinsic" relaxation time ($T_1$ and $T_2$).
   • If it decays as a power law, you've found that the damping is zero at the Fermi surface, and the exponent of that power law tells you how the damping behaves elsewhere.

Summary in a Nutshell

Think of the superconductor as a giant choir.

• The Light Pulse is the conductor shouting "Start!"
• The Higgs Mode is the choir's collective hum.
• The Relaxation Rates are how quickly the choir gets tired and stops humming.

This paper teaches us how to listen to the choir's hum, even when it's very quiet, by using a special microphone (nonlinear light) and by asking the choir to sing from different angles (polarization control). By doing this, we can figure out exactly why they are getting tired—whether it's because the room is hot, the singers are sick, or the music is too hard. This helps us understand the fundamental limits of superconductivity.

Technical Summary

1. Problem Statement

The paper addresses the challenge of measuring intrinsic relaxation rates in superconductors, specifically the Higgs mode decay rate, the quasiparticle redistribution rate ($1/T_1$), and the quasiparticle dephasing rate ($1/T_2$).

• Context: In the pure BCS model (clean limit), these intrinsic rates are zero. Observed decay in experiments is often attributed to "Landau damping" (dephasing due to the inhomogeneity of the pseudomagnetic field across the Fermi surface), which produces a power-law decay ($t^{-1/2}$ or $t^{-1}$) rather than exponential decay.
• Gap: Distinguishing between this intrinsic inhomogeneous dephasing and true phenomenological damping (energy loss to the environment, e.g., phonons) is difficult. Furthermore, existing methods often struggle to probe specific symmetry channels (irreducible representations) in unconventional superconductors.
• Goal: To establish a theoretical framework for extracting these intrinsic rates ($T_1, T_2$) and their energy dependence from nonlinear optical (terahertz) responses, specifically the time-dependent gap function and the third-harmonic generation (nonlinear current).

2. Methodology

The authors employ the Anderson pseudospin formalism to model the dynamics of superconductors under optical excitation.

• Formalism:

  • The superconducting state is mapped to a system of pseudospins $\vec{s}_k$ precessing around a pseudomagnetic field $\vec{b}_k$.
  • Light-Matter Coupling: Modeled via minimal coupling in the particle-hole basis. In the long-wavelength limit ($q=0$), the coupling is quadratic in the vector potential $\vec{A}(t)$, acting as a time-dependent shift in the chemical potential and a momentum-dependent twist in the pseudomagnetic field.
  • Damping: Phenomenological damping is introduced into the equations of motion. The authors consider three models:
    1. Homogeneous damping ($\gamma = \text{const}$).
    2. Inhomogeneous damping scaling with frequency ($\gamma \propto \omega_k^p$).
    3. Inhomogeneous damping scaling with energy distance from the Fermi level ($\gamma \propto |\epsilon_k|^p$).
  • Relaxation Times: The dynamics are decomposed into longitudinal ($s_\parallel$, related to $T_1$) and transverse ($s_\perp$, related to $T_2$) components relative to the pseudomagnetic field.

• Simulation Setup:

  • Lattices: 2D square lattice with nearest-neighbor ($J$) and further-neighbor hoppings ($J', J''$) to model both s-wave (isotropic gap) and d-wave (anisotropic gap, $\Delta_k \propto \cos k_x - \cos k_y$) superconductors.
  • Pump-Probe Scheme: A strong pump pulse (Gaussian vector potential) quenches the gap and induces pseudospin precession. A weak probe pulse measures the nonlinear current at a delay time $\tau$.
  • Polarization Control: The authors systematically vary the polarization of the pump ($\alpha$) and the direction of current measurement to selectively excite and detect different irreducible representations (irreps) of the $D_{4h}$ point group ($A_{1g}, B_{1g}, B_{2g}$).

3. Key Contributions

1. Extension of Nonlinear Response Theory: The work extends previous studies (e.g., Cui et al.) by including third-harmonic generation (nonlinear current) alongside gap dynamics, allowing for the probing of irreps that do not couple to the gap itself.
2. Energy-Dependent Damping: The paper explicitly models damping rates as functions of energy/momentum, moving beyond constant phenomenological damping.
3. Polarization as a Control Knob: It demonstrates that by rotating light polarization (e.g., from $0^\circ$ to $45^\circ$), one can selectively excite different symmetry modes ($B_{1g}$ vs. $B_{2g}$) even in the $q=0$ limit, enabling the isolation of specific relaxation channels.
4. Decoupling Intrinsic vs. Phenomenological Decay: The authors provide analytical and numerical criteria to distinguish between power-law decay (caused by intrinsic inhomogeneity) and exponential decay (caused by true relaxation).

4. Key Results

A. s-Wave Superconductors

• Collisionless Limit: In the absence of phenomenological damping, the Higgs mode and nonlinear current decay as $t^{-1/2}$ due to Landau damping (dephasing of pseudospins with different $\omega_k$).
• With Damping:
  • Homogeneous Damping ($\gamma = \text{const}$): The decay becomes exponential with a $t^{-1/2}$ envelope: $e^{-t/T_2}/\sqrt{t}$ for oscillations and $e^{-t/T_1}/\sqrt{t}$ for recovery. $T_1$ and $T_2$ can be extracted from the exponential envelope.
  • Energy-Dependent Damping ($\gamma \propto |\epsilon|^p$): If damping vanishes at the Fermi surface, the decay reverts to a pure power law ($t^{-1/p}$ for recovery). The exponent $p$ reveals the energy dependence of the damping mechanism.
• Extraction: By correcting for the intrinsic $1/\sqrt{t}$ decay, one can extract the intrinsic $T_1$ (from gap recovery) and $T_2$ (from oscillation decay) at the Fermi surface.

B. d-Wave Superconductors

• Collisionless Limit: The decay is faster than in s-wave, following $t^{-1}$ for the nonlinear current and $t^{-2.5}$ (numerically) for the Higgs amplitude. This is due to the additional dephasing from the angular dependence of the gap (nodal vs. antinodal points).
• Symmetry Selectivity:
  • x-polarized pump: Primarily excites the $B_{1g}$ mode (gap symmetry). The $B_{2g}$ mode is not excited.
  • x'-polarized pump ($45^\circ$): Excites both $A_{1g}$ and $B_{2g}$ modes.
  • Nonlinear Current: Allows probing of $B_{2g}$ via cross-polarized measurements ($J_{xy}$), which is impossible via the gap function alone.
• Relaxation Dynamics:
  • Oscillations: Dominated by antinodal regions. If damping exists at antinodes, oscillations decay exponentially ($e^{-t/T_{2, \text{antinode}}}$).
  • Recovery: Dominated by nodal regions. If damping vanishes at nodes (e.g., $\gamma \propto \omega^p$), the recovery follows a power law ($t^{-4/p}$).
  • Conclusion: One can separately extract the damping characteristics at the nodes (via recovery) and antinodes (via oscillation decay).

5. Significance

• Experimental Roadmap: The paper provides a concrete protocol for experimentalists using terahertz pump-probe spectroscopy to measure intrinsic quasiparticle lifetimes ($T_1, T_2$) in both conventional and unconventional superconductors.
• Symmetry Resolution: It highlights polarization control as a critical tool for disentangling complex relaxation mechanisms in materials with anisotropic gaps (like cuprates), where different parts of the Fermi surface may relax at different rates.
• Fundamental Insight: By distinguishing between "pure" dephasing (intrinsic to the band structure) and energy relaxation (coupling to a bath), the work clarifies the interpretation of nonlinear optical signals in quantum materials, offering a path to study fractionalization and spin liquids in future applications.

In summary, the authors demonstrate that nonlinear optical spectroscopy, combined with polarization control and careful analysis of decay exponents, serves as a powerful microscope for the microscopic relaxation dynamics of superconductors.