⚛️ quantum physics

Josephson energy of superconducting junctions: amorphous versus crystalline tunnel barriers

第一原理に基づくNEGF-DFTモデリングを用いた本研究は、超伝導接合におけるジョセフソンエネルギーがトンネル障壁の原子構造に対して指数関数的な感度を持つことを示しており、アモルファスAl $_2$ O $_3$ 障壁は、化学量論的な不均一性がパーコレーション的なトンネリング経路を生み出すことにより、結晶質のものよりも著しく高く、かつ変動の大きいエネルギー値を示すことを明らかにしている。

原著者： Wanting Zhang, Aldilene Saraiva-Souza, Félix Beaudoin, Xianghua Kong, Hong Guo, Yu Zhu

公開日 2026-09-10

📖 1 分で読めます🧠 じっくり読む

CC BY 4.0

原著者： Wanting Zhang, Aldilene Saraiva-Souza, Félix Beaudoin, Xianghua Kong, Hong Guo, Yu Zhu

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

電気抵抗なく電流が流れる超伝導回路の静かな世界において、ジョセフソン接合と呼ばれる極めて小さな部品が、不可欠なスイッチやバルブとして機能しています。これらの接合は、最も高度な量子コンピュータ、特にトランスモンとして知られる種類の量子ビットの心臓部です。これらの量子ビットの挙動は特定のエネルギー値によって支配されており、その値がコンピュータの計算速度や、異なる情報の状態をどれほど明確に区別できるかを決定します。このエネルギーは、わずか1、2ナノメートルという極めて薄い酸化アルミニウムの層を、電子がいかに容易にトンネルできるかに完全に依存しています。この障壁は非常に薄いため、その原子配列――原子が完璧な格子状に並んでいるのか、あるいはランダムでガラスのような無秩序な構造をしているのか――が、電子の動きに多大な影響を与えます。この微視的な無秩序が回路のエネルギーにどのように影響するかを正確に理解することは、信頼性の高い量子機械を構築する上で極めて重要ですが、障壁の混沌とした原子構造とデバイスのマクロなエネルギーとの間のつながりを特定することは、これまで困難でした。

研究チームは、これらの接合の詳細なコンピュータモデルを構築し、完全に秩序だった結晶性の障壁と、10種類の異なる無定形（アモルファス）の障壁を比較することで、このつながりを解明しました。彼らは、量子力学の法則と個々の原子の具体的な配置を組み合わせた強力な手法を用い、これらの構造を電子がどのように通過するかをシミュレートしました。研究は、同じ厚さでありながら内部構造が異なる酸化アルミニウム障壁に焦点を当てました。結晶モデルでは、原子は整然とした繰り返しのパターンを形成していますが、無定形モデルでは、溶融と急速冷却を模倣したプロセスによって生成された原子が、ランダムで凍結された無秩序な状態となっています。研究者たちは、これらの障壁を通る電子の透過率を計算し、その結果として得られるジョセフソンエネルギー、すなわち量子ビットの動作周波数を設定する鍵となるパラメータを算出しました。

結果は、これら2種類の障壁の間に劇的な違いがあることを明らかにしました。単一の結晶性接合の場合、計算されたエネルギーは安定して0.73ギガヘルツでした。しかし、10個の無定形接合は激しい値のばらつきを示しました。無定形グループの平均エネルギーは2.78ギガヘルツでしたが、個々の結果は、0.18ギガヘルツという低値から15.1ギガヘルツという高値まで劇的に変動していました。これは、無定形のサンプルが単に少し異なる結果を生み出したのではなく、ほぼ2桁に及ぶ範囲の結果を生み出したことを意味します。分布は大きく偏っており、ほとんどのサンプルは結晶値の近くに集まっていましたが、いくつかの外れ値が極端に高いエネルギーを持ち、平均値を押し上げていました。この変動性は、無定形障壁における原子のランダムな配置が、電子の経路が均一ではなく、各サンプルにおける偶然の原子配置に依存するような景観を作り出していることを示唆しています。

なぜこのようなことが起こるのかを理解するために、研究者たちは障壁を通る電子の動きを詳しく調査しました。彼らは、電子が障壁を乗り越えるためのエネルギーを持っていないにもかかわらず、障壁を通り抜ける「量子トンネル効果」が輸送メカニズムであることを発見しました。無定形サンプルでは、酸化物の化学組成が完全に均一ではなく、アルミニウムが豊富な領域もあれば、酸素が豊富な領域もあります。アルミニウムに富む領域は、電子にとってより低く通りやすい障壁として機能し、酸素に富む領域はより高く通りにくい障壁となります。一部の無定形サンプルでは、これらの低い障壁の領域がたまたま連なり、電子が均一な障壁よりもはるかに容易に流れることができる、連続的なパーコレーションのような経路を作り出していました。他のサンプルでは、これらの「易しい経路」が接続されなかったため、電子は高い障壁を突破するために苦労することを余儀なくされました。これが、一部の無定形接合が結晶の参照値よりもはるかに高いエネルギー値を持ち、他のものがそれよりも低くなった理由です。

本研究は、これらの微視的な知見が実験室で使用される実際のデバイスとどのように関連しているかについても検討しました。コンピュータモデルは、およそ1.44 × 1.25ナノメートルの極めて小さな材料片を表しています。量子コンピュータで使用される実際のジョセフソン接合はもっと大きく、通常は200 × 200ナノメートル程度の大きさがあります。これは、実際のデバイスには、これら微視的な領域が2万個以上含まれていることを意味します。実際のデバイスはこれほど大きいため、個々のシミュレーションで見られた極端な変動を実質的に平均化します。実際のデバイスの最終的なエネルギーは、多くの異なる微視的な経路の混合物となり、激しい変動を滑らかにしていると考えられます。これは、微小スケールでは原子の無秩序が著しい変動を生み出す一方で、マクロなデバイスは自己平均化の効果によって性能が安定することを示唆しています。

研究者たちは、ジョセフソンエネルギーの変動は、無定形障壁内における化学量論的な不均一性（アルミニウムと酸素の不均一な混合）の直接的な結果であると結論付けました。アルミニウムに富む低障壁領域が存在し、それが酸化物を横断して接続されることで、導電性を大幅に高める経路が形成され、高エネルギーの外れ値が生じます。逆に、これらの接続された経路を持たないサンプルは、完全な結晶よりも透過性が低くなることがあります。本研究はシミュレーションに依拠しており、使用された特定の数学的手法は酸化物のエネルギーギャップをわずかに過小評価する傾向がありますが、観察された傾向は、酸化物の微視的構造から超伝導回路のエネルギースケールへと至る、明確で定量的な経路を提供しています。この研究は、障壁のランダムな原子構造が単なる些細な詳細ではなく、量子コンピュータのエネルギー景観を決定づける根本的な要因であることを確立し、これらのデバイスの構築方法や性能の最適化について考えるための新しい視点を提供しています。

技術要約：超伝導接合におけるジョセフソン・エネルギー：アモルファス対結晶トンネル障壁

問題提起
ジョセフソン・エネルギー（ $E_J$ ）は、トランスモン超伝導量子ビットの周波数と非調和性を決定する基本的なパラメータである。Al/AlO $_x$ /Al接合において、 $E_J$ は極薄の酸化物障壁を通じた電子のトンネリングによって支配されており、障壁の原子構造に対して指数関数的な依存性を示す。マクロな回路モデルでは接合を集中定数要素として扱うが、実際の障壁は構造的および化学的に不均一であり、厚さの変動、酸素欠損、および化学量論的な不均一性を包含している。既存の原子論的モデルの多くは、結晶質障壁または特定の界面終端に焦点を当ててきた。しかし、現実のAlO $_x$ 障壁はアモルファス（非晶質）であり、単一の構造では表現できない確率論的な原子構成の変動をもたらす。課題は、アモルファス酸化物障壁の微細構造から、マクロなエネルギー尺度である $E_J$ への定量的な原子論的経路を確立することにある。

手法
著者らは、Al/Al $_2$ O $_3$ /Al接合における量子輸送をシミュレートするために、密度汎関数理論（DFT）と非平衡グリーン関数（NEGF）形式を組み合わせた、第一原理デバイスモデリングのワークフローを採用している。

原子構造の生成:
- 結晶参照モデル: Al(111)– $\alpha$ -Al $_2$ O $_3$ –Al(111)の整合超セル（420原子）を構築した。酸化物は、一方の側にAl、もう一方の側にOが配置される非対称な終端とした。
- アモルファス・アンサンブル: LAMMPS（Matsui力場を使用）を用いたメルト・クエンチ（溶融・急冷）分子動力学プロトコルにより、10個の独立したアモルファスAl $_2$ O $_3$ 障壁を生成した。これらはその後、不均質な金属–酸化物界面を扱うために、機械学習原子間ポテンシャルであるCHGNetを用いて、結晶質Al(111)スラブの間に挿入され、緩和された。各アモルファス接合は1080個の原子を含む。
- 幾何学的形状: 結晶モデルとアモルファスモデルは異なる横方向の周期性と断面積（ $A_0$ ）を持つが、これらは輸送解析において明示的に考慮されている。
量子輸送計算:
- 輸送計算は、線形結合原子軌道（LCAO）基底とノルム保存擬ポテンシャル（GGAレベル）を用いたNanoDCALパッケージを使用して実行された。
- 緩和された接合に対して、エネルギー依存の透過係数 $T(E)$ を計算した。
- 常伝導状態のコンダクタンス（ $G_N$ ）は、シミュレーションセルの面積 $A_0$ とデバイス面積 $A$ （比較のために $200 \times 200$ nm $^2$ に設定）でスケーリングされた、フェルミ準位における透過 $T(E_F)$ からランドゥアーの公式を用いて導出した。
- 超伝導ギャップ $\Delta = 0.18$ meVを仮定し、ゼロ温度のアンベガ・バラトフの関係式および定義 $E_J = \hbar I_c / 2e$ を用いて、臨界電流（ $I_c$ ）およびジョセフソン・エネルギー（ $E_J$ ）を抽出した。
散乱状態解析:
- 輸送メカニズムを理解するため、著者らはフェルミ準位における散乱状態の確率、特にエバネッセント減衰パターンと実空間での波動関数の局在化を分析した。

主な結果

$E_J$ の変動性: 10個のアモルファス接合は、ジョセフソン・エネルギーの広い分布を示した。平均 $E_J/h$ は 2.78 GHz で、標準偏差は 4.67 GHz である。対照的に、結晶参照モデルは 0.73 GHz を示した。
歪んだ分布: 分布は強く歪んでいる。中央値（0.86 GHz）は結晶参照値に近いが、個々の値は0.18 GHz（サンプル5）から15.1 GHz（サンプル6）まで、2桁近くに及ぶ。
透過特性: 面積正規化された透過曲線（ $T(E)/A_0$ ）は、すべてのサンプルで同様の絶縁体ギャップ構造を示している。しかし、透過の谷に対するフェルミ準位の位置は、サンプル間で大きく異なる。谷の位置が $E_F$ に対してわずかにシフトするだけで、トンネリングの指数関数的な性質により $T(E_F)$ に大きな変動が生じる。
輸送メカニズム: 散乱状態解析により、輸送は量子トンネリングを介して行われることが確認された。この変動は、アモルファス酸化物内の化学量論的不均一性に起因する。
- パーコレーション経路: 高い $E_J$ を持つサンプルでは、Alに富む領域（有効な障壁高さが低い領域）が酸化物を横切って接続し、パーコレーションのようなトンネリング経路を形成している。これにより、電子が障壁をより容易に通過できるようになり、コンダクタンスが大幅に増強される。
- 局在化: 低い $E_J$ を持つサンプルでは、このような接続された低障壁パスが存在せず、散乱状態は入射電極付近に局在しており、結晶質接合よりも低い透過率を示す。

意義と主張
本論文は、アモルファス酸化物障壁の微視的な原子構造と、超伝導回路のマクロなジョセフソン・エネルギーとの間の定量的な関連性を確立している。著者らは、アモルファス接合で見られる $E_J$ の大きな分散はランダムなノイズではなく、局所的な化学量論的揺らぎによって駆動される、パーコレーション様のトンネリングチャネルの確率論的な形成の直接的な結果であると主張している。

重要な点として、個々の原子論的モデルは極端な変動を示すが、現実的なマクロな接合（例： $200 \times 200$ nm $^2$ ）は $2 \times 10^4$ を超える微小な領域を包含している。したがって、デバイスレベルの $E_J$ は、これらの構成に対して「自己平均化」され、小さなシミュレーションセルで見られるような極端な外れ値ではなく、アンサンブル平均を反映することが期待される。

限界と今後の方向性
使用されたGGA汎関数は、酸化物のバンドギャップを過小評価する可能性があり、その結果、絶対的な透過値がわずかに過大評価されている可能性があるが、サンプル間の傾向はロバストであると著者らは認めている。彼らは、より定量的な精度を得るために、将来的な研究ではGGA+Uやハイブリッド汎関数（HSEなど）を用いるべきであると示唆している。さらに、より大きなアモルファス・アンサンブルを計算することで、実験的な接合間のばらつきと直接比較可能な $E_J$ の適切な統計分布を得られる可能性があると提案している。

関連論文