Hydrogen-like atom

In physics and chemistry, a hydrogen-like atom (or hydrogenic atom) is an atom with one electron. Except for the hydrogen atom itself (which is neutral) these atoms carry positive charge e(Z-1), where Z is the atomic number of the atom and e is the elementary charge. Because hydrogen-like atoms are two-particle systems with an interaction depending only on the distance between the two particles, their non-relativistic Schrödinger equation can be solved in analytic form. The solutions are one-electron functions and are referred to as hydrogen-like atomic orbitals. These orbitals differ from one another in one respect only: the nuclear charge eZ appears in the radial part of the wave function.

Quantum numbers

Hydrogen-like atomic orbitals are eigenfunctions of a Hamiltonian (energy operator) with eigenvalues proportional to 1/n², where n is a positive integer. Further the obitals are usually chosen such that they are also eigenfunctions of the square of the one-electron angular momentum vector operator l ≡ (l_x, l_y, l_z) and of l_z. The operator l ² ≡ l_x² + l_y² + l_z² has an eigenvalue proportional to l(l+1), where l is a non-negative integer. The operator l_z has an eigenvalue proportional to an integer usually denoted by m.

Note that the energies of the hydrogen-like orbitals do not depend on the angular momentum quantum numbers l and m, but solely on the principal quantum number n. The degeneracy (maximum number of linearly independent eigenfunctions of same energy) of energy level n is equal to n². This is the dimension of the irreducible representations of the symmetry group of hydrogen-like atoms, which is SO(4).

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. These quantum numbers are integers and we summarize their ranges:

${\begin{aligned}n&=1,2,3,4,\ldots ,\\l&=0,1,2,\ldots ,n-1,\\m&=-l,-l+1,\ldots ,l-1,l.\end{aligned}}$

To this must be added the two-valued spin quantum number ms = ±½ in application of the exclusion principle. This principle restricts the allowed values of the four quantum numbers in electron configurations of more-electron atoms: it is forbidden that two electrons have the same four quantum numbers.

Completeness of hydrogen-like orbitals

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.

Schrödinger equation

The atomic orbitals of hydrogen-like atoms are solutions of the time-independent Schrödinger equation in a potential given by Coulomb's law:

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

where

• ε₀ is the permittivity of the vacuum,
• Z is the atomic number (charge of the nucleus in unit e),
• e is the elementary charge (charge of an electron),
• r is the distance of the electron from the nucleus.

After writing the wave function as a product of functions:

$\psi (r,\theta ,\phi )=R(r)Y_{lm}(\theta ,\phi )\,$

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

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

where μ is, approximately, the mass of the electron. More accurately, it is the reduced mass of the system consisting of the electron and the nucleus. Because the electron mass is about 1836 smaller than the mass of the lightest nucleus (the proton), the value of μ is very close to the mass of the electron m_e for all hydrogenic atoms. In the remaining of the article we will often make the approximation μ = m_e. Since m_e will appear explicitly in the formulas it will be easy to correct for this approximation if necessary.

Wave function and energy

In addition to l and m, there arises a third integer n > 0 from the boundary conditions imposed on R(r). The final expression for the normalized wave function is:

$\psi _{nlm}=R_{nl}(r)\,Y_{lm}(\theta ,\phi )$

$R_{nl}(r)={\sqrt {{\left({\frac {2Z}{na_{\mu }}}\right)}^{3}{\frac {(n-l-1)!}{2n[(n+l)!]}}}}e^{-Zr/{na_{\mu }}}\left({\frac {2Zr}{na_{\mu }}}\right)^{l}L_{n-l-1}^{2l+1}\left({\frac {2Zr}{na_{\mu }}}\right)$

where:

• $L_{n-l-1}^{2l+1}$ are the generalized Laguerre polynomials in the definition given here.
• $a_{\mu }={{4\pi \varepsilon _{0}\hbar ^{2}} \over {\mu e^{2}}}$

Note that $a_{\mu }\,$ is approximately equal to $a_{0}\,$ (the Bohr radius). If the mass of the nucleus is infinite then $\mu =m_{e}\,$ and $a_{\mu }=a_{0}\,$.

• $E_{n}=-{\frac {\mu }{2}}\left({\frac {e^{2}}{4\pi \varepsilon _{0}\hbar n}}\right)^{2}$. (Energy eigenvalues. As we pointed out above they depend only on n, not on l or m).

• $Y_{lm}(\theta ,\phi )\,$ function is a spherical harmonic.

Solution

The kinetic energy operator in spherical polar coordinates is

${\frac {{\hat {p}}^{2}}{2m_{e}}}=-{\frac {\hbar ^{2}}{2m_{e}}}\nabla ^{2}=-{\frac {\hbar ^{2}}{2m_{e}\,r^{2}}}\left[{\frac {\partial }{\partial r}}{\Big (}r^{2}{\frac {\partial }{\partial r}}{\Big )}-{\hat {l}}^{2}\right].$

The spherical harmonics satisfy

${\hat {l}}^{2}Y_{lm}(\theta ,\phi )\equiv \left\{-{\frac {1}{\sin ^{2}\theta }}\left[\sin \theta {\frac {\partial }{\partial \theta }}{\Big (}\sin \theta {\frac {\partial }{\partial \theta }}{\Big )}+{\frac {\partial ^{2}}{\partial \phi ^{2}}}\right]\right\}Y_{lm}(\theta ,\phi )=l(l+1)Y_{lm}(\theta ,\phi )$.

Substituting this into the Schrödinger equation we get a one-dimensional eigenvalue equation,

$\left\{-{\hbar ^{2} \over 2m_{e}r^{2}}{d \over dr}(r^{2}{d \over dr})+{\hbar ^{2}l(l+1) \over 2m_{e}r^{2}}+V(r)\right\}R(r)=ER(r)$

Note that the first term in the kinetic energy can be rewritten

$T_{r}\equiv {\frac {1}{r^{2}}}{\frac {d}{dr}}r^{2}{\frac {d}{dr}}={\frac {1}{r}}{\frac {d^{2}}{dr^{2}}}r$.

This follows because both sides of this equation can be shown by application of the product rule to be equal to a third form of this operator:

$T_{r}={\frac {d^{2}}{dr^{2}}}+{\frac {2}{r}}{\frac {d}{dr}}.$

If subsequently the substitution $u(r)\ {\stackrel {\mathrm {def} }{=}}\ rR(r)$ is made into

$-{\hbar ^{2} \over 2m_{e}}{\frac {1}{r}}{d^{2} \over dr^{2}}r\,R(r)+\left(V(r)+{\hbar ^{2}l(l+1) \over 2m_{e}r^{2}}\right)\,R(r)=E\,R(r),$

the radial equation becomes

$-{\hbar ^{2} \over 2m_{e}}{d^{2}u(r) \over dr^{2}}+V_{\mathrm {eff} }(r)u(r)=Eu(r)$

which is precisely a Schrödinger equation for the function u(r) with an effective potential given by

$V_{\mathrm {eff} }(r)=V(r)+{\hbar ^{2}l(l+1) \over 2m_{e}r^{2}},$

where the radial coordinate r ranges from 0 to $\infty$. The correction to the potential V(r) is called the centrifugal barrier term.

In order to simplify the Schrödinger equation, we introduce the following constants that define the atomic unit of energy and length, respectively,

$E_{\textrm {h}}=m_{e}\left({\frac {e^{2}}{4\pi \varepsilon _{0}\hbar }}\right)^{2}\quad {\hbox{and}}\quad a_{0}={{4\pi \varepsilon _{0}\hbar ^{2}} \over {m_{e}e^{2}}}$.

Substitute $y=Zr/a_{0}\,$ and $W=E/(Z^{2}E_{\textrm {h}})\,$ into the radial Schrödinger equation given above. This gives an equation in which all natural constants are hidden,

$\left[-{\frac {1}{2}}{\frac {d^{2}}{dy^{2}}}+{\frac {1}{2}}{\frac {l(l+1)}{y^{2}}}-{\frac {1}{y}}\right]u_{l}=Wu_{l}.$

Two classes of solutions of this equation exist: (i) W is negative, the corresponding eigenfunctions are square integrable and the values of W are quantized (discrete spectrum). (ii) W is non-negative. Every real non-negative value of W is physically allowed (continuous spectrum), the corresponding eigenfunctions are non-square integrable. In the remaining part of this article only class (i) solutions will be considered. The wavefunctions are known as bound states, in contrast to the class (ii) solutions that are known as scattering states.

For negative W the quantity $\alpha \equiv 2{\sqrt {-2W}}$ is real and positive. The scaling of y, i.e., substitution of $x\equiv \alpha y$ gives the Schrödinger equation:

$\left[{\frac {d^{2}}{dx^{2}}}-{\frac {l(l+1)}{x^{2}}}+{\frac {2}{\alpha x}}-{\frac {1}{4}}\right]u_{l}=0,\quad {\hbox{with}}\quad x\geq 0.$

For $x\rightarrow \infty$ the inverse powers of x are negligible and a solution for large x is $\exp[-x/2]$. The other solution, $\exp[x/2]$, is physically non-acceptable. For $x\rightarrow 0$ the inverse square power dominates and a solution for small x is x^l+1. The other solution, x^-l, is physically non-acceptable. Hence, to obtain a full range solution we substitute

$u_{l}(x)=x^{l+1}e^{-x/2}f_{l}(x).\,$

The equation for f_l(x) becomes,

$\left[x{\frac {d^{2}}{dx^{2}}}+(2l+2-x){\frac {d}{dx}}+(\nu -l-1)\right]f_{l}(x)=0\quad {\hbox{with}}\quad \nu =(-2W)^{-{\frac {1}{2}}}.$

Provided $\nu -l-1$ is a non-negative integer, say k, this equation has polynomial solutions written as

$L_{k}^{(2l+1)}(x),\qquad k=0,1,\ldots ,$

which are generalized Laguerre polynomials of order k. We will take the convention for generalized Laguerre polynomials of Abramowitz and Stegun. Note that the Laguerre polynomials given in many quantum mechanical textbooks, for instance the book of Messiah, are those of Abramowitz and Stegun multiplied by a factor (2l+1+k)! The definition given in this article coincides with the one of Abramowitz and Stegun.

The energy becomes

$W=-{\frac {1}{2n^{2}}}\quad {\hbox{with}}\quad n\equiv k+l+1.$

The principal quantum number n satisfies $n\geq l+1$, or $l\leq n-1$. Since $\alpha =2/n$, the total radial wavefunction is

$R_{nl}(r)=N_{nl}\left({\tfrac {2Zr}{na_{0}}}\right)^{l}\;e^{-{\textstyle {\frac {Zr}{na_{0}}}}}\;L_{n-l-1}^{(2l+1)}\left({\tfrac {2Zr}{na_{0}}}\right),$

with normalization constant

$N_{nl}=\left[\left({\frac {2Z}{na_{0}}}\right)^{3}\cdot {\frac {(n-l-1)!}{2n[(n+l)!]}}\right]^{1 \over 2}$

which belongs to the energy

$E_{n}=-{\frac {Z^{2}}{2n^{2}}}E_{\textrm {h}},\qquad n=1,2,\ldots .$

In the computation of the normalization constant use was made of the integral

$\int _{0}^{\infty }x^{2l+2}e^{-x}\left[L_{n-l-1}^{(2l+1)}(x)\right]^{2}dx={\frac {2n(n+l)!}{(n-l-1)!}}.$