---
title: Frustrated Josephson Tunneling
url: https://www.emergentmind.com/topics/frustrated-josephson-tunneling
type: topic
---

# Frustrated Josephson Tunneling

Frustrated Josephson tunneling denotes Josephson transport in which the superconducting phases, or more generally the internal order-parameter variables that enter the Josephson energy, cannot simultaneously minimize all coupling terms. In the literature this occurs in several distinct but structurally related settings: magnetic-flux frustration in Josephson junction arrays (JJAs), sign competition between \(0\)- and \(\pi\)-junction couplings, fluxoid-quantization constraints in loop geometries, positive interband Josephson couplings in multiband superconductors, and anisotropic couplings between triplet \(d\)-vectors. The resulting phenomenology includes Berezinskii–Kosterlitz–Thouless (BKT) criticality, emergent \(\mathrm U(1)\times \mathbb Z_2\) structure, flat-band crossovers, \(\varphi\)-contacts, spontaneous currents, edge-current phases, pinning–depinning critical currents, irrational-flux superconducting dips, and nonintegral flux trapping [2001.07877], [1112.1265], [1410.2021], [1405.3269], [2604.24734].

## 1. Microscopic origin of frustration

In a JJA, frustration enters directly through the gauge-invariant phase difference. The phase-only Hamiltonian is
\[
H_{\rm JJA} = \frac{(2e)^2}{2} \sum_{\langle i,j\rangle} n_i C^{-1}_{ij} n_j - E_J\sum_{\langle i,j\rangle} \cos(\phi_i-\phi_j-A_{ij}),
\]
with
\[
A_{ij} = \frac{2\pi}{\Phi_0}\int_i^j \vec A\cdot d\vec l,
\]
so magnetic flux through a plaquette shifts the preferred bond phases. In the classical regime this maps to an XY model with \(\vec S_i=(\cos\phi_i,\sin\phi_i)\), making the Josephson phase the XY spin angle and the flux the source of frustration [2001.07877].

A second mechanism is alternating-sign Josephson coupling. In one-dimensional diamond and sawtooth arrays, the couplings \(\alpha_{ij}\) alternate in sign and are parametrized by \(\alpha=1-2f\). The frustration parameter \(f\) controls a crossover from non-frustrated to frustrated regimes, with the lowest linear band becoming completely flat at a critical \(f_c\); in the diamond chain \(f_c^{D}=2/3\), and in the sawtooth chain \(f_c^{ST}=3/4\) [1801.03842]. Closely related sign competition appears in quasi-one-dimensional ladders with \(J_1<0\) inter-rung couplings and frustrating \(J_2>0\) intra-rung couplings, where nontrivial phase differences produce chiral edge currents [1607.08201].

A third mechanism is topological rather than local. In long Josephson tunnel junctions with doubly connected electrodes, fluxoid quantization imposes a global constraint on the phase profile, so the gauge-invariant phase difference is restricted not only by the sine-Gordon equation and boundary conditions but also by loop topology. The paper identifies this as a phase-frustration constraint and attributes the reduction of critical current to the impossibility of satisfying the local static equations, the fluxoid constraint, and the minimum-energy condition simultaneously [1112.1265]. An analogous topological incompatibility appears in Pb/Ru/Sr\(_2\)RuO\(_4\) junctions when the bulk state of Sr\(_2\)RuO\(_4\) is taken to be chiral \(p\)-wave: the coupling term \(\cos[\phi(\theta,z)-N\theta]\) favors a winding phase for \(N=\pm1\), while the gradient term favors \(\phi\) constant, creating a topologically frustrated junction [1405.3269].

Internal order-parameter structure can itself be the source of frustration. In multiband London models with positive interband Josephson couplings, each Josephson term prefers a phase difference of \(\pi\), but with three components the pairwise preferences cannot all be satisfied. The ground state therefore acquires two opposite chiralities and an effective \(\mathrm U(1)\times\mathbb Z_2\) symmetry structure [1410.2021]. In triplet systems, the Josephson energy depends not only on the gauge phase difference but also on the relative orientation of the \(d\)-vectors. The effective bond energy can be written as
\[
F_J = -\sum_{\langle n,m\rangle} |\tilde{\mathcal J}_{nm}| \cos(\Delta\varphi_{nm}+\mathcal A_{nm}),
\]
where \(\mathcal A_{nm}\) is an emergent geometric phase generated by the internal spin texture of the Cooper pair; nonzero loop sums of \(\mathcal A_{nm}\) create frustrated circulating-current states [2604.24734]. Related work derives Heisenberg-like, Dzyaloshinskii–Moriya-like, and \(\Gamma\)-type anisotropic couplings between \(d\)-vectors from frustrated spin textures in the barrier, again turning Josephson tunneling into a competition among incompatible preferred configurations [2506.15661].

## 2. Lattice realizations and effective models

The canonical two-dimensional realization is the square-lattice JJA implementing the XY model. At zero flux the couplings are unfrustrated, \(J_{ij}=J>0\), while at half flux per plaquette the fully frustrated XY model has \(J_{ij}=\pm J\) with one effectively antiferromagnetic bond per plaquette and a doubly degenerate ground state apart from global \(\mathrm U(1)\) rotation [2001.07877]. This realization directly connects frustrated Josephson tunneling to the longstanding problem of coupled \(\mathrm U(1)\) and \(\mathbb Z_2\) ordering.

One-dimensional geometries display a different but related structure. In diamond and sawtooth chains the linearized Josephson dynamics has multiband dispersions \(\omega_m(k)\), and the lowest band becomes flat precisely at the crossover into the frustrated regime. In the same parameter range the ordinary phase correlator \(C_1\) crosses from long-range to short-range behavior, while higher harmonics such as \(C_2\) in the diamond chain and \(C_3\) in the sawtooth chain recover long-range order deep in the frustrated regime [1801.03842]. This establishes a direct correspondence between spectral flattening and the loss of ordinary phase coherence.

Open-boundary ladders introduce boundary-specific frustrated states. In a quasi-one-dimensional Josephson ladder with diagonal couplings, the Josephson energy maps to the expectation value of a two-band tight-binding Hamiltonian subject to a uniform-density constraint. Above \(t=0.5\), where \(t\) is the ratio of intra-rung to inter-rung coupling strength in the mapped model, circulating currents nucleate at the edges and then propagate inward as \(t\) increases; for an infinite ladder the edge-current phase disappears at \(t=1\), coincident with the opening of a gap between the flat and dispersive bands [1607.08201].

More constrained frustrated networks generate effective Ising descriptions. In the vertex-sharing Kagome lattice with periodically arranged \(0\)- and \(\pi\)-junctions, each frustrated triangle supports either a vortex or an antivortex, so the low-energy sector is represented by an Ising variable \(\sigma=\pm1\). Flux quantization around the hexagons then produces highly anisotropic long-range interactions between these effective spins [2308.07143]. In frustrated sawtooth arrays of small quantum junctions, the frustrated regime \(f>f_{cr}=0.75\) produces \(2^N\) degenerate classical minima and low-energy vortex/antivortex doublets per triangle, again motivating effective spin models for the coherent regime [2412.00391]. A related transmission-line embedding yields an effective long-range \(XX\) chain,
\[
\hat H = \sum_{n=1}^N \Delta\,\hat\sigma_{x,n} + \frac{1}{2}\sum_{n\neq m} J(|n-m|)\big[\hat\sigma_{x,n}\hat\sigma_{x,m}+\hat\sigma_{y,n}\hat\sigma_{y,m}\big],
\]
with \(J(x)\sim x^{-\beta}\) at long distance [2407.03928].

Aperiodicity offers a further variant. In a one-dimensional ladder with two incommensurate plaquette areas, the ground-state energy \(E(f)\) becomes quasiperiodic in frustration and its deep minima are well captured by an independent-plaquette model. The dominant frequencies in the power spectrum occur at \(\omega=1\), \(\omega=\alpha\), and \(\omega=\alpha\pm1\), where \(\alpha\) is the area ratio [1105.2861].

## 3. Ordered states, degeneracy, and phase transitions

In the unfrustrated square JJA, the finite-temperature transition is the standard BKT transition driven by vortex–antivortex unbinding. In the circuit-QED realization, the transition occurs near \(k_B T/E_J \simeq 0.91\), and the correlation-function exponent reaches the BKT value \(\eta=\frac14\) at the transition [2001.07877]. At full frustration \(\Phi/\Phi_0=1/2\), the same platform identifies a transition near \(k_B T/E_J \simeq 0.43\), consistent with prior numerical studies for the fully frustrated XY model. The same work emphasizes that the exact relationship between the \(\mathrm U(1)\) and \(\mathbb Z_2\) transitions remains debated.

The distinction between ordinary phase order and hidden higher-harmonic order is central in one-dimensional frustrated arrays. For the diamond and sawtooth chains, low-temperature ordinary correlations decay slowly for \(0<f<f_c\), but for \(f>f_c\) the factor \(\{\cos(pu_0)\}^n\) suppresses or oscillates the correlators and drives \(C_1\) short-ranged. Near \(f\simeq1\), however, \(C_2\) in the diamond chain and \(C_3\) in the sawtooth chain regain long-range behavior. The paper interprets the diamond-chain result as a \(2e\)-to-\(4e\) type effect: ordinary phase coherence is short-ranged while pair-like coherence persists [1801.03842]. A common misconception is therefore that frustration simply destroys order; in these arrays it can instead transfer coherence into higher harmonics.

Checkerboard frustration from four-terminal junctions modifies the usual commensuration logic. In the Nb-Pt-Nb four-terminal array, alternating plaquettes carry frustrations \(f\) and \(f'\), and superconducting resistance dips occur when \(f+f'\in\mathbb Z\), equivalently \(f=-f' \pmod{1}\). Because \(f'=\beta f\) with \(\beta=A'/A\), the dip positions
\[
f_m = \frac{m}{1+\beta}
\]
are generally irrational in the single-plaquette frustration variable. The paper interprets this as stabilization of the BKT superconducting phase even when the flux per individual plaquette is irrational [2503.14423].

Quantum frustrated arrays exhibit further ordered phases. In the sawtooth arrays of small junctions, the non-frustrated regime supports disordered (insulating) and ordered (superconducting) phases controlled by the competition between Josephson coupling and charging energies, while the frustrated regime supports a spatially disordered vortex–antivortex superposition phase and, for long arrays with strong effective inter-cell coupling, a stripe phase of vortices and antivortices [2412.00391]. In the multiband \(\mathrm U(1)\times\mathbb Z_2\) problem, the field couples mainly to the \(\mathrm U(1)\) sector, so vortex-lattice melting and chirality restoration need not coincide; the resulting chiral metallic phase is a resistive vortex liquid that still breaks time-reversal symmetry [1410.2021].

## 4. Dynamical response and nonequilibrium phase evolution

A major development has been the use of weak microwave spectroscopy to detect frustrated Josephson criticality without driving the array far from equilibrium. In the \(100\times100\) plaquette JJA inside a three-dimensional cavity, the measured observable is the complex reflection coefficient \(S_{11}\), whose linewidth is governed by the dissipative susceptibility \(\mathrm{Im}\,\chi(\omega_c)\). The effective total linewidth is
\[
\kappa_{\rm tot}=\kappa_{\rm ext}+\kappa_{\rm int}+g\frac{\mathrm{Im}\,\chi(\omega_c)}{2e}.
\]
As temperature increases, the linewidth broadens, peaks near the transition, and narrows again. The paper interprets this nonmonotonic behavior as a signature of vortex dynamics: bound pairs at low temperature, strong dissipation from unbinding near the transition, and screened free-vortex response above it [2001.07877]. The same work states explicitly that this probe is primarily sensitive to the dynamics of the \(\mathrm U(1)\) part, not directly to the \(\mathbb Z_2\) chirality sector.

Driven frustrated junctions can generate oscillatory phase dynamics even without harmonic forcing. For a multichannel junction between an \(s\)-wave superconductor and a three-band superconductor, the energy
\[
E= -\sum_{z=i,j,k}J_{zs}\cos\phi_{zs} -J_{ij}\cos(\phi_{is}-\phi_{js}) -J_{ik}\cos(\phi_{is}-\phi_{ks}) -J_{jk}\cos(\phi_{js}-\phi_{ks})
\]
has fourteen candidate ground states, including trivial \(0/\pi\) states and nontrivial \(\phi_0\) states. When the couplings are driven across boundaries in parameter space, the phase variables begin to oscillate and the voltage oscillates accordingly. The paper distinguishes coherent few-harmonic spectra from broader multi-harmonic spectra, depending on whether the trajectory settles into one minimum or moves among several nearby frustrated minima [2201.07697].

In transmission-line-embedded frustrated sawtooth chains, each triangular cell is a two-level system consisting of clockwise and anticlockwise persistent-current states. Coherent tunneling between them is described by the local \(\hat\sigma_x\) term, while the transmission line induces nonlocal exchange couplings that can dominate when the array length exceeds the screening scale \(\ell_0\simeq \sqrt{C/C_0}\). Exact diagonalization identifies a local-tunneling-dominated paramagnet, a compressible superfluid, and a weakly compressible superfluid [2407.03928].

Constraint-induced quantum dynamics also appears in the Kagome array. There the quantum regime reduces to a transverse-field Ising model in which \(\hat\sigma^z\) distinguishes vortex from antivortex and \(\hat\sigma^x\) generates macroscopic tunneling between them. In the strong-interaction limit the classical degeneracy is lifted and the ground state becomes an entangled superposition of classical vortex/antivortex patterns [2308.07143].

## 5. Multiband, topological, and triplet frustrated junctions

Multiband point contacts realize frustrated tunneling through interference among channel currents. In the dirty point contact between a single-band superconductor and a three-band superconductor, the total Josephson current is the sum of three band-resolved contributions with phase shifts \(x\), \(x+\phi\), and \(x+\theta\). When the three-band bank has broken time-reversal symmetry, the bulk equilibrium phases \((\phi,\theta)\) are nontrivial and occur in conjugate pairs, and the point contact can become a \(\varphi\)-contact with energy minima at finite phase shift rather than at \(0\) or \(\pi\) [1406.5693]. Near \(T_c\) the current-phase relation becomes close to sinusoidal, but the phase offsets remain controlled by the internal multiband state.

A related but more explicitly frustrated formulation treats positive interband Josephson couplings in a three-component London model. Because the field couples only to the overall \(\mathrm U(1)\) phase while the chiral \(\mathbb Z_2\) sector remains indirect, one may observe a vortex lattice or vortex liquid together with a separate Ising-like chirality transition [1410.2021]. This separates superconducting phase rigidity from time-reversal breaking.

Topological frustration in extended junctions produces a different critical-current mechanism. In the cylindrical Pb/Ru/Sr\(_2\)RuO\(_4\) geometry, the chiral phase of Sr\(_2\)RuO\(_4\) generates a spontaneous flux pattern on the interface. In a perfectly rotationally symmetric frustrated junction, any finite applied current shifts this flux pattern, so no static dc supercurrent can be sustained. When angular inhomogeneity pins the spontaneous flux line, a finite critical current reappears, but the limiting mechanism is then pinning–depinning of spontaneous flux rather than conventional vortex entry [1405.3269].

Triplet Josephson networks extend frustration to the internal spin geometry of Cooper pairs. For unitary triplet pairing, anisotropic Josephson coupling can be written in terms of Heisenberg-like, antisymmetric Dzyaloshinskii–Moriya-like, and \(\Gamma\)-type terms coupling neighboring \(d\)-vectors. In a three-grain ring this yields a critical antisymmetric coupling \(D_*=\frac{|J_0|}{\sqrt3}\), above which the minimum-energy texture becomes a chiral \(120^\circ\) configuration and can trap \(\Phi_0/2\) for the appropriate sign of \(J_0\) [2604.24734]. More generally, the accumulated geometric phase need not be quantized, so the trapped flux can be a nonintegral multiple of \(\Phi_0\). Parallel work derives the anisotropic \(d\)-vector couplings microscopically from frustrated noncollinear spin textures at the barrier and connects them to spatially modulated triplet pairing, anomalous vortices, and a Josephson diode effect [2506.15661].

Atomic-scale frustrated coupling has also been reported in a conventional-to-multiband setting. In scanned Josephson tunneling microscopy on Nb/FeSe, the net Josephson current is a quantum-coherent sum of two tunneling channels into two FeSe condensates with opposite signs, so the junction approaches a \(0\)-to-\(\pi\) transition as the relative transparency of the channels changes. The measured inequality
\[
I_JR_N < \frac{\pi}{2e}\Delta_{\mathrm{eff,min}}
\]
is presented as a direct signature of destructive interference from opposite-sign channels, and the extracted condensate-resolved superfluid modulations are spatially anti-correlated with cross-correlation coefficient about \(-0.72\) [2507.16758].

## 6. Experimental signatures and conceptual significance

The experimental signatures of frustrated Josephson tunneling are unusually diverse because the frustration can reside in geometry, topology, band structure, or internal pairing degrees of freedom. In JJAs and XY simulators, the central observables are linewidth anomalies in microwave reflection, vortex-sensitive damping, and transition temperatures inferred from \(\mathrm{Im}\,\chi\) [2001.07877]. In four-terminal arrays, dc resistance oscillations show a beating envelope and superconducting dips at checkerboard cancellation points \(f+f'\in\mathbb Z\), permitting an estimate of the weak-link area through \(\beta=A'/A\approx0.092\) and \(A'\approx(165\ \text{nm})^2\) [2503.14423]. In long junctions on superconducting loops, discrete critical-current shifts track trapped fluxoids, while the offset current \(I_c^{off}=I_c^+ + I_c^-\) improves sensitivity [1112.1265].

In boundary-sensitive and topological settings, the critical current itself becomes diagnostic. The Pb/Ru/Sr\(_2\)RuO\(_4\) analysis attributes the anomalous temperature dependence of \(I_c\) near the bulk \(T_c\) to a crossover from an unfrustrated topologically trivial regime to a topologically frustrated chiral regime [1405.3269]. In edge-current ladders, the onset, propagation, and eventual disappearance of chiral edge currents track the change in the underlying two-band spectrum [1607.08201].

Phase-sensitive spectroscopy is especially powerful in multiband systems. In dirty point contacts to three-band superconductors, persistent \(\varphi\)-contact behavior from low temperature to near \(T_c\) is proposed as a diagnostic of broken time-reversal symmetry [1406.5693]. In atomic-scale Nb/FeSe junctions, sublinear \(I_J\) versus \(1/R_N\), violation of the ordinary multiband lower bound, and local reconstruction of band-resolved transparencies provide direct evidence for frustrated Josephson interference in an \(s^\pm\) superconductor [2507.16758].

Across these realizations, one recurring theme is that frustration does not define a single phase but a constraint structure. It can generate BKT transitions, chirality order, flat bands, quasiperiodic energy landscapes, spontaneous supercurrents, entangled vortex–antivortex states, or nonintegral flux trapping, depending on which variables are overconstrained and which degrees of freedom remain soft. A second recurring theme is sector selectivity: some probes couple mainly to the global \(\mathrm U(1)\) dynamics, while others expose chirality, fluxoid number, channel interference, or \(d\)-vector texture. This suggests that the principal open problems are not only to identify new frustrated Josephson platforms, but also to disentangle which broken symmetry or constrained variable is actually being measured in each experiment [2001.07877], [1410.2021], [2604.24734].

Source: https://www.emergentmind.com/topics/frustrated-josephson-tunneling