409. Mesoscopic Josephson Effect in a Graphene Disk at Magnetic Field

Introduction

The Josephson effect is one of the most distinctive manifestations of phase-coherent superconductivity. In its simplest form, a supercurrent flows through a weak link between two superconductors with a current-phase relation (CPR) close to \(I(\phi)=I_c\sin\phi\), where \(I_c\) is the critical current and \(\phi\) is the superconducting phase difference. In conventional tunneling junctions, the product \(I_cR_N\) is set by the superconducting gap \(\Delta_0\) and follows the Ambegaokar-Baratoff relation \(I_cR_N=(\pi/2)\Delta_0/e\), assuming a low-transparency barrier and equilibrium conditions. This framework is foundational in superconducting electronics, metrology, and qubit engineering.

Mesoscopic Josephson junctions depart from this textbook limit. When the weak link is short compared with phase-coherence and thermal lengths, and when conduction occurs through a small number of modes with high transmission probabilities, the CPR becomes non-sinusoidal. Its asymmetry is often quantified by the skewness parameter \(S\), which is positive when the supercurrent rises more sharply than a sine wave near \(\phi=\pi/2\). Such skewness is a direct signature of channels with transmission approaching unity, and it is especially relevant in systems where ballistic transport and quantum confinement dominate.

Graphene offers a particularly rich platform for mesoscopic Josephson physics. As a two-dimensional Dirac material with linear dispersion near the charge neutrality point, high carrier mobility, tunable density, and exceptional mechanical robustness, graphene allows one to build superconducting weak links with unusual transport properties. In a superconductor-graphene-superconductor (S-g-S) junction, Andreev reflection is modified by the relativistic-like band structure, and the interference of electron and hole trajectories can be controlled by geometry, doping, and magnetic field. A disk-shaped or Corbino geometry is especially appealing because it eliminates edge-dominated transport along the perimeter and instead probes radial mode propagation between concentric circular interfaces.

This article examines the mesoscopic Josephson effect in a graphene disk under magnetic field, focusing on the regime where the critical current tends to zero and the normal-state resistance diverges, yet the product \(I_cR_N\) remains finite and anomalously enhanced. In this limit, one finds \(I_cR_N\approx1.85\,\Delta_0/e\) and skewness \(S\approx0.14\), indicating a nontrivial CPR even when the junction appears strongly suppressed in linear-response transport. These results are obtained from a quantum-mechanical mode-matching analysis of the Dirac-Bogoliubov-de Gennes (DBdG) equation and compared with a simpler incoherent-scattering model involving two circular interfaces.

Graphene as a superconducting weak-link material

Graphene is not superconducting in its pristine form under ordinary conditions, but it is an outstanding proximitized conductor. When placed between superconducting leads, Cooper pairs can leak into graphene via the proximity effect, inducing phase-coherent transport over mesoscopic distances. Several aspects make graphene uniquely valuable for Josephson devices.

First, its electronic structure is described near the Dirac points by a massless Dirac Hamiltonian. This creates electron-hole symmetry that strongly influences Andreev reflection. Second, the carrier density can be tuned electrostatically over a wide range, allowing one to move from heavily doped regimes to the vicinity of the charge neutrality point. Third, graphene’s two-dimensional nature makes it ideal for planar superconducting circuitry and for hybrid devices integrated with standard nanofabrication. Finally, its high mobility and long mean free path can support ballistic transport, which is the key ingredient for mesoscopic CPR distortion.

From an industrial perspective, graphene-based Josephson junctions are relevant to superconducting quantum circuits, low-noise detectors, cryogenic interconnects, and magnetic sensors. The possibility of tailoring supercurrent with gate voltage and magnetic flux suggests applications in reconfigurable superconducting elements. In addition, the compatibility of graphene with van der Waals heterostructures opens a pathway toward scalable, atomically sharp interfaces and hybrid devices that combine superconductivity with topological, spintronic, or optoelectronic functionality.

The Corbino geometry is particularly significant in graphene. Unlike a rectangular strip, where transport can be strongly affected by edges, a graphene disk with inner and outer contacts supports radial current flow. This geometry is useful for isolating bulk transport and for studying quantum interference in a rotationally symmetric setting. In the presence of a magnetic field, the Corbino disk exhibits quantized angular momentum channels and Landau-level physics that can couple in a highly nontrivial way to superconducting proximity effects.

Mesoscopic Josephson physics and current-phase relations

The CPR of a Josephson junction encodes the microscopic transmission properties of the weak link. In the tunneling limit, the relation is nearly sinusoidal because Cooper-pair transfer is perturbative. In contrast, when the weak link contains highly transmitting channels, the CPR becomes skewed. This skewness arises because the Andreev bound states that carry the supercurrent are strongly phase dependent and acquire a more abrupt phase evolution as transmission increases.

For a short junction with transmission eigenvalue \(T_n\), the Andreev bound-state energy is approximately
\[
E_n(\phi)=\Delta_0\sqrt{1-T_n\sin^2(\phi/2)}.
\]
The associated supercurrent is obtained by differentiating the total free energy with respect to \(\phi\). High-transparency channels produce pronounced deviations from \(\sin\phi\), and the resulting skewness is often used as an experimental proxy for ballisticity or interface transparency.

Graphene junctions are especially susceptible to such behavior because conduction can be dominated by a small set of modes whose transmission depends sensitively on incidence angle, doping, and magnetic field. In a Corbino disk, the radial geometry quantizes angular momentum and creates mode-dependent transmission through circular interfaces. This can generate a CPR that is neither purely sinusoidal nor trivially related to the normal-state conductance. The magnetic field adds another layer of control by modifying the orbital phases and the effective matching between electron and hole trajectories.

The key observation is that mesoscopic Josephson behavior can persist even in a regime where the linear normal-state conductance is strongly suppressed. When the magnetic field is tuned so that \(I_c\to0\) and \(R_N\to\infty\), the product \(I_cR_N\) need not vanish. Instead, it can approach a finite value larger than the tunneling-junction benchmark. For the graphene disk considered here, the limit yields \(I_cR_N\approx1.85\,\Delta_0/e\), which exceeds the Ambegaokar-Baratoff value and indicates that the surviving supercurrent is carried by a nontrivial distribution of transmission channels.

Dirac-Bogoliubov-de Gennes description of the graphene disk

The microscopic theory of superconducting proximity in graphene is formulated using the Dirac-Bogoliubov-de Gennes equation. In this framework, electron and hole excitations are combined into a Nambu spinor, and the low-energy Hamiltonian near each valley includes the Dirac kinetic term, electrostatic potential, magnetic vector potential, and superconducting pair potential in the leads.

For a Corbino disk, the superconducting regions occupy the inner and outer circular contacts, while the graphene annulus forms the weak link. The magnetic field enters through the vector potential in a gauge convenient for circular symmetry, and the solutions are expressed in terms of radial functions with angular-momentum quantum numbers. The matching conditions at the interfaces enforce continuity of the wave function and determine the Andreev bound states.

Mode matching in the disk geometry is particularly powerful because it respects the exact symmetries of the problem. Each angular momentum channel contributes a transmission probability that depends on the magnetic field and the radial structure of the interfaces. The DBdG approach captures coherent electron-hole conversion and phase accumulation across the annulus, which are essential for predicting the CPR and the skewness. It also reveals how graphene’s pseudospin structure modifies the interface scattering compared with ordinary Schrödinger electrons.

A notable feature of graphene is specular Andreev reflection near the Dirac point, where the reflected hole can reside in the valence band. Although the detailed regime depends on doping and contact properties, the underlying Dirac nature affects the phase coherence and the distribution of transmission eigenvalues. In a magnetic field, the radial trajectories are bent and the effective quantization changes, leading to oscillatory behavior of the supercurrent and the normal resistance.

Magnetic field effects in the Corbino geometry

The magnetic field plays a dual role in a graphene Josephson disk. On one hand, it suppresses ordinary transport by introducing orbital dephasing and reducing the overlap between propagating modes. On the other hand, it can enhance the relative importance of specific channels that remain phase coherent, thereby reshaping the CPR. In a Corbino geometry, the field does not simply create edge states as in a Hall-bar geometry; instead, it modifies radial propagation and angular-momentum matching.

As the field is increased, the normal-state conductance can decrease dramatically, causing \(R_N\) to grow large. Simultaneously, the supercurrent is suppressed, and \(I_c\) can approach zero. However, because the surviving current is carried by a few channels with relatively high transmission, the ratio \(I_cR_N\) remains finite and can exceed the tunneling value. This is a hallmark of mesoscopic transport: the conductance is not determined by a classical resistor but by a sparse set of quantum paths whose coherence survives despite strong overall suppression.

In graphene, this effect is amplified by the interplay of magnetic flux, confinement, and Dirac dispersion. The Corbino disk acts as a quantum interferometer in which each angular momentum channel accumulates a magnetic phase around the annulus. The result is a field-dependent restructuring of the transmission spectrum. The CPR becomes skewed, and the skewness parameter reaches \(S\approx0.14\), indicating a measurable departure from a sinusoidal relation. Such skewness is modest compared with the most transparent ballistic junctions, but it is significant because it appears in a regime where the junction is almost pinched off in the normal state.

Comparison with an incoherent two-interface model

To interpret the full quantum-mechanical result, it is useful to compare it with a simpler model in which scattering between the two circular interfaces is treated incoherently. In this picture, the graphene annulus is regarded as a pair of independent circular boundaries separated by a region where phase coherence is partially randomized. The model neglects detailed interference between multiple radial traversals but retains the key role of interface transmission.

This simplified description is valuable because it isolates the effect of the circular geometry and the magnetic field on the transmission distribution. It can reproduce the qualitative trends of the full DBdG calculation, including the persistence of a finite \(I_cR_N\) product and the emergence of CPR skewness. However, it does not capture all coherent resonances, mode hybridization, or subtle phase cancellations. The comparison therefore clarifies which features are robust consequences of the graphene Corbino architecture and which arise from fully coherent quantum interference.

The agreement between the two approaches is physically important. It indicates that the anomalous Josephson response is not an artifact of fine-tuned resonances but instead reflects the generic transmission landscape of a disk-shaped graphene junction in field. The incoherent model also underscores the role of the two circular interfaces as the dominant scattering elements. This is especially relevant for device design, where interface quality, contact transparency, and electrostatic profile can be engineered to optimize supercurrent performance.

Scientific and technological implications

The observation of mesoscopic Josephson behavior in a graphene disk at magnetic field has implications well beyond a specific device geometry. Scientifically, it deepens our understanding of how superconductivity is mediated through Dirac materials under strong orbital quantization. It demonstrates that the standard tunneling-junction intuition is insufficient when the weak link is ballistic, mode-selective, and geometrically constrained. The finite value of \(I_cR_N\) in a regime where both quantities individually appear to vanish highlights the subtlety of quantum transport in mesoscopic superconducting systems.

From a materials perspective, graphene remains one of the best platforms for exploring the interface between superconductivity and relativistic-like electronic structure. Its clean band structure, tunability, and compatibility with nanolithography make it ideal for systematic studies of CPR distortion, Andreev spectra, and magnetic-field-controlled supercurrents. These studies also inform the broader field of van der Waals superconductivity, where graphene can be combined with other two-dimensional materials to create designer weak links.

Technologically, skewed CPRs are useful in superconducting circuits because they alter the inductive response and switching dynamics of Josephson elements. A controlled departure from sinusoidal behavior can be exploited in phase batteries, flux-sensitive detectors, and nonlinear resonators. In graphene-based devices, electrostatic gating can tune the transparency and therefore the skewness, offering a route to reconfigurable cryogenic components. The Corbino geometry may be especially valuable where edge disorder must be minimized or where rotational symmetry is advantageous for sensor design.

The magnetic-field response is also relevant to quantum metrology and magnetometry. Because the supercurrent is sensitive to flux through the disk, graphene Corbino junctions could serve as compact flux-tunable elements with unusual response functions. Their behavior in the regime of large \(R_N\) but finite \(I_cR_N\) suggests robustness of coherent transport against strong normal-state suppression, a property that may be useful in low-dissipation superconducting electronics.

Conclusion

The mesoscopic Josephson effect in a graphene disk under magnetic field reveals how superconducting transport can remain highly structured even when ordinary conductance is nearly extinguished. In a Corbino S-g-S junction, the combination of Dirac dispersion, radial quantization, and orbital field effects produces a non-sinusoidal current-phase relation with positive skewness. The full DBdG mode-matching analysis shows that, in the limit \(I_c\to0\) and \(R_N\to\infty\), the product \(I_cR_N\) approaches \(1.85\,\Delta_0/e\), while the skewness reaches \(S\approx0.14\). These values exceed the expectations of simple tunneling theory and reflect the presence of transmission channels with probabilities close to unity.

The comparison with an incoherent two-interface model confirms that the essential physics is rooted in the graphene disk geometry and the transmission structure of its circular interfaces. Graphene’s unique combination of ballistic transport, electrostatic tunability, and compatibility with superconducting contacts makes it an exceptional laboratory for mesoscopic Josephson phenomena. More broadly, these results point toward graphene-based superconducting devices whose response can be engineered through geometry and magnetic field, with promising implications for quantum circuits, sensors, and hybrid nanoscale superconducting technologies.