Device variability of Josephson junctions induced by interface roughness

技術要約：界面粗さによって誘起されるジョセフソン接合のデバイスばらつき

問題提起
量子コンピューティング（例：ショアのアルゴリズム）に必要な数千個の量子ビットへと超伝導量子プロセッサがスケールアップするにつれ、デバイス間のばらつきが重大なボトルネックとして浮上している。超伝導量子ビット、特にトランスモンは、トンネル領域で動作するAl/AlOx/Alジョセフソン接合に依存している。この領域において、臨界電流 $I_c$ — そして結果としてジョセフソンエネルギー $E_J$ — は、微視的な障壁の詳細に対して指数関数的に依存する。Al/AlOx界面の幾何学的形状（局所的な厚さの変動や汚染など）のわずかな変化であっても、量子ビットの遷移周波数やゲートフィデリティに大幅なシフトを引き起こす可能性がある。マイクロスコピーやプロセス最適化を通じて接合の再現性を向上させるための広範な実験的努力が行われてきた一方で、微細構造の無秩序（具体的には界面の粗さ）と接合レベルのエネルギー統計を結びつける定量的な理論的枠組みは不足していた。

手法
著者らは、原子スケールの界面の無秩序とマクロな接合輸送特性を橋渡しする、定量的なモデリングフレームワークを開発した。このアプローチは、主に以下の3つのコンポーネントで構成されている：

1. 界面粗さのモデリング： Al/AlOx界面を、二つのパラメータである根平均二乗（RMS）粗さ振幅 $\sigma$ と横方向相関長 $\xi$ によって特徴付けられるガウス型ランダム場 $h(x,y)$ としてモデル化する。これらのパラメータは、既存の実験文献（例：結晶粒径や厚さのヒストグラムに関するTEM研究）および新しい分子動力学（MD）シミュレーションから導出される。CHGNet機械学習ポテンシャルを用いたMDシミュレーションでは、理想的な条件下でのAl表面の酸化をモデル化し、$\sigma$ を推定している。
2. 量子輸送近似： 現実的な接合サイズ（$\sim 100 \text{ nm} \times 100 \text{ nm}$）に対してボゴリューボフ・ド・ジャネンス（BdG）方程式を解くことの計算量的な困難さを回避するため、著者らは2つの主要な近似を用いている：
   • アンベガ・バラトフ（AB）関係式： 臨界電流 $I_c$ を常伝導抵抗 $R_N$ と結びつけるAB関係式を用いて計算する。これにより、全エネルギースペクトルにわたる積分を避け、フェルミエネルギーにおける透過係数に焦点を絞ることができ、計算の複雑さが軽減される。
   • 局所厚さ近似： フェルミ波長（$\lambda_F \approx 3.6$ Å）および障壁減衰長（$\lambda_D \approx 1.9$ Å）が、粗さの相関長（$\xi \sim 10$ nm）よりも十分に小さいことを考慮し、接合を独立したトンネルチャネルのアンサンブルとして扱う。全コンダクタンスは、局所的な障壁の厚さが粗さのプロファイルに従って変化する、局所的なコンダクタンス密度の積分として近似される。
3. 数値シミュレーション： ガウス型ランダム場モデルによって生成された様々なトンネル障壁プロファイルを持つ5,000個のジョセフソン接合のアンサンブルに対して、モデルを適用する。

主な結果
数値シミュレーションは、$E_J$ の統計的挙動に関するいくつかの定量的な洞察を与えている：

• 対数正規分布： アンサンブルにおけるジョセフソンエネルギー $E_J$ の分布はガウス分布ではなく、対数正規分布に従う。これは、$E_J$ が障壁の厚さに指数関数的に依存しており、その厚さが正規分布していることに起因する。その結果、確率密度関数は大きな値の方へ歪んでいる。
• RMS粗さ（$\sigma$）の影響： RMS粗さ振幅 $\sigma$ が増加すると、平均値 $\langle E_J \rangle$ と分散 $\text{Var}(E_J)$ の両方が著しく増加する。平均値の増加は、トンネリングの指数関数的な感度、すなわち、まれに発生する異常に薄い領域が全コンダクタンスを支配することに起因する。
• 相関長（$\xi$）の影響： 横方向相関長 $\xi$ が増加すると、平均値 $\langle E_J \rangle$ はわずかに減少するが、分散 $\text{Var}(E_J)$ は著しく増加する。分散が増加するのは、$\xi$ が大きくなると統計的に独立したトンネルチャネルの数が減少するため（自己平均化の減少）であり、それによってサンプル間のゆらぎが増幅される。平均は有限サイズ効果によりわずかに減少する。$\xi$ が接合寸法に近づくと、トンネリングを支配する極端に薄い障壁領域が実現する確率が低下するためである。
• 量子ビット周波数のばらつき： 導出された $E_J$ の統計量を用い、典型的な充電エネルギー $E_C$ を仮定すると、本モデルは、テストされた特定のパラメータにおいて、平均 6.145 GHz に対して約 0.364 GHz の量子ビット遷移周波数 $\omega_{01}/2\pi$ の標準偏差を予測する。

意義と主張
本論文は、表面の粗さがどのようにジョセフソンエネルギーのばらつきを誘起するかについて、「定量的かつ直感的な描写」を提供すると主張している。微細構造パラメータ（$\sigma$ および $\xi$）を接合レベルの統計量に直接結びつけることで、本研究は大規模量子プロセッサにおける歩留まりと均一性の課題を理解するための枠組みを提供する。著者らは、自らの結果が「接合設計に直接的な関連性」を持つことを強調しており、界面の粗さパラメータを制御することがデバイスのばらつきを最小限に抑えるために極めて重要であることを示唆している。

著者らは、自身の研究の範囲について控えめな立場を維持している。彼らは、本モデルが $\sigma$ と $\xi$ が接合の両方の界面において独立かつ同一であると仮定していることを指摘しており、これは多様な成長条件下での実験的な検証を必要とする仮定である。さらに、彼らは将来の方向性として、このモデルを製造された接合の実験的特性評価へと接続すること、および、接合と回路素子の共同設計（co-design）を可能にするために、彼らの超伝導ソルバーをマイクロ波工学ソフトウェアと統合することを挙げている。本論文は、新しい実験プロトコルを提案したり、ばらつきの問題を完全に解決すると主張したりするものではなく、界面の粗さが与える影響を定量化するための理論的ベースラインを確立するものである。