第一原理計算とは?わかりやすく解説!(量子力学:密度汎関数理論:DFT:電子状態計算:分子軌道など)
材料科学・化学・物理学の研究において、実験を行わずにコンピューターだけで物質の性質を予測する手法が急速に普及しています。
その中核をなすのが「第一原理計算(First-Principles Calculation / Ab initio Calculation)」です。
第一原理計算は原子・電子レベルの量子力学の方程式を数値的に解くことで、物質の電子状態・結合エネルギー・格子定数・バンド構造・磁性・熱的性質などを予測できる計算手法です。
本記事では、第一原理計算の定義・量子力学的背景・密度汎関数理論(DFT)の仕組み・分子軌道法との関係・計算の流れ・応用分野まで、わかりやすく体系的に解説していきます。
第一原理計算とは?その本質をひとことで言えば
それではまず、第一原理計算の定義と本質的な特徴について解説していきます。
第一原理計算とは、経験的なパラメータや実験データを使わず、量子力学の基本法則(シュレーディンガー方程式・密度汎関数理論など)だけを出発点として、物質の電子状態・物性を計算する手法です。
「第一原理(First Principles)」とは「最も根本的な原理・法則から出発する」という意味であり、ラテン語の「ab initio(始めから)」とも同義で使われます。
第一原理計算の対象と目的
第一原理計算が明らかにできる物質の性質(物性)は非常に幅広く、以下のようなものが挙げられます。
| 計算できる物性の種類 | 具体例 | 応用分野 |
|---|---|---|
| 電子構造・バンド構造 | バンドギャップ・フェルミ面・状態密度 | 半導体・絶縁体・金属の分類 |
| 結晶構造・格子定数 | 安定な結晶構造・格子定数・弾性定数 | 材料設計・圧力下での構造変化 |
| 化学結合・反応性 | 結合エネルギー・反応障壁・活性化エネルギー | 触媒設計・化学反応機構 |
| 磁性 | 磁気モーメント・磁気秩序(強磁性・反強磁性) | 磁性材料・スピントロニクス |
| 光学特性 | 吸収スペクトル・誘電率・屈折率 | 太陽電池・光学材料設計 |
| 熱力学的性質 | 生成エンタルピー・格子振動・比熱 | 材料の安定性・熱電材料 |
経験的手法との違い
計算化学・計算材料科学には第一原理計算以外に「半経験的手法(Semi-Empirical Methods)」や「古典分子動力学(Classical MD)」もあります。
これらと第一原理計算の最大の違いは「経験パラメータの有無」です。
半経験的手法・古典MDは実験データから決めたパラメータ(原子間ポテンシャルなど)を使うため、計算が速い一方で対象物質・条件によっては精度が低下することがあります。
第一原理計算は経験パラメータなしに量子力学の基本法則だけから計算するため、未知の材料・極端な条件(高温高圧など)でも高い予測精度が期待できるという強みがあります。
第一原理計算の理論的基礎:量子力学とシュレーディンガー方程式
続いては、第一原理計算の理論的基礎である量子力学とシュレーディンガー方程式について確認していきます。
第一原理計算がなぜ「原子・分子のシミュレーション」を実現できるのか、その物理的な根拠を理解することが重要です。
多体シュレーディンガー方程式
第一原理計算の出発点となる基本方程式が「シュレーディンガー方程式(Schrödinger Equation)」です。
定常状態の多体シュレーディンガー方程式
Ĥ Ψ = E Ψ
Ĥ:ハミルトニアン演算子(全エネルギーの演算子)
Ψ:多体波動関数(全粒子の量子状態を記述)
E:系の全エネルギー(固有値)
ハミルトニアンの内容(原子核と電子から成る系):
Ĥ = 電子の運動エネルギー + 原子核の運動エネルギー + 電子間反発 + 原子核間反発 + 電子−原子核引力
この方程式を解けば系の全エネルギーと電子状態が得られる。
しかし現実の物質では電子の数が数十〜数百億個にも達するため、多体シュレーディンガー方程式をそのまま解くことは原理的に不可能です。
そこで様々な近似手法が発展してきました。
ボルン・オッペンハイマー近似
多体問題を扱いやすくする最初の重要な近似が「ボルン・オッペンハイマー近似(Born-Oppenheimer Approximation)」です。
原子核は電子に比べて1836倍以上重く、動きが遅いため、「電子は原子核の位置が固定された静電場の中で運動する」と近似します。
この近似により、原子核の運動と電子の運動を分離して扱えるようになり、電子の問題(電子状態計算)を独立して解くことができます。
ハートリー・フォック法と分子軌道理論
量子化学における第一原理計算の代表的な手法のひとつが「ハートリー・フォック(HF)法」です。
HF法では、多体波動関数を「各電子が独立した軌道(分子軌道)に入る」というスレーター行列式で近似します。
各電子は他の電子が作る「平均的な静電場」の中で運動するという「平均場近似(Mean-Field Approximation)」が基礎です。
HF法は「分子軌道理論」の数学的基盤であり、化学結合・反応性・分子の電子構造の解析に広く使われています。
密度汎関数理論(DFT)の仕組みと特徴
続いては、現代の固体物理・材料科学での第一原理計算の標準的手法である「密度汎関数理論(DFT:Density Functional Theory)」について確認していきます。
DFTは1998年のノーベル化学賞(ウォルター・コーン)を受賞した理論であり、第一原理計算の主流手法として世界中で使われています。
DFTの基本定理(ホーエンベルク・コーンの定理)
DFTの理論的基盤となるのが「ホーエンベルク・コーンの定理(Hohenberg-Kohn Theorem)」です。
ホーエンベルク・コーンの定理(DFTの基礎)
定理1(一意性):基底状態の電子密度 ρ(r) が決まれば、外部ポテンシャル(原子核の位置)が一意に決まる。すなわち、電子密度がすべての情報を含む。
定理2(変分原理):真の基底状態の電子密度 ρ₀(r) は全エネルギー汎関数 E[ρ] を最小にする。
意義:多体波動関数 Ψ(N電子系では3N次元の関数)の代わりに電子密度 ρ(r)(3次元の関数)を基本変数として使える。これにより計算コストが劇的に削減される。
コーン・シャム方程式
ホーエンベルク・コーンの定理を実用的な計算スキームに落とし込んだのが「コーン・シャム方程式(Kohn-Sham Equations)」です。
コーン・シャム法では、相互作用のある多体電子系を「相互作用のない仮想的な電子系」に置き換えて扱います。
この仮想的な非相互作用電子系の各電子波動関数(コーン・シャム軌道)を自己無撞着に解くことで、元の多体問題の電子密度・全エネルギーが求められます。
電子間相互作用の難しい部分は「交換相関エネルギー(Exchange-Correlation Energy)」という汎関数に押し込められており、これの近似(LDA・GGA・ハイブリッド汎関数など)が計算精度を左右します。
DFTの計算の流れ
DFT計算の自己無撞着場(SCF)サイクル
① 初期電子密度 ρ(r) を仮定する(原子の電子密度の重ね合わせなど)
② 電子密度からコーン・シャムハミルトニアンを構成する(ハートリーポテンシャル + 交換相関ポテンシャル + 外部ポテンシャル)
③ コーン・シャム方程式を解いて新しい軌道・電子密度を求める
④ 新旧の電子密度を比較し、収束しているか確認する
⑤ 収束していない場合は新しい電子密度を混合して②に戻る(SCFループ)
⑥ 収束したら全エネルギー・力・応力などを計算・出力する
このSCF(Self-Consistent Field:自己無撞着場)サイクルを繰り返すことで、物質の基底状態の電子密度・全エネルギーが数値的に得られます。
第一原理計算の応用分野と最新動向
続いては、第一原理計算が実際にどのような研究・産業分野で活用されているかについて確認していきます。
第一原理計算の応用範囲は半導体・電池・触媒・医薬品から宇宙材料まで非常に広大です。
電池材料・太陽電池の材料設計
次世代電池(リチウムイオン電池・全固体電池・ナトリウムイオン電池)の材料設計において、第一原理計算は不可欠なツールになっています。
正極・負極・電解質材料の安定性・イオン伝導度・電子構造を計算することで、実験の前に有望な材料候補を絞り込み、開発コスト・期間を大幅に削減できます。
太陽電池材料(ペロブスカイト太陽電池・薄膜太陽電池)のバンドギャップ設計・欠陥特性の解析にも第一原理計算が活用されています。
触媒設計・表面反応の解析
第一原理計算は触媒設計において特に強力なツールです。
触媒表面での吸着エネルギー・反応経路・遷移状態のエネルギー障壁を原子レベルで計算することで、活性の高い触媒材料を理論的に予測できます。
「コンピューターによる触媒探索(高速スクリーニング)」と実験の組み合わせが、水素製造・CO₂還元・燃料電池触媒などの開発を加速しています。
マテリアルズインフォマティクスとの融合
近年、第一原理計算と機械学習・AIを組み合わせた「マテリアルズインフォマティクス(Materials Informatics)」が急速に発展しています。
大規模な第一原理計算データベース(Materials Project・AFLOW・NOMAD等)に蓄積された数百万件のデータを機械学習モデルで学習することで、「第一原理計算なみの精度を持ちながら数百万倍高速な力場(機械学習ポテンシャル)」が構築されています。
これにより、従来は不可能だった大規模系・長時間シミュレーション(数万〜数十万原子の動力学計算)が実現し、材料探索の効率が飛躍的に向上しています。
まとめ
本記事では、第一原理計算の定義・特徴・量子力学的背景(シュレーディンガー方程式・ボルン・オッペンハイマー近似)・ハートリー・フォック法・密度汎関数理論(DFT)の仕組み・SCFサイクル・電池・触媒・マテリアルズインフォマティクスへの応用まで幅広く解説してきました。
第一原理計算とは経験パラメータを一切使わず量子力学の基本法則だけを出発点として物質の電子状態・物性を計算する手法であり、その中核をなすDFT(密度汎関数理論)は現代材料科学・化学研究の標準的なツールとして世界中で活用されています。
AIと第一原理計算の融合(機械学習ポテンシャル・マテリアルズインフォマティクス)により、材料探索の速度と範囲は今後さらに拡大し、次世代エネルギー・医薬品・電子材料の開発に革命的な変化をもたらすでしょう。