计算硼原子的基态能级B---库仑排斥能

 

 

 

5个电子的轨道的径向部分是，

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)  )";

其中rj5和rj6都是第二层p电子的波函数，一个是α位一个是β位，对应两种自旋。

J12，J13，J34都是s轨道电子和s轨道电子的相互作用可以直接用公式

经计算得J12=3.125，J13=0.3，J34=0.188

但J15和J35是s电子和p电子之间的相互作用，需用公式

比如计算J15，两个电子的量子数分别是（1，0，0）（2，1，0）

因为只涉及到p轨道所以k可以只计算到4        

                                                                                         

计算Fk

经验算如果r1和r2大小顺序颠倒过来

 得到的值是一样的。

因为a1到a4都是0，所以没有计算F1到F4

因此电子15的库仑排斥能J15是1*0.31+0+0+0+0=0.31，用同样的办法计算电子35的库仑排斥能为0.2.

因此B的库仑排斥能

double fj=

ak( rj1 , rj2 ,0,0,0,0)+

ak( rj1 , rj4 ,0,0,0,0)*4+

ak( rj1 , rj6 ,0,0,1,0)*2+

ak( rj3 , rj4 ,0,0,0,0)+

ak( rj3 , rj6 ,0,0,1,0)*2;

=3.125+0.3015292384291021*4+0.312367694626497*2+

0.18798828125000008+0.20263671875*2

=5.549114061719403 Hartrees

*徐光宪《量子化学》p542