?

?

?

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

總結

以上是生活随笔為你收集整理的计算硼原子的基态能级B---库仑排斥能的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。