|Mathematical plots of hydrogen atomic orbitals
Above: left, the radial wave function for a 1s (100) atomic orbital of hydrogen plotted as a function of
distance from the atomic centre. (This function has been normalised to ensure that the integral sum
of all the probabilities is equal to 1). The Greek symbol rho (p) indicates distance from the centre
along a radius in units of the Bohr radius (the atomic radius for hydrogen in the ground state as on
the left gives the square of the radial wave function, which gives us the probability density function
(PDF). The PDF tells us the probability of encountering the electron at a given distance from the
atomic centre in any single measurement (or equivalently the proportion of measurements in which
the electron is encountered at that position in an ensemble of measurements carried out on a large
number of atoms prepared in the same quantum state). Notice that the electron is most likely to be
found at one Bohr radius from the centre, in approximate agreement with the classical atomic model.
The panel on the right shows a mathematical plot of the hydrogen orbital in 2D in greyscale (top) and
rgb colour (bottom, with blue indicating a lower probability, red a higher probability). To achieve this,
the radial probability function is multiplied by a function called a spherical harmonic, which tells us how
the radial distribution has to rotated about each axis to generate the 2D and 3D plots. The 2D plot is
essentially a cross-section of the 1s orbital. Also shown is our 3D model.
Above: left, radial wave functions and 2D plots of the 2s (200) hydrogen atomic orbital. Right, top -
greyscale computed 2D probability density distribution and our 3D model, bottom right. Note that
s-orbitals (1s, 2s, 3s, ...) are spherically symmetric, so that the radial probability distribution applies
along any radius of the atom - it is simply rotated about two of the axes. The 3D model this time shows
the natural fuzziness in the PDF of the orbital (the model for the 1s orbital omitted this and constructed
a hard smooth surface). All orbitals are in reality similarly diffuse.
When we get to the 2p orbital (210, i.e. n=2,
l=1, m=0 orbital) the radial probability function
no longer tells us all we need to know about
the geometry. The spherical harmonic causes
an angular deviation from a spherical shape.
The 3s orbital is spherically symmetric like the 1s
and 2s orbitals. However, notice that the peak of
the 3s PDF (bottom left) is further out from the
atomic centre at the origin than for the 2s orbital.
This trend continues - the orbitals get larger as
the first quantum number, the principle quantum
number (n) increases. This quantum number n
corresponds to the energy of the s-orbital, so
more energetic orbitals are larger and focused
further from the atomic centre.
Notice also that the PDF of the 1s orbital has no
zero outside the centre (but gradually decays to
zero after reaching a peak) but that the 2s orbital
has one such zero and the 3s has two such
zeros. These zeros form thin spherical shells at
which the electron is never found and they are
called nodal surfaces. The 2s nodal surface is
visible as a blue ring in-between two red rings
9where the electron is most likely to be found) in
the 2D colour plot, which is a section across the
orbital. There are n-1 nodal surfaces for an
Above the 3d (320, n=3, l=2, m=0) orbital. From the 3d orbital onwards, the shapes become more complex.
This orbital has two main lobes either side of a central torus.
In general a hydrogen atomic orbital has
n-l-1 nodal surfaces (l=0 for an s-orbital
giving n-1 nodal surfaces). Since, for a
p-orbital, l=1, there are n-2 nodal surfaces.
Orbitals for which n=l-1 have no nodal
surface and just a single peak. For the 3p
orbital there is one nodal surface (n=3, l=1,
so n-l-1 = 1) giving rise to the first minimum
in the PDF and separating the main lobes
from the smaller lobes toward the centre.
The 2D plots for these atomic orbitals were generated using software written by Bot for Windows using visual
C# for the .NET framework. Below is a screen shot of this software, called OrbPlotter, which has just finished
plotting a 2s orbital in greyscale:
The shapes of the hydrogen atom atomic orbitals are given by solving Schroedinger's wave equation for
an electron trapped inside a Coulomb potential well. The Coulomb force is the force due to the electric
attraction between two electric charges, in this case the attractive force between the negatively charged
electron and the positively charged nucleus (in the hydrogen atom the nucleus contains only a single
positively charged proton). A potential energy well is a force field that keeps a particle in place, rather
like a water well in which water has to be raised against gravity to lift it out (the physical water well is also
a gravitational potential well due to the gravitational force field of the Earth). In a similar manner, work
must be done (and energy supplied) to an electron to pull it away from the nucleus.
According to the Schroedinger equation, a particle behaves like a wave (as indeed they do) and so the
solutions to the equation, which tell us the behaviour of the particle, are waves, called wave functions.
Essentially, a particle in an energy well is like a water wave trapped in a harbour - it bounces about
between opposite 'walls' and interferes with its own reflection to establish a stationary wave which
appears not to be travelling from one wall to the next but simply moves up and down. The walls in the
case of the atom are provided by our Coulomb force-field - when the electron flies too far from the
nucleus it is pulled back toward it.
Now, on Cronodon we don't normally include much in the way of maths (unless it's something unusual),
since this can be found in standard textbooks for those who want to understand such technical things.
Thus, Bot shall not explain how the Schroedinger equation is solved for the hydrogen atom, but Bot will
give you some background information:
Here, the h with the bar across it (h-bar) is Planck's constant (h) divided by 2 x pi. This
constant determines the scale of quantisation, or the graininess of energy, since energy and
momentum are always found in multiples of Planck's constant. V is the force field, in this case
the Coulomb potential, m is the mass of the particle (the electron mass in this case) and E is
the energy of the system. A particle not confined by a force field is called a free particle and is
also described by the TDSWE. The TISWE is a second-order linear partial differential equation
and can be solved by a mathematical technique called separation of variables. Essentially this
separates out three solutions, one for each variable, the variables here being spatial
coordinates. The best method uses spherical polar coordinates, in which each point is
specified by its radial distance (r) from the origin, its azimuthal angle (phi) from the positive
x-axis, and its zenith angle (theta) with the positive z-axis.
In this case, the radial PDF is plotted along the long-axis (the vertical; axis) of the 2p orbital shown in
the 3D model above and so shows the probability density along one of the two lobes of the 'dumbbell'
shape. There are two other p-orbitals - the 211 and 212 orbitals (not shown). These three 2p orbitals
have very different shapes, however they are easily perturbed by neighbouring atoms in which case
they average out to form three 2p orbitals which have the same dumbbell shape but directed along a
different perpendicular axis in each case (hence these are labelled the 2px, 2py and 2pz orbitals).
Above: the position of the point P (on the surface of a sphere of radius r) can be described by
the standard Cartesian coordinates x, y and z. However, for problems involving spherical
symmetry the maths gets much easier if one uses spherical polars (these are also more intuitive
for such problems once you are used to them).
Separation of variables then gives us a solution (wave function) consisting of three functions
multiplied together: a radial function for r and an angular function for theta and an angular
function for phi. The two angular functions together form the spherical harmonic and the radial
function is as plotted in the figures above. The 2D plots used the radial function multiplied by
the angular function for theta.
Finally the functions have to be normalised, that is they must sum (by integration) to one. This
is because the functions describe the probability of finding the electron at a particular location
in space, and since the electron must be somewhere the probabilities must add to one! Both
the radial functions and the spherical harmonics must be normalised to give accurate
probabilities. The 2D plots only illustrate the shape of the orbital and so normalisation is not
important for this, since normalisation does not alter the shape only the scale on the axes of the
graph. To obtain radial plots with the right shape again normalisation is not required, however,
to make the vertical axis read correct probabilities the radial functions have been normalised.
Here, let Bot add a note of caution! Many textbooks state the incorrect normalisation, they are
in error! Only by plotting the functions does this become obvious (which is probably why those
authour's never noticed), and other authour's doing similar plots have found the same problem.
Without resorting to calculating each one directly (quite a bit of work) Bot found two quick
methods that gave (apparently) correct results (although one of these methods, which uses
recurrence formulae, appeared to be only accurate for s and p orbitals).
Note that atomic orbitals have no
definite surface. Instead we
visualise an orbital with the
darkest regions indicating the
highest probability of finding the
electron at that position. We can
thus visualise say the total area
in which the electron is 90% or
99% likely to be found simply by
altering the intensity of the plot.
Two different intensity scales
have been used here (not
Right: the 3d (m=2) orbital.
Note the two inner nodal planes (spheres) where the rings
left). The electron density also drops off toward zero in the
centre of the atom. Note, however, that the probability of
finding the electron is only nominally zero in these regions
at mathematical 'points' as shown in the radial plot. Hence,
we should not think of the electron being absent near
these regions - raising the intensity of the plot will narrow
the white rings visualised in the plot as this area is really a
very pale shade of grey where there is a low but
nevertheless definite probability of finding the electron.