First-principles calculations of internal conversion processes in spin defects

本論文は、多配置TDDFTと解析的な非断熱結合を組み合わせた予測的な第一原理フレームワークを紹介するものであり、これにより光学活性なスピン欠陥における内部転換率を正確に計算することが可能となり、ダイヤモンドのNV$^-$センターおよびSiCのディバカンシーに関する実験データとの長年の不一致を解決することに成功した。

要約

結晶（ダイヤモンドやシリコンカーバイドなど）の中に存在する、小さく光る欠陥を想像してみてください。それは微小な量子コンピュータとして機能しています。これらの欠陥は、まるで小さな舞台俳優のようです。レーザーを照射すると、彼らは励起され、より高いエネルギーレベル（「ステージ」）へと跳ね上がります。しかし、元の静止状態に戻るためには、進むべき道を選ばなければなりません。明るく光る道（放射失活）か、あるいは音を立てずに静かに滑り落ちる道（非放射失活）かです。

長い間、科学者たちがこれら「俳優」がどれほどの速さで「滑り落ちる」か（このプロセスは内部転換と呼ばれます）を予測しようとする際、非常に大まかな地図を使用してきました。彼らの計算は、単一車線の道路を走る一台の車だけを見て交通量を予測しようとするようなものでした。彼らは速度が極めて遅いと推測し続けていましたが、実際には交通は高速で動いていました。彼らの予測は、時には千倍も遅れていました。

この論文は、これらの予測を修正するために、高精細なGPSシステムを導入するものです。著者がどのようにこれを行ったのか、簡単な比喩を用いて説明します。

1. 「多体」問題：オーケストラ全体を見る

従来の計算手法は、欠陥内の電子を、単一の音を奏でるソロミュージシャンのように扱っていました。しかし実際には、これらの電子は複雑なジャズバンドであり、全員が同時に即興演奏を行い、互いに反応し合っています。

• 旧来の方法： バンドの相互作用を無視し、電子をあたかも一人の人間であるかのように扱っていました。
• 新しい方法： 著者らは、高度な手法（ハイブリッド汎関数を用いたTDDFT）を用いて、オーケストラ全体の音を聞く手法を用いました。すべての電子がどのように共に踊るか（多構成効果）を考慮することで、ようやく真の複雑なエネルギーレベルを聴き取ることができたのです。

2. 「振動」問題：一歩一歩を数える

電子がエネルギーレベルを下降するとき、単に落下するだけではありません。電子は余分なエネルギーを周囲の原子に放出させ、それらを振動させなければなりません。結晶を、何百万ものバネで作られた巨大なトランポリンだと考えてみてください。

• 旧来の方法： 科学者たちは時間を節約するために、トランポリンには一つのバネしかない、あるいは数個の「主要な」バネしかないと仮定して計算していました。彼らは、そのわずかな数のバネに基づいたエネルギー放出を計算していました。
• 新しい方法： 著者らは、すべての単一のバネが落下に寄与していることに気づきました。彼らは、欠陥の近くにある原子だけでなく、すべての振動する原子と同時に相互作用する方法を開発しました。これは、「非断熱結合（電子が原子をどれほど強く押すかを測定する専門的な方法）」を、試行錯誤による推測ではなく、解析的な数式として算出することによって実現しました。

結果：地図の修正

著者らは、この新しいGPSを二つの有名な「俳優」でテストしました。

1. ダイヤモンドの俳優（NVセンター）：

   • 謎： 科学者たちは、この俳優が特定の励起状態で非常に短い寿命を持つことを知っていましたが、古い計算ではもっと長く生き残るはずだとされていました。
   • 解決策： 新しい手法で「滑り落ちる」速度を計算したところ、それは驚異的に速い（1秒間に約1000億回）ことが分かりました。これは、最新の超高速実験による測定値と完璧に一致しました。これにより、「滑り落ち」こそが、この俳優が励起状態を長く維持できない主な理由であることが証明されました。

2. シリコンカーバイドの俳優（ディバカンシーセンター）：

   • 謎： この俳優については、古い計算では（光る現象のみに基づくと）約37ナノ秒間は励起状態を維持できるはずでした。しかし、実験では実際には15ナノ秒しか持たないことが示されていました。何かが欠けていたのです。
   • 解決策： 新しい手法は、科学者が見落としていた「隠れたドア」を発見しました。彼らは、これまで見過ごされていた重要な「滑り落ち」の経路（非放射チャネル）を発見し、それが崩壊を加速させていることを突き止めました。この隠れた経路を数式に加えたことで、予測はついに実験結果（15ナノ秒）と一致しました。

なぜこれが重要なのか

この論文は単に数学の問題を解決しただけではありません。これはユニバーサルなツールキットを提供しています。

• 「オーケストラ全体」（電子の相互作用）や「すべてのバネ」（振動）を無視することが、いかに極端に誤った答えを導くかを証明しました。
• 科学者が、高価な実験を事前に行うことなく、これらの量子欠陥がどのように振る舞うかを正確に予測することを可能にします。
• これらの小さな「量子ビット」（欠陥の磁気状態）が、エネルギーを失う前にどれくらいの期間持続するかを正確に把握することで、より優れた量子コンピュータの設計に向けた舞台を整えます。

要約すると、著者らは、電子の複雑なダンスとあらゆる原子の振動の両方を捉えることができる顕微鏡を作り上げ、これにより、これらの量子欠陥がどれほどの速さで「オフ」になるのかを正確に予測することがついに可能になったのです。

技術概要

技術要約：スピン欠陥における内部転換プロセスの第一原理計算

問題提起
ダイヤモンド中の負電荷窒素空孔（$NV^-$）中心や、シリコンカーバイド（SiC）中の中性ダイバカンシー（$VV^0$）中心などの光学活性スピン欠陥は、量子技術にとって極めて重要なプラットフォームである。これらの系における放射遷移や系間交差（ISC）速度は、第一原理の枠組みを用いて広く研究されてきたが、電子－フォノン結合によって駆動される同一スピン間の電子状態間の非放射遷移である内部転換（IC）プロセスについては、依然として十分に記述されていない。単一モード近似や、コーン・シャモ（Kohn-Sham）軌道に基づくアクセプティングモード近似といった一般的な近似法は、歴史的にIC速度を数桁過小評価してきた。この不一致は、光学的サイクルの完全な第一原理的記述を妨げ、励起状態の寿命、温度効果、および光学検出磁気共鳴（ODMR）スペクトルの正確な予測を困難にしている。

手法
著者らは、ICプロセスにおける非放射遷移速度（NRTR）を計算するための、広範かつ予測的な第一原理フレームワークを提案している。この手法は、2つの重要な進展を統合している：

1. 多体電子構造： 本フレームワークは、ハイブリッド汎関数（具体的にはDDH汎関数）を用いた線形応答時間依存密度汎関数理論（LR-TDDFT）を利用している。このアプローチは、単一軌道のコーン・シャモ記述では欠落しがちな、励起状態の多構成的性質を捉えることができる。著者らは、特定の励起状態を扱うためにスピン反転TDDFTの枠組みを採用している。
2. 解析的非断熱結合（NAC）および全フォノンモード： 波関数の微分を有限差分法で計算する代わりに、著者らはWESTコード内での拡張ラグランジュ法に基づく、解析的なNACベクトル（$d_{IF,A}$）の定式化を実装した。これにより、数百の原子を含む大規模なスーパーセルにおけるすべてのフォノンモードの電子－フォノン行列要素（$M^k_{IF}$）を効率的に計算することが可能となり、従来の計算手法におけるボトルネックを回避している。

NRTR $\Gamma_{NR,IF}(T)$ は、フェルミの黄金律を用い、すべての $N_{modes}$ フォノンモードの寄与を合算して算出される。計算には、多フォノン過程を考慮した生成関数アプローチと、電子－フォノンスペクトル密度関数 $D(E)$ が用いられる。

主な結果
本フレームワークは、ダイヤモンドの $NV^-$ 中心および4H-SiCの $VV^0$ 中心の光学的サイクルに適用された：

• $NV^-$ 中心（ダイヤモンド）：

  • シングレット状態（$^1A_1 \to ^1E$）： 計算されたNRTRは $\approx 10^1$ GHz のオーダーである。結果は、フェムト秒過渡吸収実験（例：78 Kにおける $\tau_{^1A_1} \approx 117.4$ ps に対し、実験値は約 92–102 ps）と優れた定量的一致を示している。これは、ICが上部シングレット状態の主要な崩壊メカニズムであることを裏付けている。
  • トリプレット状態（$^3E \to ^3A_2$）： 計算されたNRTRは、放射速度（$\approx 83.3$ MHz）と比較して無視できるほど小さい（$\lesssim 5$ MHz）であり、この遷移が主に放射的であるという理解と一致している。
  • モード解析： 本研究は、$\sim 170$ meV における局在振動モードが強い非断熱的性質を持つ一方で、全速度の約半分しか占めていないことを明らかにしている。全フォノンモードを含めることは、精度を確保するために不可欠であり、単一モード近似では完全なスペクトル密度を捉えることができない。

• $VV^0$ 中心（4H-SiC）：

  • トリプレット状態（$^3E \to ^3A_2$）： 放射的寄与のみを考慮した従来の第一原理的な推定では、実験的な寿命（$\sim 15$ ns）を過大評価していた。ICプロセスを組み込むことで、著者らはこの長年の不一致を解決し、非放射的崩壊が実験的な観測値と一致するように寿命を著しく減少させることを示した。
  • シングレット状態（$^1A_1 \to ^1E$）： 実験データの不在にもかかわらず、本フレームワークは $\approx 10^1$ GHz のNRTRを予測しており、この遷移もまたICプロセスによって支配されていることを示唆している。

意義および主張
本論文は、実験値と定量的に一致する、IC速度を計算可能な、効率的かつ正確な初の ab initio フレームワークを提供すると主張している。著者らは、主に以下の2点を貢献として強調している：

1. 不一致の解決： 本フレームワークは、電子状態の多体的な性質と全フォノンモードの寄与を正しく捉えることにより、従来の理論的推定値と実験的寿命との間の数桁に及ぶ不一致を解消する。
2. 手法の進展： 解析的なNACの実装により、有限差分による微分の計算上のボトルネックを排除し、数百の原子を含む系におけるすべてのフォノンモードの包含を可能にした。

著者らは、効果的な単一モード近似およびアクセプティングモード近似は、非放射速度を系統的に数桁過小評価することを結論づけている。彼らのアプローチは、スピン欠陥の完全な光学的サイクルおよびODMRスペクトルのパラメータフリーな予測を可能にする。さらに、解析的なNACの実装は、拡張系における ab initio 非断熱分子動力学、およびポラロンやエキシトン動力学の研究へのステップとして提示されている。