ヘルマン–ファインマンの定理
ヘルマン–ファインマンの定理(ヘルマン–ファインマンのていり、Hellmann–Feynman theorem)とは、量子力学において、パラメータ依存したハミルトニアンとそのエネルギー固有値に関する定理である。量子化学および数理固体物理学において特にヘルマン–ファインマン力(ヘルマン–ファインマンりょく、Hellmann–Feynman force)の計算に応用され、重要である。定理の名は、ドイツの物理学者H. ヘルマンと米国の物理学者R. P. ファインマンに因む。定理を最初に明示的な形で表したのは、P. Güttingerであるが、W. パウリやヘルマンの論文にも記されている。特にヘルマンは分子への適用に向けて、変分形式で表現した。また、1939年に当時、マサチューセッツ工科大学の学生であったファインマンは、この定理を示すともに、化学結合した原子において、電子及び他の原子核が原子核に及ぼす力は古典的な静電力として、扱えることを示した。「ヘルマン–ファインマンの定理」の名が定着したのは、J. C. スレイターがその著書の中でその名で呼んだことによる。
系のハミルトニアンが、あるパラメータ λ に依存するとして、それを H(λ) と表現する。 H(λ) の固有状態 ψλ があって、 H(λ)ψλE(λ)ψλ 及び規格化条件 ψλ1 が満足されるとする。このとき、
\frac{\mathrm d E (\lambda)}{\mathrm d \lambda} = \left \langle \psi_\lambda \left | \frac{\mathrm d \hat{H} (\lambda)}{\mathrm d \lambda}  \right | \psi_\lambda \right \rangle
が成り立つ。これがヘルマン–ファインマンの定理の主張である。
ここで、パラメータ λ が、原子位置座標 Rα の場合、ヘルマン–ファインマン力となる。
\begin{align}
\frac{\mathrm d E (\lambda)}{\mathrm d \lambda }
&= \frac{\mathrm d}{\mathrm d\lambda} \left \langle \psi_\lambda \left| \hat{H} (\lambda) \right | \psi_\lambda \right \rangle \\
&= \left( \frac{\mathrm d}{\mathrm d\lambda} \langle \psi_\lambda | \right) \hat{H} (\lambda) | \psi_\lambda \rangle
+ \left \langle \psi_\lambda \left | \frac{\mathrm d \hat{H}(\lambda)}{\mathrm d\lambda} \right | \psi_\lambda \right \rangle
+ \langle \psi_\lambda | \hat{H} (\lambda) \left (\frac{\mathrm d}{\mathrm d\lambda} | \psi_\lambda \rangle \right) \\
&= E(\lambda) \left( \frac{\mathrm d}{\mathrm d\lambda} \langle \psi_\lambda | \right) | \psi_\lambda \rangle
+ \left \langle \psi_\lambda \left | \frac{\mathrm d \hat{H}(\lambda)}{\mathrm d\lambda} \right | \psi_\lambda \right \rangle
+ E(\lambda) \langle \psi_\lambda | \left( \frac{\mathrm d}{\mathrm d\lambda} |\psi_\lambda \rangle \right) \\
&= E(\lambda) \frac{\mathrm d}{\mathrm d\lambda} \langle \psi_\lambda | \psi_\lambda \rangle
+ \left \langle \psi_\lambda \left | \frac{\mathrm d \hat{H}(\lambda)}{\mathrm d\lambda} \right | \psi_\lambda \right \rangle
\end{align}
ここで、 ψλ の規格化を ψλ1 と選んであるので、 dψλ/dλ0 である。よって、
\frac{\mathrm d E (\lambda)}{\mathrm d \lambda} = \left \langle \psi_\lambda \left | \frac{\mathrm d \hat{H}(\lambda)}{\mathrm d \lambda} \right | \psi_\lambda \right \rangle
が得られる。
定理の応用の1つとして、分子内力(intramolecular force)の計算があり、この手法で計算された力を特にヘルマン–ファインマン力と呼ぶ。ファインマンは1939年の「分子内の力(Forces in Molecules)」と題する論文の中で、ヘルマン–ファインマンの定理の証明を与えるとともに、分子や固体原子において、原子核に働く量子論的な力は、電子雲と他の原子核からの古典的な静電力として扱えることを示した。
電子及び位置が固定された原子核からなる系において、ポテンシャル V を原子核の位置座標で微分したものは、原子核に働く力に相当する。例えば、電子の位置を ri (xi, yi, zi)  (i 1, …,n) とし、原子核の位置を Rα(Xα, Yα, Zα)  (α 1,…,m) とする。このとき、系のハミルトニアン H はボルン–オッペンハイマー近似により、運動エネルギー
\hat{T} = -\sum_{i=1}^n \frac{\hbar^2}{2m} \nabla_i^2
\begin{align}
\hat{V} &= \sum_{i>j} \frac{e^2}
+ \sum_{\alpha > \beta} \frac{Z_{\alpha}Z_{\beta}e^2}
+ \sum_{i, \alpha} \frac{-Z_{\alpha}e^2} \\
&= \sum_{i>j} \hat{V}_{ij} + \sum_{\alpha > \beta} \hat{V}_{\alpha \beta} + \sum_{i, \alpha} \hat{V}_{i\alpha}
\end{align}
の和 T+V で表される。このとき、パラメータ λ として、原子核の位置座標 Rα をとったときに、その導関数
\boldsymbol{F}_{\alpha}= -\frac{\partial \hat{V}{\partial \boldsymbol{\lambda}= -\frac{\partial \hat{H}{\partial \boldsymbol{\lambda}
で定義される Fα(FXα, FYα, FZα)が、古典論では原子核に働く力となる。一方、量子論では、 H(λ)ψλE(λ)ψλ を満たす固有状態 ψλ により、
\begin{align}
\boldsymbol{F}_{\alpha}
&= - \left \langle \psi_{\boldsymbol{\lambda} \left| \frac{\partial \hat{H}{\partial \boldsymbol{\lambda} \right| \psi_{\boldsymbol{\lambda} \right \rangle \\
&= - \int \mathrm dV  \,\psi^* (\boldsymbol{q}_1, \dotsc, \boldsymbol{q}_n) \frac{\partial \hat{H}{\partial \boldsymbol{\lambda} } \psi(\boldsymbol{q}_1, \dotsc, \boldsymbol{q}_n)
\end{align}
が対応する。ここで2行目の波動関数 ψ における座標 qi は i 番目の電子の位置座標 ri とスピン座標 σi を合わせたものであり、dVdq1…dqn である。ファインマンの論文以前は、分子の量子力学では、これを
\boldsymbol{F'}_{\alpha} = - \frac{\mathrm d E(\boldsymbol{\lambda})}{\mathrm d \boldsymbol{\lambda}
とエネルギー固有値の微分で計算するのが、一般的であった。ファインマンはヘルマン–ファインマンの定理によって、Fα とF′α が等しいことを示した。
実際の計算において、微分係数 dE(λ)/dλで与えられる Fα を求めるには、エネルギー固有値のパラメータ λRα についての依存性の傾きを計算することになり、複数のパラメータ値に対して固有値問題を解く必要がある。一方、直接的に ψλdH/dλ ψλ で与えられるFα では、計算の労力を減らすことができる。
Fα の具体形については、波動関数による存在確率から定まる各電子の電荷密度
\rho_i(\boldsymbol{r}_i) = e \int | \psi(\boldsymbol{q}_1, \cdots, \boldsymbol{q}_n)|^2 \mathrm d\boldsymbol{q}_1 \dotsb \mathrm d\boldsymbol{q}_{i-1} d \sigma_i \mathrm d\boldsymbol{q}_{i+1} \dotsb \mathrm d\boldsymbol{q}_n
の和として、全電子の電荷密度 ρ(r) を、
\rho(\boldsymbol{r})=\sum_{i=1}^n \rho_i(\boldsymbol{r})
で導入すれば、
\boldsymbol{F}_{\alpha}= Z_{\alpha}e \int \rho(\boldsymbol{r})
\frac{\boldsymbol{R}_{\alpha}-\boldsymbol{r}〔\boldsymbol{R}_{\alpha}-\boldsymbol{r}|^3}\mathrm d^3 \boldsymbol{r} - Z_{\alpha} e \sum_{\beta(\neq \alpha)} Z_{\beta}e \frac{\boldsymbol{R}_{\alpha}-\boldsymbol{R}_{\beta}〔\boldsymbol{R}_{\alpha}-\boldsymbol{R}_{\beta}|^3}
と書き表すことができる。第1項は電子の電荷密度と電子による電場の積であり、第2項は電荷 Zβ e を持つほかの原子核による電場の効果である。この結果を静電定理 (electrostatic theorem) と呼ぶ。