Sub-Sharvin conductance and Josephson effect in graphene

This paper numerically demonstrates that tuning the electrostatic potential profile in superconductor-graphene-superconductor junctions from rectangular to parabolic causes the critical current-resistance product ($I_cR_N$) to transition from a graphene-specific value to the ballistic limit in the unipolar regime, while remaining graphene-specific and suppressing conductance in the tripolar regime.

The Gist

Imagine a superhighway made of a single layer of carbon atoms called graphene. This highway is incredibly fast; electrons zip through it like race cars with no friction. Now, imagine we block off a section of this highway with two superconducting "toll booths" (materials where electricity flows with zero resistance). This setup is called a Josephson Junction.

Normally, when you put a barrier between two superconductors, the electrons can "tunnel" through it, creating a special kind of super-current. The strength of this current depends heavily on how the barrier looks.

This paper is like a traffic engineer studying what happens when we change the shape of the roadblock in the middle of this graphene highway.

The Main Characters

1. The Electrons: They are the cars.
2. The Barrier (The Gate): This is the roadblock. The researchers can change its shape using electric gates.
   • The "Square" Wall: A sudden, sharp cliff. The road drops off instantly.
   • The "Curved" Hill: A gentle, smooth slope (like a parabola). The road rises and falls gradually.
3. The Traffic Flow (Current): How much electricity can get through.
4. The "Skewness": This is a fancy word for the shape of the traffic pattern. Is the flow symmetrical (like a perfect bell curve), or does it lean to one side?

The Big Discovery

The researchers wanted to know: Does the shape of the roadblock matter?

They found two very different scenarios, depending on which "direction" the traffic is flowing:

Scenario A: The "One-Way" Street (Unipolar Regime)

Imagine all the cars are going in the same direction (all electrons).

• The Sharp Wall: When the barrier is a sudden cliff, the traffic behaves in a very specific, "graphene-only" way. It's unique to this material.
• The Smooth Hill: When they smoothed the barrier into a gentle hill, the traffic behavior started to change. It stopped acting like "special graphene traffic" and started acting like "generic, perfect highway traffic."
• The Analogy: It's like driving through a toll booth. If the booth is a sudden wall, you have to stop and pay in a very specific, awkward way. If the booth is a gentle ramp, you just glide through like you're on a normal highway. The "special graphene magic" disappears when the road gets too smooth.

Scenario B: The "Three-Lane" Street (Tripolar Regime)

Imagine the middle of the road is going one way, but the sides are going the other way (electrons and "holes" mixing).

• The Result: No matter if the barrier is a sharp cliff or a smooth hill, the traffic always behaves in that special "graphene-only" way.
• The Analogy: It's like a complex roundabout. Even if you smooth out the curves of the roundabout, the cars still get stuck in that specific, unique pattern because the traffic flow itself is too complicated to be changed by the road shape. The "graphene magic" is robust here.

Why Does This Matter?

In the world of quantum computers, we use these tiny junctions to store information (qubits). To make them work, we need to predict exactly how much current will flow.

• The "Sub-Sharvin" Surprise: The paper also talks about "Sub-Sharvin conductance." Think of this as a speed limit. In a perfect, smooth highway, cars should go at a theoretical maximum speed (Sharvin limit). But in graphene, even on a smooth road, the cars often go slightly slower (Sub-Sharvin) because of how the carbon atoms are arranged.
• The Takeaway: If you are building a quantum device with graphene, you can't just assume the road shape doesn't matter.
  • If your device is in the "Three-Lane" mode, you are safe; the special graphene rules apply no matter how you build the barrier.
  • If your device is in the "One-Way" mode, you have to be very careful. If you make the barrier too smooth, you lose those special graphene properties and your device might behave like a generic, less interesting material.

The Bottom Line

This paper is a guide for engineers building the next generation of super-fast computers. It tells them: "Don't just build a wall; build the right shape of wall."

If you want to keep the unique, super-efficient properties of graphene, you need to understand whether your electrons are flowing in a simple line or a complex mix. If they are simple, a smooth road ruins the magic. If they are complex, the magic stays strong no matter what.

It's a reminder that in the quantum world, the shape of the road is just as important as the cars driving on it.

Technical Summary

1. Problem Statement

The paper investigates the transport properties of Superconductor-Graphene-Superconductor (S-g-S) Josephson junctions, specifically focusing on how the electrostatic potential profile within the graphene barrier affects the critical current ($I_c$), the normal-state resistance ($R_N$), and the skewness of the current-phase relation.

While previous studies (e.g., Titov and Beenakker, 2006) established that for an idealized rectangular barrier with infinite height, the product $I_c R_N$ (in units of $e/\Delta_0$) falls within a specific "graphene-specific" range ($2.1 \lesssim I_c R_N \lesssim 2.4$), it was unclear how these values behave when the potential barrier is smooth (e.g., parabolic or graded) rather than abrupt. This is crucial for realistic experimental setups where gate electrodes create smooth potential variations. The study aims to determine if the unique graphene-specific signatures persist in smooth potential regimes or if they revert to generic mesoscopic limits (tunneling or ballistic).

2. Methodology

The author employs a numerical approach based on the Dirac-Bogoliubov-de Gennes (DBdG) equations to model the system.

• System Model: A graphene strip of width $W$ and length $L$ is contacted by two superconducting electrodes. The system is modeled using the Dirac Hamiltonian for massless Dirac fermions in graphene, coupled with a superconducting pair potential $\Delta(x)$.
• Potential Profile: The electrostatic potential $V(x)$ in the normal region is defined by a tunable parameter $m$:
  $$V(x) = -V_0 \times \begin{cases} 1 & |x| > L/2 \\ |2x/L|^m & |x| \leq L/2 \end{cases}$$
  • $m=2$: Parabolic potential.
  • $m \to \infty$: Rectangular (step) potential.
  • The study varies $m$ from 2 to $\infty$ to simulate the transition from smooth to sharp barriers.
• Regimes: The analysis covers two doping regimes:
  • Unipolar ($\mu > 0$): Electron doping throughout (n-n-n structure).
  • Tripolar ($\mu < 0$): Hole doping in the center, electron doping in leads (n-p-n structure).
• Calculation:
  • The scattering problem is solved numerically using a fourth-order Runge-Kutta algorithm to find transmission amplitudes $t$ for each transverse mode $k_y$.
  • Transmission probabilities $T_n$ are calculated for the normal state ($\Delta_0 = 0$).
  • The critical current $I_c$ and normal resistance $R_N$ are derived from the transmission eigenvalues using the multichannel Josephson equation and the Landauer-Büttiker formula, respectively.
  • The skewness $S$ is calculated as $S = \frac{2\theta_c}{\pi} - 1$, where $\theta_c$ is the phase difference at maximum current.

3. Key Contributions

1. Systematic Analysis of Potential Smoothness: The paper provides the first comprehensive numerical study of how smoothing the potential barrier (changing from rectangular to parabolic) affects the $I_c R_N$ product and skewness in S-g-S junctions.
2. Distinguishing Regimes: It rigorously differentiates between the unipolar and tripolar regimes, revealing that the potential profile affects them in fundamentally opposite ways.
3. Toy Model for Crossover: The author proposes a "toy model" (Eq. 28) using a single transmission eigenvalue parameterized by a crossover parameter $\Theta$ to approximate the transition between tunneling and ballistic limits, successfully reproducing the numerical data for the tripolar regime.
4. Robustness of Graphene Signatures: It identifies conditions under which the unique "graphene-specific" values of $I_c R_N$ and skewness are preserved despite non-ideal (smooth) barrier shapes.

4. Key Results

A. Normal-State Conductance ($1/R_N$)

• Unipolar Regime ($\mu > 0$): As the potential smoothens (increasing $m$), the conductance increases from the sub-Sharvin value ($1/R_N \approx \frac{\pi}{4} G_{Sharvin}$) toward the full Sharvin value ($G_{Sharvin}$).
• Tripolar Regime ($\mu < 0$): Both the normal-state conductance and critical current are suppressed as the potential is smoothed. The n-p-n structure creates a barrier that hinders transport more effectively when the barrier is smooth.

B. Critical Current Product ($I_c R_N$) and Skewness ($S$)

The most significant findings concern the dimensionless product $I_c R_N$ (in units of $e/\Delta_0$) and the skewness $S$:

• Tripolar Regime ($\mu < 0$):

  • The values of $I_c R_N$ and $S$ remain robust and stay within the graphene-specific range ($2.1 \lesssim I_c R_N \lesssim 2.4$ and $0.25 \lesssim S \lesssim 0.42$) regardless of the potential shape (from parabolic to rectangular).
  • This suggests that the n-p-n structure inherently enforces the graphene-specific transport characteristics, even with smooth barriers.

• Unipolar Regime ($\mu > 0$):

  • The behavior is highly dependent on the potential shape.
  • For smooth potentials (low $m$, e.g., parabolic), $I_c R_N$ approaches the ballistic limit ($\pi \approx 3.14$).
  • As the barrier becomes rectangular (high $m$), $I_c R_N$ evolves downward toward the graphene-specific range.
  • Unlike the tripolar case, the unipolar regime does not maintain the graphene-specific values for smooth barriers; it behaves more like a generic ballistic conductor.

• Near the Dirac Point ($\mu \approx 0$):

  • Both $I_c R_N$ and $S$ are remarkably insensitive to the potential shape ($m$). They remain close to the values derived for the ideal rectangular barrier at the Dirac point ($I_c R_N \approx 2.08$, $S \approx 0.25$).

C. Current-Phase Relation (CPR)

• The CPR in graphene is forward-skewed (non-sinusoidal).
• In the tripolar regime, this skewness is preserved across all potential profiles.
• In the unipolar regime with smooth potentials, the CPR becomes less skewed, approaching the behavior of a ballistic contact.

5. Significance and Conclusion

• Experimental Relevance: The results clarify that experimental observations of "graphene-specific" Josephson effects (values of $I_c R_N \approx 2.1-2.4$) are not artifacts of idealized rectangular barriers. They are robust features, particularly in the tripolar (n-p-n) regime and near the Dirac point, where smooth potential barriers (common in gated devices) do not destroy these signatures.
• Diagnostic Tool: The study suggests that measuring the skewness $S$ and the product $I_c R_N$ serves as a robust probe for identifying graphene-specific transport mechanisms, as these intensive quantities are less sensitive to contact resistances and barrier smoothness than $I_c$ or $R_N$ individually.
• Theoretical Advancement: By bridging the gap between idealized rectangular models and realistic smooth potentials, the paper refines the understanding of mesoscopic transport in Dirac materials, confirming that the unique physics of graphene (such as Klein tunneling and evanescent modes) dominates in specific doping regimes even in non-ideal geometries.

In summary, the paper demonstrates that while the normal-state conductance is highly sensitive to the smoothness of the potential barrier, the Josephson critical current product and skewness retain their unique graphene-specific values in the tripolar regime and near the neutrality point, validating the use of S-g-S junctions as a platform for quantum information processing and fundamental physics studies in graphene.