1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
---
title: "Hellmann-Feynman theorem"
firstLetter: "H"
publishDate: 2021-11-29
categories:
- Physics
- Quantum mechanics
date: 2021-11-29T10:30:37+01:00
draft: false
markup: pandoc
---
# Hellmann-Feynman theorem
Consider the time-independent Schrödinger equation,
where the Hamiltonian $\hat{H}$ depends on a general parameter $\lambda$,
whose meaning or type we will not specify:
$$\begin{aligned}
\hat{H}(\lambda) \ket{\psi_n(\lambda)}
= E_n(\lambda) \ket{\psi_n(\lambda)}
\end{aligned}$$
Assuming all eigenstates $\ket{\psi_n}$ are normalized,
this gives us the following basic relation:
$$\begin{aligned}
\matrixel{\psi_m}{\hat{H}}{\psi_n}
= E_n \braket{\psi_m}{\psi_n}
= \delta_{mn} E_n
\end{aligned}$$
We differentiate this with respect to $\lambda$,
which could be a scalar or a vector.
This yields:
$$\begin{aligned}
\delta_{mn} \nabla_\lambda E_n
&= \nabla_\lambda \matrixel{\psi_m}{\hat{H}}{\psi_n}
\\
&= \matrixel{\nabla_\lambda \psi_m}{\hat{H}}{\psi_n}
+ \matrixel{\psi_m}{\nabla_\lambda \hat{H}}{\psi_n}
+ \matrixel{\psi_m}{\hat{H}}{\nabla_\lambda \psi_n}
\\
&= E_m \braket{\psi_m}{\nabla_\lambda \psi_n} + E_n \braket{\nabla_\lambda \psi_m}{\psi_n} + \matrixel{\psi_m}{\nabla_\lambda \hat{H}}{\psi_n}
\end{aligned}$$
In order to simplify this,
we differentiate the orthogonality relation
$\braket{\psi_m}{\psi_n} = \delta_{mn}$,
which ends up telling us that
$\braket{\nabla_\lambda \psi_m}{\psi_n} = - \braket{\psi_m}{\nabla_\lambda \psi_n}$:
$$\begin{aligned}
0
= \nabla_\lambda \delta_{mn}
= \nabla_\lambda \braket{\psi_m}{\psi_n}
= \braket{\nabla_\lambda \psi_m}{\psi_n} + \braket{\psi_m}{\nabla_\lambda \psi_n}
\end{aligned}$$
Using this result to replace $\braket{\nabla_\lambda \psi_m}{\psi_n}$
in the previous equation leads to:
$$\begin{aligned}
\delta_{mn} \nabla_\lambda E_n
&= (E_m - E_n) \braket{\psi_m}{\nabla_\lambda \psi_n} + \matrixel{\psi_m}{\nabla_\lambda \hat{H}}{\psi_n}
\end{aligned}$$
For $m = n$, we therefore arrive at the **Hellmann-Feynman theorem**,
which is useful when doing numerical calculations
to minimize energies with respect to $\lambda$:
$$\begin{aligned}
\boxed{
\nabla_\lambda E_n
= \matrixel{\psi_n}{\nabla_\lambda \hat{H}}{\psi_n}
}
\end{aligned}$$
While for $m \neq n$, we get the **Epstein generalization**
of the Hellmann-Feynman theorem, which is for example relevant for
the [Berry phase](/know/concept/berry-phase/):
$$\begin{aligned}
\boxed{
(E_n - E_m) \braket{\psi_m}{\nabla_\lambda \psi_n}
= \matrixel{\psi_m}{\nabla_\lambda \hat{H}}{\psi_n}
}
\end{aligned}$$
## References
1. G. Grosso, G.P. Parravicini,
*Solid state physics*,
2nd edition, Elsevier.
|