Nearly-Symmetric SQUIDs: Theory, Design & Applications

• Nearly-symmetric SQUIDs are dc superconducting interferometers with nearly identical weak links and loop arms that yield interference patterns close to the ideal cosine law.
• They leverage controlled deviations in critical currents, junction parameters, and CPR harmonics to adjust modulation depth, flux sensitivities, and operational stability.
• Fabrication techniques, from single-step lithography to hybrid SNS and flexible material designs, enable high modulation depth and low noise in quantum sensing applications.

Searching arXiv for recent and foundational papers on nearly-symmetric SQUIDs. arXiv search query: "nearly symmetric SQUID current-phase relation nanoSQUID symmetry" Nearly-symmetric superconducting quantum interference devices (SQUIDs) are dc SQUIDs in which the two weak links and the two loop arms are close to the ideal symmetric limit, but not exactly identical. In the ideal symmetric case, both junctions are identical, the loop is geometrically symmetric, and the low-inductance approximation yields a $\Phi_0$-periodic critical-current modulation with full destructive interference. In practice, small deviations in critical current, shunt resistance, capacitance, arm inductance, interface transparency, or current–phase relation (CPR) leave finite residual critical current at nominal minima, flux-axis offsets, or distortions of $I_c(\Phi)$ and $V(\Phi)$. Recent work shows that near symmetry can be a design condition, a fabrication outcome, or a gate-tuned operating regime; in InSb nanoflag devices, for example, back-gate control moves the same interferometer continuously from symmetric to nearly-symmetric to strongly asymmetric behavior by modifying transparency and higher-harmonic CPR content (Martínez-Pérez et al., 2016, Chieppa et al., 26 Apr 2025).

1. Ideal reference and operational definition

A dc SQUID consists of a superconducting loop of inductance $L$ interrupted by two Josephson junctions. The standard phase constraint is

$$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$

and in the low-inductance limit $\beta_L \ll 1$ one neglects the self-induced flux $LJ$, so that

$$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$

For an ideal symmetric dc SQUID, $I_{01}=I_{02}=I_0$, the loop inductance is shared equally between the two arms, and the critical current is

$$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$

This is the reference behavior: a symmetric, $I_c(\Phi)$0-periodic interference pattern with full modulation depth from $I_c(\Phi)$1 down to $I_c(\Phi)$2 as $I_c(\Phi)$3 (Martínez-Pérez et al., 2016).

Near symmetry denotes the regime in which these symmetric formulas remain useful but acquire modest corrections. For unequal critical currents, the standard extension is

$I_c(\Phi)$4

If $I_c(\Phi)$5 and $I_c(\Phi)$6 with $I_c(\Phi)$7, then $I_c(\Phi)$8, $I_c(\Phi)$9, and even small asymmetry leaves a finite residual critical current at $V(\Phi)$0. Arm-inductance asymmetry $V(\Phi)$1 produces a flux-axis shift and unequal slopes in $V(\Phi)$2 and $V(\Phi)$3 (Martínez-Pérez et al., 2016).

A concrete realization of this regime is the hybrid InAs nanowire–vanadium proximity SQUID, in which the two SNS weak links are adjacent uncovered segments of the same nanowire. Fits yield $V(\Phi)$4 and $V(\Phi)$5, so $V(\Phi)$6, and the measured modulation depth is approximately $V(\Phi)$7. The data are well fit by the standard two-junction expression above, indicating a nearly-symmetric device with a CPR close to sinusoidal (Spathis et al., 2010).

2. Mechanisms that move a SQUID away from exact symmetry

In practice, near symmetry is limited by three broad classes of perturbation: critical-current asymmetry, dissipative-parameter asymmetry, and geometric asymmetry. The nanoSQUID review emphasizes $V(\Phi)$8, different shunt resistances and capacitances, and unequal widths and thicknesses of the two arms, which alter both geometric and kinetic inductance. These departures reduce modulation depth, shift the interference pattern, distort $V(\Phi)$9, and complicate optimum biasing, even when they are “small” rather than catastrophic (Martínez-Pérez et al., 2016).

Modern ballistic and highly transparent weak links introduce a second mechanism: symmetry of geometry does not imply symmetry of the effective interferometer response if the CPRs are skewed or contain large higher harmonics. In the InSb nanoflag SQUIDs, the Josephson current is modeled by a 2D two-band Hamiltonian with Rashba spin–orbit and Zeeman terms, and the CPR is expanded as

$L$0

At $L$1 V the CPR is significantly skewed, with maximum at $L$2, whereas at $L$3 V it is closer to sinusoidal. In the symmetric geometry, the interference at $L$4 V remains $L$5-periodic but its minima do not reach zero; at $L$6 V the same device reaches zero at the minima. The paper interprets this as a nearly-symmetric SQUID whose deviation from full modulation is a probe of higher-harmonic CPR content rather than gross geometric mismatch (Chieppa et al., 26 Apr 2025).

A common misconception is that shallow modulation necessarily diagnoses poor matching of the two weak links. Planar NbN nanobridge nanoSQUIDs provide a counterexample. Their two bridges are patterned nominally identically in a single e-beam exposure and a single etch step, with rounded bridge ends to avoid current crowding, yet the measured modulation is only about $L$7 at $L$8 K and about $L$9 in a $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$0-nm device at $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$1 mK. The same work reports an effective inductance $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$2 and $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$3, and explicitly attributes the zig-zag interference pattern to long weak-link bridges. In that case, shallow modulation is dominated by inductive screening and long-bridge physics, not necessarily by large bridge asymmetry (Holzman et al., 2019).

A third route away from exact symmetry is geometric deformation in three dimensions. Flexible amorphous MoSi and WSi SQUIDs are lithographically symmetric planar loops with two nanowire weak links patterned identically, but bending introduces a curved geometry and strain. In one representative aMo$$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$4Si$$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$5 device, the oscillation period changes from $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$6 in the flat state to $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$7 at $$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$8, while conventional superconducting parameters change only weakly. The authors argue that simple geometric area change and simple field-lensing estimates are insufficient, and identify geometry, magnetic-field inhomogeneity, and strain as the relevant perturbations of an otherwise nearly-symmetric interferometer (Suleiman et al., 2020).

3. Interference, CPR structure, and vorticity stability

In nearly-symmetric two-junction SQUIDs with negligible self-inductance, the total current can be written as

$$\delta_1 - \delta_2 + 2\pi n = \frac{2\pi}{\Phi_0}\,(\Phi + L J),$$9

so the critical current is

$\beta_L \ll 1$0

For identical sinusoidal CPRs this reduces to the textbook cosine law. For nearly-symmetric skewed CPRs, higher harmonics alter both the shape and the modulation depth. In the InSb nanoflag platform the loop inductance is estimated as $\beta_L \ll 1$1 pH, the critical current as $\beta_L \ll 1$2 nA, and the inductance parameter as $\beta_L \ll 1$3, so self-inductance is negligible and the deviations from ideality are traced to CPR structure rather than loop screening (Chieppa et al., 26 Apr 2025).

The notion of near symmetry extends beyond two-wire interferometers when the weak links themselves are superconducting nanowires with linear CPR. In the multiple-nanowire SQUID model,

$\beta_L \ll 1$4

and neighboring-wire phases satisfy

$\beta_L \ll 1$5

where $\beta_L \ll 1$6 is the loop vorticity between wires $\beta_L \ll 1$7 and $\beta_L \ll 1$8. The model predicts that the critical current is a multi-valued function of magnetic field, and introduces vorticity stability regions (VSRs), namely regions in the current–magnetic-field plane where, for a given distribution of vortices, the currents in all wires remain below their critical values. For two-wire devices the VSRs are rhombic; for three- and four-wire devices they become diamonds, flat-top or flat-bottom diamonds, large triangles, small triangles, gliders, trapezoids, tilted triangles, and kites, depending on vorticity differences and disorder. The maximum critical current curves obey $\beta_L \ll 1$9 symmetry, while each VSR obeys $LJ$0 symmetry (Sun et al., 29 May 2025).

This multivortex framework clarifies a broader point about nearly-symmetric interferometers. In a symmetric geometry, full modulation does not follow automatically from “two equal branches”; it depends on the CPR, the admissible vorticity sectors, and the threshold phase $LJ$1. The multiple-wire model identifies a condition under which VSRs become disjoint and the modulation reaches $LJ$2:

$LJ$3

For $LJ$4 this gives the familiar $LJ$5; for $LJ$6 it gives $LJ$7; for $LJ$8, $LJ$9 (Sun et al., 29 May 2025).

4. Materials platforms and fabrication strategies

The most direct route to near symmetry is a fabrication flow that defines both weak links in the same material stack, the same lithographic step, and the same local environment. Cross-type Nb/AlO$$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$0/Nb SIS nanoSQUIDs exemplify this strategy. The junctions are defined by overlap of a trilayer strip and a U-shaped Nb wiring layer, with no unintended Nb/Nb overlap and essentially no extra parasitic capacitance beyond the intrinsic junction capacitance $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$1. The junction area is $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$2, the capacitance density is $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$3, and the total junction capacitance is $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$4 per junction. For the optimized devices at $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$5 K, the per-junction parameters are $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$6, $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$7, $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$8, and $$\delta_1 - \delta_2 \approx \frac{2\pi}{\Phi_0}\,\Phi .$$9, with $I_{01}=I_{02}=I_0$0 near unity and device-to-device spreads below $I_{01}=I_{02}=I_0$1. This is an explicitly nearly symmetric implementation in both geometry and RCSJ parameters (Schmelz et al., 2017).

Hybrid SNS implementations achieve near symmetry by defining both junctions in a single normal conductor. In the InAs nanowire–vanadium SQUID, both weak links are two uncovered segments of the same nanowire contacted by Ti/V electrodes, with $I_{01}=I_{02}=I_0$2 nm per segment. Because both junctions share the same nanowire diameter, doping, and surface treatment, the device reaches $I_{01}=I_{02}=I_0$3 and approximately $I_{01}=I_{02}=I_0$4 modulation depth, with Josephson coupling visible up to $I_{01}=I_{02}=I_0$5 K and voltage modulation observed up to the vanadium $I_{01}=I_{02}=I_0$6 K (Spathis et al., 2010).

Planar nanobridge devices pursue the same goal in ultrathin films. The on-chip integrable planar NbN nanoSQUID is patterned in a single lithographic step from $I_{01}=I_{02}=I_0$7-NbN films as thin as $I_{01}=I_{02}=I_0$8 nm on Si, with a $I_{01}=I_{02}=I_0$9 washer, two $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$0 nm by $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$1 nm weak links, and rounded constriction geometry to suppress current crowding. The device operates from $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$2 mK to $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$3 K, in $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$4 T parallel field and $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$5 mT perpendicular field, while potential operation higher than $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$6 T has also been shown. The symmetry is lithographic and process-driven, while the nonideal modulation depth is dominated by kinetic inductance and long-bridge physics (Holzman et al., 2019).

The InSb nanoflag platform adds electrostatic tunability to a nominally symmetric planar SNS geometry. In its symmetric SQUID, both junctions are narrow planar Josephson junctions with

$$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$7

embedded in a Nb loop of geometrical area $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$8 and effective area $$I_{c,\rm SQUID}(\Phi) = 2 I_0 \left|\cos\left(\frac{\pi\Phi}{\Phi_0}\right)\right| .$$9. The same geometry supports both textbook-like full modulation at low transparency and nearly-symmetric incomplete modulation at high transparency because the gate modifies the CPR rather than the loop itself (Chieppa et al., 26 Apr 2025).

A more unconventional materials direction uses flexible amorphous superconductors. In the MoSi and WSi devices, both weak links are patterned identically and the loop arms are equal in length to first order, so the flat device is nearly symmetric. Bending then perturbs the interferometer through curvature and strain rather than lithographic mismatch. The control experiment in which the bent device is rotated by $I_c(\Phi)$00 so that the weak links are perpendicular to the bending direction yields an interference pattern identical to the flat case, isolating the role of weak-link orientation relative to flexure (Suleiman et al., 2020).

5. Flux-to-voltage transfer, noise, and low-back-action readout

For conventional dc SQUID operation above the critical current, the relevant quantity is the flux-to-voltage transfer function

$I_c(\Phi)$01

and in nearly-symmetric devices a large modulation depth and a symmetric, sinusoid-like $I_c(\Phi)$02 are advantageous because they enlarge the linear range and increase $I_c(\Phi)$03. In the InSb nanoflag SQUIDs, the devices are non-hysteretic at $I_c(\Phi)$04 mK and can be operated as flux-to-voltage transducers. In the asymmetric device, $I_c(\Phi)$05 oscillations have amplitude $I_c(\Phi)$06 at $I_c(\Phi)$07 nA, with $I_c(\Phi)$08 at $I_c(\Phi)$09 V and up to $I_c(\Phi)$10 at $I_c(\Phi)$11 V. In the symmetric device, a similar modulation of about $I_c(\Phi)$12 occurs at $I_c(\Phi)$13 nA, with $I_c(\Phi)$14. Using a room-temperature preamplifier noise of $I_c(\Phi)$15, the reported flux noises are $I_c(\Phi)$16 for the asymmetric device and $I_c(\Phi)$17 for the symmetric device (Chieppa et al., 26 Apr 2025).

Near symmetry is also central to pushing toward intrinsic white-noise limits. In the cross-type Nb/AlO$I_c(\Phi)$18/Nb devices, the $I_c(\Phi)$19 loop SQUID exhibits measured white flux noise $I_c(\Phi)$20, energy resolution $I_c(\Phi)$21, usable voltage swing $I_c(\Phi)$22, and transfer function up to $I_c(\Phi)$23. The paper attributes the near-ideal behavior to self-aligned symmetric junction definition, low parasitic capacitance, $I_c(\Phi)$24, and $I_c(\Phi)$25, emphasizing that asymmetry or excess parasitics would push the device away from the nearly quantum-limited regime (Schmelz et al., 2017).

A different readout principle appears in the nonlocal SNS SQUID. There, the Josephson junctions are diffusive SNS segments with additional normal-metal leads, and flux is inferred from the flux dependence of the critical current without applying a finite voltage across the superconducting loop. The relevant estimate is

$I_c(\Phi)$26

with $I_c(\Phi)$27, $I_c(\Phi)$28, and $I_c(\Phi)$29, giving a theoretical flux noise $I_c(\Phi)$30. The significance for nearly-symmetric SQUIDs is that symmetry now has to include nonequilibrium quasiparticle distributions as well as geometry, because selective quasiparticle injection can make a nominally symmetric SNS interferometer functionally asymmetric (Noh et al., 2020).

6. Arrays, differential architectures, and absolute-field operation

Near symmetry is beneficial at the single-cell level, but several advanced applications require breaking exact periodicity or enlarging linear range at the system level. A clear example is the trade-off between identical SQUID cells and absolute magnetometry in arrays. In two-dimensional SQUID arrays, existing absolute-magnetometer designs rely on incommensurate loop areas, but varying the physical area also changes the loop inductance and therefore $I_c(\Phi)$31. The synthetic-area-spread formalism resolves this by adding bare superconducting loops with no Josephson junctions. If $I_c(\Phi)$32 is the vector of physical SQUID-loop areas, $I_c(\Phi)$33 the bare-loop areas, and $I_c(\Phi)$34, then the effective SQUID areas are

$I_c(\Phi)$35

The result is an array that behaves as if it had an incommensurate physical area spread while keeping all physical SQUID loops identical, so $I_c(\Phi)$36 is identical for every DC SQUID loop (Monaghan et al., 19 Nov 2025).

A second route is to couple more than one SQUID degree of freedom into a common readout resonator. In the superconducting microwave magnetometer for absolute flux detection, two asymmetric dc SQUIDs are embedded in a common tank circuit and the two lowest eigenfrequencies are read out in reflection. The device reaches a modulation period $I_c(\Phi)$37 and yields a magnetic offset field $I_c(\Phi)$38. This is not a nearly-symmetric single SQUID in the textbook sense; rather, it uses controlled asymmetry between two SQUID loops to overcome the $I_c(\Phi)$39-periodicity that constrains symmetric single-loop readout (Günzler et al., 2021).

System-level symmetry also underlies gradiometric low-noise designs. The femto-Tesla Nb design for quantum-ready readouts uses many identical washers in parallel, a gradiometric layout, and an input-coil architecture chosen so that the total SQUID inductance remains near the optimum value $I_c(\Phi)$40 implied by $I_c(\Phi)$41. In the fractional-SQUID implementation, $I_c(\Phi)$42 washers in parallel give $I_c(\Phi)$43, while $I_c(\Phi)$44 one-turn input coils in series give $I_c(\Phi)$45, preserving inductance matching and gradiometric balance while suppressing parasitic resonances. The paper’s CAD example is a $I_c(\Phi)$46 aT-class gradiometer in a $I_c(\Phi)$47 footprint (Sochnikov et al., 2020).

Taken together, these developments indicate that nearly-symmetric SQUIDs are best understood as a hierarchy of design regimes rather than a single topology. At the junction level, near symmetry means closely matched weak links and arm inductances. At the interferometer level, it means controlled deviations from the ideal cosine law due to higher harmonics, finite kinetic inductance, or vorticity structure. At the array and resonator level, it means preserving identical cell performance while using additional loops, synthetic areas, gradiometry, or multi-SQUID coupling to expand linear range, dynamic range, or flux unambiguity without discarding the low-noise advantages of the symmetric limit.