Helium atom

Last updated
Helium atom
Systematic IUPAC name
Helium [1]
3D model (JSmol)
EC Number
  • 231-168-5
MeSH Helium
PubChem CID
RTECS number
  • MH6520000
UN number 1046
Molar mass 4.002602 g·mol−1
AppearanceColourless gas
Boiling point −269 °C (−452.20 °F; 4.15 K)
126.151-126.155 J K−1 mol−1
V03AN03 ( WHO )
S-phrases (outdated) S9
Except where otherwise noted, data are given for materials in their standard state (at 25 °C [77 °F], 100 kPa).
X mark.svgN  verify  (what is  Yes check.svgYX mark.svgN ?)
Infobox references

A helium atom is an atom of the chemical element helium. Helium is composed of two electrons bound by the electromagnetic force to a nucleus containing two protons along with either one or two neutrons, depending on the isotope, held together by the strong force. Unlike for hydrogen, a closed-form solution to the Schrödinger equation for the helium atom has not been found. However, various approximations, such as the Hartree–Fock method, can be used to estimate the ground state energy and wavefunction of the atom.



Schematic termscheme for Para- and Orthohelium with one electron in ground state 1s and one excited electron. Helium-term-scheme.svg
Schematic termscheme for Para- and Orthohelium with one electron in ground state 1s and one excited electron.

The quantum mechanical description of the helium atom is of special interest, because it is the simplest multi-electron system and can be used to understand the concept of quantum entanglement. The Hamiltonian of helium, considered as a three-body system of two electrons and a nucleus and after separating out the centre-of-mass motion, can be written as

where is the reduced mass of an electron with respect to the nucleus, and are the electron-nucleus distance vectors and . The nuclear charge, is 2 for helium. In the approximation of an infinitely heavy nucleus, we have and the mass polarization term disappears. In atomic units the Hamiltonian simplifies to

It is important to note, that it operates not in normal space, but in a 6-dimensional configuration space . In this approximation (Pauli approximation) the wave function is a second order spinor with 4 components , where the indices describe the spin projection of both electrons (z-direction up or down) in some coordinate system. [2] It has to obey the usual normalization condition . This general spinor can be written as 2x2 matrix and consequently also as linear combination of any given basis of four orthogonal (in the vector-space of 2x2 matrices) constant matrices with scalar function coefficients

as . A convenient basis consists of one anti-symmetric matrix (with total spin , corresponding to a singlet state)

and three symmetric matrices (with total spin , corresponding to a triplet state)

It is easy to show, that the singlet state is invariant under all rotations (a scalar entity), while the triplet can be mapped to an ordinary space vector , with the three components

, and .

Since all spin interaction terms between the four components of in the above (scalar) Hamiltonian are neglected (e.g. an external magnetic field, or relativistic effects, like angular momentum coupling), the four Schrödinger equations can be solved independently. [3]

The spin here only comes into play through the Pauli exclusion principle, which for fermions (like electrons) requires antisymmetry under simultaneous exchange of spin and coordinates


Parahelium is then the singlet state with a symmetric function and orthohelium is the triplet state with an antisymmetric function . If the electron-electron interaction term is ignored, both spatial functions can be written as linear combination of two arbitrary (orthogonal and normalized)

one-electron eigenfunctions :

or for the special cases

of (both electrons have identical quantum numbers, parahelium only): . The total energy (as eigenvalue of ) is then for all cases (independent of the symmetry).

This explains the absence of the state (with ) for orthohelium, where consequently (with ) is the metastable ground state. (A state with the quantum numbers: principal quantum number , total spin , angular quantum number and total angular momentum is denoted by .)

If the electron-electron interaction term is included, the Schrödinger equation is non separable. However, also if is neglected, all states described above (even with two identical quantum numbers, like with ) cannot be written as a product of one-electron wave functions: — the wave function is entangled. One cannot say, particle 1 is in state 1 and the other in state 2, and measurements cannot be made on one particle without affecting the other.

Nevertheless, quite good theoretical descriptions of helium can be obtained within the Hartree–Fock and Thomas–Fermi approximations (see below).

Hartree–Fock method

The Hartree–Fock method is used for a variety of atomic systems. However it is just an approximation, and there are more accurate and efficient methods used today to solve atomic systems. The "many-body problem" for helium and other few electron systems can be solved quite accurately. For example, the ground state of helium is known to fifteen digits. In Hartree–Fock theory, the electrons are assumed to move in a potential created by the nucleus and the other electrons. The Hamiltonian for helium with two electrons can be written as a sum of the Hamiltonians for each electron:

where the zero-order unperturbed Hamiltonian is

while the perturbation term:

is the electron-electron interaction. H0 is just the sum of the two hydrogenic Hamiltonians:


Eni, the energy eigenvalues and , the corresponding eigenfunctions of the hydrogenic Hamiltonian will denote the normalized energy eigenvalues and the normalized eigenfunctions. So:


Neglecting the electron-electron repulsion term, the Schrödinger equation for the spatial part of the two-electron wave function will reduce to the 'zero-order' equation

This equation is separable and the eigenfunctions can be written in the form of single products of hydrogenic wave functions:

The corresponding energies are (in atomic units, hereafter a.u.):

Note that the wave function

An exchange of electron labels corresponds to the same energy . This particular case of degeneracy with respect to exchange of electron labels is called exchange degeneracy. The exact spatial wave functions of two-electron atoms must either be symmetric or antisymmetric with respect to the interchange of the coordinates and of the two electrons. The proper wave function then must be composed of the symmetric (+) and antisymmetric(-) linear combinations:

This comes from Slater determinants.

The factor normalizes . In order to get this wave function into a single product of one-particle wave functions, we use the fact that this is in the ground state. So . So the will vanish, in agreement with the original formulation of the Pauli exclusion principle, in which two electrons cannot be in the same state. Therefore, the wave function for helium can be written as

Where and use the wave functions for the hydrogen Hamiltonian. [lower-alpha 1] For helium, Z = 2 from

where E = 4 a.u. which is approximately 108.8 eV, which corresponds to an ionization potential V = 2 a.u. (54.4 eV). The experimental values are E = 2.90 a.u. (79.0 eV) and V = 0.90 a.u. ( 24.6 eV).

The energy that we obtained is too low because the repulsion term between the electrons was ignored, whose effect is to raise the energy levels. As Z gets bigger, our approach should yield better results, since the electron-electron repulsion term will get smaller.

So far a very crude independent-particle approximation has been used, in which the electron-electron repulsion term is completely omitted. Splitting the Hamiltonian showed below will improve the results:



V(r) is a central potential which is chosen so that the effect of the perturbation is small. The net effect of each electron on the motion of the other one is to screen somewhat the charge of the nucleus, so a simple guess for V(r) is

where S is a screening constant and the quantity Ze is the effective charge. The potential is a Coulomb interaction, so the corresponding individual electron energies are given (in a.u.) by

and the corresponding wave function is given by

If Ze was 1.70, that would make the expression above for the ground state energy agree with the experimental value E0 = 2.903 a.u. of the ground state energy of helium. Since Z = 2 in this case, the screening constant is S = 0.30. For the ground state of helium, for the average shielding approximation, the screening effect of each electron on the other one is equivalent to about of the electronic charge. [5]

Thomas–Fermi method

Not long after Schrödinger developed the wave equation, the Thomas–Fermi model was developed. Density functional theory is used to describe the particle density , and the ground state energy E(N), where N is the number of electrons in the atom. If there are a large number of electrons, the Schrödinger equation runs into problems, because it gets very difficult to solve, even in the atom's ground states. This is where density functional theory comes in. Thomas–Fermi theory gives very good intuition of what is happening in the ground states of atoms and molecules with N electrons.

The energy functional for an atom with N electrons is given by:


The electron density needs to be greater than or equal to 0, , and is convex.

In the energy functional, each term holds a certain meaning. The first term describes the minimum quantum-mechanical kinetic energy required to create the electron density for an N number of electrons. The next term is the attractive interaction of the electrons with the nuclei through the Coulomb potential . The final term is the electron-electron repulsion potential energy. [6]

So the Hamiltonian for a system of many electrons can be written:

For helium, N = 2, so the Hamiltonian is given by:



From the Hartree–Fock method, it is known that ignoring the electron-electron repulsion term, the energy is 8E1 = 109 eV.

The variational method

To obtain a more accurate energy the variational principle can be applied to the electron-electron potential Vee using the wave function


After integrating this, the result is:

This is closer to the experimental value, but if a better trial wave function is used, an even more accurate answer could be obtained. An ideal wave function would be one that doesn't ignore the influence of the other electron. In other words, each electron represents a cloud of negative charge which somewhat shields the nucleus so that the other electron actually sees an effective nuclear charge Z that is less than 2. A wave function of this type is given by:

Treating Z as a variational parameter to minimize H. The Hamiltonian using the wave function above is given by:

After calculating the expectation value of and Vee the expectation value of the Hamiltonian becomes:

The minimum value of Z needs to be calculated, so taking a derivative with respect to Z and setting the equation to 0 will give the minimum value of Z:

This shows that the other electron somewhat shields the nucleus reducing the effective charge from 2 to 1.69. So we obtain the most accurate result yet:

Where again, E1 represents the ionization energy of hydrogen.

By using more complicated/accurate wave functions, the ground state energy of helium has been calculated closer and closer to the experimental value 78.95 eV. [7] The variational approach has been refined to very high accuracy for a comprehensive regime of quantum states by G.W.F. Drake and co-workers [8] [9] [10] as well as J.D. Morgan III, Jonathan Baker and Robert Hill [11] [12] [13] using Hylleraas or Frankowski-Pekeris basis functions. One needs to include relativistic and quantum electrodynamic corrections to get full agreement with experiment to spectroscopic accuracy. [14] [15]

Experimental value of ionization energy

Helium's first ionization energy is 24.587387936(25) eV. [16] This value was derived by experiment. [17] The theoretic value of Helium atom's second ionization energy is 54.41776311(2) eV. [16] The total ground state energy of the helium atom is 79.005151042(40) eV, [16] or 2.90338583(13) Atomic units a.u., which equals 5.80677166 (26) Ry.

See also

Related Research Articles

Dirac equation Relativistic quantum mechanical wave equation

In particle physics, the Dirac equation is a relativistic wave equation derived by British physicist Paul Dirac in 1928. In its free form, or including electromagnetic interactions, it describes all spin-1/2 massive particles such as electrons and quarks for which parity is a symmetry. It is consistent with both the principles of quantum mechanics and the theory of special relativity, and was the first theory to account fully for special relativity in the context of quantum mechanics. It was validated by accounting for the fine details of the hydrogen spectrum in a completely rigorous way.

In physics, the Navier–Stokes equations, named after Claude-Louis Navier and George Gabriel Stokes, describe the motion of viscous fluid substances.

Potential flow velocity field as the gradient of a scalar function

In fluid dynamics, potential flow describes the velocity field as the gradient of a scalar function: the velocity potential. As a result, a potential flow is characterized by an irrotational velocity field, which is a valid approximation for several applications. The irrotationality of a potential flow is due to the curl of the gradient of a scalar always being equal to zero.

In 1851, George Gabriel Stokes derived an expression, now known as Stokes' law, for the frictional force – also called drag force – exerted on spherical objects with very small Reynolds numbers in a viscous fluid. Stokes' law is derived by solving the Stokes flow limit for small Reynolds numbers of the Navier–Stokes equations.

The Klein–Gordon equation is a relativistic wave equation, related to the Schrödinger equation. It is second-order in space and time and manifestly Lorentz-covariant. It is a quantized version of the relativistic energy–momentum relation. Its solutions include a quantum scalar or pseudoscalar field, a field whose quanta are spinless particles. Its theoretical relevance is similar to that of the Dirac equation. Electromagnetic interactions can be incorporated, forming the topic of scalar electrodynamics, but because common spinless particles like the pions are unstable and also experience the strong interaction the practical utility is limited.

In mathematics, a Green's function of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions is its impulse response.

In mathematics, the covariant derivative is a way of specifying a derivative along tangent vectors of a manifold. Alternatively, the covariant derivative is a way of introducing and working with a connection on a manifold by means of a differential operator, to be contrasted with the approach given by a principal connection on the frame bundle – see affine connection. In the special case of a manifold isometrically embedded into a higher-dimensional Euclidean space, the covariant derivative can be viewed as the orthogonal projection of the Euclidean directional derivative onto the manifold's tangent space. In this case the Euclidean derivative is broken into two parts, the extrinsic normal component and the intrinsic covariant derivative component.

Bloch sphere geometrical representation of the pure state space of a two-level quantum mechanical system

In quantum mechanics, the Bloch sphere is a geometrical representation of the pure state space of a two-level quantum mechanical system (qubit), named after the physicist Felix Bloch.

In mathematics, Green's identities are a set of three identities in vector calculus relating the bulk with the boundary of a region on which differential operators act. They are named after the mathematician George Green, who discovered Green's theorem.

The rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

In atomic physics, the electron magnetic moment, or more specifically the electron magnetic dipole moment, is the magnetic moment of an electron caused by its intrinsic properties of spin and electric charge. The value of the electron magnetic moment is approximately −9.284764×10−24 J/T. The electron magnetic moment has been measured to an accuracy of 7.6 parts in 1013.

In solid-state physics, the tight-binding model is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon superposition of wave functions for isolated atoms located at each atomic site. The method is closely related to the LCAO method used in chemistry. Tight-binding models are applied to a wide variety of solids. The model gives good qualitative results in many cases and can be combined with other models that give better results where the tight-binding model fails. Though the tight-binding model is a one-electron model, the model also provides a basis for more advanced calculations like the calculation of surface states and application to various kinds of many-body problem and quasiparticle calculations.

The following are important identities involving derivatives and integrals in vector calculus.

The theoretical and experimental justification for the Schrödinger equation motivates the discovery of the Schrödinger equation, the equation that describes the dynamics of nonrelativistic particles. The motivation uses photons, which are relativistic particles with dynamics described by Maxwell's equations, as an analogue for all types of particles.

In quantum mechanics, the Pauli equation or Schrödinger–Pauli equation is the formulation of the Schrödinger equation for spin-½ particles, which takes into account the interaction of the particle's spin with an external electromagnetic field. It is the non-relativistic limit of the Dirac equation and can be used where particles are moving at speeds much less than the speed of light, so that relativistic effects can be neglected. It was formulated by Wolfgang Pauli in 1927.

The intent of this article is to highlight the important points of the derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids.

Stokes stream function

In fluid dynamics, the Stokes stream function is used to describe the streamlines and flow velocity in a three-dimensional incompressible flow with axisymmetry. A surface with a constant value of the Stokes stream function encloses a streamtube, everywhere tangential to the flow velocity vectors. Further, the volume flux within this streamtube is constant, and all the streamlines of the flow are located on this surface. The velocity field associated with the Stokes stream function is solenoidal—it has zero divergence. This stream function is named in honor of George Gabriel Stokes.

Gyrokinetics is a theoretical framework to study plasma behavior on perpendicular spatial scales comparable to the gyroradius and frequencies much lower than the particle cyclotron frequencies. These particular scales have been experimentally shown to be appropriate for modeling plasma turbulence. The trajectory of charged particles in a magnetic field is a helix that winds around the field line. This trajectory can be decomposed into a relatively slow motion of the guiding center along the field line and a fast circular motion, called gyromotion. For most plasma behavior, this gyromotion is irrelevant. Averaging over this gyromotion reduces the equations to six dimensions rather than the seven. Because of this simplification, gyrokinetics governs the evolution of charged rings with a guiding center position, instead of gyrating charged particles.

In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

Mild-slope equation The combined effects of diffraction and refraction for water waves propagating over variable depth and with lateral boundaries

In fluid dynamics, the mild-slope equation describes the combined effects of diffraction and refraction for water waves propagating over bathymetry and due to lateral boundaries—like breakwaters and coastlines. It is an approximate model, deriving its name from being originally developed for wave propagation over mild slopes of the sea floor. The mild-slope equation is often used in coastal engineering to compute the wave-field changes near harbours and coasts.


  1. For n = 1, l = 0 and m = 0, the wavefunction in a spherically symmetric potential for a hydrogen electron is . [4] In atomic units, the Bohr radius equals 1, and the wavefunction becomes .
  1. "Helium - PubChem Public Chemical Database". The PubChem Project. USA: National Center for Biotechnology Information.
  2. P. Rennert, H. Schmiedel, C. Weißmantel. "Kleine Enzyklopädie Physik", VEB Bibliographisches Institut Leipzig, 1988, 192-194.
  3. L.D. Landau, E.M. Lifschitz. "Lehrbuch der Theoretischen Physik", Bd. III (Quantenmechanik), Akademie-Verlag, Berlin 1971, Kap. IX, pp. 218
  4. "Hydrogen Wavefunctions". Hyperphysics. Archived from the original on 1 February 2014.
  5. B.H. Bransden and C.J. Joachain's Physics of Atoms and Molecules 2nd edition Pearson Education, Inc
  6. http://www.physics.nyu.edu/LarrySpruch/Lieb.pdf
  7. David I. Griffiths Introduction to Quantum Mechanics Second edition year 2005 Pearson Education, Inc
  8. G.W.F. Drake and Zong-Chao Van (1994). "Variational eigenvalues for the S states of helium", Chem. Phys. Lett.229 486–490. [ dead link ]
  9. Zong-Chao Yan and G. W. F. Drake (1995). "High Precision Calculation of Fine Structure Splittings in Helium and He-Like Ions", Phys. Rev. Lett.74, 4791–4794.
  10. G.W.F. Drake, (1999). "High precision theory of atomic helium", Phys. Scr.T83, 83–92.
  11. J.D. Baker, R.N. Hill, and J.D. Morgan III (1989), "High Precision Calculation of Helium Atom Energy Levels", in AIP ConferenceProceedings 189, Relativistic, Quantum Electrodynamic, and Weak Interaction Effects in Atoms (AIP, New York),123
  12. Jonathan D. Baker, David E. Freund, Robert Nyden Hill, and John D. Morgan III (1990). "Radius of convergence and analytic behavior of the 1/Z expansion", Physical Review A41, 1247.
  13. Scott, T. C.; Lüchow, A.; Bressanini, D.; Morgan, J. D. III (2007). "The Nodal Surfaces of Helium Atom Eigenfunctions" (PDF). Phys. Rev. A . 75 (6): 060101. Bibcode:2007PhRvA..75f0101S. doi:10.1103/PhysRevA.75.060101. hdl:11383/1679348.
  14. G.W.F. Drake and Z.-C. Yan (1992), Phys. Rev. A46,2378-2409. .
  15. G.W.F. Drake (2006). "Springer Handbook of Atomic, molecular, and Optical Physics", Edited by G.W.F. Drake (Springer, New York), 199-219.
  16. 1 2 3 Kramida, A., Ralchenko, Yu., Reader, J., and NIST ASD Team. "NIST Atomic Spectra Database Ionization Energies Data". Gaithersburg, MD: NIST.CS1 maint: uses authors parameter (link)
  17. D. Z. Kandula, C. Gohle, T. J. Pinkert, W. Ubachs, and K. S. E. Eikema (2010). "Extreme Ultraviolet Frequency Comb Metrology". Phys. Rev. Lett. 105 (6): 063001. arXiv: 1004.5110 . Bibcode:2010PhRvL.105f3001K. doi:10.1103/PhysRevLett.105.063001. PMID   20867977.CS1 maint: uses authors parameter (link)