Switching-current distribution of a Josephson junction (source code)

= Switching-current distribution of a Josephson junction
{title2=$p(I)=\Gamma(I)S(I)/r$}

A <switching-current distribution of a Josephson junction> records the bias <electric current> at the first detected transition to a sustained nonzero <electric potential difference>. Interpret it using a stipulated monotone ramp $dI/dt=r>0$, a reproducible initial trapped state, and an instantaneous escape <hazard function> $\Gamma(I)$. Neglect retrapping and readout delay. The <survival function> $S$ obeys $dS/dt=-\Gamma S$, so changing variable from time to <electric current> gives
$$
S(I)=\exp\left[-\frac1r\int_{I_0}^I\Gamma(u)\,du\right],\qquad
p(I)=-\frac{dS}{dI}=\frac{\Gamma(I)}rS(I),\qquad
\Gamma(I)=\frac{rp(I)}{1-\int_{I_0}^Ip(u)\,du}.
$$
Here $p$ is a <probability density function> per unit <electric current>, not a rate per unit time. Division by the surviving population is essential: only trials still trapped can escape. A typical switching current is consequently not identical to the noiseless <superconducting critical current>.

\Image[/josephson-escape-switching-readout.png]
{title=From trapped phase to a voltage event and an escape probability}
{description=Original functional readout model and chosen constant-rate calculation. Boxes represent functions rather than historical wiring. The plotted <survival function> and switching <probability density function> assume $\Gamma=1000\,\mathrm{s^{-1}}$ and $r=1\,\mathrm{\mu A\,s^{-1}}$ with perfect event detection and no retrapping; they are not experimental traces.}

For an independent illustrative check, choose constant $\Gamma_0=1000\,\mathrm{s^{-1}}$ and $r=1\,\mathrm{\mu A\,s^{-1}}$. Then
$$
S(I)=e^{-\Gamma_0(I-I_0)/r},\qquad
\langle I-I_0\rangle=\frac r{\Gamma_0}=1\,\mathrm{nA}.
$$
Half the trials have switched after $\ln2/\Gamma_0=0.693147\,\mathrm{ms}$, corresponding to $0.693147\,\mathrm{nA}$ of ramp. Doubling $r$ doubles the mean increment without changing the underlying escape rate. Real junction rates vary sharply with bias; this example isolates the statistical conversion.

Repeating the preparation, measuring the ramp, resolving the voltage onset, and allowing cooling after each dissipative event are prerequisites for this interpretation. Missed events, finite detection delay, changing <temperature> and retrapping alter the measured distribution. Even a correctly reconstructed rate needs the independent thermal and spectral checks before identifying its escape mechanism.