Skip to content

Atomic Units and Scales

Hartree atomic units are a unit system adapted to nonrelativistic electrons bound by Coulomb forces. They turn the characteristic electron mass, charge, action, length, and energy into quantities of order unity, so that an equation displays its atomic structure instead of repeatedly displaying SI conversion factors.

This page derives those units, builds scale intuition, and gives a safe conversion workflow. The AMO Atomic Units reference provides the versioned conversion tables and Hartree–Rydberg diagnostics; the site-wide Atomic Units translator carries the convention across volumes. Units and Constants states the site-wide reporting policy. Atomic units simplify calculations, but measured fields, frequencies, wavelengths, rates, and uncertainty budgets should usually be reported in laboratory units as well.

Consider the fixed-nucleus, nonrelativistic hydrogen Hamiltonian in SI units,

H=−ℏ22me∇r2−e24πϵ0r,H=-\frac{\hbar^2}{2m_e}\nabla_r^2 -\frac{e^2}{4\pi\epsilon_0r},

where e>0e>0 is the elementary charge magnitude. Introduce a dimensionless coordinate ρ\boldsymbol\rho through

r=a0ρ.\mathbf r=a_0\boldsymbol\rho.

Choose a0a_0 so that the kinetic and Coulomb coefficients define the same energy scale:

ℏ2mea02=e24πϵ0a0≡Eh.\frac{\hbar^2}{m_ea_0^2} =\frac{e^2}{4\pi\epsilon_0a_0} \equiv E_{\mathrm h}.

Solving these relations gives the Bohr radius and Hartree energy,

a0=4πϵ0ℏ2mee2,Eh=e24πϵ0a0,Eh=ℏ2mea02.\begin{aligned} a_0&=\frac{4\pi\epsilon_0\hbar^2}{m_ee^2},\\ E_{\mathrm h}&=\frac{e^2}{4\pi\epsilon_0a_0},\\ E_{\mathrm h}&=\frac{\hbar^2}{m_ea_0^2}. \end{aligned}

After dividing by EhE_{\mathrm h}, the Schrödinger equation becomes

HEh=−12∇ρ2−1ρ.\frac{H}{E_{\mathrm h}} =-\frac12\nabla_{\rho}^2-\frac{1}{\rho}.

Nothing physical has been discarded. The constants have been absorbed into the definitions of length and energy. The dimensionless eigenvalues and wavefunctions can be converted back after the calculation.

Hartree atomic units set

ℏ=me=e=4πϵ0=1.\hbar=m_e=e=4\pi\epsilon_0=1.

This statement is shorthand for measuring action in units of ℏ\hbar, mass in units of mem_e, charge magnitude in units of ee, and electromagnetic quantities in the corresponding rationalized atomic units. It does not mean the electron charge is positive: in this convention the electron has charge −1-1.

The convention also does not set the speed of light to one. Instead,

c=α−1c=\alpha^{-1}

in atomic units, where α\alpha is the fine-structure constant. This fact makes powers of α\alpha visible as the small parameters controlling relativistic corrections.

The base choices generate all other units by dimensional analysis. Selected 2022 CODATA values, rounded here for practical use, are:

QuantityAtomic unitApproximate SI value
lengtha0a_05.291 772 105×10−11 m5.291\,772\,105\times10^{-11}\ \mathrm m
energyEhE_{\mathrm h}4.359 744 722×10−18 J4.359\,744\,722\times10^{-18}\ \mathrm J
energyEhE_{\mathrm h}27.211 386 246 eV27.211\,386\,246\ \mathrm{eV}
timet0=ℏ/Eht_0=\hbar/E_{\mathrm h}2.418 884 327×10−17 s2.418\,884\,327\times10^{-17}\ \mathrm s
velocityv0=a0/t0=αcv_0=a_0/t_0=\alpha c2.187 691 262×106 m s−12.187\,691\,262\times10^6\ \mathrm{m\,s^{-1}}
dipole momentd0=ea0d_0=ea_08.478 353 620×10−30 C m8.478\,353\,620\times10^{-30}\ \mathrm{C\,m}
electric fieldF0=Eh/(ea0)F_0=E_{\mathrm h}/(ea_0)5.142 206 751×1011 V m−15.142\,206\,751\times10^{11}\ \mathrm{V\,m^{-1}}

Parenthesized standard uncertainties and additional derived units are available from the NIST/CODATA constants database. Do not silently combine last digits from different CODATA adjustments.

A length reported as r=3.2r=3.2 a.u. means

r=3.2a0.r=3.2a_0.

The corresponding momentum unit is

p0=ℏa0=meαc.p_0=\frac{\hbar}{a_0}=m_e\alpha c.

Because position and momentum use reciprocal scales, a wavefunction localized to order a0a_0 naturally contains momenta of order ℏ/a0\hbar/a_0.

The Hartree energy is the Coulomb energy of two elementary charge magnitudes separated by one Bohr radius. It is also twice the infinite-mass Rydberg energy. The force unit is

Fforce,0=Eha0.F_{\mathrm{force},0}=\frac{E_{\mathrm h}}{a_0}.

The symbol FF is also widely used for electric field strength in AMO physics. Context and units must distinguish force from field; this page uses F0F_0 only for the atomic unit of electric field.

For an infinitely heavy hydrogenic nucleus, the nonrelativistic bound energies are

En=−Z22n2Eh.E_n=-\frac{Z^2}{2n^2}E_{\mathrm h}.

The ground-state binding magnitude is therefore Eh/2E_{\mathrm h}/2. This quantity is the infinite-mass Rydberg energy:

Ry⁡=Eh2.\operatorname{Ry}=\frac{E_{\mathrm h}}{2}.

Two unit conventions consequently coexist:

  • In Hartree atomic units, the hydrogen ground-state energy is −1/2-1/2.
  • In Rydberg units, the same energy is −1-1.

Electronic-structure codes may report either convention, particularly for kinetic-energy operators, pseudopotentials, and plane-wave cutoffs. The label “a.u.” alone is insufficient when the factor of two matters. Record whether the energy unit is EhE_{\mathrm h} or Ry⁡\operatorname{Ry}.

The Rydberg constant R∞R_\infty is a spectroscopic inverse-length constant related by

Eh=2hcR∞.E_{\mathrm h}=2hcR_\infty.

For a real isotope, reduced mass shifts the hydrogenic scale. Do not identify hcR∞hcR_\infty with an exact measured ionization energy without applying recoil and higher-order corrections.

Energy, Frequency, Angular Frequency, and Wavenumber

Section titled “Energy, Frequency, Angular Frequency, and Wavenumber”

An energy interval may be represented in several equivalent ways:

ΔE=hν,ΔE=ℏω,ΔE=hcν~.\begin{aligned} \Delta E&=h\nu,\\ \Delta E&=\hbar\omega,\\ \Delta E&=hc\widetilde\nu. \end{aligned}

Here ν\nu is ordinary frequency in cycles per second, ω=2πν\omega=2\pi\nu is angular frequency, and ν~=1/λ\widetilde\nu=1/\lambda is spectroscopic wavenumber. In spectroscopy, ν~\widetilde\nu is commonly reported in cm−1\mathrm{cm}^{-1}.

One Hartree corresponds to

Ehh=6.579 683 9205×1015 Hz,Ehℏ=4.134 137 3335×1016 s−1,Ehhc=219 474.6314 cm−1,EhkB=315 775.025 K.\begin{aligned} \frac{E_{\mathrm h}}{h} &=6.579\,683\,9205\times10^{15}\ \mathrm{Hz},\\ \frac{E_{\mathrm h}}{\hbar} &=4.134\,137\,3335\times10^{16}\ \mathrm{s}^{-1},\\ \frac{E_{\mathrm h}}{hc} &=219\,474.6314\ \mathrm{cm}^{-1},\\ \frac{E_{\mathrm h}}{k_B} &=315\,775.025\ \mathrm K. \end{aligned}

The second line is an angular frequency even when the radian is treated as dimensionless. Writing E/ℏE/\hbar and labeling the result “Hz” introduces a factor-of-2π2\pi error. Writing E/(hc)E/(hc) as a wavenumber and then multiplying by 2π2\pi without changing convention creates the same problem in another form.

The time unit

t0=ℏEht_0=\frac{\hbar}{E_{\mathrm h}}

is the inverse Hartree angular-frequency scale. If a time-dependent calculation returns t=1000t=1000 a.u., the physical duration is about 24.19 fs24.19\ \mathrm{fs}. A phase factor has the same form in either system:

exp⁡ ⁣(−iEtℏ)=exp⁡(−iEautau).\exp\!\left(-\frac{iEt}{\hbar}\right) =\exp(-iE_{\mathrm{au}}t_{\mathrm{au}}).

The right-hand expression is valid only when energy is in Hartree and time is in ℏ/Eh\hbar/E_{\mathrm h}.

Electric Fields, Dipoles, and Polarizabilities

Section titled “Electric Fields, Dipoles, and Polarizabilities”

For a dipole in an electric field,

HE=−d⋅F.H_E=-\mathbf d\cdot\mathbf F.

Choosing d0=ea0d_0=ea_0 and F0=Eh/(ea0)F_0=E_{\mathrm h}/(ea_0) makes their product exactly one Hartree:

d0F0=Eh.d_0F_0=E_{\mathrm h}.

Thus a dipole matrix element dmathrmaud_{mathrm{au}} in a field FmathrmauF_{mathrm{au}} produces the characteristic interaction energy

ΔEau∼dauFau.\Delta E_{\mathrm{au}} \sim d_{\mathrm{au}}F_{\mathrm{au}}.

The atomic field unit is extremely large because it is the field that changes an electron’s potential energy by one Hartree across one Bohr radius. Typical laboratory fields are often small in atomic units, even when they strongly mix nearly degenerate states.

The atomic dipole unit is

ea0≃2.54174647 D,ea_0\simeq2.54174647\ \mathrm D,

where D\mathrm D denotes the debye. The polarizability unit follows from the quadratic Stark energy −αEF2/2-\alpha_EF^2/2:

αE,0=e2a02Eh=4πϵ0a03.\alpha_{E,0} =\frac{e^2a_0^2}{E_{\mathrm h}} =4\pi\epsilon_0a_0^3.

The factor 4πϵ04\pi\epsilon_0 matters when converting a polarizability reported as a volume in a03a_0^3 to SI. Atomic, Gaussian, and SI polarizability conventions should not be mixed by dimensional appearance alone.

The Fine-Structure Constant as a Scale Parameter

Section titled “The Fine-Structure Constant as a Scale Parameter”

Atomic units make the nonrelativistic Coulomb problem order unity while leaving relativity parametrically visible. Since v0=αcv_0=\alpha c,

v0c=α.\frac{v_0}{c}=\alpha.

The electron rest energy in Hartree is

mec2=α−2Eh.m_ec^2=\alpha^{-2}E_{\mathrm h}.

A nonrelativistic electronic energy is typically of order EhE_{\mathrm h}, while leading fine-structure corrections are often of relative order α2\alpha^2 for light atoms. For a hydrogenic ion, the orbital velocity and correction parameter grow with ZZ:

vc∼Zα,ΔEfineEgross∼(Zα)2.\frac{v}{c}\sim Z\alpha, \qquad \frac{\Delta E_{\mathrm{fine}}}{E_{\mathrm{gross}}} \sim(Z\alpha)^2.

These estimates organize an expansion; they do not replace a state-specific calculation. Coefficients, cancellations, finite nuclear size, and electron correlation can control the actual uncertainty.

Setting me=1m_e=1 does not set every mass to one. Nuclear masses are expressed in electron-mass units. For a nucleus of mass MM, the one-electron reduced mass is

μ=meMme+M,\mu=\frac{m_eM}{m_e+M},

or, in atomic units,

μau=Mau1+Mau.\mu_{\mathrm{au}} =\frac{M_{\mathrm{au}}}{1+M_{\mathrm{au}}}.

The correction 1−μmathrmau1-\mu_{mathrm{au}} is of order me/Mm_e/M. Molecular electronic energies remain naturally Hartree-scale, while nuclear vibration and rotation inherit the larger nuclear masses and become smaller. A common hierarchy is

Eelectronic≫Evibrational≫Erotational,E_{\mathrm{electronic}} \gg E_{\mathrm{vibrational}} \gg E_{\mathrm{rotational}},

but the ratios depend on the potential surface, equilibrium geometry, and isotopic masses. This scale separation motivates the Born–Oppenheimer approximation; it does not make nonadiabatic coupling identically zero.

When to Use SI and When to Use Atomic Units

Section titled “When to Use SI and When to Use Atomic Units”

Use atomic units when they clarify internal electronic structure:

  • deriving or solving Coulomb Hamiltonians;
  • reporting orbital energies, matrix elements, and electronic-structure convergence;
  • comparing atomic or molecular calculations that declare the same convention;
  • identifying powers of α\alpha, mass ratios, and scaled field strengths.

Use SI-compatible or laboratory units when connecting to apparatus and metrology:

  • wavelengths, ordinary frequencies, angular frequencies, and linewidths;
  • electric and magnetic fields delivered by an instrument;
  • pulse durations, intensities, powers, temperatures, and pressures;
  • transition rates, cross sections, densities, and count rates;
  • calibrated measurements and uncertainty budgets.

A strong presentation often gives both: the calculation in atomic units and the final observable in eV, cm−1\mathrm{cm}^{-1}, Hz, nm, V/m, or seconds. The unit system should be declared near the equation, not left for the reader to infer from magnitude.

For any quantity QQ, write

Q=QauQ0,Q=Q_{\mathrm{au}}Q_0,

where Q0Q_0 is the atomic unit with the same dimensions. Then:

  1. identify the physical dimensions of QQ;
  2. construct Q0Q_0 from a0a_0, EhE_{\mathrm h}, ee, mem_e, and ℏ\hbar;
  3. multiply by the numerical atomic-unit value;
  4. convert the resulting SI quantity to the desired laboratory unit;
  5. retain only precision supported by the input and chosen CODATA adjustment.

Example: energy in several spectroscopic units

Section titled “Example: energy in several spectroscopic units”

For ΔE=0.0500Eh\Delta E=0.0500E_{\mathrm h},

ΔE=1.36057 eV,ΔEhc=10 973.7 cm−1,ΔEh=328.984 THz.\begin{aligned} \Delta E&=1.36057\ \mathrm{eV},\\ \frac{\Delta E}{hc} &=10\,973.7\ \mathrm{cm}^{-1},\\ \frac{\Delta E}{h} &=328.984\ \mathrm{THz}. \end{aligned}

Each line describes the same interval. The last is ordinary frequency because the conversion used hh, not ℏ\hbar.

A field of 10 kV cm−110\ \mathrm{kV\,cm^{-1}} is 106 V m−110^6\ \mathrm{V\,m^{-1}}. Its atomic-unit value is

Fau=1065.1422×1011≃1.945×10−6.F_{\mathrm{au}} =\frac{10^6}{5.1422\times10^{11}} \simeq1.945\times10^{-6}.

This small dimensionless number can still produce strong mixing if the relevant opposite-parity level spacing is comparably small in Hartree.

A calculated dipole magnitude d=2.10d=2.10 a.u. corresponds to

d=2.10ea0,d≃5.34 D,d≃1.78×10−29 C m.\begin{aligned} d&=2.10ea_0,\\ d&\simeq5.34\ \mathrm D,\\ d&\simeq1.78\times10^{-29}\ \mathrm{C\,m}. \end{aligned}

The sign and phase of a matrix element depend on state conventions; transition probabilities use convention-invariant combinations such as squared magnitudes and angular sums.

Confusing atomic units with the atomic mass unit

Section titled “Confusing atomic units with the atomic mass unit”

The atomic unit of mass in the Hartree system is the electron mass mem_e. The unified atomic mass unit, symbol u\mathrm u, is approximately a nucleon-scale mass. They differ by a factor of about 1822.91822.9 and are not interchangeable.

Treating Hartree and Rydberg units as identical

Section titled “Treating Hartree and Rydberg units as identical”

1Eh=2Ry⁡1E_{\mathrm h}=2\operatorname{Ry}. A hydrogen ground-state energy of −1/2-1/2 and one of −1-1 can describe the same physics in different conventions. Inspect the Hamiltonian and code documentation.

E/hE/h is ordinary frequency in Hz, while E/ℏE/\hbar is angular frequency. They differ by 2π2\pi. The same distinction separates spectroscopic wavenumber 1/λ1/\lambda from angular wave number 2π/λ2\pi/\lambda.

The atomic-unit convention sets the positive magnitude e=1e=1. The electron charge remains qe=−1q_e=-1, so signs in scalar and vector potential couplings must be derived from qeq_e, not guessed from the unit convention.

Mixing atomic and SI fields inside one Hamiltonian

Section titled “Mixing atomic and SI fields inside one Hamiltonian”

A dipole in atomic units multiplied directly by a field in V/m does not produce an energy in Hartree. Convert the field by F0F_0 or convert the dipole to C m before multiplying.

Reporting polarizability without its convention

Section titled “Reporting polarizability without its convention”

An atomic polarizability quoted in a03a_0^3 usually denotes an atomic-unit value whose SI unit contains 4πϵ04\pi\epsilon_0. A bare volume and an SI polarizability do not have the same dimensions.

Assuming a dimensionless number has no units

Section titled “Assuming a dimensionless number has no units”

The printed number 0.050.05 is not enough. It could mean Hartree, Rydberg, eV, or a ratio to a problem-specific scale. Atomic-unit values are dimensionless numerical representations of dimensional physical quantities.

Keeping more digits than the calculation supports

Section titled “Keeping more digits than the calculation supports”

CODATA conversion factors may be precise to many digits, but a model energy reported as 0.0500Eh0.0500E_{\mathrm h} has only the precision justified by that input and the model. Conversion does not create information.

Exercise 1: Nondimensionalize a many-electron atom

Section titled “Exercise 1: Nondimensionalize a many-electron atom”

Starting from the fixed-nucleus Coulomb Hamiltonian for NN electrons, use ri=a0ρi\mathbf r_i=a_0\boldsymbol\rho_i and divide by EhE_{\mathrm h}. Write the resulting Hamiltonian in Hartree atomic units.

Solution

The SI Hamiltonian is

H=∑i[−ℏ22me∇ri2−Ze24πϵ0ri]+∑i<je24πϵ0rij.\begin{aligned} H={}&\sum_i\left[ -\frac{\hbar^2}{2m_e}\nabla_{r_i}^2 -\frac{Ze^2}{4\pi\epsilon_0r_i} \right]\\ &+\sum_{i<j} \frac{e^2}{4\pi\epsilon_0r_{ij}}. \end{aligned}

Because ∇ri=a0−1∇ρi\nabla_{r_i}=a_0^{-1}\nabla_{\rho_i}, each kinetic coefficient becomes Eh/2E_{\mathrm h}/2. Each Coulomb factor e2/(4πϵ0a0)e^2/(4\pi\epsilon_0a_0) becomes EhE_{\mathrm h}. Therefore

HEh=∑i[−12∇ρi2−Zρi]+∑i<j1ρij.\begin{aligned} \frac{H}{E_{\mathrm h}} ={}&\sum_i\left[ -\frac12\nabla_{\rho_i}^2 -\frac{Z}{\rho_i} \right]\\ &+\sum_{i<j}\frac{1}{\rho_{ij}}. \end{aligned}

After declaring atomic units, one normally renames ρi\boldsymbol\rho_i as ri\mathbf r_i and writes the dimensionless operator without the explicit factor H/EhH/E_{\mathrm h}.

For a vacuum wavelength λ=589 nm\lambda=589\ \mathrm{nm}, find the ordinary frequency, spectroscopic wavenumber, photon energy in eV, and photon energy in Hartree. Use enough digits to show the conversions, but do not imply that the wavelength is known more precisely than stated.

Solution

Use

ν=cλ,ν~=1λ,E=hcλ.\nu=\frac{c}{\lambda}, \qquad \widetilde\nu=\frac{1}{\lambda}, \qquad E=\frac{hc}{\lambda}.

For λ=589 nm\lambda=589\ \mathrm{nm},

ν≃5.09×1014 Hz,ν~≃1.70×104 cm−1,E≃2.10 eV,E≃0.0774Eh.\begin{aligned} \nu&\simeq5.09\times10^{14}\ \mathrm{Hz},\\ \widetilde\nu&\simeq1.70\times10^4\ \mathrm{cm}^{-1},\\ E&\simeq2.10\ \mathrm{eV},\\ E&\simeq0.0774E_{\mathrm h}. \end{aligned}

Using unrounded intermediate values gives approximately 508.99 THz508.99\ \mathrm{THz}, 16 977.9 cm−116\,977.9\ \mathrm{cm}^{-1}, 2.1050 eV2.1050\ \mathrm{eV}, and 0.077357Eh0.077357E_{\mathrm h}. Because the stated wavelength has three significant figures, the shorter rounded set is the defensible reported result.

Exercise 3: Field scale and interaction energy

Section titled “Exercise 3: Field scale and interaction energy”

A state has a dipole matrix element of one atomic unit. Estimate the interaction energy with a field of 1.00 MV m−11.00\ \mathrm{MV\,m^{-1}} in Hartree and eV. Is the field “weak” for every atomic transition?

Solution

The field in atomic units is

Fau=1.00×1065.1422×1011≃1.945×10−6.F_{\mathrm{au}} =\frac{1.00\times10^6}{5.1422\times10^{11}} \simeq1.945\times10^{-6}.

For d=1d=1 a.u., the characteristic interaction energy is therefore

∣dF∣≃1.945×10−6Eh,∣dF∣≃5.29×10−5 eV.\begin{aligned} |dF|&\simeq1.945\times10^{-6}E_{\mathrm h},\\ |dF|&\simeq5.29\times10^{-5}\ \mathrm{eV}. \end{aligned}

This is small compared with a Hartree-scale gross-structure interval, but it may be large compared with hyperfine splittings, narrow avoided crossings, or separations between Rydberg states. “Weak” must be defined relative to the manifold and observable.

Two programs solve the infinite-mass hydrogen problem. Program A reports E1s=−0.5E_{1s}=-0.5 and program B reports E1s=−1.0E_{1s}=-1.0, each labeled “a.u.” What should you inspect before concluding that one result is wrong?

Solution

Inspect the declared energy unit and the kinetic and Coulomb coefficients. In Hartree units,

H=−12∇2−1r,E1s=−12.H=-\frac12\nabla^2-\frac{1}{r}, \qquad E_{1s}=-\frac12.

In Rydberg units the Hamiltonian divided by Ry⁡=Eh/2\operatorname{Ry}=E_{\mathrm h}/2 is

HRy=−∇2−2r,E1s=−1.H_{\mathrm{Ry}}=-\nabla^2-\frac{2}{r}, \qquad E_{1s}=-1.

The outputs agree physically if program B uses Rydberg units. Also inspect whether either program reports orbital eigenvalues, total energies, binding magnitudes, or a shifted reference zero. A numerical label without its Hamiltonian convention is not enough for comparison.

Exercise 5: Estimate the fine-structure scale

Section titled “Exercise 5: Estimate the fine-structure scale”

Using α−1≃137.036\alpha^{-1}\simeq137.036, estimate α2Eh\alpha^2E_{\mathrm h} in eV, cm−1\mathrm{cm}^{-1}, and ordinary frequency. Explain why this is not a prediction for a particular fine-structure interval.

Solution

The dimensionless factor is

α2≃5.325×10−5.\alpha^2\simeq5.325\times10^{-5}.

Multiplying the Hartree conversions gives

α2Eh≃1.45×10−3 eV,α2Ehhc≃11.7 cm−1,α2Ehh≃350 GHz.\begin{aligned} \alpha^2E_{\mathrm h} &\simeq1.45\times10^{-3}\ \mathrm{eV},\\ \frac{\alpha^2E_{\mathrm h}}{hc} &\simeq11.7\ \mathrm{cm}^{-1},\\ \frac{\alpha^2E_{\mathrm h}}{h} &\simeq350\ \mathrm{GHz}. \end{aligned}

This is only an order-of-magnitude scale for a relative α2\alpha^2 correction to a Hartree-scale energy. A particular interval contains powers of ZZ, quantum-number-dependent coefficients, angular factors, cancellations, reduced-mass effects, and possibly electron correlation.

  1. P. J. Mohr, D. B. Newell, B. N. Taylor, and E. Tiesinga, “CODATA recommended values of the fundamental physical constants: 2022,” Reviews of Modern Physics 97, 025002 (2025), DOI: 10.1103/RevModPhys.97.025002.
  2. E. Tiesinga, P. J. Mohr, D. B. Newell, and B. N. Taylor, The 2022 CODATA Recommended Values of the Fundamental Physical Constants, Web Version 9.0, National Institute of Standards and Technology, 2024.
  3. Bureau International des Poids et Mesures, The International System of Units (SI), 9th ed. (2019), updated 2026, DOI: 10.59161/AUEZ1291.
  4. C. J. Foot, Atomic Physics, Oxford University Press, 2005.
  5. B. H. Bransden and C. J. Joachain, Physics of Atoms and Molecules, 2nd ed., Pearson, 2003.
  6. T. Helgaker, P. Jørgensen, and J. Olsen, Molecular Electronic-Structure Theory, Wiley, 2000.
  7. I. N. Levine, Quantum Chemistry, 7th ed., Pearson, 2014.