🔬 mesoscale physics

Evanescent-mode Casimir-Josephson force and gate-controlled resonances in ballistic graphene Josephson junctions

This paper develops a microscopic scattering theory to describe the gate- and phase-dependent Casimir-Josephson force in ballistic graphene Josephson junctions, revealing how charge neutrality yields an evanescent-mode force while gate doping induces propagating channels that cause large oscillations and sign reversals in the mechanical response.

Original authors: Shahrukh Salim

Published 2026-09-01

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Shahrukh Salim

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 quiet world of quantum physics, materials often behave as if they are made of waves rather than solid particles. When a special type of material called a superconductor is placed next to a normal conductor, a strange phenomenon occurs: electrons can cross the boundary and transform into their antimatter counterparts, holes, without losing energy. This process, known as Andreev reflection, creates a delicate bridge between the two materials. Scientists have long known that the energy levels of these electrons depend on the distance between the superconductors and the phase of their quantum waves. Just as a stretched spring exerts a force to return to its original length, these energy levels create a physical push or pull on the junction itself. This is a mechanical force born purely from the quantum nature of matter, a concept that links the invisible world of subatomic particles to the tangible world of motion and force.

A researcher has now mapped out exactly how this force behaves in a specific, highly pure form of carbon known as graphene. Unlike ordinary metals, graphene's electrons move in a way that mimics massless particles, and they possess a unique internal structure that allows them to reflect off boundaries in two distinct ways: either staying within their original energy band or jumping to a different one. The researcher built a detailed theoretical model of a graphene strip sandwiched between two superconductors to calculate the invisible force acting on this strip. They found that when the graphene is perfectly balanced, with no extra electrons or holes added, the force follows a simple, predictable rule determined by the width of the strip and the square of its length. In this state, the force is generated by electrons that do not travel freely but instead fade away exponentially as they move across the junction, a behavior the researcher describes as evanescent.

However, the story changes dramatically when the researcher introduces a gate voltage to add more charge carriers to the graphene. This doping turns the fading electrons into traveling waves that bounce back and forth between the superconducting contacts, creating a pattern of interference similar to light in an optical cavity. This shift causes the mechanical force to oscillate wildly, growing and shrinking as the voltage changes, and even flipping its direction from a push to a pull. The study reveals that these oscillations are not random; they are governed by a specific relationship between the number of charge carriers and the length of the strip. The researcher also clarified a common misconception in the field: the force is not simply a sum of separate parts from traveling and non-traveling electrons. Instead, the total force emerges from a complex, coherent interplay of all possible electron states, and trying to separate them into distinct categories can lead to misleading conclusions.

The calculations show that for a typical device with a superconducting gap of about 0.20 millielectron volts, a length of 200 nanometers, and a width of 4 micrometers, the resulting force is approximately 1.4 femtonewtons. If the materials used have a larger energy gap and the device is smaller, the force can reach about 35 femtonewtons. While these numbers are incredibly small, they are significant enough to be detected with sensitive instruments, provided the experiment is designed to filter out static electrical noise and focus on the specific rhythm of the force as the quantum phase changes. The researcher suggests that by modulating the phase of the superconductors and measuring the resulting vibration, scientists could isolate this quantum mechanical force from the background noise of the environment.

This work provides a new way to understand how quantum mechanics can drive mechanical motion in the most fundamental materials. It distinguishes between two different scales of control: one that determines whether electrons reflect by staying in their band or jumping to another, and another that controls how they interfere as they travel across the junction. By separating these effects, the study offers a clearer picture of the forces at play in future quantum devices. The findings suggest that the mechanical response of such junctions is not just a side effect but a rich, tunable signal that carries detailed information about the electronic structure of the material. As experimentalists work to build these devices, the theoretical map provided by this study will help them know exactly what to look for, turning a subtle quantum push into a measurable reality.

Technical Summary: Evanescent-mode Casimir-Josephson force and gate-controlled resonances in ballistic graphene Josephson junctions

Problem Statement
The paper addresses the calculation of the phase-dependent equilibrium mechanical force (the Casimir-Josephson force) in a ballistic superconductor-graphene-superconductor (SGS) Josephson junction. While phase-controlled mechanical forces in metallic weak links have been studied, graphene presents a distinct physical regime due to its Dirac quasiparticles, which possess conduction/valence-band and sublattice-pseudospin structures. This leads to two distinct types of Andreev reflection: retro (intraband) and specular (interband). Furthermore, unlike standard metallic point contacts where transmission is often modeled as length-independent, the transmission eigenvalues in graphene near charge neutrality are evanescent and retain explicit length dependence even in the short-junction limit. The author aims to develop a microscopic theory that captures the full energy-dependent spectrum, distinguishing between different crossover regimes and avoiding the methodological pitfalls of separating forces into "bound" and "continuum" components in a way that misrepresents the physics of mixed regimes.

Methodology
The author employs a microscopic scattering theory based on an energy-dependent Dirac scattering matrix embedded within a normalized Matsubara determinant.

Model: A monolayer graphene strip of length $L$ and width $W$ is contacted by heavily doped $s$ -wave superconductors with a phase difference $\phi$ . The system is described by the Dirac-Bogoliubov-de Gennes (BdG) equation.
Scattering Matrix: The normal-region scattering matrix $S_e(z; q)$ is derived for complex excitation energy $z$ , accounting for transverse modes $q$ . The hole scattering matrix is obtained via particle-hole symmetry.
Thermodynamic Potential: The phase-dependent grand potential $\delta\Omega$ is calculated using a normalized scattering determinant (Eq. 10), which includes contributions from both Andreev bound states and continuum quasiparticles. The mechanical force is derived as $F = -\partial \Omega / \partial L$ .
Comparative Benchmarks: The author compares the full energy-dependent result against a "frozen-scattering" short-junction approximation (where the scattering matrix is fixed at the Fermi energy). The difference is defined as a "beyond-short-junction remainder" rather than a pure continuum force, acknowledging that this remainder contains finite- $L/\xi$ corrections to subgap levels.
Spectral Diagnostics: Exact real-energy subgap poles are solved to diagnose the microscopic character (specular/interband vs. retro/intraband) of the states, using the band indices $s_e = \text{sgn}(\mu + E)$ and $s_h = \text{sgn}(\mu - E)$ .

Key Contributions and Results

Neutral Evanescent-Mode Law:
At charge neutrality ( $\mu=0$ ) and zero temperature, the full determinant approaches a closed-form analytical result proportional to $\Delta W / L^2$ . Specifically, at phase difference $\phi = \pi$ , the force correction is $\delta F = (2 \ln 2 / \pi) \Delta W / L^2$ . This result arises because the evanescent transmission eigenvalues $\tau(q) = \text{sech}^2(qL)$ possess an explicit length dependence, providing the leading mechanical sensitivity without requiring energy-dependent corrections above the gap to dominate.
Gate-Controlled Resonances and Sign Reversals:
When the graphene is gate-doped ( $\mu \neq 0$ ), propagating channels open, introducing Fabry-Pérot phase factors ( $k_x L$ ) into the scattering matrix. This leads to large oscillations in the phase-dependent force correction $\delta F(\phi) = F(\phi) - F(0)$ . Crucially, the gate voltage can induce sign reversals of $\delta F$ at fixed phase, a phenomenon driven by coherent resonances.
Distinction of Crossover Scales:
The paper explicitly distinguishes two dimensionless control parameters that are often conflated:
- Band Character Crossover: Controlled by $|\mu|/\Delta$ . This determines the transition from interband/specular Andreev reflection to intraband/retro reflection.
- Propagation Crossover: Controlled by $|\mu|L/(\hbar v_F)$ . This governs the onset of Fabry-Pérot interference and coherent propagation resonances.
  For short junctions ( $L/\xi \ll 1$ ), these scales are parametrically distinct; the system can become predominantly intraband before entering the regime of strong propagation oscillations.
Methodological Clarification on "Continuum" Forces:
The author demonstrates that subtracting the frozen-scattering short-junction result from the full calculation does not isolate a pure continuum force. The resulting "remainder" contains both corrections to the subgap spectrum and above-gap phase shifts. Consequently, they avoid assigning a unique "bound" vs. "continuum" hierarchy to the force, instead treating the full determinant as the primary thermodynamic quantity.
Thermal and Length Scaling:
Numerical results confirm that the neutral force scales as $\Delta W / L^2$ for $L/\xi \to 0$ . Thermal effects suppress the force magnitude and round the response near $\phi = \pi$ , with the characteristic scale set by the smaller of the superconducting gap or the phase-dependent Andreev minigap.

Significance and Claims
The paper claims to provide a Dirac-material extension of earlier proposals regarding superconductivity-induced mechanical forces. Its primary significance lies in:

Establishing that the mechanical response of graphene Josephson junctions is fundamentally different from metallic point contacts due to the length-dependent nature of evanescent modes at neutrality.
Providing a rigorous microscopic framework that separates the band-character crossover from the propagation crossover, preventing the misinterpretation of gate-dependent oscillations.
Offering a precise diagnostic tool (exact subgap poles) to identify the microscopic Andreev character without artificially decomposing the total force into phenomenological partial forces in mixed regimes.
Quantifying the magnitude of the force (ranging from femtonewtons to tens of femtonewtons for realistic parameters) and proposing a lock-in detection strategy based on phase modulation to distinguish the signal from static electrostatic and electromagnetic backgrounds.

The author maintains a modest stance, noting that their results apply to an ideal ballistic model with abrupt interfaces. They acknowledge that realistic factors such as smooth electrostatic profiles, disorder, and contact inhomogeneity will broaden the predicted sharp resonances, though the qualitative features of gate sensitivity and sign-changing resonances are expected to remain robust.

More like this