Crystal field operator
In Crystal-field theory the interactions of a local atom with its environment (covalent bonding) is approximated by an effective potential. These potentials do not really exist in a solid, but should be seen as effective fields describing the (anti-) bonding states of an atom. We assume that the local spin-orbitals are given by a radial equation times a spherical harmonic describing the angular part. φn,l,m(→r)=Rn,l(r)Yl,m(θ,ϕ), with r=|→r| and the spherical harmonics having the same phase as used in Mathematica. The Hamiltonian is given by a potential V(→r). In second quantization this yields the operator: HCF=∑τ1,τ2⟨φτ1(→r)|V(→r)|φτ2(→r)⟩a†τ1a†τ2, whereby the quantum number τ stands for the combination of n (principle quantum number) l (angular momentum) m (projection of angular momentum on z axis) σ (projection of the spin on the z axis).
In order to evaluate ⟨φτ1(→r)|V(→r)|φτ2(→r)⟩ we expand V(→r) on renormalized spherical Harmonics: Ck,m(θ,ϕ)=√4π2k+1Yk,m(θ,ϕ). Using a tailor series in rk we find: V(r,θ,ϕ)=∞∑k=0k∑m=−k1√(k−m)!1√(k+m)!×∂k−|m|z(−Sign[m]∂x+ı∂y)|m|V(r,θ,ϕ)|r=0×rkCk,m(θ,ϕ). The integral ⟨φτ1(→r)|V(→r)|φτ2(→r)⟩ can now be rewritten into a r dependent part and a θ and ϕ dependent part: ⟨φτ1(→r)|V(→r)|φτ2(→r)⟩=∞∑k=0k∑m=−kAk,m⟨Yl1,m1|Ck,m|Yl2,m2⟩, with: Ak,m=1√(k−m)!1√(k+m)!×⟨Rn1,l1(r)|∂k−|m|z(−Sign[m]∂x+ı∂y)|m|V(r,θ,ϕ)|r=0|Rn2,l2(r)⟩. In crystal-field theory V(→r) is an effective potential which is taken to fit the experiment, in practice one takes the values Ak,m as the fitting parameters. The crystal-field Hamiltonian thus becomes: HCF=∑τ1,τ2∑k,mAk,m⟨Yl1,m1|Ck,m|Yl2,m2⟩a†τ1a†τ2.
In Quanty one can create crystal-field operators with the function “NewOperator()” and as a first input the string “CF”. The function furthermore needs to know the effective potential expanded on renormalized spherical harmonics. This is given as a list of the form “{{k1,m1,Ak1,m1},{k2,m2,Ak2,m2},...}”.
- Example.Quanty
-- crystal field operator in Oh -- symmetry acting on a d shell NF = 10 NB = 0 IndexDn_3d = {0,2,4,6,8} IndexUp_3d = {1,3,5,7,9} tenDq = 1.1 Akm = {{4, 0,(21/10)*tenDq}, {4,-4,(3/2)*math.sqrt(7/10)*tenDq}, {4, 4,(3/2)*math.sqrt(7/10)*tenDq}} OppCF = NewOperator("CF", NF, IndexUp_3d, IndexDn_3d, Akm)
The crystal-field operator is also used to create transition operators in spectroscopy. Also here one needs to create an operator that represents an electrical field (of the photon). In this case however the field has a physical representation and one can calculate the true expansion of the field on renormalized spherical harmonics. The term rk for large k can lead to divergent or badly converging series and therefor it is often desirable to expand on spherical Bessel functions jk(r) as will be discussed in more detail below.