Skip to content

Atomic Units

Hartree atomic units are the default working units of much of atomic physics, electronic-structure theory, and quantum chemistry. They make electron-scale Coulomb Hamiltonians compact, but the printed number is interpretable only after the unit convention, quantity, and conversion release have been stated.

This page is the AMO quick-lookup entry. It gives convention checks, versioned conversion factors, and dimensional restoration rules. The Atomic Units and Scales page owns the derivation from the Coulomb problem and the hierarchy of physical scales. The site-wide Atomic Units translator gives the short cross-volume convention entry, while Constants and Conversions owns general conversions among eV, Hz, cm−1\mathrm{cm}^{-1}, K, magnetic units, and laboratory dipole units.

Use this page when a calculation, table, or code output says “a.u.” and you need to answer one of four questions:

  1. Does “a.u.” mean Hartree atomic units or a Rydberg-based convention?
  2. Which atomic unit belongs to the reported quantity?
  3. What multiplier converts the number to SI or a laboratory unit?
  4. Which metadata are needed so another researcher can reproduce the conversion?

The central rule is

X=XauX0,X=X_{\mathrm{au}}X_0,

where XX is the physical quantity, XauX_{\mathrm{au}} is its numerical value in Hartree atomic units, and X0X_0 is the atomic unit with the same physical dimensions. A bare statement such as “X=0.03X=0.03 a.u.” is incomplete unless XX is named: the unit for an energy is EhE_{\mathrm h}, the unit for a time is ℏ/Eh\hbar/E_{\mathrm h}, and the unit for an electric field is Eh/(ea0)E_{\mathrm h}/(ea_0).

Hartree atomic units set the numerical values

ℏ=me=e=4πε0=1.\hbar=m_e=e=4\pi\varepsilon_0=1.

This compact statement means:

  • action is measured in units of ℏ\hbar;
  • mass is measured in units of the electron mass mem_e;
  • charge is measured in units of the positive elementary-charge magnitude ee; and
  • the rationalized Coulomb constant 1/(4πε0)1/(4\pi\varepsilon_0) has numerical value one.

It does not mean that the physical constants cease to exist. They define the conversion between the dimensionless numerical representation and SI. Nor does it imply

ε0=1.\varepsilon_0=1.

Because the convention sets 4πε0=14\pi\varepsilon_0=1, the numerical value of ε0\varepsilon_0 itself is 1/(4π)1/(4\pi) when it appears separately.

The symbol ee is a positive magnitude. The electron charge remains

qe=−e=−1q_e=-e=-1

in Hartree atomic units. Nuclear charges are qA=ZAq_A=Z_A, and nuclear masses are large dimensionless ratios MA/meM_A/m_e, not unity.

The speed of light is also not set to one:

c=α−1≈137.036.c=\alpha^{-1}\approx137.036.

This makes the fine-structure constant α\alpha explicit as the parameter that separates nonrelativistic electron velocities from cc. An expression that simultaneously assumes c=1c=1 and Hartree atomic units is using an additional natural-unit convention and must say so.

Unit convention versus numerical conversion

Section titled “Unit convention versus numerical conversion”

There are two different kinds of equality on this page:

ℏ=me=e=4πε0=1\hbar=m_e=e=4\pi\varepsilon_0=1

is exact by unit convention, whereas

a0=5.291 772 105 44(82)×10−11 ma_0 =5.291\,772\,105\,44(82)\times10^{-11}\ \mathrm m

is a measured conversion factor from the 2022 CODATA adjustment. The parentheses give the standard uncertainty in the final quoted digits. Choosing atomic units exactly does not make their SI realizations exact.

The conversion factors below use the 2022 CODATA recommended values, NIST Web Version 9.0. NIST identifies this as the latest available CODATA adjustment as of 26 July 2026; the next regularly scheduled adjustment is the 2026 adjustment. SI status and defining constants follow version 4.01 of the ninth edition of the BIPM SI Brochure, published in June 2026.

This page uses three status labels:

StatusMeaning in this entryExample
convention exactfixed to one by the Hartree unit conventionℏ\hbar, mem_e, ee, 4πε04\pi\varepsilon_0 as atomic-unit numerical values
SI exact or derived exactfixed by the SI, perhaps through a nonterminating expressionee, hh, ℏ=h/(2π)\hbar=h/(2\pi)
CODATA measuredhas a recommended value and standard uncertaintya0a_0, EhE_{\mathrm h}, t0t_0

The first and third labels can apply to the same physical constant in different statements. For example, me=1m_e=1 is exact in the atomic-unit representation, while me=9.109 383 7139(28)×10−31 kgm_e=9.109\,383\,7139(28)\times10^{-31}\ \mathrm{kg} is a measured SI conversion.

Do not splice final digits from different adjustments. A reproducible result should name the release, such as “2022 CODATA,” or cite the source table.

The atomic unit of length is the Bohr radius,

a0=4πε0ℏ2mee2,=ℏαmec.\begin{aligned} a_0 &=\frac{4\pi\varepsilon_0\hbar^2}{m_e e^2},\\ &=\frac{\hbar}{\alpha m_e c}. \end{aligned}

Its 2022 CODATA value is

a0=5.291 772 105 44(82)×10−11 m,=0.529 177 210 544(82) A˚,=52.917 721 0544(82) pm.\begin{aligned} a_0 &=5.291\,772\,105\,44(82)\times10^{-11}\ \mathrm m,\\ &=0.529\,177\,210\,544(82)\ \text{\AA},\\ &=52.917\,721\,0544(82)\ \mathrm{pm}. \end{aligned}

The bohr is an electron-scale length, not a generic atomic radius. Many electron densities and bond lengths are of order a few a0a_0, but their values depend on state, charge, bonding, and the radius definition.

Useful reciprocal conversions are

1 A˚≈1.889 726 126 a0,1 nm≈18.897 261 26 a0.1\ \text{\AA} \approx1.889\,726\,126\,a_0, \qquad 1\ \mathrm{nm} \approx18.897\,261\,26\,a_0.

The atomic unit of energy is the Hartree energy,

Eh=ℏ2mea02,=e24πε0a0,=α2mec2,=2hcR∞.\begin{aligned} E_{\mathrm h} &=\frac{\hbar^2}{m_ea_0^2},\\ &=\frac{e^2}{4\pi\varepsilon_0a_0},\\ &=\alpha^2m_ec^2,\\ &=2hcR_\infty. \end{aligned}

The identities expose several interpretations: kinetic energy on the scale a0a_0, Coulomb energy at separation a0a_0, a small fraction α2\alpha^2 of the electron rest energy, and twice the infinite-mass Rydberg energy.

The 2022 CODATA conversion is

Eh=4.359 744 722 2060(48)×10−18 J,=27.211 386 245 981(30) eV.\begin{aligned} E_{\mathrm h} &=4.359\,744\,722\,2060(48)\times10^{-18}\ \mathrm J,\\ &=27.211\,386\,245\,981(30)\ \mathrm{eV}. \end{aligned}

Define the infinite-mass Rydberg energy by

ERy=hcR∞=Eh2.E_{\mathrm{Ry}} =hcR_\infty =\frac{E_{\mathrm h}}{2}.

Therefore

1Eh=2ERy,ERy=13.605 693 122 990(15) eV.1E_{\mathrm h} =2E_{\mathrm{Ry}}, \qquad E_{\mathrm{Ry}} =13.605\,693\,122\,990(15)\ \mathrm{eV}.

Three objects must not be collapsed into the same symbol:

ObjectQuantity typeRelation
Hartree energy EhE_{\mathrm h}energy2hcR∞2hcR_\infty
Rydberg energy ERyE_{\mathrm{Ry}}energyhcR∞=Eh/2hcR_\infty=E_{\mathrm h}/2
Rydberg constant R∞R_\inftyinverse lengthERy/(hc)E_{\mathrm{Ry}}/(hc)

ERyE_{\mathrm{Ry}} is an infinite-nuclear-mass scale. It is not exactly the measured ionization energy of ordinary hydrogen, which includes reduced-mass, recoil, relativistic, radiative, finite-size, and hyperfine effects.

The atomic unit of time is

t0=ℏEh=2.418 884 326 5864(26)×10−17 s.t_0=\frac{\hbar}{E_{\mathrm h}} =2.418\,884\,326\,5864(26)\times10^{-17}\ \mathrm s.

Equivalently,

t0=24.188 843 265 864(26) as=0.024 188 843 265 864(26) fs.t_0 =24.188\,843\,265\,864(26)\ \mathrm{as} =0.024\,188\,843\,265\,864(26)\ \mathrm{fs}.

It is the inverse angular-frequency scale:

1t0=Ehℏ.\frac{1}{t_0} =\frac{E_{\mathrm h}}{\hbar}.

The ordinary frequency corresponding to one Hartree is instead Eh/hE_{\mathrm h}/h and is smaller by 2π2\pi.

The tables give multipliers for converting a numerical Hartree-atomic-unit value to SI:

XSI=Xau×(SI value of X0).X_{\mathrm{SI}} =X_{\mathrm{au}}\times \left(\text{SI value of }X_0\right).

To convert SI to atomic units, divide by the same multiplier. Parenthetical digits are standard uncertainties; displayed trailing digits should not be treated as the accuracy of a physical model.

QuantityAtomic unit X0X_02022 CODATA SI value
chargeee1.602 176 634×10−19 C1.602\,176\,634\times10^{-19}\ \mathrm C exact
massmem_e9.109 383 7139(28)×10−31 kg9.109\,383\,7139(28)\times10^{-31}\ \mathrm{kg}
actionℏ\hbar1.054 571 817…×10−34 J s1.054\,571\,817\ldots\times10^{-34}\ \mathrm{J\,s} derived exact
lengtha0a_05.291 772 105 44(82)×10−11 m5.291\,772\,105\,44(82)\times10^{-11}\ \mathrm m
energyEhE_{\mathrm h}4.359 744 722 2060(48)×10−18 J4.359\,744\,722\,2060(48)\times10^{-18}\ \mathrm J
timet0=ℏ/Eht_0=\hbar/E_{\mathrm h}2.418 884 326 5864(26)×10−17 s2.418\,884\,326\,5864(26)\times10^{-17}\ \mathrm s
velocitya0/t0=αca_0/t_0=\alpha c2.187 691 262 16(34)×106 m s−12.187\,691\,262\,16(34)\times10^6\ \mathrm{m\,s^{-1}}
momentumℏ/a0\hbar/a_01.992 851 915 45(31)×10−24 kg m s−11.992\,851\,915\,45(31)\times10^{-24}\ \mathrm{kg\,m\,s^{-1}}
forceEh/a0E_{\mathrm h}/a_08.238 723 5038(13)×10−8 N8.238\,723\,5038(13)\times10^{-8}\ \mathrm N
QuantityAtomic unit X0X_02022 CODATA SI value
electric potentialEh/eE_{\mathrm h}/e27.211 386 245 981(30) V27.211\,386\,245\,981(30)\ \mathrm V
electric fieldEh/(ea0)E_{\mathrm h}/(ea_0)5.142 206 751 12(80)×1011 V m−15.142\,206\,751\,12(80)\times10^{11}\ \mathrm{V\,m^{-1}}
electric-field gradientEh/(ea02)E_{\mathrm h}/(ea_0^2)9.717 362 4424(30)×1021 V m−29.717\,362\,4424(30)\times10^{21}\ \mathrm{V\,m^{-2}}
electric dipole momentea0ea_08.478 353 6198(13)×10−30 C m8.478\,353\,6198(13)\times10^{-30}\ \mathrm{C\,m}
electric quadrupole momentea02ea_0^24.486 551 5185(14)×10−40 C m24.486\,551\,5185(14)\times10^{-40}\ \mathrm{C\,m^2}
electric polarizabilitye2a02/Ehe^2a_0^2/E_{\mathrm h}1.648 777 272 12(51)×10−41 C2 m2 J−11.648\,777\,272\,12(51)\times10^{-41}\ \mathrm{C^2\,m^2\,J^{-1}}
magnetic flux densityℏ/(ea02)\hbar/(ea_0^2)2.350 517 570 77(73)×105 T2.350\,517\,570\,77(73)\times10^5\ \mathrm T
magnetic dipole momentℏe/me=2μB\hbar e/m_e=2\mu_{\mathrm B}1.854 802 013 15(58)×10−23 J T−11.854\,802\,013\,15(58)\times10^{-23}\ \mathrm{J\,T^{-1}}

The atomic unit of magnetic dipole moment is 2μB2\mu_{\mathrm B}, not μB\mu_{\mathrm B}. Many atomic-structure tables deliberately use the Bohr magneton instead. A value labeled only “a.u.” is therefore not enough to identify a magnetic-moment convention.

Representation of EhE_{\mathrm h}2022 CODATA valueConversion used
joules4.359 744 722 2060(48)×10−18 J4.359\,744\,722\,2060(48)\times10^{-18}\ \mathrm Jdirect
electronvolts27.211 386 245 981(30) eV27.211\,386\,245\,981(30)\ \mathrm{eV}Eh/eE_{\mathrm h}/e
ordinary frequency6.579 683 920 5000(72)×1015 Hz6.579\,683\,920\,5000(72)\times10^{15}\ \mathrm{Hz}Eh/hE_{\mathrm h}/h
angular frequency4.134 137 333 517(5)×1016 rad s−14.134\,137\,333\,517(5)\times10^{16}\ \mathrm{rad\,s^{-1}}Eh/ℏE_{\mathrm h}/\hbar
vacuum wavenumber219 474.631 363 14(24) cm−1219\,474.631\,363\,14(24)\ \mathrm{cm^{-1}}Eh/(hc)E_{\mathrm h}/(hc)
energy-equivalent temperature315 775.024 804(35) K315\,775.024\,804(35)\ \mathrm KEh/kBE_{\mathrm h}/k_{\mathrm B}

The last displayed digits in values derived here from rounded table entries are for lookup, not a replacement for machine-readable CODATA data. For precision metrology, retrieve the value and covariance information directly from NIST.

Starting quantityEquivalent
1 A˚1\ \text{\AA}1.889 726 126 a01.889\,726\,126\,a_0
1a01a_00.529 177 210 544 A˚0.529\,177\,210\,544\ \text{\AA}
1Eh1E_{\mathrm h}27.211 386 246 eV27.211\,386\,246\ \mathrm{eV}
1Eh1E_{\mathrm h}2ERy2E_{\mathrm{Ry}}
1Eh1E_{\mathrm h}219 474.6314 cm−1219\,474.6314\ \mathrm{cm^{-1}}
1t01t_024.188 8433 as24.188\,8433\ \mathrm{as}
1ea01ea_02.541 746 47 D2.541\,746\,47\ \mathrm D
1 D1\ \mathrm D0.393 430 27 ea00.393\,430\,27\,ea_0
1 kV cm−11\ \mathrm{kV\,cm^{-1}}1.944 690 38×10−71.944\,690\,38\times10^{-7} a.u. of electric field
1a021a_0^22.800 285 20×10−21 m22.800\,285\,20\times10^{-21}\ \mathrm{m^2}

Here the debye conversion uses

1 D=3.335 640 951…×10−30 C m.1\ \mathrm D =3.335\,640\,951\ldots\times10^{-30}\ \mathrm{C\,m}.

The coefficients in a Hamiltonian are often the fastest way to identify the unit convention.

For an infinitely heavy nucleus of charge ZZ, the Hartree-unit Hamiltonian is

H=−12∇2−Zr.H =-\frac12\nabla^2-\frac{Z}{r}.

The hydrogenic bound energies are

En=−Z22n2.E_n=-\frac{Z^2}{2n^2}.

Thus hydrogen has E1s=−1/2E_{1s}=-1/2 in Hartree. The full derivation and wavefunction normalization live in Hydrogen Atom and Atomic Units and Scales.

For nuclei at fixed positions RA\mathbf R_A, the electronic Hamiltonian in Hartree units is

Hel=−12∑i∇i2−∑i,AZAriA+∑i<j1rij+∑A<BZAZBRAB.\begin{aligned} H_{\mathrm{el}} ={}&-\frac12\sum_i\nabla_i^2 -\sum_{i,A}\frac{Z_A}{r_{iA}}\\ &+\sum_{i<j}\frac{1}{r_{ij}} +\sum_{A<B}\frac{Z_AZ_B}{R_{AB}}. \end{aligned}

Some electronic-structure programs omit the final nuclear-repulsion term from the electronic eigenvalue and add it when forming the total Born–Oppenheimer energy. The unit convention does not determine this bookkeeping choice; the program documentation must.

In the electric-dipole approximation,

Hint=−d⋅E.H_{\mathrm{int}} =-\mathbf d\cdot\boldsymbol{\mathcal E}.

With d\mathbf d in ea0ea_0 and E\boldsymbol{\mathcal E} in Eh/(ea0)E_{\mathrm h}/(ea_0), their product is automatically in Hartree:

(ea0)(Ehea0)=Eh.(ea_0) \left(\frac{E_{\mathrm h}}{ea_0}\right) =E_{\mathrm h}.

The field may be a static value, a peak optical amplitude, a root-mean-square amplitude, or a complex positive-frequency amplitude. Atomic units do not remove that convention choice.

Time-dependent Schrödinger evolution illustrates why inverse atomic time is an angular-frequency unit. In SI,

∣ψ(t)⟩=e−iHt/ℏ∣ψ(0)⟩.\lvert\psi(t)\rangle =e^{-iHt/\hbar}\lvert\psi(0)\rangle.

Writing

H=EhHau,t=t0tau,H=E_{\mathrm h}H_{\mathrm{au}}, \qquad t=t_0t_{\mathrm{au}},

and using Eht0/ℏ=1E_{\mathrm h}t_0/\hbar=1 gives

∣ψ(tau)⟩=e−iHautau∣ψ(0)⟩.\lvert\psi(t_{\mathrm{au}})\rangle =e^{-iH_{\mathrm{au}}t_{\mathrm{au}}} \lvert\psi(0)\rangle.

For an energy gap ΔE\Delta E,

ω=ΔEℏ,ν=ΔEh=ω2π.\begin{aligned} \omega &=\frac{\Delta E}{\hbar},\\ \nu &=\frac{\Delta E}{h} =\frac{\omega}{2\pi}. \end{aligned}

Consequently,

1 a.u. of inverse time=Ehℏ≈4.1341×1016 rad s−1,1Ehh≈6.5797×1015 Hz.\begin{aligned} 1\ \text{a.u. of inverse time} &=\frac{E_{\mathrm h}}{\hbar} \approx4.1341\times10^{16}\ \mathrm{rad\,s^{-1}},\\ \frac{1E_{\mathrm h}}{h} &\approx6.5797\times10^{15}\ \mathrm{Hz}. \end{aligned}

Calling Eh/ℏE_{\mathrm h}/\hbar “hertz” creates a factor-of-2π2\pi error. For a reported Rabi frequency, linewidth, or decay parameter, retain the defining equation and state whether the displayed number is divided by 2π2\pi.

Electric Fields, Dipoles, and Polarizabilities

Section titled “Electric Fields, Dipoles, and Polarizabilities”

The atomic field

E0=Ehea0≈5.1422×1011 V m−1\mathcal E_0 =\frac{E_{\mathrm h}}{ea_0} \approx5.1422\times10^{11}\ \mathrm{V\,m^{-1}}

changes the potential energy of one elementary charge by one Hartree over one Bohr radius. Laboratory fields are therefore often small as atomic-unit numbers:

10 kV cm−1≈1.9447×10−6 a.u.10\ \mathrm{kV\,cm^{-1}} \approx1.9447\times10^{-6}\ \text{a.u.}

Small in atomic units does not automatically mean perturbatively weak. The relevant comparison is between the coupling matrix element and the energy separation of the states being mixed.

The atomic dipole unit is

d0=ea0.d_0=ea_0.

A calculated transition dipole dau=2.0d_{\mathrm{au}}=2.0 therefore means

d=2.0ea0,≈5.0835 D,≈1.6957×10−29 C m.\begin{aligned} d &=2.0ea_0,\\ &\approx5.0835\ \mathrm D,\\ &\approx1.6957\times10^{-29}\ \mathrm{C\,m}. \end{aligned}

The sign or phase of an off-diagonal dipole matrix element depends on state phase conventions. Observable line strengths use phase-invariant combinations, together with angular and degeneracy factors.

For a static scalar polarizability defined by

ΔE=−12αEE2,\Delta E =-\frac12\alpha_E\mathcal E^2,

the Hartree atomic unit is

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

The final identity explains a persistent notation trap. A tabulated polarizability “in a03a_0^3” usually reports the same numerical value as the Hartree atomic-unit polarizability, but the SI polarizability is not obtained by multiplying by a bare volume. It is obtained from

αE(SI)=αE,au(4πε0a03).\alpha_E^{(\mathrm{SI})} =\alpha_{E,\mathrm{au}} \left(4\pi\varepsilon_0a_0^3\right).

If a source instead tabulates the polarizability volume αE/(4πε0)\alpha_E/(4\pi\varepsilon_0), the conversion multiplier is a03a_0^3. Record which quantity the table means.

For a plane wave in vacuum whose electric-field peak amplitude is Epk\mathcal E_{\mathrm{pk}},

I=12cε0Epk2.I =\frac12c\varepsilon_0 \mathcal E_{\mathrm{pk}}^2.

If the quoted field is root-mean-square,

I=cε0Erms2.I =c\varepsilon_0 \mathcal E_{\mathrm{rms}}^2.

An “atomic unit of intensity” can therefore differ by a factor of two when the amplitude convention is omitted. Convert the explicitly defined field amplitude first, then apply the corresponding intensity formula.

Memorized tables are useful, but dimensional restoration handles quantities that are not listed.

Suppose a quantity XX has dimensions

[X]=MaLbTcQd,[X] =M^aL^bT^cQ^d,

where QQ denotes electric-charge dimension. Its Hartree atomic unit is

X0=meaa0bt0ced.X_0 =m_e^a a_0^b t_0^c e^d.

Then

X=XauX0.X=X_{\mathrm{au}}X_0.

Examples follow directly:

QuantityDimensionsAtomic unit
area or cross sectionL2L^2a02a_0^2
volumeL3L^3a03a_0^3
number densityL−3L^{-3}a0−3a_0^{-3}
rateT−1T^{-1}t0−1t_0^{-1}
currentQT−1QT^{-1}e/t0e/t_0
momentumMLT−1MLT^{-1}mea0/t0=ℏ/a0m_ea_0/t_0=\hbar/a_0
forceMLT−2MLT^{-2}mea0/t02=Eh/a0m_ea_0/t_0^2=E_{\mathrm h}/a_0
dipole momentQLQLea0ea_0

A normalized three-dimensional wavefunction satisfies

∫∣ψ(r)∣2 d3r=1.\int \lvert\psi(\mathbf r)\rvert^2\,d^3r =1.

If r=a0rau\mathbf r=a_0\mathbf r_{\mathrm{au}}, then

ψ(r)=a0−3/2ψau(rau).\psi(\mathbf r) =a_0^{-3/2} \psi_{\mathrm{au}}(\mathbf r_{\mathrm{au}}).

Therefore a probability density reported in atomic units carries the scale a0−3a_0^{-3}. A radial function may use a different normalization and power of a0a_0; inspect its defining integral before restoring units.

A decay rate Γau\Gamma_{\mathrm{au}} converts as

ΓSI=Γaut0.\Gamma_{\mathrm{SI}} =\frac{\Gamma_{\mathrm{au}}}{t_0}.

A cross section converts as

σSI=σaua02.\sigma_{\mathrm{SI}} =\sigma_{\mathrm{au}}a_0^2.

Because

a02≈2.8003×10−17 cm2,a_0^2 \approx2.8003\times10^{-17}\ \mathrm{cm^2},

a cross section of 10−410^{-4} a.u. is approximately 2.80×10−21 cm22.80\times10^{-21}\ \mathrm{cm^2}. This conversion says nothing about whether the underlying scattering approximation is accurate.

The abbreviation “a.u.” is sometimes used for two related but numerically different systems.

In Hartree units,

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

Dividing the same physical Hamiltonian by ERy=Eh/2E_{\mathrm{Ry}}=E_{\mathrm h}/2 gives Rydberg units:

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

For a many-electron Coulomb Hamiltonian, every energy coefficient doubles under the same change of energy unit. A number expressed in Rydberg converts to Hartree through

E[Ha]=12E[Ry].E[\mathrm{Ha}] =\frac12E[\mathrm{Ry}].

Do not infer the convention from a file extension or a field name alone. Check:

  1. the kinetic-energy coefficient;
  2. the documented energy unit;
  3. a known benchmark such as hydrogen 1s1s;
  4. whether a cutoff, eigenvalue, total energy, and force use the same unit;
  5. whether the output labels Hartree as Ha, Eh, or hartree.

For example, official Quantum ESPRESSO documentation states that much of PWscf uses Rydberg units while CP uses Hartree atomic units. This is a program-specific contract, not a universal rule for plane-wave codes.

The atomic unit of mass is the electron mass:

m0=me.m_0=m_e.

It is not the unified atomic mass unit u\mathrm u. A nuclear mass supplied in unified atomic mass units must be converted to electron-mass units before entering a Hartree-unit Hamiltonian:

Mau=Mme.M_{\mathrm{au}} =\frac{M}{m_e}.

For a two-body Coulomb problem, the reduced mass is

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

or numerically in electron-mass units,

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

One may introduce problem-specific reduced-mass-scaled units, but those are not the CODATA bohr and Hartree unless explicitly declared. In standard Hartree units, isotope and recoil effects belong in the Hamiltonian rather than in a silent redefinition of a0a_0.

The same care applies to molecular coordinates. Electronic coordinates are often in bohr, nuclear masses in electron masses, vibrational frequencies in Eh/ℏE_{\mathrm h}/\hbar, and final spectra in cm−1\mathrm{cm}^{-1}. Every handoff must preserve the quantity and unit, not only the number.

For a bond length R=1.40 A˚R=1.40\ \text{\AA},

Rau=1.40 A˚0.529 177 210 544 A˚,≈2.646.\begin{aligned} R_{\mathrm{au}} &=\frac{1.40\ \text{\AA}} {0.529\,177\,210\,544\ \text{\AA}},\\ &\approx2.646. \end{aligned}

The result should not be printed to twelve digits when the input was given to three significant figures.

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

ΔE≈3.4014 eV,ΔEhc≈27 434.3 cm−1,ΔEh≈822.460 THz,ΔEkB≈39 471.9 K.\begin{aligned} \Delta E &\approx3.4014\ \mathrm{eV},\\ \frac{\Delta E}{hc} &\approx27\,434.3\ \mathrm{cm^{-1}},\\ \frac{\Delta E}{h} &\approx822.460\ \mathrm{THz},\\ \frac{\Delta E}{k_{\mathrm B}} &\approx39\,471.9\ \mathrm K. \end{aligned}

The kelvin value is an energy-equivalent temperature, not a claim that the system is in thermal equilibrium at that temperature.

A real-time calculation propagated for 2500t02500t_0 covers

t=2500(2.418 884 3266×10−17 s),≈60.4721 fs.\begin{aligned} t &=2500 \left(2.418\,884\,3266\times10^{-17}\ \mathrm s\right),\\ &\approx60.4721\ \mathrm{fs}. \end{aligned}

The timestep must also be converted. Reporting only the number of steps does not specify the simulated time interval.

Take a dipole matrix element d=2.3ea0d=2.3ea_0 and a peak field E=50 kV cm−1=5.0×106 V m−1\mathcal E=50\ \mathrm{kV\,cm^{-1}}=5.0\times10^6\ \mathrm{V\,m^{-1}}. The field in atomic units is

Eau=5.0×1065.1422×1011≈9.72×10−6.\mathcal E_{\mathrm{au}} =\frac{5.0\times10^6} {5.1422\times10^{11}} \approx9.72\times10^{-6}.

The coupling scale is

dE≈2.24×10−5Eh,≈6.09×10−4 eV,dEh≈147 GHz.\begin{aligned} d\mathcal E &\approx2.24\times10^{-5}E_{\mathrm h},\\ &\approx6.09\times10^{-4}\ \mathrm{eV},\\ \frac{d\mathcal E}{h} &\approx147\ \mathrm{GHz}. \end{aligned}

This is the magnitude before angular factors, detuning, polarization projection, and any peak-versus-complex-amplitude factor are applied.

At minimum, a computational or experimental handoff should record:

FieldExample
quantitytotal electronic energy
numerical value−75.983 948-75.983\,948
unitHartree, not merely a.u.
conventionHartree atomic units
physical modelclamped nuclei, nonrelativistic electronic Hamiltonian
reference zeroseparated nuclei and electrons, or code-specific zero
conversion release2022 CODATA
output provenanceprogram, version, input, method, basis, and convergence settings

For a field or time-dependent interaction, also record:

  • peak, root-mean-square, or complex field amplitude;
  • cyclic or angular frequency;
  • gauge and interaction convention when relevant;
  • pulse-envelope definition and full-width convention; and
  • whether quoted rates are population or amplitude decay rates.

The converted number should not imply more knowledge than the original result. If

E=0.125EhE=0.125E_{\mathrm h}

is specified to three significant figures, conversion with a twelve-significant-figure CODATA multiplier does not make the physical energy twelve-digit accurate.

Treating “a.u.” as a complete unit label

Section titled “Treating “a.u.” as a complete unit label”

Atomic-unit dimensions depend on the quantity. Write Hartree for energy, bohr for length, ea0ea_0 for dipole moment, or explicitly define the unit.

1Eh=2ERy1E_{\mathrm h}=2E_{\mathrm{Ry}}. Inspect the kinetic coefficient and a known eigenvalue before comparing outputs.

Confusing atomic mass and unified atomic mass

Section titled “Confusing atomic mass and unified atomic mass”

The Hartree atomic mass unit is mem_e. The unified atomic mass unit u\mathrm u is roughly a nucleon mass and is about 1823me1823m_e.

e=1e=1 denotes a positive magnitude. The electron charge is −1-1.

t0−1=Eh/ℏt_0^{-1}=E_{\mathrm h}/\hbar is an angular-frequency scale. Divide by 2π2\pi to obtain cycles per second.

A dipole in ea0ea_0 times a field in V/m is not yet an energy. Convert both quantities to one coherent system before multiplying.

Converting polarizability as a bare volume

Section titled “Converting polarizability as a bare volume”

The SI multiplier for Hartree atomic polarizability is 4πε0a034\pi\varepsilon_0a_0^3, not just a03a_0^3.

Peak and root-mean-square fields differ by 2\sqrt{2} for a sinusoid, so intensities inferred from them differ by a factor of two.

Reduced-mass-scaled coordinates can be useful, but they are a separate declared convention. Standard a0a_0 uses mem_e.

Values are correlated products of one adjustment. Name one release and keep it throughout a precision calculation.

Reporting converted digits as model accuracy

Section titled “Reporting converted digits as model accuracy”

Unit conversion preserves uncertainty and model error; it does not improve either.

Exercise 1: Identify the energy convention

Section titled “Exercise 1: Identify the energy convention”

Two programs give the hydrogen ground-state energy as −0.5-0.5 and −1.0-1.0, respectively. Program A uses

HA=−12∇2−1r,H_A=-\frac12\nabla^2-\frac1r,

while program B uses

HB=−∇2−2r.H_B=-\nabla^2-\frac2r.

Are the physical predictions inconsistent?

Solution

No. Program A divides physical energy by EhE_{\mathrm h}, while program B divides by ERy=Eh/2E_{\mathrm{Ry}}=E_{\mathrm h}/2. Therefore

−0.5Eh=−1.0ERy.-0.5E_{\mathrm h} =-1.0E_{\mathrm{Ry}}.

The kinetic and Coulomb coefficients change together under the factor-of-two rescaling. A valid comparison first converts both outputs to the same unit and checks that their reference zeros and physical approximations also agree.

Exercise 2: Convert an electronic energy gap

Section titled “Exercise 2: Convert an electronic energy gap”

Convert ΔE=0.0200Eh\Delta E=0.0200E_{\mathrm h} to eV, ordinary frequency in THz, and vacuum wavenumber in cm−1\mathrm{cm}^{-1}. State an appropriate number of significant figures.

Solution

Use the one-Hartree table:

ΔE=0.0200(27.211 386 246 eV),=0.544 eV,ΔEh=0.0200(6579.6839 THz),=132 THz,ΔEhc=0.0200(219 474.6314 cm−1),=4.39×103 cm−1.\begin{aligned} \Delta E &=0.0200 \left(27.211\,386\,246\ \mathrm{eV}\right),\\ &=0.544\ \mathrm{eV},\\ \frac{\Delta E}{h} &=0.0200 \left(6579.6839\ \mathrm{THz}\right),\\ &=132\ \mathrm{THz},\\ \frac{\Delta E}{hc} &=0.0200 \left(219\,474.6314\ \mathrm{cm^{-1}}\right),\\ &=4.39\times10^3\ \mathrm{cm^{-1}}. \end{aligned}

The input has three significant figures, so defensible reported values are 0.544 eV0.544\ \mathrm{eV}, 132 THz132\ \mathrm{THz}, and 4.39×103 cm−14.39\times10^3\ \mathrm{cm^{-1}}. Extra CODATA digits belong in the intermediate calculation, not the final claim.

A two-state energy gap is 0.010Eh0.010E_{\mathrm h}. Find its angular frequency, ordinary frequency, and period. Explain the role of 2π2\pi.

Solution

The angular frequency is

ω=0.010Ehℏ,≈4.13×1014 rad s−1.\begin{aligned} \omega &=0.010\frac{E_{\mathrm h}}{\hbar},\\ &\approx4.13\times10^{14}\ \mathrm{rad\,s^{-1}}. \end{aligned}

The ordinary frequency is

ν=ω2π,=0.010Ehh,≈6.58×1013 Hz.\begin{aligned} \nu &=\frac{\omega}{2\pi},\\ &=0.010\frac{E_{\mathrm h}}{h},\\ &\approx6.58\times10^{13}\ \mathrm{Hz}. \end{aligned}

Hence

T=1ν≈1.52×10−14 s=15.2 fs.T=\frac1\nu\approx1.52\times10^{-14}\ \mathrm s=15.2\ \mathrm{fs}.

Eh/ℏE_{\mathrm h}/\hbar counts phase in radians per second, while Eh/hE_{\mathrm h}/h counts cycles per second. Omitting 2π2\pi would confuse these two quantities.

A transition dipole is 3.00 D3.00\ \mathrm D and a field is 1.00 kV cm−11.00\ \mathrm{kV\,cm^{-1}}. Convert each to atomic units and estimate dEd\mathcal E in Hartree.

Solution

The dipole is

dau=3.00(0.393 430 27)≈1.18.d_{\mathrm{au}} =3.00 \left(0.393\,430\,27\right) \approx1.18.

The field is

Eau=1.944 69×10−7.\mathcal E_{\mathrm{au}} =1.944\,69\times10^{-7}.

Their product is

dE≈(1.18)(1.94469×10−7)Eh,≈2.30×10−7Eh.\begin{aligned} d\mathcal E &\approx (1.18)(1.94469\times10^{-7})E_{\mathrm h},\\ &\approx2.30\times10^{-7}E_{\mathrm h}. \end{aligned}

This is a coupling scale. A transition Hamiltonian may include polarization, angular-momentum, rotating-wave, or field-amplitude factors, so one should not identify it with an observed splitting without specifying the model.

A scattering calculation reports

σ=6.0×10−3\sigma=6.0\times10^{-3}

in Hartree atomic units. Convert the cross section to cm2\mathrm{cm^2}.

Solution

A cross section has dimensions of area, so its atomic unit is a02a_0^2:

σ=6.0×10−3a02.\sigma =6.0\times10^{-3}a_0^2.

Using

a02≈2.8003×10−17 cm2,a_0^2 \approx2.8003\times10^{-17}\ \mathrm{cm^2},

gives

σ≈1.7×10−19 cm2.\sigma \approx1.7\times10^{-19}\ \mathrm{cm^2}.

The two-significant-figure result reflects the precision of the input. No special cross-section conversion rule was required; dimensional restoration was enough.

A database lists a static polarizability as 120120 a.u. Write the conversion to SI polarizability and to the polarizability volume in a03a_0^3. Why are these not the same dimensional statement?

Solution

For the Stark-shift definition

ΔE=−12αEE2,\Delta E=-\frac12\alpha_E\mathcal E^2,

the SI polarizability is

αE(SI)=120(1.648 777 272 12×10−41)×C2 m2 J−1,≈1.98×10−39 C2 m2 J−1.\begin{aligned} \alpha_E^{(\mathrm{SI})} &=120 \left( 1.648\,777\,272\,12\times10^{-41} \right)\\ &\quad\times \mathrm{C^2\,m^2\,J^{-1}},\\ &\approx1.98\times10^{-39} \ \mathrm{C^2\,m^2\,J^{-1}}. \end{aligned}

The corresponding polarizability volume is

αE4πε0=120a03.\frac{\alpha_E}{4\pi\varepsilon_0} =120a_0^3.

The numerical value 120120 is shared because

αE,0=4πε0a03,\alpha_{E,0}=4\pi\varepsilon_0a_0^3,

but SI polarizability and volume have different dimensions. A source that prints only “120a03120a_0^3” should state that it is reporting αE/(4πε0)\alpha_E/(4\pi\varepsilon_0).

A paper replaces mem_e by the electron–nucleus reduced mass μ\mu when defining its length scale and then calls the result “one bohr.” Is that standard Hartree atomic-unit usage?

Solution

No. The standard CODATA bohr is

a0=4πε0ℏ2mee2a_0 =\frac{4\pi\varepsilon_0\hbar^2}{m_ee^2}

and uses mem_e. Replacing mem_e by μ\mu defines a useful problem-specific scaled length,

aμ=4πε0ℏ2μe2,a_\mu =\frac{4\pi\varepsilon_0\hbar^2}{\mu e^2},

but aμa_\mu is isotope-dependent and is not numerically identical to a0a_0. The paper can use that convention if it declares it; a data handoff must preserve the distinction. In standard Hartree units, one instead keeps a0a_0 and inserts the reduced-mass factor in the Hamiltonian.

Exercise 8: Audit a reproducibility record

Section titled “Exercise 8: Audit a reproducibility record”

A methods section states: “The pulse amplitude was 0.0010.001 a.u., the frequency was 0.050.05 a.u., and the calculation ran for 5000 steps.” List the missing information needed to reproduce the physical field and duration.

Solution

At least the following are missing:

  1. whether the amplitude is an electric field, vector potential, or another quantity;
  2. whether the field amplitude is peak, root-mean-square, or a complex positive-frequency amplitude;
  3. whether “frequency” means angular frequency in t0−1t_0^{-1} or cyclic frequency;
  4. the timestep in t0t_0 or seconds;
  5. the pulse-envelope definition, carrier phase, and duration convention;
  6. whether the code uses Hartree or Rydberg units;
  7. the code and version implementing the convention; and
  8. the CODATA release used for any SI conversion.

If the first number is a peak electric field in Hartree units, it would mean 0.001Eh/(ea0)0.001E_{\mathrm h}/(ea_0). If the second is an angular frequency, it means 0.05Eh/ℏ0.05E_{\mathrm h}/\hbar. The total physical duration cannot be recovered from 5000 steps without the timestep.

  1. P. J. Mohr, D. B. Newell, B. N. Taylor, and E. Tiesinga, “CODATA recommended values of the fundamental physical constants: 2022,” Journal of Physical and Chemical Reference Data 54, 033105 (2025), DOI: 10.1063/5.0279860. Table XXXV is the primary source for the atomic-unit conversions used here.
  2. 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.
  3. 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).
  4. Bureau International des Poids et Mesures, The International System of Units (SI), 9th ed., version 4.01 (2026), DOI: 10.59161/AUEZ1291.
  5. C. M. A. Brett, J. G. Frey, R. Hinde, Y. Kuroda, R. Marquardt, F. Pavese, M. Quack, J. Stohner, and A. J. Thor, eds., Quantities, Units and Symbols in Physical Chemistry, abridged 4th ed., IUPAC and Royal Society of Chemistry, 2023.
  6. E. Tiesinga, “Units and constants”, in Springer Handbook of Atomic, Molecular, and Optical Physics (2023).
  7. T. Helgaker, P. Jørgensen, and J. Olsen, Molecular Electronic-Structure Theory, Wiley, 2000.
  8. P. W. Atkins and R. S. Friedman, Molecular Quantum Mechanics, 5th ed., Oxford University Press, 2010.
  9. Quantum ESPRESSO Foundation, “What are the units for quantity XYZ?”, official user documentation, accessed 26 July 2026.