時間依存密度汎関数法のソースを表示

'''時間依存密度汎関数理論'''（じかんいぞんみつどはんかんすうりろん、{{lang-en-short|Time-dependent density-functional theory}}、略称: '''TDDFT'''）は、[[電場]]や[[磁場]]といった時間依存的ポテンシャルの存在下での多体系の性質と動力学を調べるために[[物理学]]および[[化学]]において使われる[[量子力学]]理論である。こういった場の分子や固体に対する効果はTDDFTを使って研究することが可能であり、励起エネルギー、周波数依存応答特性、光吸収スペクトルのような特徴を抽出できる。

TDDFTは[[密度汎関数理論]]（DFT）の拡張であり、概念的、計算的基礎は類似している。（時間に依存する）[[波動関数]]は（時間に依存する）[[電子密度]]と等価であることを示し、次に任意の相互作用のある系と同じ密度を返す相互作用のない架空の系の有効ポテンシャルを導く。こういった系を構築するうえでの問題はTDDFTでより複雑である。これは、とりわけ全ての瞬間における時間依存有効ポテンシャルがそれより前の全ての時間における密度の値に依存するためである。その結果として、TDDFTの実装についての時間依存近似の開発はDFTに遅れた。応用ではこの記憶の必要性はいつも決まって無視されている。

== 概要 ==
TDDFTの形式的基礎は'''[[ルンゲ・グロスの定理]]'''（'''RG定理'''、1984年）である<ref>{{cite journal|last=Runge|first=Erich|author2=Gross, E. K. U. |year=1984|title=Density-Functional Theory for Time-Dependent Systems|journal=Phys. Rev. Lett.|volume=52|issue=12|pages=997–1000|doi=10.1103/PhysRevLett.52.997|bibcode=1984PhRvL..52..997R}}</ref>（1964年の[[ホーヘンベルク・コーンの定理]]<ref>{{cite journal|title=Inhomogeneous electron gas|journal=Phys. Rev.|year=1964|first=P.|last=Hohenberg|author-link= ピエール・ホーエンバーグ|author2=Kohn, W. |volume=136|issue=3B|pages=B864–B871|doi=10.1103/PhysRev.136.B864|bibcode = 1964PhRv..136..864H |url=https://doi.org/10.1103/PhysRev.136.B864}}</ref>の時間依存版）。RG定理は、所与の初期波動関数について、系の時間依存外部ポテンシャルとその時間依存密度との間に唯一の写像が存在することを示す。これは、3''N''個の変数に依存する多体波動関数がわずか3個の変数のみに依存する密度と等価であること、そして系の全ての性質は電子密度の知識だけから決定することができることを含意する。DFTとは異なり、時間に依存する量子力学において一般的な最小化原理は存在しない。その結果として、RG定理の証明はHK定理よりもややこしい。

RG定理を所与とすると、計算的に有用な手法を開発するうえでの次の段階は、興味のある物理的（相互作用のある）系と同じ電子密度を持つ架空の相互作用のない系を決定することである。DFTと同様に、これは（時間に依存する）コーン–シャム系と呼ばれる。この系は[[ケルディッシュ形式]]において定義される[[作用 (物理学)|作用]]汎関数の{{仮リンク|停留点|en|Stationary point}}として形式的に見出される<ref>{{cite journal|last=van Leeuwen|first=Robert|year=1998|title=Causality and Symmetry in Time-Dependent Density-Functional Theory|journal=Phys. Rev. Lett.|volume=80|issue=6|pages=1280–283|doi=10.1103/PhysRevLett.80.1280|bibcode=1998PhRvL..80.1280V}}</ref>。

TDDFTの最も人気のある応用は、孤立系やそれほど多くはないが固体の[[励起状態]]のエネルギーの計算である。こういった計算は、線形応答関数 —すなわち、外部ポテンシャルが変化する時に電子密度がどのように変化するか— が系の厳密な励起エネルギーで極を持つ、という事実に基づく。こういった計算は、交換-相関ポテンシャルに加えて、交換-相関[[核 (代数学)|核]] — 密度に関する交換-相関ポテンシャルの[[汎関数微分]]— を必要とする<ref>{{cite book|last=Casida |first=M. E.|author2=C. Jamorski |author3=F. Bohr |author4=J. Guan |author5=D. R. Salahub |title=Theoretical and Computational Modeling of NLO and Electronic Materials|editor=S. P. Karna and A. T. Yeates|publisher=ACS Press|location=Washington, D.C.|year=1996|url=http://dcm.ujf-grenoble.fr/PERSONNEL/CT/casida/research/chong.ps|page=145–}}</ref><ref>{{cite journal|last=Petersilka|first=M. |author2=U. J. Gossmann |author3=E.K.U. Gross|year=1996|title=Excitation Energies from Time-Dependent Density-Functional Theory|journal=Phys. Rev. Lett.|volume=76|issue=8|pages=1212–1215|doi=10.1103/PhysRevLett.76.1212|pmid=10061664|bibcode=1996PhRvL..76.1212P|arxiv=cond-mat/0001154}}</ref>。

=== 詳細 ===
[[波動関数]]を''&Psi;''(''t'')、''t'' を時間、[[ハミルトニアン]]を<math> \hat{H} </math>とする。この時、出発点としての時間を含む[[シュレーディンガー方程式]]は、
:<math> i \hbar { \partial \Psi (t) \over {\partial t}} = \hat{H} \Psi (t) </math>
となる。ここで、<math> \hbar = h / 2 \pi </math>（''h'' は[[プランク定数]]）である。時刻 ''t''<sub>0</sub> 、''t'' でのそれぞれの波動関数の関係をシュレーディンガー表示で表すと、
:<math> \Psi (t) = e^{ -i \hat{H} (t - t_0) / \hbar} \Psi(t_0) </math>
となる。少々厳密ではないが、<math> t \to t + \Delta t </math>、<math> t_0 \to \, t </math>とし、波動関数の時間発展を、&Delta;''t'' の時間刻みによる逐次的な発展として考えると上式は、
:<math> \Psi (t + \Delta t) = e^{ - i \hat{H} \Delta t / \hbar} \Psi (t) </math>
となる。問題となるのは、<math> e^{ - i \hat{H} \Delta t / \hbar} </math>の部分の処理で、<math> \hat{H} \Delta t </math> に関して冪展開したり、指数関数部分に関して分解（[[鈴木＝トロッター分解]]）を施すなどして&Delta;''t'' に関しての逐次計算が行われ方程式が解かれる。

TDDFT (TDLDA) は、ポテンシャル部分が時間に依存する場合、例として時間によって変動する動的な[[電場]]、[[磁場]]中での電子の振舞いや、[[断熱近似]]が成り立たないような[[化学反応]]を扱うような場合などに適用される。たがし、この手法は密度汎関数理論が前提であり、正しい結果を与えることが保証されるのは[[基底状態]]に対してのみである。上記のような時間依存する系は、準位の交差など[[励起状態]]を扱う計算となっている。これまでの実際の計算例などから経験的にこのような励起状態をTDDFTは良く記述できていることが分かっているが、どのような場合でも正しい結果を与える保証はない。

==形式==
===ルンゲ・グロスの定理===
{{main|ルンゲ・グロスの定理}}

ルンゲとグロスのアプローチは、[[ハミルトニアン]]が
:<math>\hat{H}(t)=\hat{T}+\hat{V}_{\mathrm{ext}}(t)+\hat{W},</math>
の形式を取る時間に依存する[[スカラー場]]の存在下での単一要素系について考える。上式において、''T''は運動エネルギー演算子、''W''は電子-電子相互作用、''V''<sub>ext</sub>(''t'') は電子の数と連動して系を定義する外部ポテンシャルである。通常、外部ポテンシャルは系の核との電子の相互作用を含む。非自明な時間依存性について、追加の明示的に時間依存的なポテンシャルが存在する。これは、例えば、時間に依存する電場あるいは磁場から生じうる。多体波動関数は単一の[[初期条件]]の下で[[時間に依存するシュレーディンガー方程式]]にしたがって発展する。
:<math>\hat{H}(t)|\Psi(t)\rangle=i\hbar\frac{\partial}{\partial t}|\Psi(t)\rangle,\ \ \  |\Psi(0)\rangle=|\Psi\rangle</math>

その出発点としてシュレーディンガー方程式を利用し、ルンゲ・グロスの定理は、いかなる時点においても、密度は外部ポテンシャルを一意的に決定することを示す。これは2つの段階で成される。
# 外部ポテンシャルが任意の時間に関する[[テイラー級数]]で展開できることを仮定すると、追加の定数よりも異なる2つの外部ポテンシャルが異なる[[電流密度]]を生成することが示される。
# [[連続の方程式]]を使うと、有限の系について、異なる電流密度が異なる電子密度に対応することが示される。

===時間に依存するコーン–シャム系===
所与の相互作用ポテンシャルについて、RG定理は外部ポテンシャルが電子密度を一意的に決定することを示す。コーン–シャム・アプローチは、相互作用のある系と等しい電子密度を形成する相互作用のない系（相互作用ポテンシャルがゼロ）を選ぶ。こうすることの利点は、相互作用のない系を容易に解くことができること —相互作用のない系の波動関数は単一粒子[[分子軌道|軌道]]の[[スレイター行列式]]として表わすことができ、個々の軌道は3つの変数をもつ単一の[[偏微分方程式]]によって決定される— そして、相互作用ない系の運動エネルギーはこれらの軌道の観点から厳密に表すことができることである。したがって、問題は相互作用のないハミルトニアン''H''<sub>s</sub>を決定するポテンシャル（''v''<sub>s</sub>('''''r''''', ''t'') または ''v''<sub>KS</sub>('''''r''''', ''t'') と表わされる）を決定することである。
:<math>\hat{H}_{s}(t) = \hat{T}+\hat{V}_{s}(t)</math>
次に、このハミルトニアンが行列式波動関数を決定する。
:<math>\hat{H}_{s}(t)|\Phi(t)\rangle=i\frac{\partial}{\partial t}|\Phi(t)\rangle,\ \ \  |\Phi(0)\rangle=|\Phi\rangle</math>
行列式波動関数は方程式
:<math>\left(-\frac{1}{2}\nabla^{2}+v_{s}(\boldsymbol{r},t)\right)\phi_{i}(\boldsymbol{r},t)=i\frac{\partial}{\partial t}\phi_{i}(\boldsymbol{r},t)\ \ \ \phi_{i}(\boldsymbol{r},0)=\phi_{i}(\boldsymbol{r}),</math>
に従う一式の''N''個の軌道の観点から構築され、''ρ''<sub>s</sub>が相互作用のある系の密度と常に等しい
:<math>\rho_{s}(\boldsymbol{r},t)=\rho(\boldsymbol{r},t)</math>
ような時間に依存する密度
:<math>\rho_{s}(\boldsymbol{r},t)=\sum_{i = 1}^{N_{\textrm{b}}}f_{i}(t)|\phi_{i}(\boldsymbol{r},t)|^{2}</math>
を生成する。

ここで留意すべきは、上記の密度の式において、総和が<math>N_{\textrm{b}}</math>個「全ての」コーン–シャム軌道にわたること、<math>f_i(t)</math>が軌道<math>i</math>についての時間に依存する占有数であることである。もしポテンシャル''v''<sub>s</sub>('''r''',''t'') が決定できる、あるいは少くともよく近似できるならば、次に元のシュレーディンガー方程式（3''N''個の変数を持つ単一粒子微分方程式）は、それぞれ初期条件のみが異なる3次元における''N''個の微分方程式に置き換えられる。

コーン–シャム・ポテンシャルに対する近似を決定する問題は難易度が高い。DFTと類似して、時間に依存するKSポテンシャルは系の外部ポテンシャルと時間に依存するクーロン相互作用''v''<sub>J</sub>を抽出するために分解される。残った要素は交換–相関ポテンシャルである。
:<math>v_{s}(\boldsymbol{r},t)=v_{\rm ext}(\boldsymbol{r},t)+v_{J}(\boldsymbol{r},t)+v_{\rm xc}(\boldsymbol{r},t)</math>
彼らの独創的な論文において、ルンゲとグロスはディラック場を出発点とした場に基づく議論を通してKSポテンシャルの定義に取り組んだ。
:<math>A[\Psi]=\int\mathrm{d}t\ \langle\Psi(t)|\hat{H}-i\frac{\partial}{\partial t}|\Psi(t)\rangle</math>
波動関数の汎関数 ''A''[Ψ] として取り扱った、波動関数の変分は停留点として多体シュレーディンガー方程式をもたらす。電子密度と波動関数との間の一意的な写像を考え、ルンゲとグロスは次に密度汎関数
:<math>A[\rho]=A[\Psi[\rho]],\,</math>
としてディラック場を扱い、場の交換–相関要素についての形式的式を導いた。これが汎関数微分によって交換–相関ポテンシャルを決定する。後に、ディラック場の基づくやり方は、それを生成する応答関数の因果律を考えた時に逆説的結論をもたらすことが見つかった<ref>{{cite book|last=Gross|first=E. K. U. |author2=C. A. Ullrich |author3=U. J. Gossman|title=Density Functional Theory|editor=E. K. U. Gross and R. M. Dreizler|publisher=Plenum Press|location=New York|year=1995|series=B|volume=337|isbn=0-387-51993-9}}</ref> 。密度応答関数（外部ポテンシャルに関する電子密度の汎関数微分）は因果的でなければならない（いかなる時間におけるポテンシャルの変化もそれより前の時間の密度に影響を与えない）。しかしながら、ディラック場から応答関数は時間について対称的であり、必要とされる因果構造を欠いている。この問題に悩まされないやり方が後に複素時間経路積分の[[ケルディッシュ形式]]に基づく場を通して導入された。「実時間」における場の原理の精緻化による因果律パラドックスの別の解決法が最近{{仮リンク|ジョヴァンニ・ヴィナーレ|en|Giovanni Vignale}}によって提唱された<ref name="Vignale2008">{{cite journal|last1=Vignale|first1=Giovanni|title=Real-time resolution of the causality paradox of time-dependent density-functional theory|journal=Physical Review A|volume=77|issue=6|year=2008|doi=10.1103/PhysRevA.77.062511}}</ref>。

== 線形応答TDDFT ==
外部摂動が系の基底状態構造を完全に破綻しないという意味で小さければ、線形応答TDDFTを使うことができる。この場合、系の線形応答を解析することができる。これは、1次まで、系の変分が基底状態波動関数のみに依存し、DFTの全ての性質を単純に使うことができるため、大きな利点である。

小さな時間に依存する外部摂動<math>\delta V^{\text{ext}}(t)</math>を考える。
:<math>\begin{align}
\hat{H}'(t) &= \hat{H} + \delta V^{\text{ext}}(t) \\
\hat{H}'_{\text{KS}}[\rho](t) &= \hat{H}_{\text{KS}}[\rho]+\delta V_H[\rho](t)+\delta V_{\text{xc}}[\rho](t)+\delta V^{\text{ext}}(t)
\end{align}</math>
そして、電子密度の線形応答からすると、
:<math>\begin{align}
\delta \rho(\boldsymbol{r}t) &= \chi(\boldsymbol{r}t,\boldsymbol{r}'t')
\delta V^{\text{ext}}(\boldsymbol{r}'t') \\
\delta \rho(\boldsymbol{r}t) &= \chi_{\text{KS}}(\boldsymbol{r}t,\boldsymbol{r}'t')
\delta V^{\text{eff}}[\rho](\boldsymbol{r}'t')
\end{align}</math>
上式において、<math>\delta V^{\text{eff}}[\rho](t)=\delta V^{\text{ext}}(t)+\delta V_H[\rho](t)+\delta V_{\text{xc}}[\rho](t)</math>である。ここで、そしてこれ以後、プライム記号付きの変数は積分されているものと見なす。

線形応答領域内において、ハートリー（H）ポテンシャルと交換-相関（xc）ポテンシャルの線形順序への変分は密度変分に関して展開できる。
:<math>\begin{align}
\delta V_H[\rho](\boldsymbol{r}) &= \frac{\delta V_H[\rho]}{\delta\rho}\delta\rho=
\frac{1}{|\boldsymbol{r}-\boldsymbol{r}'|}\delta\rho(\boldsymbol{r}') \\
\delta V_{\text{xc}}[\rho](\boldsymbol{r}) &= \frac{\delta V_{\text{xc}}[\rho]}{\delta\rho}\delta\rho=
f_{\text{xc}}(\boldsymbol{r}t,\boldsymbol{r}'t')\delta\rho(\boldsymbol{r}')
\end{align}</math>

最後に、この関係をKS系に対する応答方程式に挿入し、得られた方程式と物理的系についての応答方程式を比較すると、TDDFTのDyson方程式が得られる。
:<math>\chi(\boldsymbol{r}_1t_1,\boldsymbol{r}_2t_2)=\chi_{\text{KS}}(\boldsymbol{r_1}t_1,\boldsymbol{r}_2t_2)+
\chi_{\text{KS}}(\boldsymbol{r_1}t_1,\boldsymbol{r}_2't_2')
\left(\frac{1}{|\boldsymbol{r}_2'-\boldsymbol{r}_1'|}+f_{\text{xc}}(\boldsymbol{r}_2't_2',\boldsymbol{r}_1't_1')\right)
\chi(\boldsymbol{r}_1't_1',\boldsymbol{r}_2t_2) </math>

この最後の方程式から、系の励起エネルギーを導くことが可能である（これらは単に応答関数の極であるため）。

その他の線形応答アプローチには、Casida形式（電子孔対における展開）やSternheimer方程式（密度汎関数摂動理論）がある。

==TDDFTプログラム==
* [http://elk.sourceforge.net ELK]
* [[Firefly (プログラム)|Firefly]]
* [[GAMESS (US)]]
* [[Gaussian]]
* [[Amsterdam Density Functional]]
* [[CP2K]]
* [[Dalton (プログラム)|Dalton]]
* [[NWChem]]
* [[octopus (ソフトウェア)|Octopus]]
* [http://www.pw-teleman.org  pw-teleman library] 
* [[PARSEC]]
* [https://github.com/LLNL/qball Qbox/Qb@ll]
* [[Q-Chem]]
* [[Spartan]]
* [[TeraChem]]
* [[TURBOMOLE]]
* [[YAMBO code]]
* [[ORCA (量子化学プログラム)|ORCA]]
* [[Jaguar (ソフトウェア)|Jaguar]]
* [https://wiki.fysik.dtu.dk/gpaw/ GPAW]
* [http://www.onetep.org/ ONETEP]

== 脚注　==
{{reflist}}

== 参考文献 ==
* {{cite book|editor1=M.A.L. Marques |editor2=C.A. Ullrich |editor3=F. Nogueira |editor4=A. Rubio |editor5=K. Burke |editor6=E.K.U. Gross |title=Time-Dependent Density Functional Theory |publisher=[[Springer-Verlag]]|year=2006| isbn =978-3-540-35422-2}}
* {{cite book|author=Carsten Ullrich |title=Time-Dependent Density-Functional Theory: Concepts and Applications (Oxford Graduate Texts) |publisher=[[Oxford University Press]]|year=2012 |isbn =978-0199563029}}

== 関連項目 ==
*[[局所密度近似]]
*[[第一原理バンド計算]]

== 外部リンク ==
* [http://tddft.org tddft.org]
* [http://th.physik.uni-frankfurt.de/~engel/tddft.html Brief introduction of TD-DFT]

{{DEFAULTSORT:しかんいそんみつとはんかんすうほう}}
[[Category:固体物理学]]
[[Category:計算物理学]]
[[Category:密度汎関数理論]]
[[Category:量子化学]]