Thermal phase slips in superconducting films

Imagine a superconductor as a super-highway where cars (electrons) can zip along without ever hitting a pothole or getting stuck in traffic. This is the "dissipationless" state: perfect, frictionless flow.

However, even on this perfect highway, there are occasional glitches. Sometimes, a group of cars spontaneously decides to stop, turn around, or get confused, causing a temporary traffic jam. In the world of physics, this glitch is called a phase slip.

This paper is about understanding exactly how and why these glitches happen in thin films of superconducting material, especially when the traffic is heavy (high current).

Here is the breakdown in simple terms:

1. The Problem: The "Dark Count" Mystery

Scientists use these superconducting films to build super-sensitive cameras that can detect single photons (particles of light). They are so sensitive that they can tell when a single photon hits them.

But there's a problem: sometimes the camera "clicks" even when there is no light. This is called a dark count.

The Cause: It's not a broken camera; it's a thermal glitch. Heat causes the supercurrent to hiccup (a phase slip), creating a tiny voltage spike that the camera mistakes for a photon.
The Goal: To fix the camera, we need to understand exactly how these hiccups happen so we can predict and stop them.

2. The Old Theory vs. The New Reality

For decades, scientists had a great theory for how these hiccups happen in thin wires (like a single-lane road). They knew that as you push more current through the wire, the "energy barrier" (the hill the cars have to climb to cause a glitch) gets lower.

But modern cameras use wide strips (like a multi-lane highway), not just thin wires.

The Old Theory Failed: The old math assumed the glitch happened in a tiny, round spot. But in a wide strip, the glitch behaves differently. It's not a round spot; it's a weird, stretched-out shape.
The Challenge: Calculating the shape of this glitch in 2D is incredibly hard. It's like trying to predict the exact shape of a ripple in a pond when you throw a stone, but the pond is made of jelly and the stone is invisible. Previous attempts required supercomputers to guess the answer numerically.

3. The Breakthrough: The "Perfect Wave"

The authors of this paper found a way to solve this mathematically without needing a supercomputer. They discovered that the shape of the glitch follows a famous, elegant equation from fluid dynamics called the Boussinesq equation.

The Analogy:
Imagine you are watching a wave travel across a calm lake. Usually, waves crash and break. But in this specific mathematical world, there is a special type of wave called a soliton. A soliton is a perfect, solitary wave that travels forever without changing its shape or losing energy.

The authors realized that the "glitch" in the superconductor is exactly like this perfect soliton wave.

Because it's a "perfect wave" (mathematically integrable), they could write down the exact solution for what the glitch looks like.
They didn't have to guess; they solved it like a puzzle with a known, beautiful answer.

4. What Does the Glitch Look Like?

The paper reveals that this glitch is highly stretched out, like a long, thin oval.

Along the flow of current: It's short and compact.
Across the flow: It's very long and wide.

Think of it like a traffic jam on a highway.

In a narrow road (1D wire), the jam is a tight knot of cars.
In a wide highway (2D film), the jam stretches out for miles in the direction perpendicular to the traffic, creating a long, thin "valley" where the cars slow down.

5. The "Overheated" Zones

One of the coolest findings is that inside this glitch, the current actually gets faster than the maximum speed limit in two specific spots (on the sides of the glitch).

Analogy: Imagine a river flowing around a rock. Usually, the water slows down behind the rock. But in this superconductor glitch, the water speeds up so much on the sides that it technically breaks the "speed limit" of the river, yet it doesn't crash because the water level (the superconducting state) is changing in a very specific way to hold it together.

6. Why This Matters

The Math: They found that the energy needed to cause a glitch drops off in a very specific way as you get closer to the maximum current. It's not a simple curve; it follows a precise power law (specifically, the 3/4 power).
The Application: This helps engineers design better photon detectors. By knowing the exact shape and energy of these glitches, they can figure out how wide to make the strips and how much current to run to minimize those annoying "dark counts."
The Transition: They also predict a "phase transition." If you lower the current enough, the glitch changes shape entirely, turning from this smooth, stretched wave into a pair of swirling vortices (like two tiny tornadoes spinning in opposite directions).

Summary

This paper is a detective story where the detectives (the authors) solved a 20-year-old mystery about why superconducting cameras make false alarms. They realized the "culprit" (the thermal glitch) isn't a messy blob, but a perfectly shaped, mathematical wave (a soliton). By recognizing this wave, they unlocked the exact formula for how these glitches behave, paving the way for better, more reliable quantum sensors.

Here is a detailed technical summary of the paper "Thermal phase slips in superconducting films" by Skvortsov and Polkin.

1. Problem Statement

Superconducting nanowire single-photon detectors (SNSPDs) rely on the dissipationless flow of supercurrents. However, thermal fluctuations can trigger thermal phase slips, events where the superconducting order parameter is locally suppressed, creating a resistive state and causing "dark counts" (false detections) in the absence of photons.

The Gap: The theoretical framework for thermal phase slips in one-dimensional (1D) wires is well-established (Langer-Ambegaokar theory), where the activation barrier scales as $(1 - I/I_c)^{5/4}$ . However, modern SNSPDs often utilize two-dimensional (2D) superconducting strips with widths ( $w$ ) much larger than the coherence length ( $\xi$ ).
The Challenge: In 2D, the problem is significantly more complex due to:
1. The coexistence of different saddle-point configurations (topologically trivial vs. vortex-like).
2. Inhomogeneous current distributions.
3. The difficulty of finding an analytical solution for the saddle-point configuration (instanton) of the Ginzburg-Landau (GL) equations in 2D. Previous studies were limited to numerical approximations.

2. Methodology

The authors investigate an infinite 2D superconducting film in the vicinity of the critical current ( $I \to I_c$ ) and the critical temperature ( $T_c$ ), working within the Ginzburg-Landau regime.

Stream Function Formulation: To reduce the complexity of the coupled non-linear GL equations (which involve both the order parameter modulus $|\Delta|$ $∣Δ∣$ and phase $\phi$ $ϕ$ ), the authors introduce a scalar stream function $\psi(x,y)$ $ψ (x, y)$ .
- The supercurrent density is defined as $\mathbf{j} = (\partial_y \psi, -\partial_x \psi)$ .
- This approach assumes the flow is laminar and vortex-free (valid near $I_c$ ), allowing the phase $\phi$ to be uniquely reconstructed from $\psi$ .
Perturbation Theory: They expand the stream function around the uniform current state using a small parameter $\epsilon \approx \sqrt{1 - I/I_c}$ .
Derivation of the Effective Action: By substituting the perturbation series into the GL free energy functional, they derive an effective non-local action. Near the critical current, the anisotropy of the instanton ( $L_y \gg L_x$ ) allows them to retain only the dominant terms, resulting in a local functional involving second and third derivatives.
Connection to Integrable Systems: The saddle-point equation for the reduced free energy functional is identified as the elliptic form of the Boussinesq equation:
$u_{\bar{x}\bar{x}} + u_{\bar{y}\bar{y}} - u_{\bar{x}\bar{x}\bar{x}\bar{x}} + (u^2)_{\bar{x}\bar{x}} = 0$
where $u$ is related to the second derivative of the stream function.
Exact Solution: The authors solve this equation using Hirota's method, a technique for finding exact soliton solutions in integrable systems. They also perform numerical variational searches using Boussinesq-like ansatz functions to validate the analytical results and explore the regime further from $I_c$ .

3. Key Contributions

First Analytical Determination: This is the first analytical derivation of the phase slip rate and activation barrier for 2D superconducting films near the critical current.
Exact Integrability: The discovery that the saddle-point configuration is governed by the exactly integrable Boussinesq equation, allowing for an exact solution rather than numerical approximation.
Scaling Laws: Derivation of the specific scaling laws for the instanton size and activation energy in 2D, which differ fundamentally from the 1D case.
Topological Transition Hypothesis: Identification of a potential second-order topological transition between a vortex-free (trivial) instanton and a vortex-antivortex pair configuration at $I \approx 0.9 I_c$ .

4. Key Results

A. Instanton Geometry and Anisotropy

The saddle-point configuration (instanton) is highly anisotropic:

Longitudinal size (along current): $L_x \sim \xi (1 - I/I_c)^{-1/4}$ (similar to 1D).
Transverse size (perpendicular to current): $L_y \sim \xi (1 - I/I_c)^{-1/2}$ .
As $I \to I_c$ , $L_y \gg L_x \gg \xi$ . The instanton becomes a long, narrow strip perpendicular to the current flow.

B. Activation Energy Scaling

The activation barrier $\Delta F$ for thermal phase slips in 2D scales as:
$\Delta F_{2D} \approx c_2 \varepsilon_{cond} d \xi^2 (1 - I/I_c)^{3/4}$

Exponent: The scaling exponent is 3/4, distinct from the 1D Langer-Ambegaokar exponent of 5/4.
Physical Origin: The barrier scales as $\Delta F_{2D} \sim (L_y/\xi) \Delta F_{1D}$ . The large transverse dimension $L_y$ modifies the power law.
Constant: $c_2 = 28.55$ .

C. Current Distribution

The supercurrent density at the saddle point exhibits a unique "overheated" profile:

The current density dips at the center of the instanton.
To conserve total current, the density increases in the surrounding regions, exceeding the nominal critical current density ( $j_c$ ) in two lobes along the transverse axis. This is stabilized by a finite gradient of the order parameter.

D. Boundary Effects and Strip Width

Wide Strips ( $w \gg L_y$ ): The instanton forms at the edge of the strip. The activation energy is exactly half of the infinite film value ( $\Delta F_{2D}/2$ ).
Narrow Strips ( $w \ll L_y$ ): The system behaves as 1D, reverting to the Langer-Ambegaokar scaling.

E. Topological Transition

Numerical simulations suggest that as the current decreases from $I_c$ :

The vortex-free (Boussinesq) solution remains valid down to $I \approx 0.9 I_c$ .
At $I_{top} \approx 0.9 I_c$ , the order parameter at the center vanishes, marking a second-order transition to a topologically non-trivial state (vortex-antivortex pair). This contradicts previous hypotheses of a first-order transition where solutions might coexist.

5. Significance

SNSPD Performance: The results provide a rigorous theoretical basis for calculating dark count rates in wide-strip SNSPDs, which are critical for high-fidelity quantum photon detection. The derived scaling law ($3/4$) explains experimental observations that deviated from 1D predictions.
Theoretical Physics: The work bridges non-equilibrium superconductivity with integrable systems theory (Boussinesq equation), demonstrating that complex saddle-point problems in condensed matter can sometimes be solved exactly.
Design Guidelines: The findings clarify the transition between 1D and 2D behavior based on strip width and current bias, aiding in the optimization of superconducting detector geometries to minimize thermal noise.