🔬 mesoscale physics

Dynamical Reduction of Two Series Josephson Junctions to a Synthetic High-Transparency Josephson Element

This paper establishes the quantitative validity and limitations of modeling two series-connected Josephson junctions as a single synthetic high-transparency element under finite-frequency drive, demonstrating that the reduced single-degree-of-freedom description remains accurate at low frequencies but deviates significantly as the drive frequency approaches the plasma-frequency scale.

Original authors: Claudio Guarcello, Sergio Pagano, Carlo Barone, Alessandro Bruno, A. Mert Bozkurt, Giovanni Filatrella

Published 2026-09-07

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Claudio Guarcello, Sergio Pagano, Carlo Barone, Alessandro Bruno, A. Mert Bozkurt, Giovanni Filatrella

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the world of superconducting electronics, where electricity flows without any resistance, the behavior of tiny barriers known as Josephson junctions is the key to building powerful quantum devices. These junctions act as special gates that allow pairs of electrons to tunnel through, creating a unique relationship between the electrical current flowing across them and the phase of the superconducting wave. While the simplest version of this relationship is a smooth, predictable curve, nature often provides more complex shapes when the junctions are engineered with specific properties. These complex shapes, which contain extra wiggles or harmonics in their response, are highly desirable for creating advanced circuits that can mix signals or amplify quantum information. However, building a single junction with these precise, complex traits is difficult. A clever workaround has emerged: instead of trying to forge a perfect, complex junction from scratch, scientists can connect two ordinary, simple junctions in a line. Under steady, unchanging conditions, this pair behaves exactly like a single, sophisticated junction with a tunable character.

The question researchers set out to answer was whether this clever shortcut holds up when the system is put under pressure. In real-world applications, these circuits are rarely static; they are driven by rapidly changing electrical currents that oscillate at high speeds. When the current changes quickly, the internal components of the junctions, such as their ability to store electrical charge and their resistance to flow, begin to act independently. This raises a critical doubt: does the pair of simple junctions still mimic the single complex one when the system is moving fast, or do the internal differences between the two break the illusion? To find out, a team of physicists simulated the behavior of two junctions connected in series, comparing the full, detailed motion of both components against a simplified model that treated them as a single effective unit.

The researchers focused on a specific setup where two junctions were placed one after the other, with one being slightly smaller than the other. They varied the size difference between them and subjected the system to alternating currents of different speeds and strengths. By running thousands of computer simulations, they tracked the voltage produced by the full two-junction system and compared it directly to the voltage predicted by the simplified single-junction model. They measured the difference between the two waveforms with extreme precision, looking for even the smallest mismatch in timing, height, or shape. The results showed that the simplified model is remarkably accurate, but only within a specific range. When the driving current oscillates slowly or with moderate strength, the two-junction system and the single-junction model produce nearly identical results, with errors so small they are almost invisible.

However, the agreement begins to crumble as the conditions become more extreme. When the frequency of the driving current approaches a specific natural speed at which the junctions naturally vibrate, or when the current amplitude becomes very large, the two systems start to diverge. In these high-energy regimes, the internal dynamics of the two separate junctions become important, and the simplified model can no longer capture the full complexity of the motion. The researchers found that the breakdown is most pronounced when the junctions are very different in size, but if the two junctions are nearly identical, the simplified model remains accurate even under more demanding conditions. In the most symmetric case, where the two junctions are essentially twins, the single-junction description becomes almost perfectly exact, even as the system vibrates intensely.

This work provides a clear map for engineers designing superconducting circuits. It confirms that treating a pair of series junctions as a single, tunable element is a valid and powerful strategy, but it also draws a hard line around where that strategy stops working. The simplified description is safe to use for low-frequency operations and moderate signal strengths, offering a reliable way to design complex quantum devices without needing to model every internal detail. Yet, as the system approaches its natural resonant frequencies or is pushed into highly nonlinear states, the full complexity of the two-junction system must be taken into account. The study does not suggest that the simplified model is wrong, but rather that it has a defined domain of validity, much like a map that is perfect for a city center but loses detail as you drive toward the distant countryside. By quantifying exactly where the error remains negligible and where it grows, the researchers have given circuit designers the confidence to use these synthetic elements where they work best, while knowing precisely when to switch to a more detailed approach.

Technical Summary: Dynamical Reduction of Two Series Josephson Junctions

Problem Statement
The current-phase relation (CPR) of a Josephson junction (JJ) is a fundamental parameter in superconducting electronics. While conventional tunnel junctions exhibit a sinusoidal CPR, nonsinusoidal CPRs with higher harmonics arise in weak links with finite channel transparency. Recent work has demonstrated that two conventional (sinusoidal) JJs connected in series can synthesize the energy-phase relation of a single effective weak link with tunable transparency, controlled by the geometric asymmetry of the junctions.

However, this mapping has primarily been established in the static limit, where capacitive and resistive currents are negligible. The validity of reducing a two-junction system to a single "synthetic" element under finite-frequency drive remains an open question. When driven at finite frequencies, the individual capacitive and dissipative responses of each junction introduce internal dynamical degrees of freedom and time scales that are eliminated in static reductions. These internal dynamics can significantly affect the collective behavior of coupled JJ systems, potentially invalidating the single-element approximation in driven superconducting circuits.

Methodology
The authors investigate the dynamical validity of this reduction by comparing a full two-junction model against a derived effective single-element model.

System Definition: The study considers two conventional JJs in series with critical currents $I_{c1}$ and $I_{c2}$ , capacitances $C_1$ and $C_2$ , and resistances $R_1$ and $R_2$ . The asymmetry is defined by $\alpha = A_2/A_1 = I_{c2}/I_{c1} = C_2/C_1$ , assuming fabrication on the same chip with identical process parameters. The resistive terms are treated phenomenologically and do not necessarily scale with area.
Full Model: The dynamics are described by the Resistively and Capacitively Shunted Junction (RCSJ) equations for both junctions, driven by a common AC current $i(\tau) = i_{ac} \sin(\Omega\tau)$ .
Effective Model Derivation:
- Static CPR: The static mapping yields a nonsinusoidal effective CPR, $i_{eff}(\phi)$ , corresponding to a specific transparency $T$ .
- Dynamic Parameters: To extend this to finite frequencies, the authors derive effective capacitive ( $C_{eff}$ ) and dissipative ( $R_{eff}$ ) coefficients for the collective phase $\phi = \phi_1 + \phi_2$ . These are obtained by equating the electrostatic energy and power dissipation of the effective element to the sum of the individual junctions in the low-frequency, small-phase limit.
- Plasma Frequency: A key finding in the derivation is that due to the area scaling of both $I_c$ and $C$ , the effective plasma frequency of the synthetic element ( $\omega_{p,eff}$ ) is identical to the bare plasma frequencies of the individual junctions ( $\omega_{p1} = \omega_{p2}$ ).
Numerical Analysis: The authors perform numerical simulations comparing the full two-junction voltage waveform ( $v_{full}$ ) with the effective single-junction voltage ( $v_{eff}$ ) across a broad range of drive frequencies ( $\Omega$ ) and amplitudes ( $i_{ac}$ ) for various asymmetry values ( $\alpha \in \{0.25, 0.50, 0.75, 0.99\}$ ).
Metric: Accuracy is quantified using a normalized root-mean-square (RMS) waveform error ( $\epsilon_V$ ), which captures discrepancies in amplitude, phase, and harmonic content.

Key Results

Validity Domain: A broad region of low error ( $\epsilon_V \lesssim 10^{-4} - 10^{-5}$ ) exists at low drive frequencies and moderate amplitudes. In this regime, the effective model accurately reproduces the full system's dynamics, including the finite-frequency dynamical enhancement near resonance.
Symmetry Dependence: The accuracy of the reduction improves significantly as the junctions become more symmetric ( $\alpha \to 1$ ). In the limit $\alpha = 1$ , the two junctions are identical, and the effective description becomes dynamically exact.
Breakdown Conditions: Pronounced deviations emerge under two main conditions:
1. High Frequency: As the drive frequency approaches the plasma frequency scale ( $\Omega \approx 1$ ), internal relative-phase dynamics become relevant, causing the single-element approximation to fail.
2. Large Amplitude: At larger drive amplitudes, nonlinear effects and internal phase differences cause significant waveform mismatches.
Resonance Features: The analysis identifies a third-order superharmonic resonance condition ( $\Omega \approx \Omega_p/3$ ) where narrow structures in the error map appear.

Significance and Claims
The paper provides a quantitative dynamical criterion for the use of reduced single-element descriptions in superconducting circuits. It establishes that the static mapping of two series JJs to a synthetic high-transparency element is not merely a static approximation but remains operationally valid over a wide range of finite-frequency drives, provided the system is not driven deep into the nonlinear, near-resonant regime.

The authors claim that their work defines the specific frequency and driving-amplitude ranges where the static CPR reduction can be safely employed. This is particularly relevant for the design of driven superconducting circuits, such as those utilizing synthetic CPRs for three-wave mixing or parametric amplification, ensuring that the simplified models used for design do not overlook critical internal dynamical effects that arise at higher frequencies or amplitudes. The results confirm that the reduction is robust for adiabatic and near-adiabatic regimes but requires the full two-junction description when internal degrees of freedom become active.

More like this