@@ -708,9 +708,9 @@ def cf_trigonal_d(case, Delta1, Delta2):
708
708
- '100': :math:`z` axis vertical
709
709
710
710
Delta1: float scalar
711
- The :math:`e^\sigma_g` orbitals are Delta1 above the :math:`a^_{1g}`
711
+ The :math:`e^\\ sigma_g` orbitals are Delta1 above the :math:`a^_{1g}`
712
712
Delta1: float scalar
713
- The :math:`e^\pi_g` orbitals are Delta2 below the :math:`a^_{1g}`
713
+ The :math:`e^\\ pi_g` orbitals are Delta2 below the :math:`a^_{1g}`
714
714
715
715
Returns
716
716
-------
@@ -719,23 +719,23 @@ def cf_trigonal_d(case, Delta1, Delta2):
719
719
"""
720
720
721
721
cf_trig = np .zeros ((5 , 5 ), dtype = complex )
722
- cf_trig [0 , 0 ] = 0 # a1g
723
- cf_trig [1 , 1 ] = - Delta2 # epi_g
724
- cf_trig [2 , 2 ] = - Delta2 # epi_g
725
- cf_trig [3 , 3 ] = Delta1 # esigma_g
726
- cf_trig [4 , 4 ] = Delta1 # esigma_g
722
+ cf_trig [0 , 0 ] = 0 # a1g
723
+ cf_trig [1 , 1 ] = - Delta2 # epi_g
724
+ cf_trig [2 , 2 ] = - Delta2 # epi_g
725
+ cf_trig [3 , 3 ] = Delta1 # esigma_g
726
+ cf_trig [4 , 4 ] = Delta1 # esigma_g
727
727
728
728
cf_trig_spin = np .zeros ((10 , 10 ), dtype = complex )
729
729
cf_trig_spin [::2 , ::2 ] = cf_trig
730
730
cf_trig_spin [1 ::2 , 1 ::2 ] = cf_trig
731
731
732
732
if case == '111' :
733
733
cf_real = cb_op (cf_trig_spin , tmat_trig2r ('111' , ispin = True ))
734
- cf = cb_op (cf_real , tmat_r2c ('d' , ispin = True ))
734
+ cf = cb_op (cf_real , tmat_r2c ('d' , ispin = True ))
735
735
return cf
736
736
elif case == '100' :
737
737
cf_real = cb_op (cf_trig_spin , tmat_trig2r ('100' , ispin = True ))
738
- cf = cb_op (cf_real , tmat_r2c ('d' , ispin = True ))
738
+ cf = cb_op (cf_real , tmat_r2c ('d' , ispin = True ))
739
739
return cf
740
740
else :
741
741
raise Exception (f"case { case } not implemented."
0 commit comments