# Hydrogen-like atom

Last updated

A hydrogen-like atom (or hydrogenic atom) is any atom or ion with a single valence electron. These atoms are isoelectronic with hydrogen. Examples of hydrogen-like atoms include, but are not limited to, hydrogen itself, all alkali metals such as Rb and Cs, singly ionized alkaline earth metals such as Ca+ and Sr+ and other ions such as He+, Li2+, and Be3+ and isotopes of any of the above. A hydrogen-like atom includes a positively charged core consisting of the atomic nucleus and any core electrons as well as a single valence electron. Because helium is common in the universe, the spectroscopy of singly ionized helium is important in EUV astronomy, for example, of DO white dwarf stars.

## Contents

The non-relativistic Schrödinger equation and relativistic Dirac equation for the hydrogen atom can be solved analytically, owing to the simplicity of the two-particle physical system. The one-electron wave function solutions are referred to as hydrogen-like atomic orbitals. Hydrogen-like atoms are of importance because their corresponding orbitals bear similarity to the hydrogen atomic orbitals.

Other systems may also be referred to as "hydrogen-like atoms", such as muonium (an electron orbiting an antimuon), positronium (an electron and a positron), certain exotic atoms (formed with other particles), or Rydberg atoms (in which one electron is in such a high energy state that it sees the rest of the atom effectively as a point charge).

## Schrödinger solution

In the solution to the Schrödinger equation, which is non-relativistic, hydrogen-like atomic orbitals are eigenfunctions of the one-electron angular momentum operator L and its z component Lz. A hydrogen-like atomic orbital is uniquely identified by the values of the principal quantum number n, the angular momentum quantum number l, and the magnetic quantum number m. The energy eigenvalues do not depend on l or m, but solely on n. To these must be added the two-valued spin quantum number ms = ±12, setting the stage for the Aufbau principle. This principle restricts the allowed values of the four quantum numbers in electron configurations of more-electron atoms. In hydrogen-like atoms all degenerate orbitals of fixed n and l, m and s varying between certain values (see below) form an atomic shell.

The Schrödinger equation of atoms or ions with more than one electron has not been solved analytically, because of the computational difficulty imposed by the Coulomb interaction between the electrons. Numerical methods must be applied in order to obtain (approximate) wavefunctions or other properties from quantum mechanical calculations. Due to the spherical symmetry (of the Hamiltonian), the total angular momentum J of an atom is a conserved quantity. Many numerical procedures start from products of atomic orbitals that are eigenfunctions of the one-electron operators L and Lz. The radial parts of these atomic orbitals are sometimes numerical tables or are sometimes Slater orbitals. By angular momentum coupling many-electron eigenfunctions of J2 (and possibly S2) are constructed.

In quantum chemical calculations hydrogen-like atomic orbitals cannot serve as an expansion basis, because they are not complete. The non-square-integrable continuum (E > 0) states must be included to obtain a complete set, i.e., to span all of one-electron Hilbert space. [1]

In the simplest model, the atomic orbitals of hydrogen-like atoms/ions are solutions to the Schrödinger equation in a spherically symmetric potential. In this case, the potential term is the potential given by Coulomb's law:

${\displaystyle V(r)=-{\frac {1}{4\pi \varepsilon _{0}}}{\frac {Ze^{2}}{r}}}$

where

After writing the wave function as a product of functions:

${\displaystyle \psi (r,\theta ,\phi )=R_{nl}(r)Y_{\ell m}(\theta ,\phi )}$

(in spherical coordinates), where ${\displaystyle Y_{\ell m}}$ are spherical harmonics, we arrive at the following Schrödinger equation:

${\displaystyle -{\frac {\hbar ^{2}}{2\mu }}\left[{\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial R(r)}{\partial r}}\right)-{\frac {l(l+1)R(r)}{r^{2}}}\right]+V(r)R(r)=ER(r),}$

where ${\displaystyle \mu }$ is, approximately, the mass of the electron (more accurately, it is the reduced mass of the system consisting of the electron and the nucleus), and ${\displaystyle \hbar }$ is the reduced Planck constant.

Different values of l give solutions with different angular momentum, where l (a non-negative integer) is the quantum number of the orbital angular momentum. The magnetic quantum number m (satisfying ${\displaystyle -l\leq m\leq l}$) is the (quantized) projection of the orbital angular momentum on the z-axis. See here for the steps leading to the solution of this equation.

### Non-relativistic wavefunction and energy

In addition to l and m, a third integer n > 0, emerges from the boundary conditions placed on R. The functions R and Y that solve the equations above depend on the values of these integers, called quantum numbers . It is customary to subscript the wave functions with the values of the quantum numbers they depend on. The final expression for the normalized wave function is:

${\displaystyle \psi _{n\ell m}=R_{n\ell }(r)\,Y_{\ell m}(\theta ,\phi )}$
${\displaystyle R_{n\ell }(r)=-{\sqrt {{\left({\frac {2Z}{na_{\mu }}}\right)}^{3}{\frac {(n-\ell -1)!}{2n{(n+\ell )!}}}}}e^{-Zr/{na_{\mu }}}\left({\frac {2Zr}{na_{\mu }}}\right)^{\ell }L_{n-\ell -1}^{(2\ell +1)}\left({\frac {2Zr}{na_{\mu }}}\right)}$

where:

• ${\displaystyle L_{n-\ell -1}^{(2\ell +1)}}$ are the generalized Laguerre polynomials.
• ${\displaystyle a_{\mu }={\frac {4\pi \varepsilon _{0}\hbar ^{2}}{\mu e^{2}}}={\frac {\hbar c}{\alpha \mu c^{2}}}={\frac {m_{\mathrm {e} }}{\mu }}a_{0}}$
where ${\displaystyle \alpha }$ is the fine-structure constant. Here, ${\displaystyle \mu }$ is the reduced mass of the nucleus-electron system, that is, ${\displaystyle \mu ={{m_{\mathrm {N} }m_{\mathrm {e} }} \over {m_{\mathrm {N} }+m_{\mathrm {e} }}}}$ where ${\displaystyle m_{\mathrm {N} }}$ is the mass of the nucleus. Typically, the nucleus is much more massive than the electron, so ${\displaystyle \mu \approx m_{\mathrm {e} }.}$ (But for positronium ${\displaystyle \mu =m_{\mathrm {e} }/2}$.) ${\displaystyle a_{0}}$ is the Bohr radius.
• ${\displaystyle E_{n}=-\left({\frac {Z^{2}\mu e^{4}}{32\pi ^{2}\epsilon _{0}^{2}\hbar ^{2}}}\right){\frac {1}{n^{2}}}=-\left({\frac {Z^{2}\hbar ^{2}}{2\mu a_{\mu }^{2}}}\right){\frac {1}{n^{2}}}=-{\frac {\mu c^{2}Z^{2}\alpha ^{2}}{2n^{2}}}.}$
• ${\displaystyle Y_{\ell m}(\theta ,\phi )\,}$ function is a spherical harmonic.

parity due to angular wave function is ${\displaystyle {\left({-1}\right)}^{\ell }}$.

### Quantum numbers

The quantum numbers ${\displaystyle n}$, ${\displaystyle \ell }$ and ${\displaystyle m}$ are integers and can have the following values:

${\displaystyle n=1,2,3,4,\dots }$
${\displaystyle \ell =0,1,2,\dots ,n-1}$
${\displaystyle m=-\ell ,-\ell +1,\ldots ,0,\ldots ,\ell -1,\ell }$

For a group-theoretical interpretation of these quantum numbers, see this article. Among other things, this article gives group-theoretical reasons why ${\displaystyle \ell and ${\displaystyle -\ell \leq m\leq \,\ell }$.

### Angular momentum

Each atomic orbital is associated with an angular momentum L. It is a vector operator, and the eigenvalues of its square L2Lx2 + Ly2 + Lz2 are given by:

${\displaystyle {\hat {L}}^{2}Y_{\ell m}=\hbar ^{2}\ell (\ell +1)Y_{\ell m}}$

The projection of this vector onto an arbitrary direction is quantized. If the arbitrary direction is called z, the quantization is given by:

${\displaystyle {\hat {L}}_{z}Y_{\ell m}=\hbar mY_{\ell m},}$

where m is restricted as described above. Note that L2 and Lz commute and have a common eigenstate, which is in accordance with Heisenberg's uncertainty principle. Since Lx and Ly do not commute with Lz, it is not possible to find a state that is an eigenstate of all three components simultaneously. Hence the values of the x and y components are not sharp, but are given by a probability function of finite width. The fact that the x and y components are not well-determined, implies that the direction of the angular momentum vector is not well determined either, although its component along the z-axis is sharp.

These relations do not give the total angular momentum of the electron. For that, electron spin must be included.

This quantization of angular momentum closely parallels that proposed by Niels Bohr (see Bohr model) in 1913, with no knowledge of wavefunctions.

### Including spin–orbit interaction

In a real atom, the spin of a moving electron can interact with the electric field of the nucleus through relativistic effects, a phenomenon known as spin–orbit interaction. When one takes this coupling into account, the spin and the orbital angular momentum are no longer conserved, which can be pictured by the electron precessing. Therefore, one has to replace the quantum numbers l, m and the projection of the spin ms by quantum numbers that represent the total angular momentum (including spin), j and mj, as well as the quantum number of parity.

See the next section on the Dirac equation for a solution that includes the coupling.

## Solution to Dirac equation

In 1928 in England Paul Dirac found an equation that was fully compatible with special relativity. The equation was solved for hydrogen-like atoms the same year (assuming a simple Coulomb potential around a point charge) by the German Walter Gordon. Instead of a single (possibly complex) function as in the Schrödinger equation, one must find four complex functions that make up a bispinor. The first and second functions (or components of the spinor) correspond (in the usual basis) to spin "up" and spin "down" states, as do the third and fourth components.

The terms "spin up" and "spin down" are relative to a chosen direction, conventionally the z direction. An electron may be in a superposition of spin up and spin down, which corresponds to the spin axis pointing in some other direction. The spin state may depend on location.

An electron in the vicinity of a nucleus necessarily has non-zero amplitudes for the third and fourth components. Far from the nucleus these may be small, but near the nucleus they become large.

The eigenfunctions of the Hamiltonian, which means functions with a definite energy (and which therefore do not evolve except for a phase shift), have energies characterized not by the quantum number n only (as for the Schrödinger equation), but by n and a quantum number j, the total angular momentum quantum number. The quantum number j determines the sum of the squares of the three angular momenta to be j(j+1) (times ħ2, see Planck constant). These angular momenta include both orbital angular momentum (having to do with the angular dependence of ψ) and spin angular momentum (having to do with the spin state). The splitting of the energies of states of the same principal quantum number n due to differences in j is called fine structure. The total angular momentum quantum number j ranges from 1/2 to n−1/2.

The orbitals for a given state can be written using two radial functions and two angle functions. The radial functions depend on both the principal quantum number n and an integer k, defined as:

${\displaystyle k={\begin{cases}-j-{\tfrac {1}{2}}&{\text{if }}j=\ell +{\tfrac {1}{2}}\\j+{\tfrac {1}{2}}&{\text{if }}j=\ell -{\tfrac {1}{2}}\end{cases}}}$

where ℓ is the azimuthal quantum number that ranges from 0 to n−1. The angle functions depend on k and on a quantum number m which ranges from −j to j by steps of 1. The states are labeled using the letters S, P, D, F et cetera to stand for states with ℓ equal to 0, 1, 2, 3 et cetera (see azimuthal quantum number), with a subscript giving j. For instance, the states for n=4 are given in the following table (these would be prefaced by n, for example 4S1/2):

m = −7/2m = −5/2m = −3/2m = −1/2m = 1/2m = 3/2m = 5/2m = 7/2
k = 3, ℓ = 3F5/2F5/2F5/2F5/2F5/2F5/2
k = 2, ℓ = 2D3/2D3/2D3/2D3/2
k = 1, ℓ = 1P1/2P1/2
k = 0
k = −1, ℓ = 0S1/2S1/2
k = −2, ℓ = 1P3/2P3/2P3/2P3/2
k = −3, ℓ = 2D5/2D5/2D5/2D5/2D5/2D5/2
k = −4, ℓ = 3F7/2F7/2F7/2F7/2F7/2F7/2F7/2F7/2

These can be additionally labeled with a subscript giving m. There are 2n2 states with principal quantum number n, 4j+2 of them with any allowed j except the highest (j=n−1/2) for which there are only 2j+1. Since the orbitals having given values of n and j have the same energy according to the Dirac equation, they form a basis for the space of functions having that energy.

The energy, as a function of n and |k| (equal to j+1/2), is:

${\displaystyle {\begin{array}{rl}E_{n\,j}&=\mu c^{2}\left(1+\left[{\dfrac {Z\alpha }{n-|k|+{\sqrt {k^{2}-Z^{2}\alpha ^{2}}}}}\right]^{2}\right)^{-1/2}\\&\\&\approx \mu c^{2}\left\{1-{\dfrac {Z^{2}\alpha ^{2}}{2n^{2}}}\left[1+{\dfrac {Z^{2}\alpha ^{2}}{n}}\left({\dfrac {1}{|k|}}-{\dfrac {3}{4n}}\right)\right]\right\}\end{array}}}$

(The energy of course depends on the zero-point used.) Note that if Z were able to be more than 137 (higher than any known element) then we would have a negative value inside the square root for the S1/2 and P1/2 orbitals, which means they would not exist. The Schrödinger solution corresponds to replacing the inner bracket in the second expression by 1. The accuracy of the energy difference between the lowest two hydrogen states calculated from the Schrödinger solution is about 9 ppm (90 μeV too low, out of around 10 eV), whereas the accuracy of the Dirac equation for the same energy difference is about 3 ppm (too high). The Schrödinger solution always puts the states at slightly higher energies than the more accurate Dirac equation. The Dirac equation gives some levels of hydrogen quite accurately (for instance the 4P1/2 state is given an energy only about 2×10−10 eV too high), others less so (for instance, the 2S1/2 level is about 4×10−6 eV too low). [2] The modifications of the energy due to using the Dirac equation rather than the Schrödinger solution is of the order of α2, and for this reason α is called the fine-structure constant.

The solution to the Dirac equation for quantum numbers n, k, and m, is:

${\displaystyle \Psi ={\begin{pmatrix}g_{n,k}(r)r^{-1}\Omega _{k,m}(\theta ,\phi )\\if_{n,k}(r)r^{-1}\Omega _{-k,m}(\theta ,\phi )\end{pmatrix}}={\begin{pmatrix}g_{n,k}(r)r^{-1}{\sqrt {(k+{\tfrac {1}{2}}-m)/(2k+1)}}Y_{k,m-1/2}(\theta ,\phi )\\-g_{n,k}(r)r^{-1}\operatorname {sgn} k{\sqrt {(k+{\tfrac {1}{2}}+m)/(2k+1)}}Y_{k,m+1/2}(\theta ,\phi )\\if_{n,k}(r)r^{-1}{\sqrt {(-k+{\tfrac {1}{2}}-m)/(-2k+1)}}Y_{-k,m-1/2}(\theta ,\phi )\\-if_{n,k}(r)r^{-1}\operatorname {sgn} k{\sqrt {(-k+{\tfrac {1}{2}}-m)/(-2k+1)}}Y_{-k,m+1/2}(\theta ,\phi )\end{pmatrix}}}$

where the Ωs are columns of the two spherical harmonics functions shown to the right. ${\displaystyle Y_{a,b}(\theta ,\phi )}$ signifies a spherical harmonic function:

${\displaystyle Y_{a,b}(\theta ,\phi )={\begin{cases}(-1)^{b}{\sqrt {{\frac {2a+1}{4\pi }}{\frac {(a-b)!}{(a+b)!}}}}P_{a}^{b}(\cos \theta )e^{ib\phi }&{\text{if }}a>0\\Y_{-a-1,b}(\theta ,\phi )&{\text{if }}a<0\end{cases}}}$

in which ${\displaystyle P_{a}^{b}}$ is an associated Legendre polynomial. (Note that the definition of Ω may involve a spherical harmonic that doesn't exist, like ${\displaystyle Y_{0,1}}$, but the coefficient on it will be zero.)

Here is the behavior of some of these angular functions. The normalization factor is left out to simplify the expressions.

${\displaystyle \Omega _{-1,-1/2}\propto {\binom {0}{1}}}$
${\displaystyle \Omega _{-1,1/2}\propto {\binom {1}{0}}}$
${\displaystyle \Omega _{1,-1/2}\propto {\binom {(x-iy)/r}{z/r}}}$
${\displaystyle \Omega _{1,1/2}\propto {\binom {z/r}{(x+iy)/r}}}$

From these we see that in the S1/2 orbital (k = −1), the top two components of Ψ have zero orbital angular momentum like Schrödinger S orbitals, but the bottom two components are orbitals like the Schrödinger P orbitals. In the P1/2 solution (k = 1), the situation is reversed. In both cases, the spin of each component compensates for its orbital angular momentum around the z axis to give the right value for the total angular momentum around the z axis.

The two Ω spinors obey the relationship:

${\displaystyle \Omega _{k,m}={\begin{pmatrix}z/r&(x-iy)/r\\(x+iy)/r&-z/r\end{pmatrix}}\Omega _{-k,m}}$

To write the functions ${\displaystyle g_{n,k}(r)}$ and ${\displaystyle f_{n,k}(r)}$ let us define a scaled radius ρ:

${\displaystyle \rho \equiv 2Cr}$

with

${\displaystyle C={\frac {\sqrt {\mu ^{2}c^{4}-E^{2}}}{\hbar c}}}$

where E is the energy (${\displaystyle E_{n\,j}}$) given above. We also define γ as:

${\displaystyle \gamma \equiv {\sqrt {k^{2}-Z^{2}\alpha ^{2}}}}$

When k = −n (which corresponds to the highest j possible for a given n, such as 1S1/2, 2P3/2, 3D5/2...), then ${\displaystyle g_{n,k}(r)}$ and ${\displaystyle f_{n,k}(r)}$ are:

${\displaystyle g_{n,-n}(r)=A(n+\gamma )\rho ^{\gamma }e^{-\rho /2}}$
${\displaystyle f_{n,-n}(r)=AZ\alpha \rho ^{\gamma }e^{-\rho /2}}$

where A is a normalization constant involving the gamma function:

${\displaystyle A={\frac {1}{\sqrt {2n(n+\gamma )}}}{\sqrt {\frac {C}{\gamma \Gamma (2\gamma )}}}}$

Notice that because of the factor Zα, f(r) is small compared to g(r). Also notice that in this case, the energy is given by

${\displaystyle E_{n,n-1/2}={\frac {\gamma }{n}}\mu c^{2}={\sqrt {1-{\frac {Z^{2}\alpha ^{2}}{n^{2}}}}}\,\mu c^{2}}$

and the radial decay constant C by

${\displaystyle C={\frac {Z\alpha }{n}}{\frac {\mu c^{2}}{\hbar c}}.}$

In the general case (when k is not −n), ${\displaystyle g_{n,k}(r){\text{ and }}f_{n,k}(r)}$ are based on two generalized Laguerre polynomials of order ${\displaystyle n-|k|-1}$ and ${\displaystyle n-|k|}$:

${\displaystyle g_{n,k}(r)=A\rho ^{\gamma }e^{-\rho /2}\left(Z\alpha \rho L_{n-|k|-1}^{(2\gamma +1)}(\rho )+(\gamma -k){\frac {\gamma \mu c^{2}-kE}{\hbar cC}}L_{n-|k|}^{(2\gamma -1)}(\rho )\right)}$
${\displaystyle f_{n,k}(r)=A\rho ^{\gamma }e^{-\rho /2}\left((\gamma -k)\rho L_{n-|k|-1}^{(2\gamma +1)}(\rho )+Z\alpha {\frac {\gamma \mu c^{2}-kE}{\hbar cC}}L_{n-|k|}^{(2\gamma -1)}(\rho )\right)}$

with A now defined as

${\displaystyle A={\frac {1}{\sqrt {2k(k-\gamma )}}}{\sqrt {{\frac {C}{n-|k|+\gamma }}{\frac {(n-|k|-1)!}{\Gamma (n-|k|+2\gamma +1)}}{\frac {1}{2}}\left(\left({\frac {Ek}{\gamma \mu c^{2}}}\right)^{2}+{\frac {Ek}{\gamma \mu c^{2}}}\right)}}}$

Again f is small compared to g (except at very small r) because when k is positive the first terms dominate, and α is big compared to γ−k, whereas when k is negative the second terms dominate and α is small compared to γ−k. Note that the dominant term is quite similar to corresponding the Schrödinger solution – the upper index on the Laguerre polynomial is slightly less (2γ+1 or 2γ−1 rather than 2ℓ+1, which is the nearest integer), as is the power of ρ (γ or γ−1 instead of ℓ, the nearest integer). The exponential decay is slightly faster than in the Schrödinger solution.

The normalization factor makes the integral over all space of the square of the absolute value equal to 1.

### 1S orbital

Here is the 1S1/2 orbital, spin up, without normalization:

${\displaystyle \Psi \propto {\begin{pmatrix}(1+\gamma )r^{\gamma -1}e^{-Cr}\\0\\iZ\alpha r^{\gamma -1}e^{-Cr}z/r\\iZ\alpha r^{\gamma -1}e^{-Cr}(x+iy)/r\end{pmatrix}}}$

Note that γ is a little less than 1, so the top function is similar to an exponentially decreasing function of r except that at very small r it theoretically goes to infinity. But the value of the ${\displaystyle r^{\gamma -1}}$ only surpasses 10 at a value of r smaller than ${\displaystyle 10^{1/(\gamma -1)},}$ which is a very small number (much less than the radius of a proton) unless Z is very large.

The 1S1/2 orbital, spin down, without normalization, comes out as:

${\displaystyle \Psi \propto {\begin{pmatrix}0\\(1+\gamma )r^{\gamma -1}e^{-Cr}\\iZ\alpha r^{\gamma -1}e^{-Cr}(x-iy)/r\\-iZ\alpha r^{\gamma -1}e^{-Cr}z/r\end{pmatrix}}}$

We can mix these in order to obtain orbitals with the spin oriented in some other direction, such as:

${\displaystyle \Psi \propto {\begin{pmatrix}(1+\gamma )r^{\gamma -1}e^{-Cr}\\(1+\gamma )r^{\gamma -1}e^{-Cr}\\iZ\alpha r^{\gamma -1}e^{-Cr}(x-iy+z)/r\\iZ\alpha r^{\gamma -1}e^{-Cr}(x+iy-z)/r\end{pmatrix}}}$

which corresponds to the spin and angular momentum axis pointing in the x direction. Adding i times the "down" spin to the "up" spin gives an orbital oriented in the y direction.

### 2P1/2 and 2S1/2 orbitals

To give another example, the 2P1/2 orbital, spin up, is proportional to:

${\displaystyle \Psi \propto {\begin{pmatrix}\rho ^{\gamma -1}e^{-\rho /2}\left(Z\alpha \rho +(\gamma -1){\frac {\gamma \mu c^{2}-E}{\hbar cC}}(-\rho +2\gamma )\right)z/r\\\rho ^{\gamma -1}e^{-\rho /2}\left(Z\alpha \rho +(\gamma -1){\frac {\gamma \mu c^{2}-E}{\hbar cC}}(-\rho +2\gamma )\right)(x+iy)/r\\i\rho ^{\gamma -1}e^{-\rho /2}\left((\gamma -1)\rho +Z\alpha {\frac {\gamma \mu c^{2}-E}{\hbar cC}}(-\rho +2\gamma )\right)\\0\end{pmatrix}}}$

(Remember that ${\displaystyle \rho =2rC}$. C is about half what it is for the 1S orbital, but γ is still the same.)

Notice that when ρ is small compared to α (or r is small compared to ${\displaystyle \hbar c/(\mu c^{2})}$) the "S" type orbital dominates (the third component of the bispinor).

For the 2S1/2 spin up orbital, we have:

${\displaystyle \Psi \propto {\begin{pmatrix}\rho ^{\gamma -1}e^{-\rho /2}\left(Z\alpha \rho +(\gamma +1){\frac {\gamma \mu c^{2}+E}{\hbar cC}}(-\rho +2\gamma )\right)\\0\\i\rho ^{\gamma -1}e^{-\rho /2}\left((\gamma +1)\rho +Z\alpha {\frac {\gamma \mu c^{2}+E}{\hbar cC}}(-\rho +2\gamma )\right)z/r\\i\rho ^{\gamma -1}e^{-\rho /2}\left((\gamma +1)\rho +Z\alpha {\frac {\gamma \mu c^{2}+E}{\hbar cC}}(-\rho +2\gamma )\right)(x+iy)/r\end{pmatrix}}}$

Now the first component is S-like and there is a radius near ρ = 2 where it goes to zero, whereas the bottom two-component part is P-like.

### Negative-energy solutions

In addition to bound states, in which the energy is less than that of an electron infinitely separated from the nucleus, there are solutions to the Dirac equation at higher energy, corresponding to an unbound electron interacting with the nucleus. These solutions are not normalizable, but solutions can be found which tend toward zero as r goes to infinity (which is not possible when ${\displaystyle |E|<\mu c^{2}}$ except at the above-mentioned bound-state values of E). There are similar solutions with ${\displaystyle E<-\mu c^{2}.}$ These negative-energy solutions are just like positive-energy solutions having the opposite energy but for a case in which the nucleus repels the electron instead of attracting it, except that the solutions for the top two components switch places with those for the bottom two.

Negative-energy solutions to Dirac's equation exist even in the absence of a Coulomb force exerted by a nucleus. Dirac hypothesized that we can consider almost all of these states to be already filled. If one of these negative-energy states is not filled, this manifests itself as though there is an electron which is repelled by a positively-charged nucleus. This prompted Dirac to hypothesize the existence of positively-charged electrons, and his prediction was confirmed with the discovery of the positron.

## Beyond Gordon's solution to the Dirac equation

The Dirac equation with a simple Coulomb potential generated by a point-like non-magnetic nucleus was not the last word, and its predictions differ from experimental results as mentioned earlier. More accurate results include the Lamb shift (radiative corrections arising from quantum electrodynamics) [3] and hyperfine structure.

## Notes

1. This was observed as early as 1928 by E. A. Hylleraas, Z. f. Physik vol. 48, p. 469 (1928). English translation in H. Hettema, Quantum Chemistry, Classic Scientific Papers, p. 81, World Scientific, Singapore (2000). Later it was pointed out again by H. Shull and P.-O. Löwdin, J. Chem. Phys. vol. 23, p. 1362 (1955).
2. Calculated from Table 4.1 in Felix Nendzig. "The Quantum Theory of the Hydrogen Atom" (PDF). Archived from the original (PDF) on October 20, 2013. Retrieved Oct 20, 2013.
3. For the radiative correction, see Nendzig, opus citatum.

## Related Research Articles

A hydrogen atom is an atom of the chemical element hydrogen. The electrically neutral atom contains a single positively charged proton and a single negatively charged electron bound to the nucleus by the Coulomb force. Atomic hydrogen constitutes about 75% of the baryonic mass of the universe.

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-12 massive particles, called "Dirac 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 structure of the hydrogen spectrum in a completely rigorous way.

In the quantum mechanics description of a particle in spherical coordinates, a spherically symmetric potential, is a potential that depends only on the distance between the particle and a defined center point. In particular, if the particle in question is an electron and the potential is derived from Coulomb's law, then the problem can be used to describe a hydrogen-like (one-electron) atom.

In atomic physics, hyperfine structure is defined by small shifts in otherwise degenerate energy levels and the resulting splittings in those energy levels of atoms, molecules, and ions, due to electromagnetic multipole interaction between the nucleus and electron clouds.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics. The Hamilton–Jacobi equation is particularly useful in identifying conserved quantities for mechanical systems, which may be possible even when the mechanical problem itself cannot be solved completely.

In rotordynamics, 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.

The Curie–Weiss law describes the magnetic susceptibility χ of a ferromagnet in the paramagnetic region above the Curie point:

In atomic physics, the electron magnetic moment, or more specifically the electron magnetic dipole moment, is the magnetic moment of an electron resulting from its intrinsic properties of spin and electric charge. The value of the electron magnetic moment is −9.2847647043(28)×10−24 J⋅T−1. The electron magnetic moment has been measured to an accuracy of 1.7×10−13 relative to the Bohr magneton.

In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

In quantum physics, the spin–orbit interaction is a relativistic interaction of a particle's spin with its motion inside a potential. A key example of this phenomenon is the spin–orbit interaction leading to shifts in an electron's atomic energy levels, due to electromagnetic interaction between the electron's magnetic dipole, its orbital motion, and the electrostatic field of the positively charged nucleus. This phenomenon is detectable as a splitting of spectral lines, which can be thought of as a Zeeman effect product of two relativistic effects: the apparent magnetic field seen from the electron perspective and the magnetic moment of the electron associated with its intrinsic spin. A similar effect, due to the relationship between angular momentum and the strong nuclear force, occurs for protons and neutrons moving inside the nucleus, leading to a shift in their energy levels in the nucleus shell model. In the field of spintronics, spin–orbit effects for electrons in semiconductors and other materials are explored for technological applications. The spin–orbit interaction is at the origin of magnetocrystalline anisotropy and the spin Hall effect.

In electromagnetism, the electromagnetic tensor or electromagnetic field tensor is a mathematical object that describes the electromagnetic field in spacetime. The field tensor was first used after the four-dimensional tensor formulation of special relativity was introduced by Hermann Minkowski. The tensor allows related physical laws to be written very concisely.

In quantum mechanics, the Hellmann–Feynman theorem relates the derivative of the total energy with respect to a parameter, to the expectation value of the derivative of the Hamiltonian with respect to that same parameter. According to the theorem, once the spatial distribution of the electrons has been determined by solving the Schrödinger equation, all the forces in the system can be calculated using classical electrostatics.

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

In mathematical physics, the Dirac algebra is the Clifford algebra . This was introduced by the mathematical physicist P. A. M. Dirac in 1928 in developing the Dirac equation for spin-½ particles with a matrix representation of the gamma matrices, which represent the generators of the algebra.

In mathematics, a Coulomb wave function is a solution of the Coulomb wave equation, named after Charles-Augustin de Coulomb. They are used to describe the behavior of charged particles in a Coulomb potential and can be written in terms of confluent hypergeometric functions or Whittaker functions of imaginary argument.

In linear algebra, a raising or lowering operator is an operator that increases or decreases the eigenvalue of another operator. In quantum mechanics, the raising operator is sometimes called the creation operator, and the lowering operator the annihilation operator. Well-known applications of ladder operators in quantum mechanics are in the formalisms of the quantum harmonic oscillator and angular momentum.

An LC circuit can be quantized using the same methods as for the quantum harmonic oscillator. An LC circuit is a variety of resonant circuit, and consists of an inductor, represented by the letter L, and a capacitor, represented by the letter C. When connected together, an electric current can alternate between them at the circuit's resonant frequency:

An electric dipole transition is the dominant effect of an interaction of an electron in an atom with the electromagnetic field.

In physics, relativistic angular momentum refers to the mathematical formalisms and physical concepts that define angular momentum in special relativity (SR) and general relativity (GR). The relativistic quantity is subtly different from the three-dimensional quantity in classical mechanics.

In mathematical physics, the Gordon decomposition of the Dirac current is a splitting of the charge or particle-number current into a part that arises from the motion of the center of mass of the particles and a part that arises from gradients of the spin density. It makes explicit use of the Dirac equation and so it applies only to "on-shell" solutions of the Dirac equation.