计算碳原子系综的能级C

C原子核外6个电子的排布是

轨道的径向部分是

String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )";

String rj2 ="(z/a0)**(1.5)*2*sympy.exp(-z*r2/a0 )";

String rj3="(z/( 2*a1))**(1.5)*(2-z*r1/a1)*sympy.exp(-z*r1/(2*a1) )";

String rj4="(z/( 2*a1))**(1.5)*(2-z*r2/a1)*sympy.exp(-z*r2/(2*a1) )";

String rj5="(z/( 2*a1))**(1.5)*(z*r1/(3**0.5*a1) )*sympy.exp(-z*r1/(2*a1) )";

String rj6="(z/( 2*a1))**(1.5)*(z*r2/(3**0.5*a1) )*sympy.exp(-z*r2/(2*a1) )";

能级可以展开为

因为

所以

其中E1=-18

E3=-1.96，E5=-1.96。因此动能和势能的和为-43.875Hartrees

设两个P电子的量子数都是（2，1，0）

计算库仑能用公式

计算的J12=3.75，J13=0.36，J15=0.37，J34=0.23，J35=0.24，J56=0.29。库仑排斥能的和为8.19Hartrees

计算交换能用公式

计算的K13=0.06840270133467292，

K15=0.000749325366733499，

K35=0.0439453125，

交换能的和为0.452Hartrees

因此C的能级为

-43.875+8.188502153666878-0.4523893568056257=-36.13 Hartrees

计算程序用前面计算硼原子能级得到的方法，只需将子程序eB（）改为

public static void eC() throws  IOException, ParseException, InterruptedException {int z=6;int a0=1;int a1=4;int cn5=0;int cn6=0;String fx1="(z/a0)**(1.5)*2*sympy.exp(-z*r/a0  )*(4*pi)**(-0.5)";String fx2="(z/( 2*a1))**(1.5)*(2-z*r/a1)*sympy.exp(-z*r/(2*a1)  )*(4*pi)**(-0.5)";String fx3="(z/( 2*a1))**(1.5)*(z*r/(3**0.5*a1) )*sympy.exp(-z*r/(2*a1)  )*(3/(4*pi))**(0.5)* cos(θ)";String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0  )";String rj2 ="(z/a0)**(1.5)*2*sympy.exp(-z*r2/a0  )";String rj3="(z/( 2*a1))**(1.5)*(2-z*r1/a1)*sympy.exp(-z*r1/(2*a1)  )";String rj4="(z/( 2*a1))**(1.5)*(2-z*r2/a1)*sympy.exp(-z*r2/(2*a1)  )";String rj5="(z/( 2*a1))**(1.5)*(z*r1/(3**0.5*a1) )*sympy.exp(-z*r1/(2*a1)  )";String rj6="(z/( 2*a1))**(1.5)*(z*r2/(3**0.5*a1) )*sympy.exp(-z*r2/(2*a1)  )";double fh=hin(fx1,fx1)*2+hin(fx2,fx2)*2+hin(fx3,fx3)*2;double fj=ak( rj1 , rj2 ,0,0,0,0)+ak( rj1 , rj4 ,0,0,0,0)*4+ak( rj1 , rj6 ,0,0,1,cn6)*4+ak( rj3 , rj4 ,0,0,0,0)+ak( rj3 , rj6 ,0,0,1,cn6)*4+ak( rj5 , rj6 ,1,cn5,1,cn6);double fk=bk( rj1 ,rj4, rj2, rj3,0,0 ,0,0)*4+ bk( rj1 ,rj6, rj2, rj5,0,0 ,1,cn6)*4+bk( rj3 ,rj6, rj4, rj5,0,0 ,1,cn6)*4;        System.out.println( fh+fj-fk +"  "+fh+"  "+fj+"  "+fk +"  "+(fh+fj-fk)*27.2 +"  "+(fh+fj-fk)/37.68 );//System.out.println(  ak( rj5 , rj6 , 1 ,-1 ,1, 1 )  );}

*徐光宪《量子化学》

计算碳原子系综的能级C相关推荐

亚1纳米制程晶体管，一个碳原子栅极厚度：清华重大突破登上Nature
来源:机器之心本文约1800字,建议阅读8分钟小到晶体管栅极开关时等效尺寸都不一样了. 从 20 世纪 50 年代集成电路问世以来,硅晶体管像摩尔定律预测的那样逐渐缩小.微芯片上的晶体管数量越来越 ...
计算硼原子的基态能级的java程序
这个程序由2个java程序,1个python程序共3部分构成,具体包括 public class eneB:用于计算硼原子的动能,势能,库仑能和相互作用能 public class overlap3: ...
石墨烯是一种由碳原子紧密堆积构成的二维晶体，你了解它的特性和应用吗？
石墨烯石墨烯是一种由碳原子紧密堆积构成的二维晶体,是石墨中的单层薄片晶体.它是人类已知强度最高.韧性最好.重量最轻.透光率最高.导电性最佳的材料. 1.石墨烯的特性石墨烯是目前已知的最薄的一种材料 ...
计算Be原子基态能级
用类氢轨道计算Be原子基态能级,Be原子基态4个核外电子的轨道分别是1s和2s,1s轨道的半径是1, 2s轨道的半径是4.轨道分别是计算能级因为E1=E2,E3=E4,J13=J14=J23=J2 ...
用类氢轨道计算交换积分和Li原子2S谱项能级
Li原子核外有3个电子,让他们的轨道分别是(1,0,0)和(2,0,0),用类氢轨道代入,3个电子的轨道分别是其中z是核电荷数等于3,a0是第一层轨道半径等于1,a1是第二层轨道半径等于4. 将能量 ...
计算硼原子的基态能级B---交换能
5个电子轨道的径向部分是, String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )"; String rj2 ="(z/a0 ...
计算氦原子的基态能级
把He核外的两个电子看作有库仑斥力的类H核外电子,由H的径向函数和球谐函数得到He核外电子的波函数得到总的波函数于是能级为其中哈密顿量电子1的动能+势能+电子2的动能+势能+两个电子的库仑排 ...
计算硼原子的基态能级B---库仑排斥能
5个电子的轨道的径向部分是, String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )"; String rj2 ="(z/a ...
计算硼原子基态能级B---动能和势能
设B原子核外5个电子的轨道分别是这5个电子的主,角,磁量子数分别是(1,0,0), (1,0,0), (2,0,0), (2,0,0), (2,1,0). 硼原子的能量为其中E1=E2,E3=E4 ...

计算碳原子系综的能级C

计算碳原子系综的能级C相关推荐

最新文章

热门文章