Quaternionic superconductivity with a single-field Bogoliubov-de Gennes--Ginzburg-Landau framework and charge-4e couplings

This paper establishes a compact quaternionic field theory framework for spinful superconductivity that unifies singlet and triplet pairing into a single field, enabling the analytical derivation and numerical verification of charge-4e signatures, including vestigial order criteria, fractional flux vortices, and dominant second-harmonic Josephson effects.

The Gist

Imagine you are trying to describe a complex dance performed by a pair of electrons. In the world of superconductivity, these electrons usually pair up to form a "Cooper pair," which acts like a single, super-fluid particle that can flow without resistance.

For decades, physicists have used a complicated set of mathematical tools (like complex numbers and 2x2 matrices) to describe these dances, especially when the electrons have "spin" (a quantum property like a tiny internal compass). It's like trying to describe a ballet using a spreadsheet: accurate, but messy and hard to read.

This paper introduces a new, much cleaner way to describe this dance using Quaternions.

1. The New Language: Quaternions as a "Super-Compass"

Think of a standard number as a point on a line. A "complex number" is a point on a flat sheet (2D). A Quaternion is a point in 4D space.

The authors realized that the "spin" of an electron pair (whether they are spinning together in a specific way or opposite ways) fits perfectly into this 4D structure.

• The Old Way: You had to write down four different equations to track the different spin combinations.
• The New Way: You just write down one quaternion number, let's call it $q$. This single number holds all the information about the spin-singlet (opposite spins) and spin-triplet (aligned spins) states simultaneously.

Analogy: Imagine you are describing a weather system. The old way was to write separate reports for wind speed, wind direction, humidity, and temperature. The new way is to use a single "Weather Cube" that contains all four data points in one neat package. It makes the math compact and easier to manipulate.

2. The "Ghost" Dance: Charge-4e Superconductivity

Usually, superconductors carry a charge of 2e (because a Cooper pair is two electrons). But the paper explores a rare, exotic state where four electrons team up to form a "quartet" (Charge-4e).

Think of this like a dance troupe:

• Standard Superconductor: Two dancers hold hands and glide across the floor.
• Charge-4e Superconductor: Two pairs of dancers link up to form a square formation, gliding together as a single unit.

The authors show how to predict when this "quartet" dance happens. They found that even if the standard "pair" dance disappears, the "quartet" dance can survive as a "vestigial" order (a leftover echo of the original dance).

3. The Magic of the "Half-Step"

One of the coolest predictions of this theory is about magnetic flux (how magnetic fields pass through the superconductor).

• In a normal superconductor, magnetic fields come in "packets" of a specific size (let's call it a "Full Step").
• In this new Charge-4e state, because the dancers are holding hands in groups of four, the magnetic field packets become half the size.

Analogy: Imagine a staircase where you usually take full steps. If you suddenly start taking half-steps, you can fit twice as many steps in the same amount of space. This "half-step" is a smoking gun that tells scientists, "Hey, we have a Charge-4e superconductor here!"

4. The "Double Beat" in the Music

The paper also looks at what happens when you push electricity through these materials (the Josephson effect).

• Normal Superconductor: If you apply a voltage, the current oscillates at a specific rhythm (frequency).
• Charge-4e Superconductor: The rhythm doubles. It's like a drummer who usually hits the snare once per beat suddenly hitting it twice as fast.

The authors created computer simulations to prove this. They built a virtual "lattice" (a grid of atoms) and watched the electrons dance. They confirmed that:

1. The math works perfectly with their new Quaternion language.
2. The "half-step" magnetic flux appears exactly as predicted.
3. The "double beat" electrical signal appears when the quartet formation is strong.

Why Does This Matter?

This isn't just about fancy math. It provides a universal translator for physicists.

• It connects the microscopic world (how individual electrons pair up) to the macroscopic world (how the material behaves in a device).
• It helps scientists design new materials that might work at higher temperatures or be used in quantum computers.
• It gives a clear checklist for experimentalists: "If you see a half-step in the magnetic field and a double-beat in the electrical rhythm, you've found a Charge-4e superconductor."

Summary

The authors took a messy, multi-page description of spinning electron pairs and compressed it into a single, elegant 4D number (a Quaternion). Using this new lens, they proved that electrons can form "quartets" (groups of four), which leads to exotic behaviors like magnetic fields splitting in half and electrical rhythms doubling. It's a new, simpler map for navigating the strange and wonderful world of quantum superconductivity.

Technical Summary

1. Problem Statement

Modern superconductors often exhibit complex behaviors involving strong spin-orbit coupling (SOC), time-reversal symmetry (TRS) with $T^2 = -1$, and multicomponent pairing that mixes singlet and triplet channels. Current theoretical frameworks face three primary limitations:

1. Complexity: Describing the spin content, Kramers constraints, and symmetry relations requires multiple matrices and gauge choices in standard complex $2 \times 2$ spin language.
2. Opacity: The placement of a system into Altland-Zirnbauer (AZ) topological classes is not immediately visible in the working variables.
3. Inconsistency: Topological invariants (e.g., $Z_2$ indices for helical Majorana states) are often computed using notation different from that used to write the Hamiltonian or Ginzburg-Landau (GL) theory.

Furthermore, there is a growing interest in charge-4e (quartet) superconductivity, where a uniform charge-4e state emerges as a vestigial order even when the standard charge-2e order vanishes. A unified, symmetry-transparent formalism is needed to treat spin, topology, and multi-fermion condensation simultaneously.

2. Methodology

The authors propose a Quaternionic Field Theory framework to recast spinful superconductivity.

• Quaternionic Encoding: The singlet-triplet gap is encoded into a single quaternion field $q(\mathbf{k}) = \psi(\mathbf{k})\mathbf{1} + \mathbf{d}(\mathbf{k}) \cdot (\mathbf{e}_x, \mathbf{e}_y, \mathbf{e}_z)$, where $\psi$ is the singlet and $\mathbf{d}$ is the triplet vector. This field lives in the complexified quaternion algebra (biquaternions).
• Compact Hamiltonian: The Bogoliubov-de Gennes (BdG) Hamiltonian is rewritten in a compact form:
  $$H_{BdG} = \xi_k \tau_z + \tau_+ q(\mathbf{k}) + \tau_- q^\dagger(\mathbf{k})$$
  where $\tau$ are Nambu matrices and $q^\dagger$ is the biquaternion adjoint. This formulation preserves time-reversal symmetry ($T: q(\mathbf{k}) \to q^\dagger(-\mathbf{k})$) and spin rotation covariance ($Sp(1)$) naturally.
• Symmetry and Topology: The framework maps the AZ classification (Classes C, CI, DIII) directly to algebraic constraints on $q(\mathbf{k})$. It connects naturally to quaternionic Berry connections and the FKMM invariant for topological diagnostics.
• Charge-4e Coupling: A quartet field $Q$ is defined as the symmetric product of two pairs: $Q \propto \text{Sc}(q^2)$. A coupled Ginzburg-Landau functional is constructed with covariant derivatives $(\nabla - 2ie\mathbf{A})q$ and $(\nabla - 4ie\mathbf{A})Q$, explicitly assigning charge $4e$ to the quartet field.
• Analytical and Numerical Verification:
  • Analytical: A one-loop evaluation of the fluctuation bubble $\Pi(0)$ is performed to derive a quantitative criterion for the vestigial charge-4e phase ($\mu_{\text{eff}} = \mu - g^2\Pi(0) < 0$).
  • Numerical: The authors simulate a 2D Class-DIII lattice model, a GL simulation of a charge-4e vortex, and a microscopic 1D two-orbital chain using Density Matrix Renormalization Group (DMRG).

3. Key Contributions

1. Unified Formalism: The paper establishes a single quaternionic variable $q$ that simultaneously encodes the gap structure, gauge charge, time-reversal constraints, and topological classification, simplifying the bookkeeping for spinful superconductors.
2. Charge-4e Phenomenology: It provides a compact derivation of charge-4e signatures, including:
   • Phase locking: $\phi_Q \approx 2\phi_q$.
   • Halved flux quantum: $\Phi_0/2 = h/4e$.
   • Doubled Josephson frequency: $f_{4e} = 2f_J$.
3. Quantitative Vestigial Criterion: The authors derive an explicit analytical expression for the effective quartet mass $\mu_{\text{eff}}$, determining the conditions under which a charge-4e condensate forms as a vestigial order above the $2e$ transition temperature.
4. Distinction of 4e Signatures: The work clarifies that a dominant second harmonic in the current-phase relation ($I_2 \gg I_1$) is not sufficient proof of charge-4e transport (as it can arise from high transparency or 0-$\pi$ transitions). Instead, it emphasizes the need for consistency between flux periodicity, frequency doubling, and a microscopic mechanism (like the coherent quartet channel).

4. Key Results

• Topological Validation: In a 2D Class-DIII lattice model, the $Z_2$ index computed from the occupied BdG eigenvectors (via the sewing-matrix/Pfaffian method) correctly predicts the existence of helical Majorana edge modes, confirming the framework's topological fidelity.
• Vortex Simulation: GL simulations of a pure quartet vortex demonstrate:
  • Flux quantization of exactly $h/4e$ (within $\sim 2\%$ numerical error).
  • A core size scaling as $\xi_Q \propto \sqrt{\eta/|\mu_{\text{eff}}|}$, consistent with analytical predictions.
• Microscopic Correlations: DMRG simulations of a 1D two-orbital chain show that in specific parameter regimes, the charge-4e correlation function decays more slowly than the charge-2e function, indicating a quartet-dominant regime.
• Josephson Effects: The study models a short junction where the second harmonic dominates ($I_2 \gg I_1$) due to a coherent 4e channel. This leads to:
  • Doubled AC Josephson emission ($2f_J$).
  • Shapiro steps at $V_n = n h f / 4e$ (even-only ladder).
  • SQUID interference patterns with a period of $\Phi_0/2$.

5. Significance

This paper provides a symmetry-faithful and computationally efficient route from microscopic pairing mechanisms to device-level observables in exotic superconductors.

• Theoretical Impact: It resolves the fragmentation between spin, topology, and GL theory by unifying them under quaternion algebra, making the derivation of topological invariants and symmetry constraints more transparent.
• Experimental Relevance: The framework offers clear, quantitative criteria to distinguish true charge-4e superconductivity from superficially similar phenomena (like $4\pi$ Josephson effects in topological platforms). This is crucial for interpreting recent experimental reports of 4e supercurrents in engineered Josephson structures and high-pressure hydrides.
• Future Applications: The formalism is readily extendable to multi-band materials, non-centrosymmetric crystals, and mesoscopic devices, providing a unified language for the next generation of topological and high-temperature superconductor research.