Impact of spin-orbit coupling and Zeeman interaction on the multiple Andreev reflections subharmonic gap structure in nanoscopic Josephson junctions

This paper theoretically demonstrates that spin-orbit coupling induces avoided crossings in the dispersion relation of spinful Josephson junctions, leading to pronounced multiple Andreev reflection features in the subharmonic gap structure whose visibility under an external magnetic field is governed by the spin polarization of the bands.

The Gist

The Big Picture: A Traffic Jam of Electrons

Imagine a superhighway where cars (electrons) are trying to drive through a tunnel (a nanowire) to get from one city (a superconductor) to another. Usually, in a superconductor, these cars travel in perfect pairs called "Cooper pairs."

When you apply a voltage (push), these cars try to cross the tunnel. Sometimes, they hit a wall and bounce back. In a special quantum world, when they bounce off the superconductor, they don't just bounce back as a car; they turn into a "hole" (a missing car) and bounce back again. This process is called Andreev reflection.

If the voltage is low, the cars get stuck in a loop, bouncing back and forth many times before they finally have enough energy to escape. This is called Multiple Andreev Reflection (MAR). The paper studies how this "bouncing game" changes when you add two specific ingredients: Spin-Orbit Coupling (a twist in the road) and Zeeman Interaction (a magnetic wind).

The Setup: The Experiment

The researchers built a theoretical model of a tiny wire connecting two superconducting cities.

• The Wire: It's a semiconductor nanowire.
• The Push: They apply a voltage to force electrons through.
• The Twist (Spin-Orbit Coupling): Imagine the road itself is twisted like a corkscrew. As a car drives, its "spin" (a quantum property like a tiny internal compass) gets rotated.
• The Wind (Zeeman Interaction): They apply a magnetic field, which acts like a strong wind pushing the cars with "North" spins one way and "South" spins the other way.

What Happens Without the Twist? (No Spin-Orbit Coupling)

First, the researchers looked at what happens if the road is straight (no twist), but the magnetic wind is blowing.

• The Result: As long as the wind isn't too strong, the traffic pattern (the conductance) looks exactly the same. The cars bounce back and forth in the same predictable rhythm.
• The Change: Once the wind gets very strong, it starts pushing the "North" and "South" lanes apart so much that the energy gap (the tunnel entrance) changes shape. Suddenly, the traffic jams shift to different locations. The "bouncing" happens at different voltages because the lanes have moved.

What Happens With the Twist? (With Spin-Orbit Coupling)

This is where the paper gets interesting. They added the "corkscrew road" (Spin-Orbit Coupling).

• The Magic Gap: Even with a weak magnetic wind, the twist in the road creates a new, stable gap right in the middle of the energy spectrum. It's like a new, permanent speed bump that stays there even when the wind changes.
• New Bouncing Patterns: Because of this new gap, the cars can now bounce in new ways. The researchers found that the "traffic jams" (peaks in conductance) split and move.
  • Some peaks move to lower voltages (easier to cross).
  • Some peaks move to higher voltages (harder to cross).
  • New peaks appear that didn't exist before.

Think of it like a game of pinball. Without the twist, the flippers (magnetic field) just push the ball slightly. With the twist, the whole table tilts, creating new bumpers and new paths for the ball to take. The ball hits new spots, creating new "dings" (peaks in the data).

The Secret Rule: Spin Selection

The most important discovery in the paper is a "rule of the road" regarding the cars' internal compasses (spin).

• The Rule: Cars can only successfully cross the tunnel and be absorbed if their compasses are pointing in compatible directions.
• The Blockage: If a car has a "North" spin and tries to enter a lane that only accepts "South" spins, it gets rejected. It's like trying to plug a square peg into a round hole.
• The Result: Even though the math says a certain energy level should create a traffic jam (a peak), if the spins don't match, the peak disappears. The researchers showed that the visibility of these peaks is entirely controlled by whether the spins align or clash.

Summary of Findings

1. Twist + Wind = New Roads: Combining the spin-orbit twist and the magnetic wind creates a complex map of energy gaps that didn't exist before.
2. Traffic Shifts: This new map causes the "bouncing" peaks in the electrical current to split, move, and change shape as you turn up the magnetic field.
3. Spin is the Gatekeeper: You can't just see these new peaks; they only appear if the electron spins are allowed to mix. If the spins are too different (orthogonal), the signal vanishes.

Why This Matters (According to the Paper)

The paper doesn't claim to build a new device or cure a disease. Instead, it provides a diagnostic tool.

• Scientists use these "bouncing peaks" (MAR) to measure the properties of tiny wires.
• This paper explains that if you see these peaks shifting or splitting in a specific way, it tells you exactly how the "twist" (spin-orbit coupling) and the "wind" (magnetic field) are interacting inside the wire.
• It helps researchers understand the "energy map" of these materials, which is crucial for studying exotic states of matter (like Majorana states) that scientists are trying to create in the lab.

In short: The paper maps out how twisting the road and blowing the wind changes the traffic flow of electrons, revealing that the "compass" of the electron is the ultimate traffic cop.

Technical Summary

Technical Summary: Impact of Spin-Orbit Coupling and Zeeman Interaction on Multiple Andreev Reflections

Problem Statement
Multiple Andreev reflections (MAR) in voltage-biased Josephson junctions generate a characteristic subharmonic gap structure in conductance spectra, typically appearing at voltages $V_b = 2\Delta/en$. This structure is widely used to characterize transport properties and estimate induced superconducting gaps in semiconductor-superconductor hybrids. However, in systems with strong spin-orbit interaction (SOI) and large g-factors—crucial for topological superconductivity and Majorana bound state research—the interplay between SOI, Zeeman interaction, and superconductivity creates a complex energy landscape. Previous experiments have shown that magnetic fields can induce additional subgap peaks and modify MAR spectra, yet a comprehensive theoretical understanding linking the specific evolution of the lead band structure to these conductance features, particularly regarding spin selection rules, remains to be fully established.

Methodology
The authors employ a direct, coherent, microscopic description of time-dependent current transport through a quasi-one-dimensional Josephson junction. The system consists of a semiconducting nanowire with proximity-induced superconductivity, modeled by a Hamiltonian that explicitly includes:

• Rashba spin-orbit coupling (strength $\alpha$).
• Zeeman interaction (strength $E_z$) from a perpendicular magnetic field.
• An opaque scattering region (normal segment) with tunable transparency $D$.

The superconducting phase difference is treated as time-dependent due to the voltage bias, $\phi(t) = (2e/\hbar)\int V_b dt$. The authors utilize the Tkwant package to perform time-evolution simulations, calculating the dissipative quasiparticle current without resorting to approximations like the short-junction limit or the Andreev limit. The model uses material parameters typical for InSb nanowires proximitized by high-field superconductors (e.g., Nb). Theoretical predictions are validated by comparing numerical time-dependent results with analytical solutions for the dispersion relations of the superconducting leads.

Key Contributions and Results

1. Decoupling of Spin Sectors without SOI:
   In the absence of spin-orbit coupling, the system Hamiltonian is block-diagonal, separating electron-hole pairs with opposite spins. The authors demonstrate that under these conditions, the Zeeman effect shifts the two spin sectors in opposite directions but does not alter the effective distance between the gap edges within each sector. Consequently, the MAR conductance peaks remain fixed at integer fractions of the superconducting gap ($2\Delta/en$) regardless of the magnetic field strength, until the bands cross the Fermi energy. Once bands cross the Fermi level, the relevant energy scale shifts from the full gap to the distance between the Fermi energy and the gap edge, causing the MAR peaks to shift to higher biases.

2. Spin Mixing and Gap Evolution with SOI:
   The introduction of Rashba SOI fundamentally alters the dispersion relation by opening avoided crossings. The authors analytically derive and numerically confirm that SOI induces a stable superconducting gap at zero energy (which persists up to large magnetic fields) and modifies higher-energy gaps. As the magnetic field increases, the energy differences between these various gap edges evolve non-monotonically: outer gaps widen while inner gaps near zero energy narrow.

3. Emergence of New MAR Features:
   The complex evolution of the band structure leads to a rich modification of the MAR conductance spectra. The authors show that new conductance peaks emerge corresponding to the energy differences between the newly opened gap edges and their subharmonics ($1/n$).

   • As the magnetic field increases, peaks associated with the narrowing zero-energy gap shift to lower voltages.
   • Peaks associated with the widening outer gaps shift to higher voltages.
   • At specific field strengths, first-order and second-order transitions can coalesce, merging peaks in the conductance trace.

4. Spin Selection Rules and Visibility:
   A critical finding is that the visibility of specific MAR features is governed by spin selection rules. The authors calculate the spin polarization ($\langle \sigma_z \tau_0 \rangle$) of the bands. They demonstrate that transitions between bands with opposite spin polarizations (e.g., the outermost gap edges at high fields) are suppressed due to spin orthogonality, rendering the corresponding conductance features invisible. Conversely, transitions between bands with mixed or aligned spin polarizations are allowed and manifest as distinct peaks. This explains why certain theoretically possible energy transitions do not appear in the conductance spectra.

5. Robustness to Transparency:
   The study confirms that while the amplitude and line shape of the conductance depend on the junction transparency ($D$), the qualitative evolution of the MAR subharmonic features driven by the band structure changes remains robust across low to intermediate transparencies ($D=0.05$ to $0.1$). At high transparencies ($D \gtrsim 0.8$), the features broaden and may appear as dips, consistent with established theory.

Significance
The paper establishes a direct theoretical link between the microscopic band structure of proximitized semiconductors and the macroscopic conductance response in Josephson junctions. It clarifies that the complex MAR spectra observed in experiments with magnetic fields are not merely artifacts but are direct signatures of the interplay between SOI and Zeeman splitting. By identifying spin selection rules as the governing factor for peak visibility, the work provides a refined tool for interpreting transport data in hybrid semiconductor-superconductor devices. This understanding is essential for accurately characterizing the induced gap and transport properties in systems engineered for topological quantum computing, where distinguishing between trivial and topological subgap states is critical.