[論文レビュー] 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 および歪度の傾向を明らかにした。

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.

研究の動機と目的

• グラフェンベースのジョセフソン接合における電気的障壁形状が超伝導輸送に与える影響を理解する動機づけ。
• 障壁プロファイルとキャリアドーピングに応じて I_c R_N と歪度 S がどのように変化するかを解明する。
• グラフェン特有の輸送特性が現れる領域（ディラック点、ユニポーラ、トリポーラ）を特定する。
• 障壁の滑らかさと有限障壁高さに対するグラフェン特有の特徴の頑健性を評価する。

提案手法

• V(x) が位置依存ポテンシャル Δ(x) の超伝導対ポテンシャルを持つグラフェンに対して、Dirac-Bogoliubov-De-Gennes 方程式を解く。
• モード整合と数値積分（ランゲ-クッタ法）を用いて、各横断モードとエネルギーに対する透過率 T_n を計算し、多モードジョセフソン公式から I(θ) を得て、ランドーア=ビュティカー理論から R_N を求める。
• |x| ≤ L/2 に対して V(x) = -V0 * |2x/L|^m によって V0 を一定に保ちつつ m の指標でパラボラ型から矩形型へ調整する。
• 短尺接合極限では、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）ドーピング領域の両方を考慮し、θ_c（最大電流の位相）を用いて歪度 S を解析する。

実験結果

リサーチクエスチョン

• RQ1 障壁プロファイル形状（放物線型から矩形型） は S-g-S グラフェン接合の I_c R_N にどのような影響を与えるか？
• RQ2 障壁が滑らかになるにつれて、ユニポーラとトリポーラのドーピング領域における常態伝導 1/R_N および電流位相関係の挙動はどうなるか？
• RQ3 ディラック点付近および障壁平滑化の下で、グラフェン特有の I_c R_N および歪度 S の値は残るか？
• RQ4 単純なカ crossover（トンネル結合からバリスティック）説明が、グラフェン接合の I_c R_N–S 関係を捉えられるか？

主な発見

• ユニポーラ領域では、障壁の滑らかさの増加に伴い（m の成長）グラフェン特有の値から ballistic 極限へ I_c R_N が進化する。
• ユニポーラドーピングでは、障壁が増大するにつれて 1/R_N がサーウィン値付近の Sharvin 値へ向かって上昇する（サーウィン値は G_Sharvin の近く）。
• トリポーラ領域では、障壁の滑らかさが進むと 1/R_N および I_c がともに抑制されるが、矩形ポテンシャルでも I_c R_N はグラフェン特有の値に近い。
• ディラック点近傍では、障壁形状の変化に対して 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 もグラフェン特有の値へ傾向を示す。
• おもちゃのクロスオーバーモデル T_k_y^(Θ) はトンネル結合からバリスティックへの進化を捉え、I_c R_N と S を単一パラメータ Θ にリンクさせ、μ<0 かつ滑らかなポテンシャル領域でグラフェン特有の領域と整合する。