[论文解读] Sub-Sharvin conductance and Josephson effect in graphene

Adam Rycerz|arXiv (Cornell University)|Mar 1, 2026

Graphene research and applications被引用 0

一句话总结

该论文在数值层面分析了在将介电势垒从抛物线形调到矩形形时，石英相对效应与石墨烯基 S-g-S 结的常态导电性如何演变，揭示了石墨烯特定的 I_c R_N 与不同掺杂 regimes（单极 vs 三极）中的偏斜趋势。

ABSTRACT

Titov and Beenakker [Phys. Rev. B 74, 041401(R) (2006)] found, by solving the Dirac-Bogoliubov-De-Gennes equation, that the product of critical current and normal-state resistance for superconductor-graphene-superconductor (S-g-S) Josephson junction takes values (for a short junction and zero temperature) between $I_cR_N\approx{}2.1$ and $I_cR_N\approx{}2.4$ in units of $e/Δ_0$, where $Δ_0$ is the superconducting gap. These values are notably higher than the tunnelling bound ($π/2$), but lower than the ballistic bound ($π$). Here we analyze numerically the tunneling of Cooper pairs through S-g-S junctions in which the longitudinal electrostatic potential profile is tuned, within gates electrodes, from a rectangular to a parabolic one. In the unipolar regime (i.e., when the chemical potential is above the top of a barrier, $μ>0$), it is found that $I_cR_N$ gradually evolves from the graphene-specific to the ballistic value. At the same time, the normal-state conductance increases from the sub-Sharvin value of $1/R_N\approx(π/4)\,G_{ m Sharvin}$ towards to the Sharvin value $G_{ m Sharvin}=g_0|μ|W/(π\hbar{}v_F)$, with the conductance quantum $g_0=4e^2/h$, the junction width $W$, and the Fermi velocity in graphene $v_F$. In contrast, in the tripolar regime ($μ<0$), both normal-state conductance and the critical current are suppressed when smoothing the potential; however, $I_c{}R_N$ remains close to the graphene-specific range, even for a parabolic potential. The skewness of the current-phase relation is also discussed.

研究动机与目标

• 理解电静垒形状如何影响石墨烯基约瑟夫森结中的超导传输的动机。
• 确定 barrier profile 与载流子掺杂对 I_c R_N 与偏斜度 S 的演化。
• 识别在 Dirac 点、单极、三极等特定 regimes 下，石墨烯特定传输特征的出现。
• 评估 barrier 平滑与有限势垒高度对石墨烯特征签名的鲁棒性。

提出的方法

• 对含有位置依赖势垒 V(x) 与超导对偶势 Δ(x) 的石墨烯，求解 Dirac-Bogoliubov-De-Gennes 方程。
• 使用模式匹配与数值积分（Runge-Kutta）计算每个横向模和能量的透射率 T_n，然后通过多模约瑟夫森公式获得 I(θ)，并通过朗道-布特克关系得到 R_N。
• 通过指数 m 将 V(x) 从抛物线形调到矩形：V(x) = -V0 * |2x/L|^m（|x| ≤ L/2），V0 固定。
• 在短结近似下，将 I(θ) 与 R_N 与透射本征值 T_n 联系起来：I(θ) = (eΔ0/ħ) Σ_n T_n sin θ / sqrt(1 − T_n sin^2(θ/2))，R_N^−1 = (4e^2/h) Σ_n T_n。
• 考虑单极（μ > 0）与三极（μ < 0）掺杂 regime，并通过 θ_c（最大电流时的相位）分析偏斜度 S。

实验结果

研究问题

• RQ1 势垒轮廓形状（从抛物线到矩形）如何影响 S-g-S 石墨烯结的 I_c R_N？
• RQ2 在单极与三极掺杂 regime 下，随着势垒平滑化，常态导电 1/R_N 与电流相位关系的行为如何？
• RQ3 在 Dirac 点附近以及势垒平滑化条件下，石墨烯特定的 I_c R_N 与偏斜度 S 是否仍然存在？
• RQ4 是否存在简单的跨越描述（隧穿到准经典/炮弹极限），能捕捉石墨烯结的 I_c R_N–S 关系？

主要发现

• 在单极区域，I_c R_N 随着势垒平滑（m 增大）而从石墨烯特定值向准 ballistic 极限演变。
• 在单极掺杂时，常态导电 1/R_N 从次 Sharvin 向 Sharvin 值提升，随着 m 增大。
• 在三极区域，随着势垒平滑，1/R_N 与 I_c 都被抑制，但即使是抛物线势也能使 I_c R_N 接近石墨烯特定值。
• 在 Dirac 点附近，I_c R_N 与 S 对势垒形状变化具有鲁棒性，仍接近石墨烯特定值（Dirac 点在矩形/无势垒情况下 I_c R_N e/Δ0 ≈ 2.08，S ≈ 0.255）。
• 对于高掺杂 |μ| ≫ ħ v_F/L，在势垒平滑增加时，I_c R_N 倾向石墨烯特定值而非 ballistic 极限；偏斜度 S 亦趋向石墨烯特定值。
• 一个 toy 跨越模型 T_k_y^(Θ) 捕捉了隧穿到 ballistic 的演化，将 I_c R_N 与 S 联系到单一参数 Θ，并在 μ<0 且势垒平滑时与石墨烯特定区间对齐。