Cooper quartets and fractional vortices in frustrated Josephson junction dice arrays

This paper demonstrates through numerical simulations and tensor network techniques that frustrated Josephson junction dice arrays at one-third flux quantum frustration exhibit a superconductor-insulator transition characterized by half-vortex deconfinement and the emergence of a topologically protected 4e superconducting phase mediated by Cooper quartets.

The Gist

The Big Picture: A Dance of Superconducting Pairs

Imagine a superconductor as a crowded dance floor where people (electrons) usually pair up to dance in twos. These "Cooper pairs" (charge 2e) move in perfect sync, allowing electricity to flow without any resistance. This is standard superconductivity.

However, this paper explores a weird, exotic dance floor where, under specific conditions, the dancers don't just pair up; they form groups of four (Cooper quartets, charge 4e). The researchers are trying to figure out if they can build a machine that forces these groups of four to form and stay together.

The Stage: The "Dice" Lattice

To get these groups of four, the scientists are looking at a specific shape for their dance floor. Instead of a square grid (like a chessboard), they are using a Dice Lattice.

• The Shape: Imagine a honeycomb, but with extra connections. It looks like a grid of diamonds (rhombuses) packed together.
• The Setup: They are building this out of tiny islands of superconducting material connected by "Josephson junctions" (tiny bridges).
• The Frustration: They apply a magnetic field to the whole thing. But they don't apply just any amount of field. They apply a very specific amount: one-third of a magnetic "quantum" per diamond shape.

In physics, this is called "frustration." It's like trying to seat three people at a round table with only two chairs; they can't all be comfortable at once. This "frustration" forces the electrons to behave in unusual ways.

The Main Discovery: The "Quarter" Dance

When the researchers crunched the numbers and ran simulations on this frustrated Dice Lattice, they found something amazing at that specific "one-third" magnetic setting:

1. The Switch: The system stops acting like a normal superconductor (where pairs of two dance) and starts acting like a 4e superconductor (where groups of four dance).
2. The Evidence:
   • The Current: When they measured the electric current flowing through the system, the rhythm changed. Instead of a beat that repeated every time a single pair passed, the beat repeated only when four charges passed. It's like a drumbeat that only happens on the 4th count.
   • The Vortices (The Whirlpools): In a normal superconductor, magnetic fields create tiny whirlpools (vortices) that act like single units. In this "frustrated" state, the whirlpools split in half. These are called half-vortices.
   • The Tether: These half-vortices are tied together in pairs by invisible strings (domain walls). They can't run away alone; they are stuck in a group of two. Because they are stuck in pairs, the system effectively behaves as if the charge carriers are groups of four.

The "Half-Vortex" Analogy

Think of the magnetic field as a crowd of people trying to walk through a hallway.

• Normal Superconductor: The crowd moves in orderly lines. If someone gets stuck, the whole line stops.
• This Exotic State: The magnetic field is so "frustrated" that the crowd splits into two smaller, chaotic groups (half-vortices). These two groups are tied together by a rope. They can wiggle around, but they can't separate. Because they are bound together, the whole system moves as a single, larger unit (the quartet).

What About Disorder and Temperature?

Real-world experiments aren't perfect. The paper checked if this "group of four" dance would survive if the dance floor was slightly bumpy (disorder) or if the room got hot (temperature).

• Disorder: They found that even if the magnetic field isn't perfectly uniform or the bridges aren't identical, the "group of four" state is surprisingly robust. It survives the bumps.
• Temperature: As the system gets hotter, the "strings" tying the half-vortices together eventually snap. Once they snap, the groups of four fall apart, and the system returns to normal or stops conducting electricity entirely. The researchers calculated exactly when this "snap" happens (a phase transition).

The "Order by Disorder" Twist

The paper also looked at what happens at extremely cold temperatures (near absolute zero) when you add a tiny bit of electrical repulsion (charging energy).

• The Paradox: Usually, adding disorder (like repulsion) makes things messy. But here, the quantum rules say that the "messy" state of the groups of four is actually so crowded that the system gets confused.
• The Result: To resolve this confusion, the system suddenly snaps back into a rigid, ordered pattern (like a crystal) at ultra-low temperatures. It's as if the dancers, overwhelmed by the chaos of the group dance, decide to stand in a perfect, rigid line just to calm down. This is called "Order by Disorder."

Summary of Claims

The paper claims that:

1. Dice Lattices with a specific magnetic field (1/3 flux) are a perfect playground for creating 4e superconductivity (groups of four).
2. This state is characterized by half-vortices that are confined in pairs.
3. This state is stable against the imperfections found in real-world experiments.
4. At extremely low temperatures, quantum effects might force the system to abandon the "group of four" dance and return to a rigid, ordered state, but for a wide range of temperatures, the exotic "group of four" phase dominates.

The authors conclude that these setups are a promising way to build the hardware for future quantum computers that are protected by the laws of topology (meaning they are naturally resistant to errors), but they stop short of claiming this is ready for immediate commercial use. They are describing the physics of the phenomenon, not a finished product.

Technical Summary

Technical Summary: Cooper quartets and fractional vortices in frustrated Josephson junction dice arrays

Problem Statement
The paper investigates the emergence of superconductivity mediated by Cooper quartets (charge $4e$) rather than standard Cooper pairs ($2e$). While the condensation of $4e$ carriers is theoretically significant for understanding exotic superconductors and engineering topologically protected quantum memories, realizing this phase experimentally remains challenging. The authors focus on Josephson junction arrays (JJAs) arranged in a "dice lattice" geometry. Specifically, they analyze the low-energy physics of these arrays when subjected to magnetic frustration, defined as the magnetic flux $\Phi$ per rhombic plaquette normalized by the flux quantum $\Phi_0$ ($f = \Phi/\Phi_0$). Previous theoretical work suggested that at specific frustration values, particularly $f=1/3$, the system might support a phase characterized by fractional vortices (half-vortices) and $4e$ superconductivity, but a comprehensive numerical characterization of this regime in the presence of realistic disorder and thermal fluctuations was lacking.

Methodology
The authors employ a multi-faceted numerical approach to study the classical 2D XY model describing the dice lattice JJA:

1. Relaxation Dynamics: To determine critical currents ($I_c$), the authors simulate the relaxation dynamics of the system under adiabatic changes in external magnetic flux or current bias. This involves local updates of phases to mimic realistic phase slips and vortex motion, allowing for the estimation of switching currents in both phase-biased and current-biased scenarios.
2. Fourier Analysis: To identify the charge of the current carriers, the authors perform a harmonic decomposition of the supercurrents flowing between source and drain islands. They analyze the amplitudes of harmonics $I_n$, where $n$ corresponds to the charge $2ne$.
3. Monte Carlo Simulations: Used to calculate the helicity modulus and determine the Berezinskii-Kosterlitz-Thouless (BKT) transition temperature associated with the unbinding of half-vortices.
4. Infinite Tensor Networks: To characterize the system at finite temperatures and access the thermodynamic limit, the authors utilize 2D infinite tensor network methods (specifically Corner Transfer Matrix Renormalization Group, CTMRG). This allows for the calculation of correlation functions and entropy density, distinguishing between power-law and exponential decay behaviors.
5. Disorder Analysis: The robustness of the observed features is tested by introducing Gaussian disorder in both the magnetic flux per plaquette ($\Delta f$) and the Josephson energies ($\Delta E_J$).
6. Quantum Fluctuation Analysis: The authors extend the classical model to include charging energy ($E_C$) to assess the impact of quantum fluctuations and the "order-by-disorder" mechanism on the ground state degeneracy.

Key Results

• Critical Current Signatures: The numerical simulations reveal a distinct peak in the critical current at frustration $f=1/3$. This peak is asymmetric and robust against moderate disorder in both flux and Josephson energy. In contrast, the critical current at $f=1/2$ (full frustration) is suppressed due to Aharonov-Bohm caging.
• Dominance of Cooper Quartets: Fourier analysis of the supercurrents across the lattice shows that at $f=1/3$, the second harmonic ($I_2$, corresponding to $4e$ charge) dominates the first harmonic ($I_1$, corresponding to $2e$ charge) by a factor of approximately five. This dominance is observed for transport between islands with sixfold connectivity (hubs) and persists over a range of distances, suggesting a bulk $4e$ phase.
• Half-Vortex Dynamics: The system at $f=1/3$ exhibits an extensive degeneracy of ground states analogous to the antiferromagnetic Ising model on a triangular lattice. The lowest-energy excitations are identified as half-vortices (flux $\Phi_0/2$) formed at the junction of three zero-energy domain walls. The dynamics of these half-vortices (creation, winding, and annihilation of pairs) explain the halved periodicity of the current-phase relation.
• Phase Transition and Correlations: The helicity modulus analysis indicates a BKT transition at $T_{BKT} \approx 0.227 E_J$, consistent with the unbinding of half-vortices. Below this temperature, correlation functions reveal a unique signature of the $4e$ phase: the correlation of the $2e$ phase operator ($e^{i\phi}$) decays exponentially due to the proliferation of domain walls, while the $4e$ phase operator ($e^{i2\phi}$) retains quasi-long-range order (algebraic decay).
• Role of Charging Energy: When a small charging energy $E_C$ is introduced, quantum fluctuations lift the extensive ground-state degeneracy via an "order-by-disorder" mechanism, favoring specific vortex lattice configurations (e.g., honeycomb or stripe patterns) at $T=0$. However, the energy splitting between these ordered states is extremely small ($\sim 10$ mK). Consequently, at intermediate temperatures ($T^* < T < T_{BKT}$), the entropic contribution of the degenerate manifold overcomes the energy splitting, restoring the disordered $4e$ phase.

Significance and Claims
The paper claims to provide strong numerical evidence that frustrated Josephson junction dice arrays at $f=1/3$ realize a superconducting phase mediated by Cooper quartets. The authors assert that this phase is a classical precursor to states with topological order, as the half-vortices present in the system belong to the same topological class as magnetic excitations in Kitaev's surface codes and $Z_2$ gauge theories.

The study confirms Korshunov's conjecture that the low-energy physics of the dice lattice at $f=1/3$ is characterized by the proliferation of zero-energy domain walls and confined half-vortices. The results suggest that hybrid superconductor-semiconductor dice arrays offer a scalable and tunable platform for observing these exotic phenomena. The authors conclude that while charging energy and mutual inductance may induce ordering at extremely low temperatures ($<10$ mK), the $4e$ superfluid phase is expected to dominate the physics over a broad, experimentally accessible temperature range, making these systems suitable candidates for studying fractional vortices and potentially realizing topologically protected quantum memories.