⚡ Josephson Junction — AC & DC Josephson Effects
Explore DC and AC Josephson effects: a supercurrent tunnels through a thin insulating barrier; under applied voltage the phase oscillates at frequency 2eV/h.
About this simulation
This simulation integrates the RCSJ (resistively and capacitively shunted junction) model of a Josephson junction, in which the supercurrent I = Ic·sin(φ) depends only on the quantum phase difference φ between two weakly coupled superconductors. In DC mode a fixed drive current is applied and the phase equation of motion is solved with 4th-order Runge–Kutta, revealing whether the junction stays dissipationless or switches to a resistive state. In AC mode a fixed voltage forces the phase to wind at a constant rate, reproducing the Josephson voltage-frequency relation dφ/dt = 2eV/ħ. Three linked panels — phase portrait, time-domain supercurrent and I-V characteristic — let you watch these two regimes and their crossover in real time.
🔬 What it shows
The phase portrait plots φ mod 2π against dφ/dt, revealing closed orbits when the junction is superconducting and drifting trajectories once it switches to the resistive branch. The time-domain panel traces the supercurrent Is(t) against the ±Ic bounds, while the I-V panel sweeps drive current to build the average-voltage characteristic, showing the sharp knee at the critical current.
🎮 How to use
Switch the Mode selector between DC Josephson (current-biased) and AC Josephson (voltage-biased). Set the critical current Ic (0.1–2.0), the drive current ratio I/Ic (0–2.0) for DC mode, and the voltage bias V (0.05–3.0) for AC mode. The damping slider sets βc−½ (0.05–2.0), the inverse square root of the McCumber parameter, tuning the junction from overdamped to underdamped. Pause, Reset or open the Info panel with the on-screen buttons.
💡 Did you know?
Because the Josephson frequency fJ = 2eV/h links voltage directly to frequency through fundamental constants, arrays of junctions driven by microwaves produce quantised voltage steps (Shapiro steps) precise enough to define the international volt standard, while loops of junctions called SQUIDs measure magnetic fields down to the femtotesla scale.
Frequently asked questions
What is the difference between DC mode and AC mode here?
DC mode current-biases the junction: you set I/Ic and the simulation integrates the phase equation to find whether a static phase (zero voltage) or a rotating phase (finite voltage) results. AC mode voltage-biases it instead: you set V directly and the phase is forced to advance at the fixed rate dφ/dt = 2eV/ħ, always producing an oscillating supercurrent.
What does the damping slider actually control?
The damping control sets βc−½ , the inverse square root of the McCumber parameter that appears in the RCSJ equation of motion. Low values correspond to an underdamped, high-Q junction whose I-V curve is hysteretic, while high values correspond to an overdamped junction that switches smoothly between the superconducting and resistive branches with no hysteresis.
Why does the phase portrait sometimes trace closed loops and sometimes drift sideways?
A closed loop means the phase φ oscillates around a fixed point without ever completing a full rotation — the junction is superconducting and carries a static (or bounded) supercurrent. A trajectory that drifts continuously across the panel means φ is winding indefinitely, which corresponds to a running, time-averaged voltage: the junction has switched to its resistive branch.
What is the I-V panel plotting and why is there a knee?
It sweeps the DC drive current from zero up to 2.5·Ic, settles the phase dynamics at each value, and measures the resulting time-averaged voltage ⟨dφ/dt⟩. Below Ic the average voltage stays at zero because the phase can sit still; once the drive exceeds Ic no static phase can balance the current, so the phase starts rotating and the average voltage rises sharply — the knee marks the critical current.
What happens physically when the drive current exceeds Ic?
Once I/Ic passes 1, the junction can no longer sustain a constant phase difference, so φ begins rotating continuously instead of settling at a fixed value. Each full rotation of φ by 2π corresponds to one flux quantum passing through the junction, and the time-averaged rate of rotation shows up as a genuine voltage across it, moving the state indicator from "SC" to "Resistive".
Explore DC and AC Josephson effects: a supercurrent tunnels through a thin insulating barrier; under applied voltage the phase oscillates at frequency 2eV/h.
3D · Three.js / WebGL renderer · 60 FPS target · runs fully client-side, no install