Hellmann–Feynman theorem

Hellmann–Feynman theorem. In quantum mechanics and molecular physics, the Hellmann–Feynman theorem provides a powerful method for calculating the forces on nuclei within a molecule or crystal directly from the electronic wavefunction. It states that the derivative of the total energy with respect to a parameter, such as a nuclear coordinate, equals the expectation value of the derivative of the Hamiltonian with respect to that same parameter. This result simplifies complex calculations, particularly for geometries and lattice constants, by avoiding explicit differentiation of the many-electron wavefunction itself.

Statement of the theorem

For a non-degenerate, normalized eigenstate \(|\psi_\lambda\rangle\) of a Hamiltonian \(\hat{H}_\lambda\) that depends smoothly on a real parameter \(\lambda\), the theorem states: \[ \frac{d E_\lambda}{d \lambda} = \left\langle \psi_\lambda \left| \frac{\partial \hat{H}_\lambda}{\partial \lambda} \right| \psi_\lambda \right\rangle, \] where \(E_\lambda\) is the corresponding eigenvalue. The parameter \(\lambda\) can represent a nuclear coordinate in the Born–Oppenheimer approximation, an external electric field or magnetic field, or a coupling constant in the Hamiltonian. The theorem is foundational in quantum chemistry for computing properties like force constants and dipole moments. Its validity relies on the wavefunction being an exact eigenstate of the Hamiltonian, a condition strictly met within the variational principle.

Derivation

The derivation begins by considering the expectation value of the Hamiltonian, \(E_\lambda = \langle \psi_\lambda | \hat{H}_\lambda | \psi_\lambda \rangle\). Taking the total derivative with respect to \(\lambda\) and applying the product rule yields terms from the derivative of the wavefunction. Using the fact that \(|\psi_\lambda\rangle\) is an eigenstate, \(\hat{H}_\lambda |\psi_\lambda\rangle = E_\lambda |\psi_\lambda\rangle\), and the normalization condition \(\langle \psi_\lambda | \psi_\lambda \rangle = 1\), the terms involving derivatives of the wavefunction cancel. This cancellation, which relies on the Hermitian nature of the Hamiltonian, leaves only the expectation value of the partial derivative of the operator. This elegant proof is a standard exercise in textbooks on quantum mechanics and quantum chemistry.

Example applications

A primary application is in calculating interatomic forces within the Born–Oppenheimer approximation, where \(\lambda\) is a nuclear coordinate. The force on a nucleus is the negative gradient of the total energy, which the theorem equates to the expectation value of the gradient of the Hamiltonian, essentially the electrostatic force from the electron density. This is exploited in computational methods like density functional theory and Hartree–Fock for geometry optimization. The theorem also simplifies computation of properties like polarizability (derivative with respect to an external electric field) and magnetizability (derivative with respect to a magnetic field). In solid-state physics, it aids in calculating stress and elastic constants for materials like silicon or graphene.

Connections to other theorems

The Hellmann–Feynman theorem is closely related to the Ehrenfest theorem, which deals with the time evolution of expectation values and provides quantum analogs of Newton's laws of motion. Both theorems connect derivatives of expectation values to expectation values of derivatives. It is also a specific instance of the more general quantum mechanical variational principle, as it holds exactly for eigenstates that satisfy the stationary condition. In quantum field theory, analogous relations appear in the form of Schwinger–Dyson equations. Furthermore, the theorem underpins the force theorem in density functional theory, which is used for efficient ab initio molecular dynamics simulations.

Limitations and extensions

The standard theorem fails for degenerate states and for approximate wavefunctions that are not exact eigenstates of the Hamiltonian, such as those from truncated configuration interaction or finite basis set calculations. This leads to so-called "Pulay forces" in quantum chemistry, requiring correction terms. Extensions, like the generalized Hellmann–Feynman theorem, account for wavefunction parameter dependence, as seen in coupled cluster or multiconfigurational self-consistent field methods. Research into its applicability within time-dependent density functional theory and for open quantum systems described by Lindbladian dynamics continues to be an active area in theoretical chemistry and condensed matter physics. Category:Quantum chemistry Category:Quantum mechanics Category:Physics theorems Category:Theoretical chemistry