浅い欠陥状態に対する高精度な DFT-1/2 手法：シリコン中のドナー結合エネルギーを効率的に計算する方法

シリコンの小さな変化が重要な理由

すべてのコンピュータチップや太陽電池は、電気の流れを精密に制御するために添加された不純物（ドーパント）に依存しています。超高効率トランジスタからドナーを基盤とする量子ビットに至るまで、半導体結晶内のドーパント原子に余分な電子がどれだけ強く結び付けられているかを正確に知る必要がある技術が増えています。本論文は、その結合エネルギーを高精度で、しかもより高速かつ実用的に算出する手法を紹介します。特に現代のエレクトロニクスの基盤であるシリコンのドーパントに対して有効です。

余分な電子を与える原子たち

純粋なシリコンでは、原子は規則的な結合で電子を共有しており、室温では導電性は低めです。リン、ヒ素、アンチモン、ビスマスなどの第V族元素を微量加えると、各ドーパントは余分な1個の電子をもたらします。その余分な電子は自由に飛び回るわけではなく、ドーパントと周囲のシリコンに緩く束縛された水素様の雲として存在します。この結合の強さ、すなわちドナー結合エネルギーは、電子がどれだけ容易に解放されて電流を運ぶか、あるいは量子操作に参加できるかを決めます。実験での測定は確立されていますが、第一原理計算から信頼性高く予測することは困難で計算コストも高くなりがちです。

標準的な計算が及ばない点

材料設計の実務で広く使われる密度汎関数理論（DFT）にも欠点があり、電子の局在化を過小評価し、半導体におけるバンド端の位置を誤る傾向があります。浅いドナーに対しては、電子雲が何十個もの原子に広がるため、DFTは通常、結合エネルギーを著しく小さく予測してしまいます。ハイブリッド汎関数やGW計算のようなより高度な手法はこれらの問題を改善できますが、特に数千原子を含む大きなシミュレーションボックスが必要な場合、計算コストが極めて大きくなります。従来の“タンデム”手法は、異なるセルサイズで異なる理論レベルを混在させ結果を継ぎ合わせる必要があり、ワークフローが複雑で系依存になりがちでした。

小さな補正で大きな効果

著者らは、標準DFTに近似的な自己エネルギー補正を直接導入するDFT-1/2と呼ばれる手法を活用しています。実務的には、特定の原子のある原子軌道から概念的に半電子を取り除くことで、その原子の有効ポテンシャルをわずかに修正します。まず、バルクシリコンにこの補正を適用して計算上のバンドギャップを実験値に近づけ、伝導帯の基準を信頼できるものにします。次にドナー状態の電子的性質を調べると、第V族ドーパント全てにおいてドナー状態はドーパント原子のs軌道が主導していることが分かります。そこでその軌道に対して特化した半電子補正を適用し、ドナー準位と最も近い空の伝導状態との分離を最大化するように単一のカットオフ半径を微調整します。重要なことに、この最適化された補正はシミュレーションボックスを拡大しても有効なままで、数千原子を含むスーパーセルでも再利用できます。

手法の精度

この二段階の補正（まずホストのシリコン、その後ドナー自身に対する補正）により、得られるドナー結合エネルギーは実験値と極めてよく一致します。シリコン中のヒ素については、予測値は実験値とわずか0.3ミリ電子ボルトしか異ならず、事実上完全な一致であり、はるかに計算負荷の大きいハイブリッド計算と同等の精度です。アンチモンとリンでは誤差はそれぞれ約5および8ミリ電子ボルトで、未補正のDFTに比べて大幅な改善です。非常に重いドーパントであるビスマスについては、エネルギーレベルをわずかに変える相対論効果であるスピン軌道相互作用も含めています。これにより計算された結合エネルギーは実験値から約5ミリ電子ボルト以内に収まり、従来のより負荷の高い手法で見落とされがちだった物理を明らかにします。シリコンに限らないことを示すため、同じワークフローを酸化亜鉛中の水素ドナーにも適用し、再び測定値を数ミリ電子ボルトの精度で再現しています。

将来のチップ設計に向けた実用的な道具

専門外の読者に向けた要点は、著者らが標準DFT計算の低コストと単純さを維持しながら、はるかに重い手法と同等の精度に到達するレシピを作り上げたことです。ホスト材料の全体的なバンド構造とドーパント周辺の局所環境の両方を系統的に補正することで、彼らのDFT-1/2プロトコルは非常に大きなシミュレーションセルでも信頼できるドナー結合エネルギーを提供します。これにより、日常のエレクトロニクスや新興の量子デバイスを制御するドーパントを研究するための強力かつ汎用性の高いツールとなり、個々の不純物が意図した通りに振る舞う材料設計を支援します。

引用: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

キーワード: 浅いドナー, シリコンドーパント, 密度汎関数理論, 量子材料, 半導体欠陥