Loading [MathJax]/jax/output/CommonHTML/jax.js

Kubic Harmonics (K)

Download a notebook that generates this page

The kubic harmonics (also known as cubic harmonics) are linear combinations of the spherical harmonics and irreducible representations of the cubic (Oh) point group. For s, p and d wave-functions the kubic harmonics are besides a different order the same as the tesseral harmonics. For higher angular momentum they are different.

The relation between the kubic harmonics and the tesseral harmonics is given by the matrix TZtoK, whereby the vector {Z(m=l)l,K(m=l+1)l,,K(m=l)l} is given as TZtoK{Y(m=l)l,Y(m=l+1)l,,Y(m=l)l}. Explicit forms of TZtoK for l equal or smaller to 6 are:

Tl=0ZtoK=(1)Tl=1ZtoK=(001100010)Tl=2ZtoK=(0000100100010000001010000)Tl=3ZtoK=(01000000000322052252203220000000100000005220322322052200000000010)Tl=4ZtoK=(000073200053200000010000005320007320122072200000000007220122010000000007220122000000000012207220001000000)Tl=5ZtoK=(00010000000010000000000000001580352803728372803528015800000000000100000000000218098205285280982021800000000000000010000000740324015241524032407400000000000001000)Tl=6ZtoK=(000000122000722000000000011400054000000005400011400000072200012200011240152403400000000000000340152401124000100000000000165160916052800000000000000528091601651600000100000000031605516031128000000000000003112805516031601000000000000) ### The following table shows the kubic harmonics up to l=6. We list the explicit function in terms of the directional cosines x, y and z. The plots show the surface defined by the equation r=K(m)lK(m)l. The color of the surface is according to the phase with red for positive and cyan for negative. We show a 3D image as well as a projection along the x, y and z direction. ###



l=0

ml=0 - a1g

K(0)0=12πK(0)0=12π



l=1

ml=1 - t1u

K(1)1=123πsin(θ)cos(ϕ)K(1)1=123πx

ml=0 - t1u

K(0)1=123πsin(θ)sin(ϕ)K(0)1=123πy

ml=1 - t1u

K(1)1=123πcos(θ)K(1)1=123πz



l=2

ml=2 - eg

K(2)2=1415πsin2(θ)cos(2ϕ)K(2)2=1415π(xy)(x+y)

ml=1 - eg

K(1)2=185π(3cos(2θ)+1)K(1)2=145π(x2+y22z2)

ml=0 - t2g

K(0)2=1215πsin(θ)cos(θ)sin(ϕ)K(0)2=1215πyz

ml=1 - t2g

K(1)2=1215πsin(θ)cos(θ)cos(ϕ)K(1)2=1215πxz

ml=2 - t2g

K(2)2=1415πsin2(θ)sin(2ϕ)K(2)2=1215πxy



l=3

ml=3 - a2u

K(3)3=14105πsin2(θ)cos(θ)sin(2ϕ)K(3)3=12105πxyz

ml=2 - t1u

K(2)3=1647π(20sin(3θ)cos3(ϕ)+3sin(θ)(cos(ϕ)5cos(3ϕ)))K(2)3=147πx(2x23(y2+z2))

ml=1 - t1u

K(1)3=1647π(20sin3(θ)sin(3ϕ)+3(sin(θ)+5sin(3θ))sin(ϕ))K(1)3=147πy(3x2+2y23z2)

ml=0 - t1u

K(0)3=187πcos(θ)(5cos(2θ)1)K(0)3=147πz(2z23(x2+y2))

ml=1 - t2u

K(1)3=116105πsin(θ)cos(ϕ)(2sin2(θ)cos(2ϕ)+3cos(2θ)+1)K(1)3=14105πx(yz)(y+z)

ml=2 - t2u

K(2)3=132105πsin(θ)sin(ϕ)(4sin2(θ)cos(2ϕ)+6cos(2θ)+2)K(2)3=14105πy(z2x2)

ml=3 - t2u

K(3)3=14105πsin2(θ)cos(θ)cos(2ϕ)K(3)3=14105πz(xy)(x+y)



l=4

ml=4 - a1g

K(4)4=13221π(5sin4(θ)cos(4ϕ)+35cos4(θ)30cos2(θ)+3)K(4)4=1421π(x43x2(y2+z2)+y43y2z2+z4)

ml=3 - eg

K(3)4=3165πsin2(θ)(7cos(2θ)+5)cos(2ϕ)K(3)4=385π(xy)(x+y)(x2+y26z2)

ml=2 - eg

K(2)4=13215π(7sin4(θ)cos(4ϕ)35cos4(θ)+30cos2(θ)3)K(2)4=1815π(x4+6z2(x2+y2)12x2y2+y42z4)

ml=1 - t1g

K(1)4=33235πsin(2θ)sin(ϕ)(2sin2(θ)cos(2ϕ)+3cos(2θ)+1)K(1)4=3435πyz(yz)(y+z)

ml=0 - t1g

K(0)4=33235πsin(2θ)cos(ϕ)(2sin2(θ)cos(2ϕ)+3cos(2θ)+1)K(0)4=3435πxz(z2x2)

ml=1 - t1g

K(1)4=31635πsin4(θ)sin(4ϕ)K(1)4=3435πxy(xy)(x+y)

ml=2 - t2g

K(2)4=3325πsin(2θ)sin(ϕ)(14sin2(θ)cos(2ϕ)+7cos(2θ)3)K(2)4=345πyz(6x2+y2+z2)

ml=3 - t2g

K(3)4=3645π(7sin(4θ)sin(ϕ)sin(2ϕ)+sin(2θ)(cos(ϕ)+7cos(3ϕ)))K(3)4=345πxz(x26y2+z2)

ml=4 - t2g

K(4)4=3165πsin2(θ)(7cos(2θ)+5)sin(2ϕ)K(4)4=345πxy(x2+y26z2)



l=5

ml=5 - eu

K(5)5=1161155πsin2(θ)cos(θ)(3cos(2θ)+1)sin(2ϕ)K(5)5=141155πxyz(x2+y22z2)

ml=4 - eu

K(4)5=316385πsin4(θ)cos(θ)sin(4ϕ)K(4)5=34385πxyz(xy)(x+y)

ml=3 - t(1)1u

K(3)5=125611πsin(θ)(30(21cos4(θ)14cos2(θ)+1)cos(ϕ)+35sin2(θ)(19cos2(θ))cos(3ϕ)+63sin4(θ)cos(5ϕ))K(3)5=11611πx(8x440x2(y2+z2)+15(y2+z2)2)

ml=2 - t(1)1u

K(2)5=125611πsin(θ)(63sin4(θ)sin(5ϕ)+35sin2(θ)(9cos2(θ)1)sin(3ϕ)+30(21cos4(θ)14cos2(θ)+1)sin(ϕ))K(2)5=11611πy(15x4+x2(30z240y2)+8y440y2z2+15z4)

ml=1 - t(1)1u

K(1)5=125611π(30cos(θ)+35cos(3θ)+63cos(5θ))K(1)5=11611πz(40z2(x2+y2)+15(x2+y2)2+8z4)

ml=0 - t(2)1u

K(0)5=3385π(16sin(θ)(5cos(2ϕ)3)cos3(ϕ)+sin(3θ)(14cos(ϕ)+39cos(3ϕ)5cos(5ϕ))+sin(5θ)(42cos(ϕ)27cos(3ϕ)+cos(5ϕ)))4096K(0)5=316385πx(y46y2z2+z4)

ml=1 - t(2)1u

K(1)5=3385π(8sin5(θ)sin(5ϕ)+2sin(θ)sin(ϕ)+7(sin(3θ)+3sin(5θ))sin(ϕ)12sin3(θ)(9cos(2θ)+7)sin(3ϕ))2048K(1)5=316385πy(x46x2z2+z4)

ml=2 - t(2)1u

K(2)5=316385πsin4(θ)cos(θ)cos(4ϕ)K(2)5=316385πz(x46x2y2+y4)

ml=3 - t2u

K(3)5=11281155πsin(θ)(2(21cos4(θ)14cos2(θ)+1)cos(ϕ)+sin2(θ)(19cos2(θ))cos(3ϕ)3sin4(θ)cos(5ϕ))K(3)5=181155πx(2x2(yz)(y+z)y4+z4)

ml=4 - t2u

K(4)5=11281155πsin(θ)(sin2(θ)(3sin2(θ)sin(5ϕ)+(19cos2(θ))sin(3ϕ))+(42cos4(θ)+28cos2(θ)2)sin(ϕ))K(4)5=181155πy(xz)(x+z)(x22y2+z2)

ml=5 - t2u

K(5)5=1161155πsin2(θ)cos(θ)(3cos(2θ)+1)cos(2ϕ)K(5)5=181155πz(xy)(x+y)(x2+y22z2)



l=6

ml=6 - a1g

K(6)6=164132π(21sin4(θ)(111cos2(θ))cos(4ϕ)+231cos6(θ)315cos4(θ)+105cos2(θ)5)K(6)6=18132π(2x615x4(y2+z2)15x2(y412y2z2+z4)+(y2+z2)(2y417y2z2+2z4))

ml=5 - a2g

K(5)6=1128150152πsin2(θ)((33cos4(θ)+18cos2(θ)1)cos(2ϕ)+sin4(θ)cos(6ϕ))K(5)6=18150152π(xy)(x+y)(xz)(x+z)(yz)(y+z)

ml=4 - eg

K(4)6=11282732πsin2(θ)(5(33cos4(θ)18cos2(θ)+1)cos(2ϕ)+11sin4(θ)cos(6ϕ))K(4)6=182732π(xy)(x+y)(x45z2(x2+y2)9x2y2+y4+5z4)

ml=3 - eg

K(3)6=912π(48sin4(θ)(11cos(2θ)+9)cos(4ϕ)+105cos(2θ)+126cos(4θ)+231cos(6θ)+50)2048K(3)6=18912π(x615z2(x4+y4)+15z4(x2+y2)+y62z6)

ml=2 - t1g

K(2)6=312891πsin(θ)cos(θ)(11sin4(θ)sin(5ϕ)+5sin2(θ)(311cos2(θ))sin(3ϕ)+2(33cos4(θ)30cos2(θ)+5)sin(ϕ))K(2)6=3891πyz(yz)(y+z)(10x2+y2+z2)

ml=1 - t1g

K(1)6=312891πsin(θ)cos(θ)(2(33cos4(θ)30cos2(θ)+5)cos(ϕ)+5sin2(θ)(311cos2(θ))cos(3ϕ)+11sin4(θ)cos(5ϕ))K(1)6=3891πxz(xz)(x+z)(x210y2+z2)

ml=0 - t1g

K(0)6=36491πsin4(θ)(11cos(2θ)+9)sin(4ϕ)K(0)6=3891πxy(xy)(x+y)(x2+y210z2)

ml=1 - t(1)2g

K(1)6=125613652πsin(θ)cos(θ)(33sin4(θ)sin(5ϕ)+9sin2(θ)(311cos2(θ))sin(3ϕ)+2(33cos4(θ)30cos2(θ)+5)sin(ϕ))K(1)6=11613652πyz(16x416x2(y2+z2)+(y2+z2)2)

ml=2 - t(1)2g

K(2)6=125613652πsin(θ)cos(θ)(2(33cos4(θ)30cos2(θ)+5)cos(ϕ)+9sin2(θ)(11cos2(θ)3)cos(3ϕ)+33sin4(θ)cos(5ϕ))K(2)6=11613652πxz(x4+2x2(z28y2)+16y416y2z2+z4)

ml=3 - t(1)2g

K(3)6=125613652πsin2(θ)(60cos(2θ)+33cos(4θ)+35)sin(2ϕ)K(3)6=11613652πxy(16z2(x2+y2)+(x2+y2)2+16z4)

ml=4 - t(2)2g

K(4)6=30032π(3(5sin(2θ)+12sin(4θ)+33sin(6θ))sin(ϕ)+48sin5(θ)cos(θ)sin(5ϕ)+20sin3(θ)(21cos(θ)+11cos(3θ))sin(3ϕ))4096K(4)6=11630032πyz(3y410y2z2+3z4)

ml=5 - t(2)2g

K(5)6=30032π(240sin(2θ)sin4(ϕ)cos(ϕ)+48sin(4θ)sin2(ϕ)(7cos(ϕ)+cos(3ϕ))+sin(6θ)(198cos(ϕ)+55cos(3ϕ)+3cos(5ϕ)))8192K(5)6=11630032πxz(3x410x2z2+3z4)

ml=6 - t(2)2g

K(6)6=13230032πsin6(θ)sin(6ϕ)K(6)6=11630032πxy(3x410x2y2+3y4)

Different orbital basis sets used

Print/export