ヘルマン-ファインマンの定理

ヘルマン-ファインマンの定理（-のていり、英: Hellmann-Feynman theorem）とは、量子力学において、パラメータ依存したハミルトニアンとそのエネルギー固有値に関する定理である。定理の名は、ドイツの物理学者H. ヘルマンと米国の物理学者R. P. ファインマンに因む。定理を最初に明示的な形で表したのは、P. Güttingerであるが、W. パウリやヘルマンの論文にも記されている。特にヘルマンは分子への適用に向けて、変分形式で表現した。また、1939年に当時、マサチューセッツ工科大学の学生であったファインマンは、この定理を示すともに、化学結合した原子において、電子及び他の原子核が原子核に及ぼす力は古典的な静電力として、扱えることを示した。「ヘルマン-ファインマンの定理」の名が定着したのは、J. C. スレイターがその著書の中でその名で呼んだことによる。

定理

系のハミルトニアンが、あるパラメータλに依存するとして、それをH (λ)と表現する。H (λ)の固有状態|ψ_λ>があって、 H (λ) |ψ_λ> = E (λ) |ψ_λ>及び規格化条件<ψ_λ|ψ_λ>= 1 が満足されるとする。このとき、

${ d E (\lambda) \over { d \lambda } } = \left \langle \psi_\lambda \left | {d H(\lambda) \over {d \lambda} } \right | \psi_\lambda \right \rangle$

が成り立つ。これがヘルマン-ファインマンの定理の主張である。

ここで、パラメータλが、原子位置座標R_αの場合、ヘルマン-ファインマン力となる。

証明

$\begin{align} { d E (\lambda) \over { d \lambda } } &= \frac{d}{d\lambda} \left \langle \psi_\lambda \left| H(\lambda) \right |\psi_\lambda \right \rangle \\ &= \bigg ( \frac{d}{d\lambda} \bigl \langle \psi_\lambda \bigl |\bigg )H(\lambda) \bigr |\psi_\lambda \bigr \rangle +\bigg \langle \psi_\lambda \bigg | \frac{dH(\lambda)}{d\lambda} \bigg | \psi_\lambda \bigg \rangle +\bigl \langle \psi_\lambda \bigl |H(\lambda)\bigg ( \frac{d}{d\lambda} \bigr | \psi_\lambda \bigr \rangle \bigg ) \\ &= E(\lambda) \bigg ( \frac{d}{d\lambda}\bigl \langle \psi_\lambda \bigl | \bigg ) \bigr | \psi_\lambda \bigr \rangle +\bigg \langle \psi_\lambda \bigg|\frac{dH(\lambda)}{d\lambda} \bigg | \psi_\lambda \bigg \rangle +E(\lambda) \bigl \langle \psi_\lambda \bigl | \bigg ( \frac{d}{d\lambda}\bigr |\psi_\lambda \bigr \rangle \bigg ) \\ &= E(\lambda)\frac{d}{d\lambda}\langle \psi_\lambda | \psi_\lambda \rangle +\bigg \langle \psi_\lambda \bigg | \frac{dH(\lambda)}{d\lambda} \bigg | \psi_\lambda \bigg \rangle. \end{align}$

ここで、|ψ_λ>の規格化を<ψ_λ|ψ_λ>= 1 と選んであるので、d<ψ_λ|ψ_λ>/dλ=0である。よって、

${ d E (\lambda) \over { d \lambda } } = \left \langle \psi_\lambda \left | {d H(\lambda) \over {d \lambda} } \right | \psi_\lambda \right \rangle$

が得られる。

応用

分子内力と静電定理

定理の応用の１つとして、分子内力（intramolecular force）の計算がある。ファインマンは1939年の「分子内の力（Forces in Molecules）」と題する論文の中で、ヘルマン-ファインマンの定理の証明を与えるとともに、分子や固体原子において、原子核に働く量子論的な力は、電子雲と他の原子核からの古典的な静電力として扱えることを示した。

電子及び位置が固定された原子核からなる系において、ポテンシャルV を原子核の位置座標で微分したものは、原子核に働く力に相当する。例えば、電子の位置をr_i=（x_i, y_i, z_i）（i=1,…,n）とし、原子核の位置をR_α=（X_α, Y_α, Z_α）（α=1,…,m）とする。このとき、 系のハミルトニアンH はボルン-オッペンハイマー近似により、運動エネルギー

$T=\sum_{i=1}^n -\frac{\hbar^2}{2m} \nabla_i^2$

とポテンシャルエネルギー

$\begin{align} V &= \sum_{i>j} \frac{e^2}{|\mathbf{r}_i-\mathbf{r}_j|} +\sum_{\alpha > \beta} \frac{Z_{\alpha}Z_{\beta}e^2}{|\mathbf{R}_{\alpha}-\mathbf{R}_{\beta}|} +\sum_{i, \alpha}\frac{-Z_{\alpha}e^2}{|\mathbf{r}_i-\mathbf{R}_{\alpha}|} \\ &= \sum_{i>j} V_{ij}+\sum_{\alpha > \beta} V_{\alpha \beta}+\sum_{i, \alpha}V_{i\alpha} \end{align}$

の和T+V で表される。このとき、パラメータλとして、原子核の位置座標R_αをとったときに、その導関数

$\mathbf{F}_{\alpha}= -\frac{\partial V}{\partial \mathbf{\lambda}}= -\frac{\partial H}{\partial \mathbf{\lambda}}$

で定義されるF_α=（F_α^X, F_α^Y, F_α^Z）が、古典論では原子核に働く力となる。一方、量子論では、H_λ|ψ_λ>=E(λ)|ψ_λ>を満たす固有状態|ψ_λ>により、

$\begin{align} \mathbf{F}_{\alpha} &= - \langle \psi_{\mathbf{\lambda}}|\frac{\partial H}{\partial \mathbf{\lambda}}|\psi_{\mathbf{\lambda}} \rangle \\ &= - \int \psi^* (\mathbf{q}_1, \cdots, \mathbf{q}_n) \frac{\partial H}{\partial \mathbf{\lambda} } \psi(\mathbf{q}_1, \cdots, \mathbf{q}_n) \, dV \end{align}$

が対応する。ここで2行目の波動関数ψにおける座標q_iはi番目の電子の位置座標r_iとスピン座標ζ_iを合わせたものであり、 dV =dq₁…dq_nである。ファインマンの論文以前は、分子の量子力学では、これを

$\mathbf{F'}_{\alpha}=- \frac{d E(\mathbf{\lambda})}{d \mathbf{\lambda}}$

とエネルギー固有値の微分で計算するのが、一般的であった。ファインマンはヘルマン-ファインマンの定理によって、F_αとF_α'が等しいことを示した。

実際の計算において、微分係数dE(λ)/dλで与えられるF_α'を求めるには、エネルギー固有値のパラメータλ=R_αについての依存性の傾きを計算することになり、複数のパラメータ値に対して固有値問題を解く必要がある。一方、直接的に<ψ_λ|dH/dλ|ψ_λ>で与えられるF_αでは、計算の労力を減らすことができる。

F_αの具体形については、波動関数による存在確率から定まる各電子の電荷密度

$\rho_i(\mathbf{r}_i) = e \int | \psi(\mathbf{q}_1, \cdots, \mathbf{q}_n)|^2 d\mathbf{q}_1 \cdots d\mathbf{q}_{i-1} d \zeta_i d\mathbf{q}_{i+1} \cdots d\mathbf{q}_n$

の和として、全電子の電荷密度ρ(r)を、

$\rho(\mathbf{r})=\sum_{i=1}^n \rho_i(\mathbf{r})$

で導入すれば、

$\mathbf{F}_{\alpha}= Z_{\alpha}e \int \rho(\mathbf{r}) \frac{\mathbf{R}_{\alpha}-\mathbf{r}}{|\mathbf{R}_{\alpha}-\mathbf{r}|^3}d^3 \mathbf{r} - Z_{\alpha} e \sum_{\beta(\neq \alpha)} Z_{\beta}e \frac{\mathbf{R}_{\alpha}-\mathbf{R}_{\beta}}{|\mathbf{R}_{\alpha}-\mathbf{R}_{\beta}|^3}$

と書き表すことができる。第１項は電子の電荷密度と電子による電場の積であり、第2項は電荷Z_βeを持つほかの原子核による電場の効果である。この結果を静電定理（electrostatic theorem）と呼ぶ。

参考文献

• Richard M. Martin, Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press (2004) ISBN 978-0521782852 ；R. M. マーチン (著)、寺倉清之、寺倉郁子、善甫康成(翻訳) 『物質の電子状態 上』 : シュプリンガー・ジャパン株式会社

関連項目

• 量子力学
• 量子化学 - ヘルマン-ファインマン力