The hydrogen atom is the simplest possible atom, comprising a single proton and a single electron.
To determine the wavefunction of the atom's electron, we must solve the 3D Schrödinger wave equation with the Coulomb potential of the single proton:
−
ℏ
2
m
(
d
2
d
x
2
+
d
2
d
y
2
+
d
2
d
z
2
)
ψ
(
x
,
y
,
z
)
+
V
(
x
,
y
,
z
)
ψ
(
x
,
y
,
z
)
=
E
ψ
(
x
,
y
,
z
)
{\displaystyle -{\frac {\hbar }{2m}}\left({\frac {d^{2}}{dx^{2}}}+{\frac {d^{2}}{dy^{2}}}+{\frac {d^{2}}{dz^{2}}}\right)\psi \left(x,y,z\right)+V\left(x,y,z\right)\psi \left(x,y,z\right)=E\psi \left(x,y,z\right)}
V
(
x
,
y
,
z
)
=
V
(
r
)
=
−
e
2
4
π
ϵ
0
r
{\displaystyle V\left(x,y,z\right)=V\left(r\right)=-{\frac {e^{2}}{4\pi \epsilon _{0}r}}}
Energy of an electron near a proton, taking both kinetic and potential energy into account
T
(
r
)
+
V
(
r
)
=
E
(
r
)
{\displaystyle \,\quad T(r)+V(r)=E(r)\,}
We re-express the wavefunction in spherical co-ordinates:
ψ
(
x
,
y
,
z
)
=
ψ
(
r
,
θ
,
ϕ
)
{\displaystyle \psi (x,y,z)=\psi (r,\theta ,\phi )\,}
If there is no angular dependence, we can simple solve for the radial component, r .
d
2
d
x
2
+
d
2
d
y
2
+
d
2
d
z
2
=
d
2
d
r
2
+
2
r
d
d
r
{\displaystyle {\frac {d^{2}}{dx^{2}}}+{\frac {d^{2}}{dy^{2}}}+{\frac {d^{2}}{dz^{2}}}={\frac {d^{2}}{dr^{2}}}+{\frac {2}{r}}{\frac {d}{dr}}}
The Schrödinger equation for the radial solution of the hydrogen atom is then:
−
ℏ
2
m
(
d
2
d
r
2
+
2
r
d
d
r
)
ψ
(
r
)
−
e
2
4
π
ϵ
0
r
ψ
(
r
)
=
E
ψ
(
r
)
{\displaystyle -{\frac {\hbar }{2m}}\left({\frac {d^{2}}{dr^{2}}}+{\frac {2}{r}}{\frac {d}{dr}}\right)\psi \left(r\right)-{\frac {e^{2}}{4\pi \epsilon _{0}r}}\psi \left(r\right)=E\psi \left(r\right)}
For the lowest energy state, this differential equation can be solved by the following trial solution:
ψ
(
r
)
=
ψ
0
e
−
a
r
{\displaystyle \psi (r)=\psi _{0}e^{-ar}\,}
Substituting in, we get:
−
ℏ
2
m
(
a
2
e
−
a
r
−
2
a
e
−
a
r
r
)
−
e
2
4
π
ϵ
0
r
e
−
a
r
=
E
e
−
a
r
{\displaystyle -{\frac {\hbar }{2m}}\left(a^{2}e^{-ar}-{\frac {2ae^{-ar}}{r}}\right)-{\frac {e^{2}}{4\pi \epsilon _{0}r}}e^{-ar}=Ee^{-ar}}
−
ℏ
2
m
(
a
2
−
2
a
r
)
−
e
2
4
π
ϵ
0
r
=
E
{\displaystyle -{\frac {\hbar }{2m}}\left(a^{2}-{\frac {2a}{r}}\right)-{\frac {e^{2}}{4\pi \epsilon _{0}r}}=E}
The coefficents of r 0 and r −1 must be separately satisfied:
−
ℏ
2
a
2
2
m
=
E
{\displaystyle -{\frac {\hbar ^{2}a^{2}}{2m}}=E}
−
ℏ
2
a
m
−
e
2
4
π
ϵ
0
=
0
{\displaystyle -{\frac {\hbar ^{2}a}{m}}-{\frac {e^{2}}{4\pi \epsilon _{0}}}=0}
We therefore have values for a and E :
a
=
e
2
m
4
π
ℏ
2
ϵ
0
≈
1
0.05
nm
{\displaystyle a={\frac {e^{2}m}{4\pi \hbar ^{2}\epsilon _{0}}}\approx {\frac {1}{0.05{\mbox{ nm}}}}}
E
=
−
m
e
4
32
π
2
ϵ
0
2
ℏ
2
≈
13.8
eV
{\displaystyle E=-{\frac {me^{4}}{32\pi ^{2}\epsilon _{0}^{2}\hbar ^{2}}}\approx 13.8{\mbox{ eV}}}
Higher energy levels have a more complex wavefunction:
ψ
n
(
r
)
=
ψ
0
e
−
a
n
r
L
n
(
r
)
{\displaystyle \psi _{n}(r)=\psi _{0}e^{-a_{n}r}L_{n}(r)\,}
where L n (r) is the Laguerre polynomial of order n . It is no important here how this result was reached - for a more thorough explanation, see the Quantum Mechanics Wikibook. The associated energies are:
E
n
=
−
m
e
2
32
n
2
π
2
ϵ
0
2
ℏ
2
≈
−
13.6
n
2
eV
,
n
=
1
,
2
,
3
,
…
{\displaystyle E_{n}=-{\frac {me^{2}}{32n^{2}\pi ^{2}\epsilon _{0}^{2}\hbar ^{2}}}\approx -{\frac {13.6}{n^{2}}}{\mbox{ eV}},\quad \quad n=1,2,3,\ldots }
The energy levels become closer and closer as n increases and the energy approaches zero. n = ∞ corresponds to the vacuum level at which the electron is free of the proton. Electrons can move between the states by absorbing (to move up) or emitting (to move down) the correct amount of energy.
For every n , a spherical solution exists, in which the wavefunction is invariant with angle. These spherical orbitals are known as the s -orbitals.
The complete solution for the wavefunction has four quantum numbers:
ψ
n
,
l
,
m
,
s
(
r
,
θ
,
ϕ
)
,
l
=
0
,
1
,
…
,
n
−
1
,
m
=
−
l
,
…
,
+
l
,
s
=
−
1
2
,
+
1
2
{\displaystyle \psi _{n,l,m,s}(r,\theta ,\phi ),\quad l=0,1,\ldots ,n-1,\quad m=-l,\ldots ,+l,\quad s=-{\frac {1}{2}},+{\frac {1}{2}}}
For example, when n = 2 and l = 1 there are three solutions corresponding to m = −1,0,1. These are known as the p -orbitals. Up to six electrons can be held in the p -orbitals - two for each solution, one spin-up and one spin-down.