Processing math: 86%

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

physics_chemistry:point_groups:d4h:orientation_zxy [2018/03/21 18:32] – created Stefano Agrestiniphysics_chemistry:point_groups:d4h:orientation_zxy [2018/04/06 09:03] (current) Maurits W. Haverkort
Line 1: Line 1:
 +~~CLOSETOC~~
 +
 ====== Orientation Zxy ====== ====== Orientation Zxy ======
 +
 +===== Symmetry Operations =====
  
 ### ###
-alligned paragraph text+ 
 +In the D4h Point Group, with orientation Zxy there are the following symmetry operations 
 ### ###
  
-===== Example =====+### 
 + 
 +{{:physics_chemistry:pointgroup:d4h_zxy.png}}
  
 ### ###
-description text+
 ### ###
  
-==== Input ==== +^ Operator ^ Orientation ^ 
-<code Quanty Example.Quanty> +^ E | {0,0,0} , | 
--- some example code+^ C4 | {0,0,1} , {0,0,1} , | 
 +^ C2 | {0,0,1} , | 
 +^ C2 | {0,1,0} , {1,0,0} , | 
 +^ C2 | {1,1,0} , {1,1,0} , | 
 +^ i | {0,0,0} , | 
 +^ S4 | {0,0,1} , {0,0,1} , | 
 +^ σh | {0,0,1} , | 
 +^ σv | {1,0,0} , {0,1,0} , | 
 +^ σd | {1,1,0} , {1,1,0} , | 
 + 
 +### 
 + 
 +===== Different Settings ===== 
 + 
 +### 
 + 
 +  * [[physics_chemistry:point_groups:d4h:orientation_zxy|Point Group D4h with orientation Zxy]] 
 + 
 +### 
 + 
 +===== Character Table ===== 
 + 
 +### 
 + 
 +  ^  E(1)  ^  C4(2)  ^  C2(1)  ^  C2(2)  ^  C2(2)  ^  i(1)  ^  S4(2)  ^  σh(1)  ^  σv(2)  ^  σd(2)  ^ 
 +^ A1g1111111111
 +^ A2g1111111111
 +^ B1g1111111111
 +^ B2g1111111111
 +^ Eg2020020200
 +^ A1u1111111111
 +^ A2u1111111111
 +^ B1u1111111111
 +^ B2u1111111111
 +^ Eu2020020200
 + 
 +### 
 + 
 +===== Product Table ===== 
 + 
 +### 
 + 
 +  ^  A1g  ^  A2g  ^  B1g  ^  B2g  ^  Eg  ^  A1u  ^  A2u  ^  B1u  ^  B2u  ^  Eu  ^ 
 +^ A1g  | A1g  | A2g  | B1g  | B2g  | Eg  | A1u  | A2u  | B1u  | B2u  | Eu  | 
 +^ A2g  | A2g  | A1g  | B2g  | B1g  | Eg  | A2u  | A1u  | B2u  | B1u  | Eu  | 
 +^ B1g  | B1g  | B2g  | A1g  | A2g  | Eg  | B1u  | B2u  | A1u  | A2u  | Eu  | 
 +^ B2g  | B2g  | B1g  | A2g  | A1g  | Eg  | B2u  | B1u  | A2u  | A1u  | Eu  | 
 +^ Eg  | Eg  | Eg  | Eg  | Eg  | A1g+A2g+B1g+B2g  | Eu  | Eu  | Eu  | Eu  | A1u+A2u+B1u+B2u  | 
 +^ A1u  | A1u  | A2u  | B1u  | B2u  | Eu  | A1g  | A2g  | B1g  | B2g  | Eg  | 
 +^ A2u  | A2u  | A1u  | B2u  | B1u  | Eu  | A2g  | A1g  | B2g  | B1g  | Eg  | 
 +^ B1u  | B1u  | B2u  | A1u  | A2u  | Eu  | B1g  | B2g  | A1g  | A2g  | Eg  | 
 +^ B2u  | B2u  | B1u  | A2u  | A1u  | Eu  | B2g  | B1g  | A2g  | A1g  | Eg  | 
 +^ Eu  | Eu  | Eu  | Eu  | Eu  | A1u+A2u+B1u+B2u  | Eg  | Eg  | Eg  | Eg  | A1g+A2g+B1g+B2g  | 
 + 
 +### 
 + 
 +===== Sub Groups with compatible settings ===== 
 + 
 +### 
 + 
 +  * [[physics_chemistry:point_groups:c1:orientation_1|Point Group C1 with orientation 1]] 
 +  * [[physics_chemistry:point_groups:c2h:orientation_z|Point Group C2h with orientation Z]] 
 +  * [[physics_chemistry:point_groups:c2v:orientation_zxy|Point Group C2v with orientation Zxy]] 
 +  * [[physics_chemistry:point_groups:c2:orientation_x|Point Group C2 with orientation X]] 
 +  * [[physics_chemistry:point_groups:c2:orientation_y|Point Group C2 with orientation Y]] 
 +  * [[physics_chemistry:point_groups:c2:orientation_z|Point Group C2 with orientation Z]] 
 +  * [[physics_chemistry:point_groups:c4h:orientation_z|Point Group C4h with orientation Z]] 
 +  * [[physics_chemistry:point_groups:c4v:orientation_zxy|Point Group C4v with orientation Zxy]] 
 +  * [[physics_chemistry:point_groups:c4:orientation_z|Point Group C4 with orientation Z]] 
 +  * [[physics_chemistry:point_groups:ci:orientation_|Point Group Ci with orientation ]] 
 +  * [[physics_chemistry:point_groups:cs:orientation_x|Point Group Cs with orientation X]] 
 +  * [[physics_chemistry:point_groups:cs:orientation_y|Point Group Cs with orientation Y]] 
 +  * [[physics_chemistry:point_groups:cs:orientation_z|Point Group Cs with orientation Z]] 
 +  * [[physics_chemistry:point_groups:d2d:orientation_zxy|Point Group D2d with orientation Zxy]] 
 +  * [[physics_chemistry:point_groups:d2h:orientation_xyz|Point Group D2h with orientation XYZ]] 
 +  * [[physics_chemistry:point_groups:d2:orientation_xyz|Point Group D2 with orientation XYZ]] 
 +  * [[physics_chemistry:point_groups:d4:orientation_zxy|Point Group D4 with orientation Zxy]] 
 +  * [[physics_chemistry:point_groups:s4:orientation_z|Point Group S4 with orientation Z]] 
 + 
 +### 
 + 
 +===== Super Groups with compatible settings ===== 
 + 
 +### 
 + 
 +  * [[physics_chemistry:point_groups:oh:orientation_xyz|Point Group Oh with orientation XYZ]] 
 + 
 +### 
 + 
 +===== Invariant Potential expanded on renormalized spherical Harmonics ===== 
 + 
 +### 
 + 
 +Any potential (function) can be written as a sum over spherical harmonics. 
 +V(r,θ,ϕ)=k=0km=kAk,m(r)C(m)k(θ,ϕ) 
 +Here Ak,m(r) is a radial function and C(m)k(θ,ϕ) a renormalised spherical harmonics. C(m)k(θ,ϕ)=4π2k+1Y(m)k(θ,ϕ) 
 +The presence of symmetry induces relations between the expansion coefficients such that V(r,θ,ϕ) is invariant under all symmetry operations. For the D4h Point group with orientation Zxy the form of the expansion coefficients is: 
 + 
 +### 
 + 
 +==== Expansion ==== 
 + 
 +### 
 + 
 + $$A_{k,m} = \begin{cases} 
 + A(0,0) & k=0\land m=0 \\ 
 + A(2,0) & k=2\land m=0 \\ 
 + A(4,4) & k=4\land (m=-4\lor m=4) \\ 
 + A(4,0) & k=4\land m=0 \\ 
 + A(6,4) & k=6\land (m=-4\lor m=4) \\ 
 + A(6,0) & k=6\land m=0 
 +\end{cases}$$ 
 + 
 +### 
 + 
 +==== Input format suitable for Mathematica (Quanty.nb) ==== 
 + 
 +### 
 + 
 +<code Quanty Akm_D4h_Zxy.Quanty.nb
 + 
 +Akm[k_,m_]:=Piecewise[{{A[0, 0], k == 0 && m == 0}, {A[2, 0], k == 2 && m == 0}, {A[4, 4], k == 4 && (m == -4 || m == 4)}, {A[4, 0], k == 4 && m == 0}, {A[6, 4], k == 6 && (m == -4 || m == 4)}, {A[6, 0], k == 6 && m == 0}}, 0] 
 </code> </code>
  
-==== Result ==== +###
-<WRAP center box 100%> +
-text produced as output +
-</WRAP>+
  
-===== Table of contents ===== +==== Input format suitable for Quanty ====
-{{indexmenu>.#1}}+
  
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{0, 0, A(0,0)} , 
 +       {2, 0, A(2,0)} , 
 +       {4, 0, A(4,0)} , 
 +       {4,-4, A(4,4)} , 
 +       {4, 4, A(4,4)} , 
 +       {6, 0, A(6,0)} , 
 +       {6,-4, A(6,4)} , 
 +       {6, 4, A(6,4)} }
 +
 +</code>
 +
 +###
 +
 +==== One particle coupling on a basis of spherical harmonics ====
 +
 +###
 +
 +The operator representing the potential in second quantisation is given as:
 +O=n,l,m,n,l,mψn,l,m(r,θ,ϕ)|V(r,θ,ϕ)|ψn,l,m(r,θ,ϕ)an,l,man,l,m
 +For the quantisation of the wave-function (physical meaning of the indices n,l,m) we can choose a basis of spherical harmonics times some radial function, i.e. ψn,l,m(r,θ,ϕ)=Rn,l(r)Y(l)m(θ,ϕ). With this choice the integral for the expectation value in front of the creation and annihilation operators separates into a radial part and angular part. The angular part has an analytical solution, the radial integral is cast int a parameter.
 +Anl,nl(k,m)=Rn,l|Ak,m(r)|Rn,l
 +Note the difference between the function Ak,m and the parameter Anl,nl(k,m)
 +
 +
 +###
 +
 +
 +
 +###
 +
 +
 +we can express the operator as 
 +O=n,l,m,n,l,m,k,mAnl,nl(k,m)Y(m)l(θ,ϕ)|C(m)k(θ,ϕ)|Y(m)l(θ,ϕ)an,l,man,l,m
 +
 +
 +###
 +
 +
 +
 +###
 +
 +
 +The table below shows the expectation value of O on a basis of spherical harmonics. We suppressed the principle quantum number indices. Note that in principle Al,l(k,m) can be complex. Instead of allowing complex parameters we took Al,l(k,m)+IBl,l(k,m) (with both A and B real) as the expansion parameter.
 +
 +###
 +
 +
 +
 +###
 +
 +  ^  Y(0)0  ^  Y(1)1  ^  Y(1)0  ^  Y(1)1  ^  Y(2)2  ^  Y(2)1  ^  Y(2)0  ^  Y(2)1  ^  Y(2)2  ^  Y(3)3  ^  Y(3)2  ^  Y(3)1  ^  Y(3)0  ^  Y(3)1  ^  Y(3)2  ^  Y(3)3  ^
 +^Y(0)0|Ass(0,0)|0|0|0|0|0|Asd(2,0)5|0|0|0|0|0|0|0|0|0|
 +^Y(1)1|0|App(0,0)15App(2,0)|0|0|0|0|0|0|0|0|0|3527Apf(2,0)1327Apf(4,0)|0|0|0|2Apf(4,4)33|
 +^Y(1)0|0|0|App(0,0)+25App(2,0)|0|0|0|0|0|0|0|0|0|3537Apf(2,0)+4Apf(4,0)321|0|0|0|
 +^Y(1)1|0|0|0|App(0,0)15App(2,0)|0|0|0|0|0|2Apf(4,4)33|0|0|0|3527Apf(2,0)1327Apf(4,0)|0|0|
 +^Y(2)2|0|0|0|0|Add(0,0)27Add(2,0)+121Add(4,0)|0|0|0|13107Add(4,4)|0|0|0|0|0|0|0|
 +^Y(2)1|0|0|0|0|0|Add(0,0)+17Add(2,0)421Add(4,0)|0|0|0|0|0|0|0|0|0|0|
 +^Y(2)0|Asd(2,0)5|0|0|0|0|0|Add(0,0)+27Add(2,0)+27Add(4,0)|0|0|0|0|0|0|0|0|0|
 +^Y(2)1|0|0|0|0|0|0|0|Add(0,0)+17Add(2,0)421Add(4,0)|0|0|0|0|0|0|0|0|
 +^Y(2)2|0|0|0|0|13107Add(4,4)|0|0|0|Add(0,0)27Add(2,0)+121Add(4,0)|0|0|0|0|0|0|0|
 +^Y(3)3|0|0|0|2Apf(4,4)33|0|0|0|0|0|Aff(0,0)13Aff(2,0)+111Aff(4,0)5429Aff(6,0)|0|0|0|111143Aff(4,4)5143703Aff(6,4)|0|0|
 +^Y(3)2|0|0|0|0|0|0|0|0|0|0|Aff(0,0)733Aff(4,0)+10143Aff(6,0)|0|0|0|13370Aff(4,4)+1014314Aff(6,4)|0|
 +^Y(3)1|0|3527Apf(2,0)1327Apf(4,0)|0|0|0|0|0|0|0|0|0|Aff(0,0)+15Aff(2,0)+133Aff(4,0)25143Aff(6,0)|0|0|0|111143Aff(4,4)5143703Aff(6,4)|
 +^Y(3)0|0|0|3537Apf(2,0)+4Apf(4,0)321|0|0|0|0|0|0|0|0|0|Aff(0,0)+415Aff(2,0)+211Aff(4,0)+100429Aff(6,0)|0|0|0|
 +^Y(3)1|0|0|0|3527Apf(2,0)1327Apf(4,0)|0|0|0|0|0|111143Aff(4,4)5143703Aff(6,4)|0|0|0|Aff(0,0)+15Aff(2,0)+133Aff(4,0)25143Aff(6,0)|0|0|
 +^Y(3)2|0|0|0|0|0|0|0|0|0|0|13370Aff(4,4)+1014314Aff(6,4)|0|0|0|Aff(0,0)733Aff(4,0)+10143Aff(6,0)|0|
 +^Y(3)3|0|2Apf(4,4)33|0|0|0|0|0|0|0|0|0|111143Aff(4,4)5143703Aff(6,4)|0|0|0|Aff(0,0)13Aff(2,0)+111Aff(4,0)5429Aff(6,0)|
 +
 +
 +###
 +
 +==== Rotation matrix to symmetry adapted functions (choice is not unique) ====
 +
 +###
 +
 +
 +Instead of a basis of spherical harmonics one can chose any other basis, which is given by a unitary transformation. Here we choose a rotation that simplifies the representation of the crystal field
 +
 +###
 +
 +
 +
 +###
 +
 +  ^  Y(0)0  ^  Y(1)1  ^  Y(1)0  ^  Y(1)1  ^  Y(2)2  ^  Y(2)1  ^  Y(2)0  ^  Y(2)1  ^  Y(2)2  ^  Y(3)3  ^  Y(3)2  ^  Y(3)1  ^  Y(3)0  ^  Y(3)1  ^  Y(3)2  ^  Y(3)3  ^
 +^s|1|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|
 +^px|0|12|0|12|0|0|0|0|0|0|0|0|0|0|0|0|
 +^py|0|i2|0|i2|0|0|0|0|0|0|0|0|0|0|0|0|
 +^pz|0|0|1|0|0|0|0|0|0|0|0|0|0|0|0|0|
 +^dx2y2|0|0|0|0|12|0|0|0|12|0|0|0|0|0|0|0|
 +^d3z2r2|0|0|0|0|0|0|1|0|0|0|0|0|0|0|0|0|
 +^dyz|0|0|0|0|0|i2|0|i2|0|0|0|0|0|0|0|0|
 +^dxz|0|0|0|0|0|12|0|12|0|0|0|0|0|0|0|0|
 +^dxy|0|0|0|0|i2|0|0|0|i2|0|0|0|0|0|0|0|
 +^fxyz|0|0|0|0|0|0|0|0|0|0|i2|0|0|0|i2|0|
 +^fx(5x2r2)|0|0|0|0|0|0|0|0|0|54|0|34|0|34|0|54|
 +^fy(5y2r2)|0|0|0|0|0|0|0|0|0|i54|0|i34|0|i34|0|i54|
 +^fz(5z2r2)|0|0|0|0|0|0|0|0|0|0|0|0|1|0|0|0|
 +^fx(y2z2)|0|0|0|0|0|0|0|0|0|34|0|54|0|54|0|34|
 +^fy(z2x2)|0|0|0|0|0|0|0|0|0|i34|0|i54|0|i54|0|i34|
 +^fz(x2y2)|0|0|0|0|0|0|0|0|0|0|12|0|0|0|12|0|
 +
 +
 +###
 +
 +==== One particle coupling on a basis of symmetry adapted functions ====
 +
 +###
 +
 +After rotation we find
 +
 +###
 +
 +
 +
 +###
 +
 +  ^  s  ^  px  ^  py  ^  pz  ^  dx2y2  ^  d3z2r2  ^  dyz  ^  dxz  ^  dxy  ^  fxyz  ^  fx(5x2r2)  ^  fy(5y2r2)  ^  fz(5z2r2)  ^  fx(y2z2)  ^  fy(z2x2)  ^  fz(x2y2)  ^
 +^s|Ass(0,0)|0|0|0|0|Asd(2,0)5|0|0|0|0|0|0|0|0|0|0|
 +^px|0|App(0,0)15App(2,0)|0|0|0|0|0|0|0|0|31037Apf(2,0)+Apf(4,0)221+1356Apf(4,4)|0|0|3Apf(2,0)235+1657Apf(4,0)Apf(4,4)32|0|0|
 +^py|0|0|App(0,0)15App(2,0)|0|0|0|0|0|0|0|0|31037Apf(2,0)+Apf(4,0)221+1356Apf(4,4)|0|0|3Apf(2,0)2351657Apf(4,0)+Apf(4,4)32|0|
 +^pz|0|0|0|App(0,0)+25App(2,0)|0|0|0|0|0|0|0|0|3537Apf(2,0)+4Apf(4,0)321|0|0|0|
 +^dx2y2|0|0|0|0|Add(0,0)27Add(2,0)+121Add(4,0)+13107Add(4,4)|0|0|0|0|0|0|0|0|0|0|0|
 +^d3z2r2|Asd(2,0)5|0|0|0|0|Add(0,0)+27Add(2,0)+27Add(4,0)|0|0|0|0|0|0|0|0|0|0|
 +^dyz|0|0|0|0|0|0|Add(0,0)+17Add(2,0)421Add(4,0)|0|0|0|0|0|0|0|0|0|
 +^dxz|0|0|0|0|0|0|0|Add(0,0)+17Add(2,0)421Add(4,0)|0|0|0|0|0|0|0|0|
 +^dxy|0|0|0|0|0|0|0|0|Add(0,0)27Add(2,0)+121Add(4,0)13107Add(4,4)|0|0|0|0|0|0|0|
 +^fxyz|0|0|0|0|0|0|0|0|0|Aff(0,0)733Aff(4,0)13370Aff(4,4)+10143Aff(6,0)1014314Aff(6,4)|0|0|0|0|0|0|
 +^fx(5x2r2)|0|31037Apf(2,0)+Apf(4,0)221+1356Apf(4,4)|0|0|0|0|0|0|0|0|Aff(0,0)215Aff(2,0)+344Aff(4,0)+122352Aff(4,4)125Aff(6,0)17162528672Aff(6,4)|0|0|Aff(2,0)1514453Aff(4,0)+12276Aff(4,4)3557253Aff(6,0)5286356Aff(6,4)|0|0|
 +^fy(5y2r2)|0|0|31037Apf(2,0)+Apf(4,0)221+1356Apf(4,4)|0|0|0|0|0|0|0|0|Aff(0,0)215Aff(2,0)+344Aff(4,0)+122352Aff(4,4)125Aff(6,0)17162528672Aff(6,4)|0|0|Aff(2,0)15+14453Aff(4,0)12276Aff(4,4)+3557253Aff(6,0)+5286356Aff(6,4)|0|
 +^fz(5z2r2)|0|0|0|3537Apf(2,0)+4Apf(4,0)321|0|0|0|0|0|0|0|0|Aff(0,0)+415Aff(2,0)+211Aff(4,0)+100429Aff(6,0)|0|0|0|
 +^fx(y2z2)|0|3Apf(2,0)235+1657Apf(4,0)Apf(4,4)32|0|0|0|0|0|0|0|0|Aff(2,0)1514453Aff(4,0)+12276Aff(4,4)3557253Aff(6,0)5286356Aff(6,4)|0|0|Aff(0,0)+7132Aff(4,0)122352Aff(4,4)544Aff(6,0)+2528672Aff(6,4)|0|0|
 +^fy(z2x2)|0|0|3Apf(2,0)2351657Apf(4,0)+Apf(4,4)32|0|0|0|0|0|0|0|0|Aff(2,0)15+14453Aff(4,0)12276Aff(4,4)+3557253Aff(6,0)+5286356Aff(6,4)|0|0|Aff(0,0)+7132Aff(4,0)122352Aff(4,4)544Aff(6,0)+2528672Aff(6,4)|0|
 +^fz(x2y2)|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|Aff(0,0)733Aff(4,0)+13370Aff(4,4)+10143Aff(6,0)+1014314Aff(6,4)|
 +
 +
 +###
 +
 +===== Coupling for a single shell =====
 +
 +
 +
 +###
 +
 +Although the parameters Al,l(k,m) uniquely define the potential, there is no simple relation between these paramters and the eigenstates of the potential. In this section we replace the parameters Al,l(k,m) by paramters that relate to the eigen energies of the potential acting on or between two shells with angular momentum l and l.
 +
 +###
 +
 +
 +
 +###
 +
 +Click on one of the subsections to expand it or <hiddenSwitch expand all> 
 +
 +###
 +
 +==== Potential for s orbitals ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + \text{Ea1g} & k=0\land m=0 \\
 + 0 & \text{True}
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{Ea1g, k == 0 && m == 0}}, 0]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{0, 0, Ea1g} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +  ^  Y(0)0  ^
 +^Y(0)0|Ea1g|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +  ^  s  ^
 +^s|Ea1g|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Rotation matrix used** >
 +
 +###
 +
 +  ^  Y(0)0  ^
 +^s|1|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Irriducible representations and their onsite energy** >
 +
 +###
 +
 +^ ^Ea1g | {{:physics_chemistry:pointgroup:d4h_zxy_orb_0_1.png?150}} |
 +|ψ(θ,ϕ)=11 |12π | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |12π | ::: |
 +
 +
 +###
 +
 +</hidden>
 +==== Potential for p orbitals ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + \frac{1}{3} (\text{Ea2u}+2 \text{Eeu}) & k=0\land m=0 \\
 + \frac{5 (\text{Ea2u}-\text{Eeu})}{3} & k=2\land m=0
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{(Ea2u + 2*Eeu)/3, k == 0 && m == 0}, {(5*(Ea2u - Eeu))/3, k == 2 && m == 0}}, 0]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{0, 0, (1/3)*(Ea2u + (2)*(Eeu))} , 
 +       {2, 0, (5/3)*(Ea2u + (-1)*(Eeu))} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +  ^  Y(1)1  ^  Y(1)0  ^  Y(1)1  ^
 +^Y(1)1|Eeu|0|0|
 +^Y(1)0|0|Ea2u|0|
 +^Y(1)1|0|0|Eeu|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +  ^  px  ^  py  ^  pz  ^
 +^px|Eeu|0|0|
 +^py|0|Eeu|0|
 +^pz|0|0|Ea2u|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Rotation matrix used** >
 +
 +###
 +
 +  ^  Y(1)1  ^  Y(1)0  ^  Y(1)1  ^
 +^px|12|0|12|
 +^py|i2|0|i2|
 +^pz|0|1|0|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Irriducible representations and their onsite energy** >
 +
 +###
 +
 +^ ^Eeu | {{:physics_chemistry:pointgroup:d4h_zxy_orb_1_1.png?150}} |
 +|ψ(θ,ϕ)=11 |123πsin(θ)cos(ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |123πx | ::: |
 +^ ^Eeu | {{:physics_chemistry:pointgroup:d4h_zxy_orb_1_2.png?150}} |
 +|ψ(θ,ϕ)=11 |123πsin(θ)sin(ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |123πy | ::: |
 +^ ^Ea2u | {{:physics_chemistry:pointgroup:d4h_zxy_orb_1_3.png?150}} |
 +|ψ(θ,ϕ)=11 |123πcos(θ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |123πz | ::: |
 +
 +
 +###
 +
 +</hidden>
 +==== Potential for d orbitals ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + \frac{1}{5} (\text{Ea1g}+\text{Eb1g}+\text{Eb2g}+2 \text{Eeg}) & k=0\land m=0 \\
 + \text{Ea1g}-\text{Eb1g}-\text{Eb2g}+\text{Eeg} & k=2\land m=0 \\
 + \frac{3}{2} \sqrt{\frac{7}{10}} (\text{Eb1g}-\text{Eb2g}) & k=4\land (m=-4\lor m=4) \\
 + \frac{3}{10} (6 \text{Ea1g}+\text{Eb1g}+\text{Eb2g}-8 \text{Eeg}) & k=4\land m=0
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{(Ea1g + Eb1g + Eb2g + 2*Eeg)/5, k == 0 && m == 0}, {Ea1g - Eb1g - Eb2g + Eeg, k == 2 && m == 0}, {(3*Sqrt[7/10]*(Eb1g - Eb2g))/2, k == 4 && (m == -4 || m == 4)}, {(3*(6*Ea1g + Eb1g + Eb2g - 8*Eeg))/10, k == 4 && m == 0}}, 0]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{0, 0, (1/5)*(Ea1g + Eb1g + Eb2g + (2)*(Eeg))} , 
 +       {2, 0, Ea1g + (-1)*(Eb1g) + (-1)*(Eb2g) + Eeg} , 
 +       {4, 0, (3/10)*((6)*(Ea1g) + Eb1g + Eb2g + (-8)*(Eeg))} , 
 +       {4,-4, (3/2)*((sqrt(7/10))*(Eb1g + (-1)*(Eb2g)))} , 
 +       {4, 4, (3/2)*((sqrt(7/10))*(Eb1g + (-1)*(Eb2g)))} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +  ^  Y(2)2  ^  Y(2)1  ^  Y(2)0  ^  Y(2)1  ^  Y(2)2  ^
 +^Y(2)2|Eb1g+Eb2g2|0|0|0|Eb1gEb2g2|
 +^Y(2)1|0|Eeg|0|0|0|
 +^Y(2)0|0|0|Ea1g|0|0|
 +^Y(2)1|0|0|0|Eeg|0|
 +^Y(2)2|Eb1gEb2g2|0|0|0|Eb1g+Eb2g2|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +  ^  dx2y2  ^  d3z2r2  ^  dyz  ^  dxz  ^  dxy  ^
 +^dx2y2|Eb1g|0|0|0|0|
 +^d3z2r2|0|Ea1g|0|0|0|
 +^dyz|0|0|Eeg|0|0|
 +^dxz|0|0|0|Eeg|0|
 +^dxy|0|0|0|0|Eb2g|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Rotation matrix used** >
 +
 +###
 +
 +  ^  Y(2)2  ^  Y(2)1  ^  Y(2)0  ^  Y(2)1  ^  Y(2)2  ^
 +^dx2y2|12|0|0|0|12|
 +^d3z2r2|0|0|1|0|0|
 +^dyz|0|i2|0|i2|0|
 +^dxz|0|12|0|12|0|
 +^dxy|i2|0|0|0|i2|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Irriducible representations and their onsite energy** >
 +
 +###
 +
 +^ ^Eb1g | {{:physics_chemistry:pointgroup:d4h_zxy_orb_2_1.png?150}} |
 +|ψ(θ,ϕ)=11 |1415πsin2(θ)cos(2ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |1415π(x2y2) | ::: |
 +^ ^Ea1g | {{:physics_chemistry:pointgroup:d4h_zxy_orb_2_2.png?150}} |
 +|ψ(θ,ϕ)=11 |185π(3cos(2θ)+1) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |145π(3z21) | ::: |
 +^ ^Eeg | {{:physics_chemistry:pointgroup:d4h_zxy_orb_2_3.png?150}} |
 +|ψ(θ,ϕ)=11 |1415πsin(2θ)sin(ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |1215πyz | ::: |
 +^ ^Eeg | {{:physics_chemistry:pointgroup:d4h_zxy_orb_2_4.png?150}} |
 +|ψ(θ,ϕ)=11 |1415πsin(2θ)cos(ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |1215πxz | ::: |
 +^ ^Eb2g | {{:physics_chemistry:pointgroup:d4h_zxy_orb_2_5.png?150}} |
 +|ψ(θ,ϕ)=11 |1415πsin2(θ)sin(2ϕ) | ::: |
 +|ψ(ˆx,ˆy,ˆz)=11 |1215πxy | ::: |
 +
 +
 +###
 +
 +</hidden>
 +==== Potential for f orbitals ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + \frac{1}{7} (\text{Ea2u}+\text{Eb1u}+\text{Eb2u}+2 \text{Eeu1}+2 \text{Eeu2}) & k=0\land m=0 \\
 + \frac{5}{7} \left(\text{Ea2u}-\text{Eeu1}+\sqrt{15} \text{Meu}\right) & k=2\land m=0 \\
 + -\frac{3}{56} \left(2 \sqrt{70} \text{Eb1u}-2 \sqrt{70} \text{Eb2u}-3 \sqrt{70} \text{Eeu1}+3 \sqrt{70} \text{Eeu2}-2 \sqrt{42} \text{Meu}\right) & k=4\land (m=-4\lor m=4) \\
 + \frac{3}{28} \left(12 \text{Ea2u}-14 \text{Eb1u}-14 \text{Eb2u}+9 \text{Eeu1}+7 \text{Eeu2}-2 \sqrt{15} \text{Meu}\right) & k=4\land m=0 \\
 + -\frac{13}{560} \sqrt{3} \left(4 \sqrt{42} \text{Eb1u}-4 \sqrt{42} \text{Eb2u}+5 \sqrt{42} \text{Eeu1}-5 \sqrt{42} \text{Eeu2}+2 \sqrt{70} \text{Meu}\right) & k=6\land (m=-4\lor m=4) \\
 + \frac{13}{280} \left(40 \text{Ea2u}+12 \text{Eb1u}+12 \text{Eb2u}-25 \text{Eeu1}-39 \text{Eeu2}-14 \sqrt{15} \text{Meu}\right) & k=6\land m=0
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{(Ea2u + Eb1u + Eb2u + 2*Eeu1 + 2*Eeu2)/7, k == 0 && m == 0}, {(5*(Ea2u - Eeu1 + Sqrt[15]*Meu))/7, k == 2 && m == 0}, {(-3*(2*Sqrt[70]*Eb1u - 2*Sqrt[70]*Eb2u - 3*Sqrt[70]*Eeu1 + 3*Sqrt[70]*Eeu2 - 2*Sqrt[42]*Meu))/56, k == 4 && (m == -4 || m == 4)}, {(3*(12*Ea2u - 14*Eb1u - 14*Eb2u + 9*Eeu1 + 7*Eeu2 - 2*Sqrt[15]*Meu))/28, k == 4 && m == 0}, {(-13*Sqrt[3]*(4*Sqrt[42]*Eb1u - 4*Sqrt[42]*Eb2u + 5*Sqrt[42]*Eeu1 - 5*Sqrt[42]*Eeu2 + 2*Sqrt[70]*Meu))/560, k == 6 && (m == -4 || m == 4)}, {(13*(40*Ea2u + 12*Eb1u + 12*Eb2u - 25*Eeu1 - 39*Eeu2 - 14*Sqrt[15]*Meu))/280, k == 6 && m == 0}}, 0]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{0, 0, (1/7)*(Ea2u + Eb1u + Eb2u + (2)*(Eeu1) + (2)*(Eeu2))} , 
 +       {2, 0, (5/7)*(Ea2u + (-1)*(Eeu1) + (sqrt(15))*(Meu))} , 
 +       {4, 0, (3/28)*((12)*(Ea2u) + (-14)*(Eb1u) + (-14)*(Eb2u) + (9)*(Eeu1) + (7)*(Eeu2) + (-2)*((sqrt(15))*(Meu)))} , 
 +       {4,-4, (-3/56)*((2)*((sqrt(70))*(Eb1u)) + (-2)*((sqrt(70))*(Eb2u)) + (-3)*((sqrt(70))*(Eeu1)) + (3)*((sqrt(70))*(Eeu2)) + (-2)*((sqrt(42))*(Meu)))} , 
 +       {4, 4, (-3/56)*((2)*((sqrt(70))*(Eb1u)) + (-2)*((sqrt(70))*(Eb2u)) + (-3)*((sqrt(70))*(Eeu1)) + (3)*((sqrt(70))*(Eeu2)) + (-2)*((sqrt(42))*(Meu)))} , 
 +       {6, 0, (13/280)*((40)*(Ea2u) + (12)*(Eb1u) + (12)*(Eb2u) + (-25)*(Eeu1) + (-39)*(Eeu2) + (-14)*((sqrt(15))*(Meu)))} , 
 +       {6,-4, (-13/560)*((sqrt(3))*((4)*((sqrt(42))*(Eb1u)) + (-4)*((sqrt(42))*(Eb2u)) + (5)*((sqrt(42))*(Eeu1)) + (-5)*((sqrt(42))*(Eeu2)) + (2)*((sqrt(70))*(Meu))))} , 
 +       {6, 4, (-13/560)*((sqrt(3))*((4)*((sqrt(42))*(Eb1u)) + (-4)*((sqrt(42))*(Eb2u)) + (5)*((sqrt(42))*(Eeu1)) + (-5)*((sqrt(42))*(Eeu2)) + (2)*((sqrt(70))*(Meu))))} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +  ^  Y(3)3  ^  Y(3)2  ^  Y(3)1  ^  Y(3)0  ^  Y(3)1  ^  Y(3)2  ^  Y(3)3  ^
 +^Y(3)3|18(5Eeu1+3Eeu2215Meu)|0|0|0|18(15Eeu115Eeu2+2Meu)|0|0|
 +^Y(3)2|0|Eb1u+Eb2u2|0|0|0|Eb2uEb1u2|0|
 +^Y(3)1|0|0|18(3Eeu1+5Eeu2+215Meu)|0|0|0|18(15Eeu115Eeu2+2Meu)|
 +^Y(3)0|0|0|0|Ea2u|0|0|0|
 +^Y(3)1|18(15Eeu115Eeu2+2Meu)|0|0|0|18(3Eeu1+5Eeu2+215Meu)|0|0|
 +^Y(3)2|0|Eb2uEb1u2|0|0|0|Eb1u+Eb2u2|0|
 +^Y(3)3|0|0|18(15Eeu115Eeu2+2Meu)|0|0|0|18(5Eeu1+3Eeu2215Meu)|
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +  ^  fxyz  ^  fx(5x2r2)  ^  fy(5y2r2)  ^  fz(5z2r2)  ^  fx(y2z2)  ^  fy(z2x2)  ^  fz(x2y2)  ^
 +^fxyz|Eb1u|0|0|0|0|0|0|
 +^fx(5x2r2)|0|Eeu1|0|0|Meu|0|0|
 +^fy(5y2r2)| 0 | 0 | \text{Eeu1} | 0 | 0 | -\text{Meu} | 0 |
 +^ f_{z\left(5z^2-r^2\right)} | 0 | 0 | 0 | \text{Ea2u} | 0 | 0 | 0 |
 +^ f_{x\left(y^2-z^2\right)} | 0 | \text{Meu} | 0 | 0 | \text{Eeu2} | 0 | 0 |
 +^ f_{y\left(z^2-x^2\right)} | 0 | 0 | -\text{Meu} | 0 | 0 | \text{Eeu2} | 0 |
 +^ f_{z\left(x^2-y^2\right)} | 0 | 0 | 0 | 0 | 0 | 0 | \text{Eb2u} |
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Rotation matrix used** >
 +
 +###
 +
 +    ^  {Y_{-3}^{(3)}}   ^  {Y_{-2}^{(3)}}   ^  {Y_{-1}^{(3)}}   ^  {Y_{0}^{(3)}}   ^  {Y_{1}^{(3)}}   ^  {Y_{2}^{(3)}}   ^  {Y_{3}^{(3)}}   ^
 +^ f_{\text{xyz}} | 0 | \frac{i}{\sqrt{2}} | 0 | 0 | 0 | -\frac{i}{\sqrt{2}} | 0 |
 +^ f_{x\left(5x^2-r^2\right)} | \frac{\sqrt{5}}{4} | 0 | -\frac{\sqrt{3}}{4} | 0 | \frac{\sqrt{3}}{4} | 0 | -\frac{\sqrt{5}}{4} |
 +^ f_{y\left(5y^2-r^2\right)} | -\frac{i \sqrt{5}}{4} | 0 | -\frac{i \sqrt{3}}{4} | 0 | -\frac{i \sqrt{3}}{4} | 0 | -\frac{i \sqrt{5}}{4} |
 +^ f_{z\left(5z^2-r^2\right)} | 0 | 0 | 0 | 1 | 0 | 0 | 0 |
 +^ f_{x\left(y^2-z^2\right)} | -\frac{\sqrt{3}}{4} | 0 | -\frac{\sqrt{5}}{4} | 0 | \frac{\sqrt{5}}{4} | 0 | \frac{\sqrt{3}}{4} |
 +^ f_{y\left(z^2-x^2\right)} | -\frac{i \sqrt{3}}{4} | 0 | \frac{i \sqrt{5}}{4} | 0 | \frac{i \sqrt{5}}{4} | 0 | -\frac{i \sqrt{3}}{4} |
 +^ f_{z\left(x^2-y^2\right)} | 0 | \frac{1}{\sqrt{2}} | 0 | 0 | 0 | \frac{1}{\sqrt{2}} | 0 |
 +
 +
 +###
 +
 +</hidden>
 +<hidden **Irriducible representations and their onsite energy** >
 +
 +###
 +
 +^ ^\text{Eb1u} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_1.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{4} \sqrt{\frac{105}{\pi }} \sin ^2(\theta ) \cos (\theta ) \sin (2 \phi ) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{2} \sqrt{\frac{105}{\pi }} x y z | ::: |
 +^ ^\text{Eeu1} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_2.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{16} \sqrt{\frac{7}{\pi }} \sin (\theta ) \cos (\phi ) \left(10 \sin ^2(\theta ) \cos (2 \phi )-5 \cos (2 \theta )-7\right) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{16} \sqrt{\frac{7}{\pi }} x \left(5 x^2-15 y^2-15 z^2+3\right) | ::: |
 +^ ^\text{Eeu1} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_3.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |-\frac{1}{16} \sqrt{\frac{7}{\pi }} \sin (\theta ) \sin (\phi ) \left(10 \sin ^2(\theta ) \cos (2 \phi )+5 \cos (2 \theta )+7\right) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{16} \sqrt{\frac{7}{\pi }} y \left(-15 x^2+5 y^2-15 z^2+3\right) | ::: |
 +^ ^\text{Ea2u} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_4.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{16} \sqrt{\frac{7}{\pi }} (3 \cos (\theta )+5 \cos (3 \theta )) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{4} \sqrt{\frac{7}{\pi }} z \left(5 z^2-3\right) | ::: |
 +^ ^\text{Eeu2} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_5.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |-\frac{1}{16} \sqrt{\frac{105}{\pi }} \sin (\theta ) \cos (\phi ) \left(2 \sin ^2(\theta ) \cos (2 \phi )+3 \cos (2 \theta )+1\right) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |-\frac{1}{16} \sqrt{\frac{105}{\pi }} x \left(x^2-3 y^2+5 z^2-1\right) | ::: |
 +^ ^\text{Eeu2} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_6.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{32} \sqrt{\frac{105}{\pi }} \sin (\theta ) \sin (\phi ) \left(-4 \sin ^2(\theta ) \cos (2 \phi )+6 \cos (2 \theta )+2\right) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{16} \sqrt{\frac{105}{\pi }} y \left(-3 x^2+y^2+5 z^2-1\right) | ::: |
 +^ ^\text{Eb2u} | {{:physics_chemistry:pointgroup:d4h_zxy_orb_3_7.png?150}} |
 +|\psi(\theta,\phi)=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{4} \sqrt{\frac{105}{\pi }} \sin ^2(\theta ) \cos (\theta ) \cos (2 \phi ) | ::: |
 +|\psi(\hat{x},\hat{y},\hat{z})=\phantom{\sqrt{\frac{1}{1}}} |\frac{1}{4} \sqrt{\frac{105}{\pi }} z \left(x^2-y^2\right) | ::: |
 +
 +
 +###
 +
 +</hidden>
 +===== Coupling between two shells =====
 +
 +
 +
 +###
 +
 +Click on one of the subsections to expand it or <hiddenSwitch expand all> 
 +
 +###
 +
 +==== Potential for s-d orbital mixing ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + 0 & k\neq 2\lor m\neq 0 \\
 + A(2,0) & \text{True}
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{0, k != 2 || m != 0}}, A[2, 0]]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{2, 0, A(2,0)} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +    ^  {Y_{-2}^{(2)}}   ^  {Y_{-1}^{(2)}}   ^  {Y_{0}^{(2)}}   ^  {Y_{1}^{(2)}}   ^  {Y_{2}^{(2)}}   ^
 +^ {Y_{0}^{(0)}} | 0 | 0 | \frac{A(2,0)}{\sqrt{5}} | 0 | 0 |
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +    ^  d_{x^2-y^2}   ^  d_{3z^2-r^2}   ^  d_{\text{yz}}   ^  d_{\text{xz}}   ^  d_{\text{xy}}   ^
 +^ \text{s} | 0 | \frac{A(2,0)}{\sqrt{5}} | 0 | 0 | 0 |
 +
 +
 +###
 +
 +</hidden>
 +==== Potential for p-f orbital mixing ====
 +
 +<hidden **Potential parameterized with onsite energies of irriducible representations** >
 +
 +###
 +
 + $$A_{k,m} = \begin{cases}
 + 0 & k=0\land m=0 \\
 + A(2,0) & k=2\land m=0 \\
 + A(4,4) & k=4\land (m=-4\lor m=4) \\
 + A(4,0) & k=4\land m=0
 +\end{cases}$$
 +
 +###
 +
 +</hidden>
 +<hidden **Input format suitable for Mathematica (Quanty.nb)** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty.nb>
 +
 +Akm[k_,m_]:=Piecewise[{{0, k == 0 && m == 0}, {A[2, 0], k == 2 && m == 0}, {A[4, 4], k == 4 && (m == -4 || m == 4)}, {A[4, 0], k == 4 && m == 0}}, 0]
 +
 +</code>
 +
 +###
 +
 +</hidden><hidden **Input format suitable for Quanty** >
 +
 +###
 +
 +<code Quanty Akm_D4h_Zxy.Quanty>
 +
 +Akm = {{2, 0, A(2,0)} , 
 +       {4, 0, A(4,0)} , 
 +       {4,-4, A(4,4)} , 
 +       {4, 4, A(4,4)} }
 +
 +</code>
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of spherical Harmonics** >
 +
 +###
 +
 +    ^  {Y_{-3}^{(3)}}   ^  {Y_{-2}^{(3)}}   ^  {Y_{-1}^{(3)}}   ^  {Y_{0}^{(3)}}   ^  {Y_{1}^{(3)}}   ^  {Y_{2}^{(3)}}   ^  {Y_{3}^{(3)}}   ^
 +^ {Y_{-1}^{(1)}} | 0 | 0 | \frac{1}{15} \sqrt{\frac{2}{7}} (9 A(2,0)-5 A(4,0)) | 0 | 0 | 0 | -\frac{2 A(4,4)}{3 \sqrt{3}} |
 +^ {Y_{0}^{(1)}} | 0 | 0 | 0 | \frac{27 A(2,0)+20 A(4,0)}{15 \sqrt{21}} | 0 | 0 | 0 |
 +^ {Y_{1}^{(1)}} | -\frac{2 A(4,4)}{3 \sqrt{3}} | 0 | 0 | 0 | \frac{1}{15} \sqrt{\frac{2}{7}} (9 A(2,0)-5 A(4,0)) | 0 | 0 |
 +
 +
 +###
 +
 +</hidden>
 +<hidden **The Hamiltonian on a basis of symmetric functions** >
 +
 +###
 +
 +    ^  f_{\text{xyz}}   ^  f_{x\left(5x^2-r^2\right)}   ^  f_{y\left(5y^2-r^2\right)}   ^  f_{z\left(5z^2-r^2\right)}   ^  f_{x\left(y^2-z^2\right)}   ^  f_{y\left(z^2-x^2\right)}   ^  f_{z\left(x^2-y^2\right)}   ^
 +^ p_x | 0 | \frac{5 A(4,0)-9 A(2,0)}{10 \sqrt{21}}+\frac{1}{3} \sqrt{\frac{5}{6}} A(4,4) | 0 | 0 | \frac{5 A(4,0)-9 A(2,0)}{6 \sqrt{35}}-\frac{A(4,4)}{3 \sqrt{2}} | 0 | 0 |
 +^ p_y | 0 | 0 | \frac{5 A(4,0)-9 A(2,0)}{10 \sqrt{21}}+\frac{1}{3} \sqrt{\frac{5}{6}} A(4,4) | 0 | 0 | \frac{9 A(2,0)-5 A(4,0)}{6 \sqrt{35}}+\frac{A(4,4)}{3 \sqrt{2}} | 0 |
 +^ p_z | 0 | 0 | 0 | \frac{27 A(2,0)+20 A(4,0)}{15 \sqrt{21}} | 0 | 0 | 0 |
 +
 +
 +###
 +
 +</hidden>
 +
 +===== Table of several point groups =====
 +
 +###
 +
 +[[physics_chemistry:point_groups|Return to Main page on Point Groups]]
 +
 +###
 +
 +###
 +
 +^Nonaxial groups      | [[physics_chemistry:point_groups:c1|C]]<sub>[[physics_chemistry:point_groups:c1|1]]</sub> | [[physics_chemistry:point_groups:cs|C]]<sub>[[physics_chemistry:point_groups:cs|s]]</sub> | [[physics_chemistry:point_groups:ci|C]]<sub>[[physics_chemistry:point_groups:ci|i]]</sub> | | | | |
 +^C<sub>n</sub> groups | [[physics_chemistry:point_groups:c2|C]]<sub>[[physics_chemistry:point_groups:c2|2]]</sub> | [[physics_chemistry:point_groups:c3|C]]<sub>[[physics_chemistry:point_groups:c3|3]]</sub> | [[physics_chemistry:point_groups:c4|C]]<sub>[[physics_chemistry:point_groups:c4|4]]</sub> | [[physics_chemistry:point_groups:c5|C]]<sub>[[physics_chemistry:point_groups:c5|5]]</sub> | [[physics_chemistry:point_groups:c6|C]]<sub>[[physics_chemistry:point_groups:c6|6]]</sub> | [[physics_chemistry:point_groups:c7|C]]<sub>[[physics_chemistry:point_groups:c7|7]]</sub> | [[physics_chemistry:point_groups:c8|C]]<sub>[[physics_chemistry:point_groups:c8|8]]</sub>
 +^D<sub>n</sub> groups | [[physics_chemistry:point_groups:d2|D]]<sub>[[physics_chemistry:point_groups:d2|2]]</sub> | [[physics_chemistry:point_groups:d3|D]]<sub>[[physics_chemistry:point_groups:d3|3]]</sub> | [[physics_chemistry:point_groups:d4|D]]<sub>[[physics_chemistry:point_groups:d4|4]]</sub> | [[physics_chemistry:point_groups:d5|D]]<sub>[[physics_chemistry:point_groups:d5|5]]</sub> | [[physics_chemistry:point_groups:d6|D]]<sub>[[physics_chemistry:point_groups:d6|6]]</sub> | [[physics_chemistry:point_groups:d7|D]]<sub>[[physics_chemistry:point_groups:d7|7]]</sub> | [[physics_chemistry:point_groups:d8|D]]<sub>[[physics_chemistry:point_groups:d8|8]]</sub>
 +^C<sub>nv</sub> groups | [[physics_chemistry:point_groups:c2v|C]]<sub>[[physics_chemistry:point_groups:c2v|2v]]</sub> | [[physics_chemistry:point_groups:c3v|C]]<sub>[[physics_chemistry:point_groups:c3v|3v]]</sub> | [[physics_chemistry:point_groups:c4v|C]]<sub>[[physics_chemistry:point_groups:c4v|4v]]</sub> | [[physics_chemistry:point_groups:c5v|C]]<sub>[[physics_chemistry:point_groups:c5v|5v]]</sub> | [[physics_chemistry:point_groups:c6v|C]]<sub>[[physics_chemistry:point_groups:c6v|6v]]</sub> | [[physics_chemistry:point_groups:c7v|C]]<sub>[[physics_chemistry:point_groups:c7v|7v]]</sub> | [[physics_chemistry:point_groups:c8v|C]]<sub>[[physics_chemistry:point_groups:c8v|8v]]</sub>
 +^C<sub>nh</sub> groups | [[physics_chemistry:point_groups:c2h|C]]<sub>[[physics_chemistry:point_groups:c2h|2h]]</sub> | [[physics_chemistry:point_groups:c3h|C]]<sub>[[physics_chemistry:point_groups:c3h|3h]]</sub> | [[physics_chemistry:point_groups:c4h|C]]<sub>[[physics_chemistry:point_groups:c4h|4h]]</sub> | [[physics_chemistry:point_groups:c5h|C]]<sub>[[physics_chemistry:point_groups:c5h|5h]]</sub> | [[physics_chemistry:point_groups:c6h|C]]<sub>[[physics_chemistry:point_groups:c6h|6h]]</sub> | | | 
 +^D<sub>nh</sub> groups | [[physics_chemistry:point_groups:d2h|D]]<sub>[[physics_chemistry:point_groups:d2h|2h]]</sub> | [[physics_chemistry:point_groups:d3h|D]]<sub>[[physics_chemistry:point_groups:d3h|3h]]</sub> | [[physics_chemistry:point_groups:d4h|D]]<sub>[[physics_chemistry:point_groups:d4h|4h]]</sub> | [[physics_chemistry:point_groups:d5h|D]]<sub>[[physics_chemistry:point_groups:d5h|5h]]</sub> | [[physics_chemistry:point_groups:d6h|D]]<sub>[[physics_chemistry:point_groups:d6h|6h]]</sub> | [[physics_chemistry:point_groups:d7h|D]]<sub>[[physics_chemistry:point_groups:d7h|7h]]</sub> | [[physics_chemistry:point_groups:d8h|D]]<sub>[[physics_chemistry:point_groups:d8h|8h]]</sub>
 +^D<sub>nd</sub> groups | [[physics_chemistry:point_groups:d2d|D]]<sub>[[physics_chemistry:point_groups:d2d|2d]]</sub> | [[physics_chemistry:point_groups:d3d|D]]<sub>[[physics_chemistry:point_groups:d3d|3d]]</sub> | [[physics_chemistry:point_groups:d4d|D]]<sub>[[physics_chemistry:point_groups:d4d|4d]]</sub> | [[physics_chemistry:point_groups:d5d|D]]<sub>[[physics_chemistry:point_groups:d5d|5d]]</sub> | [[physics_chemistry:point_groups:d6d|D]]<sub>[[physics_chemistry:point_groups:d6d|6d]]</sub> | [[physics_chemistry:point_groups:d7d|D]]<sub>[[physics_chemistry:point_groups:d7d|7d]]</sub> | [[physics_chemistry:point_groups:d8d|D]]<sub>[[physics_chemistry:point_groups:d8d|8d]]</sub>
 +^S<sub>n</sub> groups | [[physics_chemistry:point_groups:S2|S]]<sub>[[physics_chemistry:point_groups:S2|2]]</sub> | [[physics_chemistry:point_groups:S4|S]]<sub>[[physics_chemistry:point_groups:S4|4]]</sub> | [[physics_chemistry:point_groups:S6|S]]<sub>[[physics_chemistry:point_groups:S6|6]]</sub> | [[physics_chemistry:point_groups:S8|S]]<sub>[[physics_chemistry:point_groups:S8|8]]</sub> | [[physics_chemistry:point_groups:S10|S]]<sub>[[physics_chemistry:point_groups:S10|10]]</sub> | [[physics_chemistry:point_groups:S12|S]]<sub>[[physics_chemistry:point_groups:S12|12]]</sub> |  | 
 +^Cubic groups | [[physics_chemistry:point_groups:T|T]] | [[physics_chemistry:point_groups:Th|T]]<sub>[[physics_chemistry:point_groups:Th|h]]</sub> | [[physics_chemistry:point_groups:Td|T]]<sub>[[physics_chemistry:point_groups:Td|d]]</sub> | [[physics_chemistry:point_groups:O|O]] | [[physics_chemistry:point_groups:Oh|O]]<sub>[[physics_chemistry:point_groups:Oh|h]]</sub> | [[physics_chemistry:point_groups:I|I]] | [[physics_chemistry:point_groups:Ih|I]]<sub>[[physics_chemistry:point_groups:Ih|h]]</sub>
 +^Linear groups      | [[physics_chemistry:point_groups:cinfv|C]]<sub>[[physics_chemistry:point_groups:cinfv|\inftyv]]</sub> | [[physics_chemistry:point_groups:cinfv|D]]<sub>[[physics_chemistry:point_groups:dinfh|\inftyh]]</sub> | | | | | |
 +
 +###
Print/export