In fields like computational chemistry and solid-state and condensed matter physics the so-called atomic orbitals, or spin-orbitals, as they appear in textbooks[1] [2] [3] on quantum physics, are often partially replaced by cubic harmonics for a number of reasons. These harmonics are usually named tesseral harmonics in the field of condensed matter physics in which the name kubic harmonics rather refers to the irreducible representations in the cubic point-group. [4]
The
2l+1
n
l
\psinlm(r)=Rnl(r)
m(\theta,\varphi) | |
Y | |
l |
in which the
Rnl(r)
m(\theta,\varphi) | |
Y | |
l |
m(\theta,\varphi) | |
Y | |
l |
In many cases, especially in chemistry and solid-state and condensed-matter physics, the system under investigation doesn't have rotational symmetry. Often it has some kind of lower symmetry, with a special point group representation, or it has no spatial symmetry at all. Biological and biochemical systems, like amino acids and enzymes often belong to low molecular symmetry point groups. The solid crystals of the elements often belong to the space groups and point groups with high symmetry. (Cubic harmonics representations are often listed and referenced in point group tables.) The system has at least a fixed orientation in three-dimensional Euclidean space. Therefore, the coordinate system that is used in such cases is most often a Cartesian coordinate system instead of a spherical coordinate system. In a Cartesian coordinate system the atomic orbitals are often expressed as
\psinlc(r)=Rnl(r)Xlc(r)
with the cubic harmonics,[6] [7] [8]
Xlc(r)
For the representations of the spherical harmonics a spherical coordinate system is chosen with a principal axis in the z-direction. For the cubic harmonics this axis is also the most convenient choice. For states of higher angular momentum quantum number
l
l(l+1)
l(l+1)
Xlc(r)=
0 | |
Y | |
l |
Xlc'(r)=
1 | ||||||
|
Xlc''(r)=
1 | ||||||
|
A substantial number of the spherical harmonics are listed in the Table of spherical harmonics.
First of all, the cubic harmonics are real functions, while spherical harmonics are complex functions. The complex numbers are two-dimensional with a real part and an imaginary part. Complex numbers offer very handsome and effective tools to tackle mathematical problems analytically but they are not very effective when they are used for numerical calculations. Skipping the imaginary part saves half the calculational effort in summations, a factor of four in multiplications and often factors of eight or even more when it comes to computations involving matrices.
The cubic harmonics often fit the symmetry of the potential or surrounding of an atom. A common surrounding of atoms in solids and chemical complexes is an octahedral surrounding with an octahedral cubic point group symmetry. The representations of the cubic harmonics often have a high symmetry and multiplicity so operations like integrations can be reduced to a limited, or irreducible, part of the domain of the function that has to be evaluated. A problem with the 48-fold octahedral Oh symmetry can be calculated much faster if one limits a calculation, like an integration, to the irreducible part of the domain of the function.
The s-orbitals only have a radial part.
\psin00(r)=Rn0(r)
0 | |
Y | |
0 |
s=X00=
0 | |
Y | |
0 |
=
1 | |
\sqrt{4\pi |
The three p-orbitals are atomic orbitals with an angular momentum quantum number ℓ = 1. The cubic harmonic expression of the p-orbitals
pz=
c | |
N | |
1 |
z | |
r |
=
0 | |
Y | |
1 |
px=
c | |
N | |
1 |
x | |
r |
=
1 | |
\sqrt{2 |
py=
c | |
N | |
1 |
y | |
r |
=
i | |
\sqrt{2 |
c | |
N | |
1 |
=\left(
3 | |
4\pi |
\right)1/2
The five d-orbitals are atomic orbitals with an angular momentum quantum number ℓ = 2. The angular part of the d-orbitals are often expressed like
\psin2c(r)=Rn2(r)X2c(r)
The angular part of the d-orbitals are the cubic harmonics
X2c(r)
d | |
z2 |
=
c | |
N | |
2 |
3z2-r2 | |
2r2\sqrt{3 |
dxz=
c | |
N | |
2 |
xz | |
r2 |
=
1 | |
\sqrt{2 |
dyz=
c | |
N | |
2 |
yz | |
r2 |
=
i | |
\sqrt{2 |
dxy=
c | |
N | |
2 |
xy | |
r2 |
=
i | |
\sqrt{2 |
d | |
x2-y2 |
=
c | |
N | |
2 |
x2-y2 | |
2r2 |
=
1 | |
\sqrt{2 |
c | |
N | |
2 |
=\left(
15 | |
4\pi |
\right)1/2
The seven f-orbitals are atomic orbitals with an angular momentum quantum number ℓ = 3. often expressed like
\psin3c(r)=Rn3(r)X3c(r)
The angular part of the f-orbitals are the cubic harmonics
X3c(r)
f | |
z3 |
=
c | |
N | |
3 |
z(2z2-3x2-3y2) | |
2r3\sqrt{15 |
f | |
xz2 |
=
c | |
N | |
3 |
x(4z2-x2-y2) | |
2r3\sqrt{10 |
f | |
yz2 |
=
c | |
N | |
3 |
y(4z2-x2-y2) | |
2r3\sqrt{10 |
fxyz=
c | |
N | |
3 |
xyz | |
r3 |
=
i | |
\sqrt{2 |
f | |
z(x2-y2) |
=
c | |
N | |
3 |
z\left(x2-y2\right) | |
2r3 |
=
1 | |
\sqrt{2 |
f | |
x(x2-3y2) |
=
c | |
N | |
3 |
x\left(x2-3y2\right) | |
2r3\sqrt{6 |
f | |
y(3x2-y2) |
=
c | |
N | |
3 |
y\left(3x2-y2\right) | |
2r3\sqrt{6 |
c | |
N | |
3 |
=\left(
105 | |
4\pi |
\right)1/2