In solid-state physics and quantum chemistry, the Hartree equations or the self-consistent field approximation are a set of non-linear equations used to study many-electron systems inside a metal. These quantum mechanical equations are self-consistent, meaning that the solutions can be found by iteration. These approximation is the result of a mean-field theory that describes one electron interacting with a field that is the result of averaging the position of the rest of electrons. The equations are named after Douglas Hartree, who introduced them in 1927. Hartree method is one of the main ingredients of Hartree–Fock method, which improves on Hartree equations by including the exchange interaction.
History In 1927, a year after the publication of the Schrödinger equation, Hartree formulated what are now known as the Hartree equations for atoms, using the concept of self-consistency that Robert Bruce Lindsay had introduced in his study of many electron systems in the context of Bohr theory. Hartree assumed that the nucleus together with the electrons formed a spherically symmetric field. The charge distribution of each electron was the solution of the Schrödinger equation for an electron in a potential v ( r ) {\displaystyle v(r)} , derived from the field. Self-consistency required that the final field, computed from the solutions, was self-consistent with the initial field, and he thus called his method the self-consistent field method. In order to solve the equation of an electron in a spherical potential, Hartree first introduced atomic units to eliminate physical constants. Then he converted the Laplacian from Cartesian to spherical coordinates to show that the solution was a product of a radial function P ( r ) / r {\displaystyle P(r)/r} and a spherical harmonic with an angular quantum number ℓ {\displaystyle \ell } , namely ψ = ( 1 / r ) P ( r ) S ℓ ( θ , ϕ ) {\displaystyle \psi =(1/r)P(r)S_{\ell }(\theta ,\phi )} . The equation for the radial function was
d 2 P ( r ) d r 2 + { 2 [ E − v ( r ) ] − ℓ ( ℓ + 1 ) r 2 } P ( r ) = 0. {\displaystyle {\frac {\mathrm {d} ^{2}P(r)}{\mathrm {d} r^{2}}}+\left\{2[E-v(r)]-{\frac {\ell (\ell +1)}{r^{2}}}\right\}P(r)=0.}
Hartree product The wavefunction which describes all of the electrons, Ψ {\displaystyle \Psi } , is almost always too complex to calculate directly. Hartree's original method was to first calculate the solutions to Schrödinger's equation for individual electrons 1, 2, 3, . . . {\displaystyle ...} , p, in the states α , β , γ , . . . , π {\displaystyle \alpha ,\beta ,\gamma ,...,\pi } , which yields individual solutions: ψ α ( x 1 ) , ψ β ( x 2 ) , ψ γ ( x 3 ) , . . . , ψ π ( x p ) {\displaystyle \psi _{\alpha }(\mathbf {x} _{1}),\psi _{\beta }(\mathbf {x} _{2}),\psi _{\gamma }(\mathbf {x} _{3}),...,\psi _{\pi }(\mathbf {x} _{p})} . Since each ψ {\displaystyle \psi } is a solution to the Schrödinger equation by itself, their product should at least approximate a solution. This simple method of combining the wavefunctions of the individual electrons is known as the Hartree product:
… excerpt ends here. Continue reading the full article.
