現代の計算化学・材料科学において、分子や固体の電子状態を理論的に計算する手法としてTDDFT(Time-Dependent Density Functional Theory:時間依存密度汎関数理論)が重要な役割を果たしています。
光の吸収・発光現象、化学反応のダイナミクス、新材料の設計など、幅広い応用分野で活用されるこの理論は、ノーベル化学賞にも関連する現代量子化学の中核をなす手法です。
本記事では、TDDFTとは何か、DFT(密度汎関数理論)との違い、基本的な仕組み、そして応用分野について詳しく解説していきます。
TDDFTの基本概念とDFTとの関係
それではまず、TDDFTの基本概念と、基礎となるDFT(密度汎関数理論)との関係について解説していきます。
DFT(密度汎関数理論)の概要
DFT(Density Functional Theory:密度汎関数理論)は、量子力学的な多体電子系の基底状態(最もエネルギーが低い安定状態)を、電子密度を基本変数として計算する理論です。
Walter Kohnがこの理論の基礎を確立した功績により、1998年にノーベル化学賞が授与されました。
DFTは計算コストを大幅に削減しながら精度の高い計算が可能なため、現在の計算化学・計算物理学で最も広く使われている電子状態計算法のひとつです。
ただし、DFTは基底状態の計算に特化しており、励起状態(光を吸収して電子が高エネルギー状態になった状態)を直接計算することはできません。
TDDFTの誕生とDFTとの違い
TDDFT(Time-Dependent Density Functional Theory)は、DFTを時間依存の状況(励起状態・動的過程)に拡張した理論です。
1984年にErich RungeとEberhard Gross(Runge-Gross定理)によって理論的な基礎が確立されました。
TDDFTでは、時間に依存した外場(光などの電磁場)に対する電子密度の変化を追跡することで、励起エネルギー・吸収スペクトル・光化学反応などを計算できます。
DFTが基底状態(安定状態)の計算に特化しているのに対し、TDDFTは時間依存の外場に対する電子系の応答を計算できます。これにより光の吸収・発光・励起状態の計算が可能になります。
Runge-Gross定理の意義
TDDFTの基礎となるRunge-Gross定理は、時間依存の外部ポテンシャルと時間依存の電子密度の間に一対一の対応関係があることを示した定理です。
これにより、「時間依存の電子密度がわかれば、系のすべての物理量が原理的に決定できる」という強力な基礎が与えられました。
この定理が成立することで、DFTのKohn-Sham方程式を時間依存版に拡張した計算スキームが正当化されます。
TDDFTの計算の仕組み
続いては、TDDFTが実際にどのような手順で計算を行うかについて確認していきましょう。
時間依存Kohn-Sham方程式
TDDFTの実際の計算は、時間依存Kohn-Sham方程式を解くことで行われます。
各電子が独立した粒子として動く仮想的な系を設定し、その一粒子の波動関数が時間とともにどのように変化するかを追跡します。
この仮想系における電子密度が実際の多体電子系の密度と一致するように、交換相関汎関数と呼ばれる項で多体効果を取り込む設計です。
線形応答TDDFTと励起エネルギー計算
実用的な計算で最もよく使われるのが線形応答TDDFT(LR-TDDFT)です。
外部電磁場(光)が微小な場合、電子系の応答が外場に線形比例するという近似のもとで、励起エネルギーと振動子強度(光の吸収強度)を効率よく計算できます。
線形応答TDDFTで計算できる主な量
・光吸収スペクトル(UV-Visスペクトル)
・励起エネルギー(S1・S2・T1など)
・振動子強度(各遷移の明るさ)
・円二色性スペクトル(CD)
交換相関汎関数の選択
TDDFTの計算精度に大きく影響するのが交換相関汎関数の選択です。
一般的によく使われる汎関数(B3LYP・PBE0・CAM-B3LYP・ωB97X-Dなど)はそれぞれ特性があり、計算対象の系や求める物性によって適切な汎関数を選ぶことが重要です。
特に電荷移動励起状態(CT状態)の計算には、範囲分離ハイブリッド汎関数(CAM-B3LYP等)が有効です。
TDDFTの主な応用分野
続いては、TDDFTが実際にどのような分野で活用されているかを確認していきましょう。
光化学・発光材料の設計
TDDFTの最も重要な応用分野のひとつが、光化学と発光材料の設計です。
有機EL(OLED)材料の発光波長・量子収率の予測、蛍光色素の吸収・発光スペクトルの計算、光触媒の電子励起過程の解析などにTDDFTが活用されています。
これらの計算により、実験前に分子設計の方向性を絞り込み、効率的な新材料開発が実現できます。
生化学・薬学への応用
生体分子の光吸収・蛍光特性の計算にもTDDFTは重要なツールです。
タンパク質内の色素分子の吸収波長のシフト(ソルバトクロミズム)の解析、蛍光プローブの設計、光線力学療法(PDT)に使う光感受性物質の電子状態解析などに活用されています。
材料科学・ナノテクノロジーへの応用
半導体・金属ナノ粒子の光学特性の理論計算にもTDDFTは広く用いられます。
金や銀のナノ粒子の表面プラズモン共鳴、量子ドットの発光波長の設計、ペロブスカイト型太陽電池材料の光吸収計算など、最先端の材料科学研究に不可欠なツールとなっています。
TDDFTの限界と発展
続いては、TDDFTの現在の限界と、それを克服するための発展的な手法について確認していきましょう。
TDDFTの既知の限界
TDDFTは非常に強力な手法ですが、いくつかの既知の限界があります。
電荷移動励起状態に対する精度の問題、コニカルインターセクション(2つのポテンシャルエネルギー面が交差する特殊な点)の記述の困難さ、強い相関をもつ多励起配置状態への対応限界などが課題として知られています。
TDDFTの改良手法
これらの限界を克服するための研究が進んでいます。
範囲分離ハイブリッド汎関数の開発、TDDFTとGW近似・多体摂動論を組み合わせた手法、機械学習を活用した新しい汎関数の開発などが活発に研究されています。
TDDFTに代わる高精度の励起状態計算手法(EOM-CCSD、CASPT2、MRCI等)との使い分けも実践的に重要です。
計算ツールとソフトウェアの利用
TDDFTの計算には、GAUSSIAN、ORCA、Q-Chem、Turboモルフォなどの量子化学計算ソフトウェアが広く使われています。
これらのソフトウェアはTDDFTの実装を提供しており、専門的な量子化学の知識があれば比較的容易に利用できます。
クラウドHPC・スパコンの利用により、実用的なサイズの分子(数十〜数百原子)のTDDFT計算が日常的に行われています。
まとめ
本記事では、TDDFTの基本概念、DFTとの違い、Runge-Gross定理の意義、線形応答TDDFTの仕組み、光化学・材料科学・薬学などへの応用、そして限界と発展について解説しました。
TDDFTはDFTを時間依存の励起状態計算に拡張した手法であり、光の吸収・発光スペクトルの計算を中心に、有機EL材料から生体分子まで幅広い応用を持つ現代量子化学の重要ツールです。
計算コストと精度のバランスが優れているため、大学・企業の研究現場での日常的な計算ツールとして広く活用されています。
量子化学・材料科学への理解を深める上で、TDDFTへの知識は今後さらに重要性を増すでしょう。