
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

52271
10.1038/s41467-024-52271-z
Article
Current induced hidden states in Josephson junctions
http://orcid.org/0000-0001-7254-1023
Chen Shaowen shaowenchen@g.harvard.edu

1
http://orcid.org/0009-0003-2851-9250
Park Seunghyun 1
http://orcid.org/0000-0001-6024-8241
Vool Uri 12
Maksimovic Nikola 1
Broadway David A. 3
Flaks Mykhailo 3
http://orcid.org/0000-0002-8451-3385
Zhou Tony X. 15
http://orcid.org/0000-0003-1699-388X
Maletinsky Patrick 3
http://orcid.org/0000-0002-9493-268X
Stern Ady 4
http://orcid.org/0000-0002-6999-1039
Halperin Bertrand I. 1
http://orcid.org/0000-0002-5737-7963
Yacoby Amir yacoby@g.harvard.edu

1
1 https://ror.org/03vek6s52 grid.38142.3c 0000 0004 1936 754X Department of Physics, Harvard University, Cambridge, MA 02138 USA
2 https://ror.org/01c997669 grid.419507.e 0000 0004 0491 351X Max Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany
3 https://ror.org/02s6k3f65 grid.6612.3 0000 0004 1937 0642 Department of Physics, University of Basel, Klingelbergstrasse 82, Basel, CH-4056 Switzerland
4 https://ror.org/0316ej306 grid.13992.30 0000 0004 0604 7563 Weizmann Institute of Science, Rehovot, 76100 Israel
5 Present Address: Northrop Grumman Mission Systems, Linthicum, MD 21090 USA
14 9 2024
14 9 2024
2024
15 80591 7 2024
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Josephson junctions enable dissipation-less electrical current through metals and insulators below a critical current. Despite being central to quantum technology based on superconducting quantum bits and fundamental research into self-conjugate quasiparticles, the spatial distribution of super current flow at the junction and its predicted evolution with current bias and external magnetic field remain experimentally elusive. Revealing the hidden current flow, featureless in electrical resistance, helps understanding unconventional phenomena such as the nonreciprocal critical current, i.e., Josephson diode effect. Here we introduce a platform to visualize super current flow at the nanoscale. Utilizing a scanning magnetometer based on nitrogen vacancy centers in diamond, we uncover competing ground states electrically switchable within the zero-resistance regime. The competition results from the superconducting phase re-configuration induced by the Josephson current and kinetic inductance of thin-film superconductors. We further identify a new mechanism for the Josephson diode effect involving the Josephson current-induced phase. The nanoscale super current flow emerges as a new experimental observable for elucidating unconventional superconductivity, and optimizing quantum computation and energy-efficient devices.

Diract imaging of supercurrent flow at a Josephson junction has been inaccessible in experiment. Here, using nanoscale magnetometry, the authors find large kinetic inductance of thin film superconductors can lead to competing Josephson vortex states hidden below the critical current, and also provide a new route towards the Josephson diode effect.

Subject terms

Superconducting properties and materials
Characterization and analytical techniques
501100001711 Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (Swiss National Science Foundation) 188521 100010661 EC | Horizon 2020 Framework Programme (EU Framework Programme for Research and Innovation H2020) LEGOTOP No. 788715 501100001659 Deutsche Forschungsgemeinschaft (German Research Foundation) CRC/Transregio 183, EI 519/7-1 501100003977 Israel Science Foundation (ISF) 2074/19 100000001 National Science Foundation (NSF) DMR 1231319 100000183 United States Department of Defense | United States Army | U.S. Army Research, Development and Engineering Command | Army Research Office (ARO) W911NF-22-1-0248, W911NF-21-2-0147 100000936 Gordon and Betty Moore Foundation (Gordon E. and Betty I. Moore Foundation) GBMF 9468 issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Characterization and control over the super current flow is critical for Josephson junctions (JJs)1–3, which have become a building block in quantum and classical technology4–11 while remained a rich area of exploration into fundamental particles12–14 and unconventional superconductivity15–17. Compared to spectroscopic probes that measures the amplitude of the superconducting (SC) wave function18, the super current flow encodes the SC phase. Mapping the spatial distribution of super current has revealed the pairing symmetry of unconventional superconductors19,20, and recently identified screening current as the source of SC diode effect in SC/ferromagnet structures21. In addition, the local super current flow affects device parameters such as the impedance of SC circuits and anharmonicity of SC qubits due to the change in kinetic inductance22. Despite the scientific and technological relevance, direct visualization of the Josephson current flow and its response to external tuning knobs such as bias current and magnetic field remains experimentally beyond reach18,23–27. This is mostly due to the sensitive nature of the JJ, which responds to small perturbations and the nanoscale spatial resolution needed to resolve the evolution of the super current flow. To date, JJ characterization has primarily relied on indirect measurements such as the critical current that separates the dissipation-less (zero electrical resistance) and resistive states. However, this only provides insight into the resistive state while the ground state below the critical current stays hidden.

Here we quantitatively visualize the current flow in a JJ device with nanoscale resolution. The spatial distribution of Josephson current flow can be modulated by varying the SC phase difference between two sides of the junction. In any JJ, the SC phase difference is governed by three factors: (i) external magnetic field; (ii) external bias current; (iii) self-field or SC phase gradient induced by the finite Josephson current density. Our measurements reveal the evolution of Josephson current flow with all three factors, including features associated with the change of the number of current loops at the junction known as the Josephson vortex (JV). In particular, factors (i) and (ii) can affect (iii), altering the super current flow even without detectable transport features. We find two previously unidentified effects of the Josephson current-induced phase from factor (iii). First, hidden ground states with different numbers of JVs are found within the zero-resistance state, which can be electrically switched below the critical current. Second, a new mechanism for the Josephson diode effect is established based on the second harmonic phase terms induced by the Josephson current when time-reversal and inversion symmetry are broken.

The measurement setup is shown in Fig. 1a. We employ a diamond tip containing a single nitrogen vacancy (NV) center to map the local magnetic field generated by the current flow28. The results are obtained from two devices with junction width W = 0.15 and 0.2 μm, length L = 1.5 μm and thickness t = 35 nm. The SC electrodes are measured to be in the thin-film limit L ≪ λp, where λp is the Pearl length (Supplementary Fig. 1). This suggests the factor (iii) contribution in our device comes from the Josephson current-induced phase associated with the kinetic inductance of the SC film, instead of the self-field effect. The junctions are diffusive (electron mean free path lmfp < W) and over-damped (no hysteresis during bias current sweeps). Throughout the paper, we refer to the transverse (longitudinal) direction as x(y), and the direction perpendicular to the plane as z. The origin x = y = 0 is set to the center of JJ.Fig. 1 Measurement setup and expected Josephson current flow.

a Schematics showing SC-normal-SC junction measured by scanning NV center embedded in a diamond tip. The SC wave function can be described by an amplitude and phase Ψ = ∣Ψ∣eiθ. Under external magnetic field Bz, the screening current near the JJ (red lines) induces a phase difference ϕe(x). The bias current causes a phase difference between the SC electrodes ϕbias. b Measured differential resistance dV/dI versus perpendicular magnetic field Bz and bias current Idc, at T = 7 K. Dashed lines are the expected critical current (see Supplementary Note 2), where red is ϕbias = π/2, blue is ϕbias = −π/2. The Ic nodes are denoted as  ±Bn. c, d Calculated Josephson current normalized by critical current density for 0- and 1-JV states, at external Bz = 1.10 mT in (c), and Bz = 1.91 mT in (d). The current flow is sine-like at zero bias current (bold lines), shifts along x direction at a finite bias current, and becomes cosine-like at the critical current. e Simulations showing the Josephson current flow (top) and local SC phase (bottom) of the 0- and 1-JV states. The screening current near the junction Jx (red arrows) is reduced by the Josephson current Jy (cyan arrows) in the 0-JV state, and enhanced in the 1-JV state. This causes the Josephson current-induced phase. Ibias = 0, Bz ≈ 1.2 mT in this simulation. f, Simulated ϕe(x) for 0- and 1-JV states, at the same Bz as (e). ϕe(x) is the difference of θ taken along the two dashed lines in each sub-panel of (e). g NV control and current bias sequence based on the ac magnetometry protocol. The X(Y) microwave (MW) pulses rotate the qubit around the X(Y) axis by π2 or π. NV qubit is put on the equator of the Bloch sphere and rotated by the magnetic field generated by current flow. Pulses of different bias current are synced with the MW pulses such that the final signal is the difference between the two current flow patterns.

Results

current-induced phase

The Josephson current density can be modeled by the sinusoidal current-phase relation221 Jy(x)=Jcsin[ϕ(x)],

where Jc is the Josephson critical current density (assumed constant for now), and ϕ(x) is the phase difference across the JJ at position x. ϕ(x) = ϕe(x) + ϕbias, where ϕe(x) arises from the external magnetic field Bz (factor i), and ϕbias is the additional phase difference due to the injected bias current (factor ii). The strength of Josephson current-induced phase (factor iii) is regulated by Jc. For small Jc, the Josephson penetration length λJ≈Φ0Lt/4πμ0JcλL2≫L29, the Josephson current-induced phase can be neglected ("weak junction” limit). Φ0 is the flux quantum, λL is the London penetration length.

In the weak-junction limit, external Bz controls the number of JV. The transport critical current Ic oscillates and reaches zero at nodes Bz = ±Bn (n is integer). It is known as the “Fraunhofer map”30,31. In each lobe where Bn−1 < ∣Bz∣ < Bn, there are n JVs at the junction; in the central lobe there is 0-JV (Fig. 1b); the only way to change the number of JV is by sweeping Bz through the nodes Bn (Supplementary Fig. 2). In weak junctions the first term ϕe(x) is induced by the screening current Jx in the SC electrodes, and scaled by ϕe0 ≡ ϕe∣x=L/2 ≈ 1.7BzL2/Φ031 (Supplementary Eqn. 9). The second term ϕbias changes from  −π/2 + nπ to π/2 + nπ when Ibias sweeps from  −∣Ic∣ to  +∣Ic∣ (Fig. 1c, d), which can be viewed as moving the JV along the xx direction.

In “strong junctions” (λJ ≪ L), ϕ(x) is altered by the Josephson current-induced phase and lacks analytical solutions. Qualitatively, the screening current Jx deviates from the weak-junction limit by an amount proportional to the Josephson current Jy (Fig. 1e), due to the continuity of current. This leads to additional phase gradient ∂θ/∂x∝LkJx, which changes ϕe(x) = θ(x)∣y = W/2 − θ(x)∣y = −W/2. Here Lk is the kinetic inductance, which is inversely proportional to the superfluid stiffness. θ is the SC phase. Specifically, the larger Jx in the 1-JV state leads to enhanced ϕe0, compared to the 0-JV state at the same external magnetic field (Fig. 1f). Ibias can further change the Josephson current and its induced phase. Thus in strong junctions, the total ϕ(x) and current flow need to be solved self-consistently. We find that our device is close to the weak-junction limit at T = 7 K, but is in an intermediate range at T = 4 K.

Visualizing Josephson current flow

To optimize magnetic field sensitivity, we utilize the “ac magnetometry" protocol, synchronizing NV control pulses with the signal (Fig. 1g). In this protocol, different bias currents I1 and I2 is applied to the junction during two halves of each cycle. The magnetic field generated by Ik rotates the prepared NV spin superposition state along the equator of the Bloch sphere by an angle φk = 2πγebnvτ/2, where γe is the gyromagnetic ratio of electron spin, bnv is the magnetic field projected along NV axis, and τ is cycle duration. After the sequence, the accumulated angle is φ1 − φ2 so each measurement records the difference between two selected scenarios. bnv is converted to vector components of the magnetic field bx,y,z using a Fourier method (see Methods). The current vector (jx, jy) is then reconstructed with the Fourier method using in-plane components of the magnetic field. We find similar results with regularization and machine learning methods (see Supplementary Note 5).

We start at T = 7 K where the transport result suggests a weak junction (Fig. 1b). To visualize the evolution of Josephson current, we use two types of sequences. In the first sequence, we highlight the effect of ϕbias by taking the difference between finite and zero Ibias (schematics in Fig. 2a). The expected current profiles jy(x) are shown in Fig. 2a, by subtracting the relevant lines in Fig. 1c. The sign of Ibias determines the direction of the profile shift and the amplitude determines the amount of the shift. Fig. 2b, c show measurements using Ibias ≈ ±∣Ic∣ in the sequence while no JV is in the junction. As expected, current features are seen at the opposite side of the junction for  ±Ibias. jy(x) at the junction also shows the lateral movement when sweeping Ibias (Supplementary Fig. 3). ϕbias can be acquired by fitting jy(x) to the calculated profiles from Fig. 2a (Supplementary Eqn. 10). The result agrees with the sinusoidal current-phase relation (Fig. 2d).Fig. 2 Visualizing Josephson current response to bias current and magnetic field.

a Top sketch shows the measurement sequence that takes the difference between zero and finite current bias. Expected current flow show changing line shapes at various bias current Ibias, with the opposite bias current leading to inverted current flow pattern around the center of the junction. b, c Color map shows z-component of the current-generated magnetic field bz. The arrows show the reconstructed current flow vector. Results are measured using the sequence shown in (a). External magnetic field is Bz = 0.95 mT. The SC electrodes are marked by solid lines and normal metal part is marked by dashed lines. d Phase difference caused by bias current ϕbias extracted from the current flow profile at the junction. The result agrees with the sinusoidal current-phase relation. Inset shows the measured jy(x) at Ibias = Ic, and the fitting results using Supplementary Eqn. 10. The data in the gray area is excluded in the fitting. e Top sketch shows the measurement sequence that takes the difference between symmetric positive and negative Ibias. The expected signals change sign when measuring the 0- and 1-JV states at their respective ∣Ic∣. All signals show cosine-like shape, with amplitude growing with Ibias. f, g bz and current flow vector maps measured using the sequence shown in e). The Josephson current changes sign between 0- and 1-JV states. External magnetic field is Bz = 1.10 mT in (f), and Bz = 1.91 mT in (g). h Effective phase difference across the junction ϕeff = ϕe∣x=W/2 deviates from the external field contribution ϕext (indicated by the shaded area), as a result of induced phase from the Josephson current. Scale bar is 0.5 μm for (b, c, f, g). The measurements are taken at T = 7 K. Error bars in d and h represent standard deviation from fitting.

Next we show the effect of magnetic flux on the Josephson current by taking the difference between  ±Ibias (schematics in Fig. 2e). Here the expected signals are cosine-like and only grow in amplitude with Ibias (Fig. 2e). For the same Ibias direction, the signal flips sign for 0- and 1-JV states (switching from red to blue branch in Fig. 1b). Measurements using this sequence are shown in Fig. 2f, g, where we use Ibias ≈ ∣Ic∣ in both cases. Results measured at Ibias ≈ 0.5∣Ic∣ show the same shape with half the amplitude (Supplementary Fig. 5a, b). We repeat the measurement at various magnetic fields, and the current profile reversal can be seen when external Bz crosses the node B0 ≈ 1.5 mT, i.e., when JV number changes by 1 (Supplementary Fig. 4). Notably, the current flow at x = 0 is parallel (anti-parallel) to Ibias when the junction contains even (odd) number of JVs. In the 2-JV state, the current flow at x = 0 and bias current are both positive (Supplementary Fig. 6), like the 0-JV state (Fig. 2f).

Our measurements provide quantitative details of the current flow compared to previous methods24,32–34. First, we confirm the SC is in the thin-film limit (L ≪ λp). The absolute value of λp ≈ 13 μm at T = 7 K is directly measured from the stray field (Supplementary Fig. 1). In comparison, indirect measurements of λp range from 1 to 5 μm (see Methods), and thus are unable to determine whether the SC is in the thin-film regime. Second, the JV extends into the SC electrode on both sides by ~350 nm (Supplementary Fig. 5c), consistent with the effective area expected from the nodes Bn. However, the measured jy(x) profile does not match the expectation under external Bz. For example, at Bz = 1.1 mT the jy(x) in Fig. 2f is expected to be negative at x = ±L/2, but stays positive in the experiment, suggesting a smaller than expected ϕe0 (Supplementary Fig. 4f).

We introduce an effective phase difference ϕeff as a fitting parameter for the measured jy(x), replacing the theoretically predicted ϕe0 (Supplementary Eqn. 11). Fig. 2h shows ϕeff is lower (higher) than the phase induced by the external field ϕext when Bz < B0 (Bz > B0). The discrepancy is a direct consequence of the Josephson current-induced phase, as expected in strong junctions. While we find λJ is comparable to L when calculated using experimentally measured parameters, no strong junction feature is observed in the “Fraunhofer map" (Fig. 1b). The Josephson current-induced phase only causes small changes to Ic, and the effect cancels out when tracking Ic over large range of Bz32. However, such an effect is still pertinent to designing SC devices such as JJ arrays because the inductance of each junction is affected by the supercurrent flow35.

Electrically switching JV below ∣Ic∣

The “Fraunhofer map” changes when measured at T = 4 K. The nodes of ∣Ic∣ at Bn are lifted although sharp kinks remain (Fig. 3a). This could be caused by a combination of reasons, including an asymmetric critical current density, Jc(x) ≠ Jc(−x) and the strong junction effect due to the increased ∣Ic∣33,36,37. However, the precise mechanism remains difficult to dissect owing to the zero-resistance below Ic. Here the origin of non-zero local minima of ∣Ic∣ is revealed by mapping the current flow. A differential magnetic field b~nv is measured using a sequence that senses the small ac bias (ĩac) response while fixing the dc bias (Idc), shown by the schematic drawing above Fig. 3b. The b~nv is measured around the kinks of ∣Ic∣ while the NV is fixed over the center of the junction (Fig. 3b main panel). Abrupt changes of b~nv at large ∣Idc∣ match the transport Ic (circles in Fig. 3b). This suggests that the junction is minimally perturbed by the NV magnetometer.Fig. 3 Electric control of JV ground state below critical current.

a Differential resistance dV/dI measured in the same device as Fig. 1c but at T = 4 K, showing order-of-magnitude increase and oscillation of Ic that does not reach zero. The dashed box shows the measurement range in (b). b Top sketch shows measurement sequence. Main panel shows the differential magnetic field projected along NV axis b~nv, generated by the current flow response to the small ac current ĩac ≈ 0.8 μA. b~nv is shown versus dc bias current Idc and perpendicular magnetic field Bz. The NV is fixed at (x, y) ≈ (−500, 0) nm, where x = y = 0 is the center of JJ. Circles are critical current extracted from transport result in (a). c Difference of Gibbs free energy between the 0- and 1-JV states below Ic, Δε = ε0V − ε1V, from TDGL simulations. Blue (red) circles show simulated critical current of state with 0- (1-) JV. The color in the non-overlap region is saturated to indicate only 0- and 1-JV state is present. d, e Spatial maps of differential b~z and current flow vector measured with sequence shown in (b), showing both (d) 0-JV, and (e) 1-JV states at the same external field Bz = 1.4 mT. The features in e are positioned asymmetrically because the JV is at right-of-center of the junction due to finite Idc. f Full vortex profile is observed when using a sequence that takes the difference between 0- and 1-JV states.

An additional sharp boundary of b~nv below Ic separates the 0- and 1-JV states. This is confirmed by the spatial maps in Fig. 3d, e. In Fig. 3d (0-JV state), the current is almost uniform despite being measured at higher external Bz than Fig. 2f, suggesting ϕe0 is strongly suppressed by the Josephson current-induced phase. Intriguingly, the 0- and 1-JV states can be reached at the same Bz but different Ibias (Fig. 3d, e). Measuring the difference between two such Ibias, the current profile that corresponds to the JV number changing by 1 is observed (Fig. 3f). This demonstrates precise JV number control using pure electric means while staying below Ic, which could be useful for low-dissipation memory and logic devices based on SC hybrid structures.

The phase boundary below Ic supports Josephson current-induced phase as the primary reason for the node-lifting in our device. The induced phase enables the co-stability of the 0- and 1-JV states originating from the existence of two local minima of the energy as a function of the order parameter Ψ(r) at the same external magnetic field, as shown by the time-dependent Ginzburg-Landau (TDGL) simulation (Fig. 3c, for details see Supplementary Note 8). The overlapping states with different number of JVs were predicted to arise from the self-field effect previously38–42. However, we do not expect self-field to be the main effect here. The measured self-field of the current is insignificant in our device (less than 5% of the external field), and the SC is still in the thin-film limit at T = 4 K (L ≪ λp, see Methods). Furthermore, the boundary of the 0- and 1-JV phase diagram only extends from the 1-JV region (Fig. 3b) in the experiment. This is independent of sweeping directions of Bz or Idc (Supplementary Fig. 9a), and similar behavior is observed between the 1- and 2-JV states (Supplementary Fig. 9b). The lack of hysteresis is quite unexpected. One possibility is that the JJ relaxes to the ground state with lower energy due to the elevated temperature and small perturbations of the measurement, although the detailed mechanism is an open question for future work. The simulated Gibbs free energy difference Δε = ε0V−ε1V shows the 1-JV state has lower energy than the 0-JV state in most, but not all of the overlap region (Fig. 3c). In fact, including the self-field effect energetically favors the 0-JV over the 1-JV state, further deviating from the experimental results (Supplementary Fig. 20).

Inversion asymmetry and Josephson diode effect

The transport Ic is non-reciprocal when Ibias is applied in opposite directions at T = 4 K, i.e. ∣Ic+∣≠∣Ic−∣ (Fig. 4a). This is referred to as the “Josephson diode effect”43. The asymmetry parameter, η=∣Ic+∣−∣Ic−∣∣Ic+∣+∣Ic−∣, exceeds 10% in our device. We identify a new mechanism for the diode effect comprising three ingredients, (i) time-reversal symmetry breaking (by Bz), (ii) inversion symmetry breaking, and (iii) Josephson current-induced phase. The first two conditions are required by symmetry17,44, while the third provides a mechanism whereby the Josephson current is not a simple sinusoidal function of ϕbias. As a result, the critical current density is reached on opposite sides of the junction at  ±Ibias, which leads to asymmetric Ic± (Fig. 4b). Interestingly, it can be shown theoretically that first two ingredients alone are not sufficient to generate the diode effect Supplementary Note 7. Combining the symmetry breaking with the Josephson current-induced phase introduces higher harmonic terms with a phase offset in the current-phase relation Supplementary Note 7. Our results confirm all three components are necessary. For example, the current profile measured at T = 7 K is asymmetric, jy(x) ≠ jy(−x), suggesting broken inversion symmetry (Supplementary Fig. 16). However, the weaker current-induced phase due to smaller Jc is insufficient to generate a non-reciprocal global critical current response (Supplementary Note 6).Fig. 4 A new mechanism for the Josephson diode effect.

a Left axis, forward and backward critical current extracted from results measured in Fig. 3a. Right axis, asymmetry parameter η showing JDE when time-reversal symmetry is broken by Bz. b Schematics showing inversion symmetry breaking (non-uniform Jc) at the junction can lead to different forward and backward critical current. c, d Differential magnetic field b~z and current flow vector measured at symmetric  ±Ibias show inversion asymmetric patterns. Measurements are taken at Bz = 0.5 mT, using the same sequence as described in Fig. 3b. e, f Simulated results corresponding to (c) and (e), when inversion symmetry at the junction is broken, modeled as non-uniform junction width W1 > W2. The measurements in (a, c, d) are taken at T = 4 K.

The current flow measurement directly reveals the broken inversion symmetry even when the global Ic± is almost symmetric. The b~nv map at  ±Ibias is measured at Bz = 0.5 mT. Although η is only about 2%, the current flow pattern is clearly asymmetric for  ±Ibias; a loop appears near the left edge for  −Ibias (Fig. 4c) but not for  + Ibias (Fig. 4d). We model the non-uniform Jc(x) with an uneven junction width W (W1 > W2) in the TDGL simulation, confirming the role of broken inversion symmetry (Fig. 4e, f); if the junction is inversion symmetric, the ac current flow for  ±Ibias should be mirrored along the x direction (Supplementary Fig. 21f–h). In fact, the broken inversion symmetry is also responsible for the skewed phase boundary for  ±Ibias in Fig. 3b, c. In reality, the non-uniform Jc(x) could be due to variations in junction width, SC/normal barrier transparency, or normal metal resistivity.

Discussion

The JV discussed in our work should be distinguished from the Abrikosov vortex in type-II superconductors45. While both move in the same direction with Ibias and exhibit normal cores, observed in spectroscopic studies18,46, only the JV configuration can be controlled by a small change of Ibias below Ic. Even when the current-induced phase is weak, the JV can be precisely moved side-to-side by the small change of Ibias from  − ∣Ic∣ to ∣Ic∣ at Bz = Bn ± ε, (ε ≪ 1). In particular, ∣Ic∣ should vanish at Bz = Bn, if the junction is symmetric about its midpoint. This control over the JV position enables us to observe the large alternating magnetic field signal at the JJ with minimal changes in Ibias (Fig. 2g). Finally, the minimal energy cost associated with JV movement (IcΔϕ) as Ic → 0 supports JV control as an energy-efficient way of communication between qubits47.

The new mechanism of the Josephson diode effect offers a blueprint to realize a scalable SC rectifier with any thin-film SC. Conventional Josephson diodes that are driven by the self-field effect require a large operating current because the geometric inductance is usually small, especially at the nanoscale48,49. However, the kinetic inductance can dominate in SC with small superfluid stiffness (e.g., low superfluid density), making it possible to reduce the device size. This also enables electric tuning of the Josephson diode by injecting a small current Jx to control the Josephson current-induced phase.

Finally, spatial mapping of the current flow J(x, y) presents an alternative observable to electrical transport in SC hybrid structures. By accurately measuring J(x, y) with high sensitivity and spatial resolution, we pinpoint the origin of the Josephson diode effect in our device, which is otherwise hidden. Our approach opens up further avenues to unseal the mechanisms for the non-reciprocity in a broad range of SC systems, and symmetry breaking in gate-tunable superconductors based on van der Waals and moiré materials. The measured current flow could be directly compared with simulations based on TDGL or quantum transport to diagnose SC circuits, such as the local transparency of the JJ barrier.

Methods

Variable temperature scanning setup

Measurements were performed in a home-built variable temperature system with optical access. There are multiple nano-pillars containing NV centers on each diamond probe, and a goniometer with both pitch and yawn control is used to set the stand-off distance between the NV and the sample28, which ranges between 130 to 180 nm throughout the study. The NV center is excited with 532 nm green laser (Coherent Sapphire), and read out with standard optical detected magnetic resonance (ODMR) technique using a 600 nm long-pass optical filter. The time-averaged power of the green laser pulses is less than 50 μW. The microwave (MW) drive is provided via on-chip transmission line next to the sample. MW is sourced from SGS-100A (Rohde & Schwarz) and modulated with the built-in IQ mixer. MW pulses are then amplified by +40 db using 30S1G6C (AR Inc), and routed through another switch (RF lambda) to reduce noise from the amplifier. MW and bias current control sequences are generated by arbitrary wave generator AWG5014C (Tektronic).

Device fabrication and electrical characterization

The SC and normal parts of the JJ are made of niobium nitride (NbN) and gold (Au) thin films, respectively. The JJs are fabricated on undoped Silicon substrate with 285 nm SiO2 on top. Standard electron beam lithography method is used to define the device geometry using double-layer e-beam resist. The normal part of the junction is first formed with thermal evaporation (2 nm Ti/ 35 nm Au). A short Argon milling process is used right before sputtering SC electrodes (2 nm Ti/6 nm Nb/30 nm NbN). The MW strip line is formed with 2 nm Ti/ 60 nm Au. Four terminal resistance result was first measured with dc bias from Keithley 2400 and dc voltage with Keithley 2100, and then taken numerical derivative to acquire differential resistance shown in the main text.

Probe fabrication and typical characteristics

The diamond fabrication process follows ref. 50. Specifically, ultra pure diamond with natural 13C abundance and [100] facet (electronic grade from Element Six) is diced into thin slabs ~50 μm thick. One side of the slab is etched by Argon/Chloride plasma to relieve surface strain, then implanted with 15N ions at a dose of 5 × 1010/cm2 and acceleration energy of 6 keV (Innovion). Then the diamond is annealed in ultrahigh vacuum (<3 × 10−8 Torr) at 800 ∘C for 24 hours to form NV centers. The diamond nano-pillars are defined with standard e-beam lithography and etched with O2 plasma. On average we get 1 NV center per diamond pillar with this process. The NV depth from the surface is about 15-20 nm. Typical ODMR red photon count is 100 k/s, contrast in pulsed measurement is 20 to 30%, and the coherence time is T2*≈1μ s and T2 ≈ 30 μs at 4 K and the small magnetic field used in this study.

Detail about NV magnetometry

NV is a spin-1 system with low energy states s=0,±1. The 0 is split in energy from ±1 by the zero field splitting (2.87 GHz) and ±1 are further split by the Zeeman energy EZ = gμBBnv, here g = 2 is the Landé g-factor for electron, μB is the Bohr magneton, and Bnv is the magnetic field along NV axis. In practice, we apply a external field of less than 50 G along the NV axis, and drive the 0 and −1 states as a qubit using MW. As mentioned in the main text, “ac” magnetometry is used to filter out low frequency noise and maximize sensitivity by utilizing the longer T2 coherence time. Specifically, the NV qubit is first prepared in the 0 state using a green laser pulse, and then driven into the superposition state 12(0+i−1) by a Xπ2 MW pulse. We use two types of dynamic decoupling sequences, the spin echo (Hahn echo) with one π-pulse, and the Carr-Purcell-Meiboom-Gill (CPMG) with n-Yπ pulses in the experiment51,52. Between neighboring π-pulses, the qubit rotates by an angle φ = 2πγebnvτn, where γe = 28 GHz/T is the gyromagnetic ratio of the electron spin, bnv is the magnetic field generated by the current projected along NV axis, and τn is the evolution time between neighboring MW pulses. The π-pulses reverse the qubit rotation direction, and the total angle is the difference of the accumulation in each half of the sequence. The frequency of NV control sequence and bias current modulation is f = 100 to 500 kHz, corresponding to <1 nA bias current due to the AC Josephson effect I = hf/2eRN22 (h is Planck’s constant, e is electron charge, RN ≈ 1Ω is the normal state resistance of the JJ). This is 3 to 4 orders of magnitude smaller than the bias current applied to the JJ.

To extract the phase accumulation angle, the NV spin is projected to the 0 and −1 states using four π2-pulses π2±X/Y and record the ODMR signal. The angle is then calculated from2 φ=arctanCπ2X−Cπ2−XCπ2Y−Cπ2−Y

where Cπ2±X/Y are the photon counts from π2±X/Y projections. The measurement sequences are averaged up to 100k times (about 10 seconds) at each point to extract the Bnv.

Reconstructing current flow from magnetic field

In this section, we discuss the three methods used to reconstruct current flow jx,y from bnv. For all methods, we first convert the magnetic field projected along NV axis bnv, to Cartesian vector magnetic field bx,y,z using the source-free constraint for the stray field53–55,3 ∇×b=0

Thus in the Fourier space the vector components are,4 bz(k)=bnv(k)unv⋅ubx(k)=−ikxkbz(k)by(k)=−ikykbz(k)

here k = (kx, ky) is the 2D wavevector, k=kx2+ky2, unv is the unit vector of the NV axis, u = (−ikx/k, −iky/k, 1). The singularity point at k = 0 is discarded during the reconstruction. Because the SC electrode is much longer than our measurement window in the y direction, the jy outside the window on the top and bottom sides also contribute to the measured bnv. In practice, we extend the measured bnv with the top and bottom lines in the y direction, and linearly extrapolates bnv to zero in the x direction. The padding size in each direction is 10 times of the measurement window, at which point increasing the size does not change the reconstruction result. More detail of the reconstruction process is discussed in Supplementary Note 5.

Supplementary information

Supplementary Information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52271-z.

Acknowledgements

We thank Y. Xie, J. Cremer, A. Hamo, T. Werkmeister for inspiring discussions. A.Y. acknowledges support from the Army Research Office under Grant numbers W911NF-22-1-0248 and W911NF-21-2-0147, the Gordon and Betty Moore Foundation through Grant GBMF 9468, and the Quantum Science Center (QSC), a National Quantum Information Science Research Center of the U.S. Department of Energy (DOE). S.C. and S.P. acknowledge partial support from the Harvard Quantum Initiative in Science and Engineering. A.S. acknowledges support by the European Union’s Horizon 2020 research and innovation program (Grant Agreement LEGOTOP No. 788715), the German Research Foundation DFG (CRC/Transregio 183, EI 519/7-1), and the Israel Science Foundation Quantum Science and Technology (2074/19). P.M. acknowledges financial support through SNSF project No. 188521 and through ERC consolidator grant QS2DM. B.I.H. acknowledges support from NSF grant DMR 1231319.

Author contributions

S.C., S.P., U.V., and A.Y. conceived and designed the experiments; S.C. and S.P. prepared the devices, performed the electrical transport and magnetometry measurements, analyzed and visualized the data with input from U.V., N.M., and A.Y.; S.C., and S.P. carried out the simulation and analysis with input from A.S., B.I.H., and A.Y.; D.A.B., M.F., and P.M. carried out the current reconstruction using machine learning method; T.Z. fabricated the diamond probes; S.C. and S.P. wrote the manuscript with input from U.V., A.S., B.I.H., A.Y., and contributions from all co-authors.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

All experimental and numerically simulated data included in this work are available in the Zenodo database56.

Competing interests

A.Y., S.C., E.P., U.V., N.M., A.S., and B.I.H. have applied for a patent partially based on this work. The other authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Shaowen Chen, Seunghyun Park.
==== Refs
References

1. Josephson BD Possible new effects in superconductive tunnelling Phys. Lett. 1962 1 251 253 10.1016/0031-9163(62)91369-0
Josephson, B. D. Possible new effects in superconductive tunnelling. Phys. Lett. 1, 251–253 (1962).10.1016/0031-9163(62)91369-0
2. de Gennes P Boundary effects in superconductors Rev. Mod. Phys. 1964 36 225 10.1103/RevModPhys.36.225
de Gennes, P. Boundary effects in superconductors. Rev. Mod. Phys. 36, 225 (1964).10.1103/RevModPhys.36.225
3. Josephson B Coupled superconductors Rev. Mod. Phys. 1964 36 216 10.1103/RevModPhys.36.216
Josephson, B. Coupled superconductors. Rev. Mod. Phys. 36, 216 (1964).10.1103/RevModPhys.36.216
4. Devoret MH Schoelkopf RJ Superconducting circuits for quantum information: an outlook Science 2013 339 1169 1174 10.1126/science.1231930 23471399
Devoret, M. H. & Schoelkopf, R. J. Superconducting circuits for quantum information: an outlook. Science 339, 1169–1174 (2013).23471399 10.1126/science.1231930
5. Fedorov KG Shcherbakova AV Wolf MJ Beckmann D Ustinov AV Fluxon readout of a superconducting qubit Phys. Rev. Lett. 2014 112 160502 10.1103/PhysRevLett.112.160502 24815629
Fedorov, K. G., Shcherbakova, A. V., Wolf, M. J., Beckmann, D. & Ustinov, A. V. Fluxon readout of a superconducting qubit. Phys. Rev. Lett. 112, 160502 (2014).24815629 10.1103/PhysRevLett.112.160502
6. Kjaergaard M Superconducting qubits: current state of play Annu. Rev. Cond. Matter Phys. 2020 11 369 395 10.1146/annurev-conmatphys-031119-050605
Kjaergaard, M. et al. Superconducting qubits: current state of play. Annu. Rev. Cond. Matter Phys. 11, 369–395 (2020).10.1146/annurev-conmatphys-031119-050605
7. Wang H Coherent terahertz emission of intrinsic Josephson junction stacks in the hot spot regime Phys. Rev. Lett. 2010 105 057002 10.1103/PhysRevLett.105.057002 20867948
Wang, H. et al. Coherent terahertz emission of intrinsic Josephson junction stacks in the hot spot regime. Phys. Rev. Lett. 105, 057002 (2010).20867948 10.1103/PhysRevLett.105.057002
8. Welp U Kadowaki K Kleiner R Superconducting emitters of THz radiation Nat. Photon. 2013 7 702 710 10.1038/nphoton.2013.216
Welp, U., Kadowaki, K. & Kleiner, R. Superconducting emitters of THz radiation. Nat. Photon. 7, 702–710 (2013).10.1038/nphoton.2013.216
9. Kirtley JR Wikswo Jr JP Scanning squid microscopy Annu. Rev. Mater Sci. 1999 29 117 148 10.1146/annurev.matsci.29.1.117
Kirtley, J. R. & Wikswo Jr, J. P. Scanning squid microscopy. Annu. Rev. Mater Sci. 29, 117–148 (1999).10.1146/annurev.matsci.29.1.117
10. Kirtley J Fundamental studies of superconductors using scanning magnetic imaging Rep. Prog. Phys. 2010 73 126501 10.1088/0034-4885/73/12/126501
Kirtley, J. Fundamental studies of superconductors using scanning magnetic imaging. Rep. Prog. Phys. 73, 126501 (2010).10.1088/0034-4885/73/12/126501
11. Holmes DS Ripple AL Manheimer MA Energy-efficient superconducting computing - power budgets and requirements IEEE Trans. Appl. Supercond. 2013 23 1701610 10.1109/TASC.2013.2244634
Holmes, D. S., Ripple, A. L. & Manheimer, M. A. Energy-efficient superconducting computing - power budgets and requirements. IEEE Trans. Appl. Supercond. 23, 1701610 (2013).10.1109/TASC.2013.2244634
12. Grosfeld E Stern A Observing Majorana bound states of Josephson vortices in topological superconductors PNAS 2011 108 11810 11814 10.1073/pnas.1101469108 21730165
Grosfeld, E. & Stern, A. Observing Majorana bound states of Josephson vortices in topological superconductors. PNAS 108, 11810–11814 (2011).21730165 10.1073/pnas.1101469108
13. Beenakker C Search for Majorana fermions in superconductors Annu. Rev. Condens. Matter Phys. 2013 4 113 136 10.1146/annurev-conmatphys-030212-184337
Beenakker, C. Search for Majorana fermions in superconductors. Annu. Rev. Condens. Matter Phys. 4, 113–136 (2013).10.1146/annurev-conmatphys-030212-184337
14. Lutchyn RM Majorana zero modes in superconductor–semiconductor heterostructures Nat. Rev. Mater. 2018 3 52 68 10.1038/s41578-018-0003-1
Lutchyn, R. M. et al. Majorana zero modes in superconductor–semiconductor heterostructures. Nat. Rev. Mater. 3, 52–68 (2018).10.1038/s41578-018-0003-1
15. Pal B Josephson diode effect from cooper pair momentum in a topological semimetal Nat. Phys. 2022 18 1228 1233 10.1038/s41567-022-01699-5 36217362
Pal, B. et al. Josephson diode effect from cooper pair momentum in a topological semimetal. Nat. Phys. 18, 1228–1233 (2022).36217362 10.1038/s41567-022-01699-5
16. Jeon K-R Zero-field polarity-reversible Josephson supercurrent diodes enabled by a proximity-magnetized pt barrier Nat. Mater. 2022 21 1008 1013 10.1038/s41563-022-01300-7 35798947
Jeon, K.-R. et al. Zero-field polarity-reversible Josephson supercurrent diodes enabled by a proximity-magnetized pt barrier. Nat. Mater. 21, 1008–1013 (2022).35798947 10.1038/s41563-022-01300-7
17. Nadeem, M., Fuhrer, M. S. & Wang, X. The superconducting diode effect. Nat. Rev. Phys. 5, 558–577 (2023).
18. Roditchev D Direct observation of Josephson vortex cores Nat. Phys. 2015 11 332 337 10.1038/nphys3240
Roditchev, D. et al. Direct observation of Josephson vortex cores. Nat. Phys. 11, 332–337 (2015).10.1038/nphys3240
19. Tsuei C Pairing symmetry and flux quantization in a tricrystal superconducting ring of YBa2Cu3O7−δ Phys. Rev. Lett. 1994 73 593 10.1103/PhysRevLett.73.593 10057486
Tsuei, C. et al. Pairing symmetry and flux quantization in a tricrystal superconducting ring of YBa2Cu3O7−δ. Phys. Rev. Lett. 73, 593 (1994).10057486 10.1103/PhysRevLett.73.593
20. Hilgenkamp H Ordering and manipulation of the magnetic moments in large-scale superconducting π-loop arrays Nature 2003 422 50 53 10.1038/nature01442 12621428
Hilgenkamp, H. et al. Ordering and manipulation of the magnetic moments in large-scale superconducting π-loop arrays. Nature 422, 50–53 (2003).12621428 10.1038/nature01442
21. Gutfreund A Direct observation of a superconducting vortex diode Nat. Commun. 2023 14 1630 10.1038/s41467-023-37294-2 36959184
Gutfreund, A. et al. Direct observation of a superconducting vortex diode. Nat. Commun. 14, 1630 (2023).36959184 10.1038/s41467-023-37294-2
22. Tinkham, M. Introduction to superconductivity (McGraw-Hill, 1996).
23. Gross R Koelle D Low temperature scanning electron microscopy of superconducting thin films and Josephson junctions Rep. Prog. Phys. 1994 57 651 10.1088/0034-4885/57/7/001
Gross, R. & Koelle, D. Low temperature scanning electron microscopy of superconducting thin films and Josephson junctions. Rep. Prog. Phys. 57, 651 (1994).10.1088/0034-4885/57/7/001
24. Dremov VV Local Josephson vortex generation and manipulation with a magnetic force microscope Nat. Commun. 2019 10 4009 10.1038/s41467-019-11924-0 31488813
Dremov, V. V. et al. Local Josephson vortex generation and manipulation with a magnetic force microscope. Nat. Commun. 10, 4009 (2019).31488813 10.1038/s41467-019-11924-0
25. Hovhannisyan RA Grebenchuk SY Baranov DS Roditchev D Stolyarov VS Lateral Josephson junctions as sensors for magnetic microscopy at nanoscale J. Phys. Chem., Lett. 2021 12 12196 12201 10.1021/acs.jpclett.1c03556 34918928
Hovhannisyan, R. A., Grebenchuk, S. Y., Baranov, D. S., Roditchev, D. & Stolyarov, V. S. Lateral Josephson junctions as sensors for magnetic microscopy at nanoscale. J. Phys. Chem., Lett. 12, 12196–12201 (2021).34918928 10.1021/acs.jpclett.1c03556
26. Grebenchuk SY Observation of interacting Josephson vortex chains by magnetic force microscopy Phys. Rev. Res. 2020 2 023105 10.1103/PhysRevResearch.2.023105
Grebenchuk, S. Y. et al. Observation of interacting Josephson vortex chains by magnetic force microscopy. Phys. Rev. Res. 2, 023105 (2020).10.1103/PhysRevResearch.2.023105
27. Stolyarov VS Revealing Josephson vortex dynamics in proximity junctions below critical current Nano Lett. 2022 22 5715 5722 10.1021/acs.nanolett.2c00647 35820103
Stolyarov, V. S. et al. Revealing Josephson vortex dynamics in proximity junctions below critical current. Nano Lett. 22, 5715–5722 (2022).35820103 10.1021/acs.nanolett.2c00647
28. Vool U Imaging phonon-mediated hydrodynamic flow in WTe2 Nat. Phys. 2021 17 1216 1220 10.1038/s41567-021-01341-w
Vool, U. et al. Imaging phonon-mediated hydrodynamic flow in WTe2. Nat. Phys. 17, 1216–1220 (2021).10.1038/s41567-021-01341-w
29. Tolpygo SK Gurvitch M Critical currents and Josephson penetration depth in planar thin-film high-tc Josephson junctions Appl. Phys. Lett. 1996 69 3914 3916 10.1063/1.117568
Tolpygo, S. K. & Gurvitch, M. Critical currents and Josephson penetration depth in planar thin-film high-t c Josephson junctions. Appl. Phys. Lett. 69, 3914–3916 (1996).10.1063/1.117568
30. Rowell J Magnetic field dependence of the Josephson tunnel current Phys. Rev. Lett. 1963 11 200 10.1103/PhysRevLett.11.200
Rowell, J. Magnetic field dependence of the Josephson tunnel current. Phys. Rev. Lett. 11, 200 (1963).10.1103/PhysRevLett.11.200
31. Clem JR Josephson junctions in thin and narrow rectangular superconducting strips Phys. Rev. B 2010 81 144515 10.1103/PhysRevB.81.144515
Clem, J. R. Josephson junctions in thin and narrow rectangular superconducting strips. Phys. Rev. B 81, 144515 (2010).10.1103/PhysRevB.81.144515
32. Dynes R Fulton T Supercurrent density distribution in Josephson junctions Phys. Rev. B 1971 3 3015 10.1103/PhysRevB.3.3015
Dynes, R. & Fulton, T. Supercurrent density distribution in Josephson junctions. Phys. Rev. B 3, 3015 (1971).10.1103/PhysRevB.3.3015
33. Mayer B Schuster S Beck A Alff L Gross R Magnetic field dependence of the critical current in YBa2Cu3O7-δ bicrystal grain boundary junctions Appl. Phys. Lett. 1993 62 783 785 10.1063/1.108578
Mayer, B., Schuster, S., Beck, A., Alff, L. & Gross, R. Magnetic field dependence of the critical current in bicrystal grain boundary junctions. Appl. Phys. Lett. 62, 783–785 (1993).10.1063/1.108578
34. Holm J Mygind J A novel cryogenic scanning laser microscope tested on Josephson tunnel junctions Rev. Sci. Inst. 1995 66 4547 4551 10.1063/1.1145288
Holm, J. & Mygind, J. A novel cryogenic scanning laser microscope tested on Josephson tunnel junctions. Rev. Sci. Inst. 66, 4547–4551 (1995).10.1063/1.1145288
35. Kuzmin R Mehta N Grabon N Manucharyan VE Tuning the inductance of Josephson junction arrays without squids Appl. Phys. Lett. 2023 123 182602 10.1063/5.0171047
Kuzmin, R., Mehta, N., Grabon, N. & Manucharyan, V. E. Tuning the inductance of Josephson junction arrays without squids. Appl. Phys. Lett. 123, 182602 (2023).10.1063/5.0171047
36. Hilgenkamp H Mannhart J Grain boundaries in high-t c superconductors Rev. Mod. Phys. 2002 74 485 10.1103/RevModPhys.74.485
Hilgenkamp, H. & Mannhart, J. Grain boundaries in high-t c superconductors. Rev. Mod. Phys. 74, 485 (2002).10.1103/RevModPhys.74.485
37. Boris AA Evidence for nonlocal electrodynamics in planar Josephson junctions Phys. Rev. Lett. 2013 111 117002 10.1103/PhysRevLett.111.117002 24074115
Boris, A. A. et al. Evidence for nonlocal electrodynamics in planar Josephson junctions. Phys. Rev. Lett. 111, 117002 (2013).24074115 10.1103/PhysRevLett.111.117002
38. Owen C Scalapino D Vortex structure and critical currents in Josephson junctions Phys. Rev. 1967 164 538 10.1103/PhysRev.164.538
Owen, C. & Scalapino, D. Vortex structure and critical currents in Josephson junctions. Phys. Rev. 164, 538 (1967).10.1103/PhysRev.164.538
39. Pagano S Ruggiero B Sarnelli E Magnetic-field dependence of the critical current in long Josephson junctions Phys. Rev. B 1991 43 5364 10.1103/PhysRevB.43.5364
Pagano, S., Ruggiero, B. & Sarnelli, E. Magnetic-field dependence of the critical current in long Josephson junctions. Phys. Rev. B 43, 5364 (1991).10.1103/PhysRevB.43.5364
40. Kuplevakhsky S Glukhov A Static solitons of the sine-gordon equation and equilibrium vortex structure in Josephson junctions Phys. Rev. B 2006 73 024513 10.1103/PhysRevB.73.024513
Kuplevakhsky, S. & Glukhov, A. Static solitons of the sine-gordon equation and equilibrium vortex structure in Josephson junctions. Phys. Rev. B 73, 024513 (2006).10.1103/PhysRevB.73.024513
41. Kuplevakhsky S Glukhov A Exact analytical solution of the problem of current-carrying states of the Josephson junction in external magnetic fields Phys. Rev. B 2007 76 174515 10.1103/PhysRevB.76.174515
Kuplevakhsky, S. & Glukhov, A. Exact analytical solution of the problem of current-carrying states of the Josephson junction in external magnetic fields. Phys. Rev. B 76, 174515 (2007).10.1103/PhysRevB.76.174515
42. Kuplevakhsky S Glukhov A Exact analytical solution of a classical Josephson tunnel junction problem Low. Temp. Phys. 2010 36 1012 1021 10.1063/1.3521573
Kuplevakhsky, S. & Glukhov, A. Exact analytical solution of a classical Josephson tunnel junction problem. Low. Temp. Phys. 36, 1012–1021 (2010).10.1063/1.3521573
43. Moll PJ Geshkenbein VB Evolution of superconducting diodes Nat. Phys. 2023 19 1379 1380 10.1038/s41567-023-02229-7
Moll, P. J. & Geshkenbein, V. B. Evolution of superconducting diodes. Nat. Phys. 19, 1379–1380 (2023).10.1038/s41567-023-02229-7
44. Zhang Y Gu Y Li P Hu J Jiang K General theory of Josephson diodes Phys. Rev. X 2022 12 041013
Zhang, Y., Gu, Y., Li, P., Hu, J. & Jiang, K. General theory of Josephson diodes. Phys. Rev. X 12, 041013 (2022).
45. Abrikosov AA Magnetic properties of superconductors of the second group Sov. Phys. 1957 5 1174 1182
Abrikosov, A. A. Magnetic properties of superconductors of the second group. Sov. Phys. 5, 1174–1182 (1957).
46. Hess H Robinson R Dynes R Valles Jr J Waszczak J Scanning-tunneling-microscope observation of the Abrikosov flux lattice and the density of states near and inside a fluxoid Phys. Rev. Lett. 1989 62 214 10.1103/PhysRevLett.62.214 10039952
Hess, H., Robinson, R., Dynes, R., Valles Jr, J. & Waszczak, J. Scanning-tunneling-microscope observation of the Abrikosov flux lattice and the density of states near and inside a fluxoid. Phys. Rev. Lett. 62, 214 (1989).10039952 10.1103/PhysRevLett.62.214
47. Wallraff A Quantum dynamics of a single vortex Nature 2003 425 155 158 10.1038/nature01826 12968173
Wallraff, A. et al. Quantum dynamics of a single vortex. Nature 425, 155–158 (2003).12968173 10.1038/nature01826
48. Goldman A Kreisman P Meissner effect and vortex penetration in Josephson junctions Phys. Rev. 1967 164 544 10.1103/PhysRev.164.544
Goldman, A. & Kreisman, P. Meissner effect and vortex penetration in Josephson junctions. Phys. Rev. 164, 544 (1967).10.1103/PhysRev.164.544
49. Golod T Krasnov VM Demonstration of a superconducting diode-with-memory, operational at zero magnetic field with switchable nonreciprocity Nat. Commun. 2022 13 3658 10.1038/s41467-022-31256-w 35760801
Golod, T. & Krasnov, V. M. Demonstration of a superconducting diode-with-memory, operational at zero magnetic field with switchable nonreciprocity. Nat. Commun. 13, 3658 (2022).35760801 10.1038/s41467-022-31256-w
50. Zhou TX Stöhr RJ Yacoby A Scanning diamond NV center probes compatible with conventional AFM technology Appl. Phys. Lett. 2017 111 163106 10.1063/1.4995813
Zhou, T. X., Stöhr, R. J. & Yacoby, A. Scanning diamond NV center probes compatible with conventional AFM technology. Appl. Phys. Lett. 111, 163106 (2017).10.1063/1.4995813
51. Biercuk M Doherty A Uys H Dynamical decoupling sequence construction as a filter-design problem J. Phys. B: At. Mol. Opt. Phys. 2011 44 154002 10.1088/0953-4075/44/15/154002
Biercuk, M., Doherty, A. & Uys, H. Dynamical decoupling sequence construction as a filter-design problem. J. Phys. B: At. Mol. Opt. Phys. 44, 154002 (2011).10.1088/0953-4075/44/15/154002
52. Pham LM Enhanced solid-state multispin metrology using dynamical decoupling Phys. Rev. B 2012 86 045214 10.1103/PhysRevB.86.045214
Pham, L. M. et al. Enhanced solid-state multispin metrology using dynamical decoupling. Phys. Rev. B 86, 045214 (2012).10.1103/PhysRevB.86.045214
53. Blakely, R. J. Potential theory in gravity and magnetic applications (Cambridge University Press, 1996).
54. Lima, E. A. & Weiss, B. P. Obtaining vector magnetic field maps from single-component measurements of geological samples. J. Geophys. Res. Solid Earth 114, B06102 (2009).
55. Casola F Van Der Sar T Yacoby A Probing condensed matter physics with magnetometry based on nitrogen-vacancy centres in diamond Nat. Rev. Mater. 2018 3 1 13 10.1038/natrevmats.2017.88
Casola, F., Van Der Sar, T. & Yacoby, A. Probing condensed matter physics with magnetometry based on nitrogen-vacancy centres in diamond. Nat. Rev. Mater. 3, 1–13 (2018).10.1038/natrevmats.2017.88
56. Chen, S. & Park, S. Current induced hidden states in Josephson junctions. Zenodo. 10.5281/zenodo.13256436 (2024).
