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電子的波函數，一個是α位一個是β位，對應兩種自旋。

計算相互作用能用公式

比如計算k15，兩個電子的量子數分別是（1，0，0）（2，1，0）

計算k=0-4的bk值

因為只有b1不為0，因此只計算G1

因此

用同樣的辦法計算K13=0.057，k35=0.036

因此相互作用能

double fk=

bk( rj1 ,rj4, rj2, rj3,0,0 ,0,0)*4+

bk( rj1 ,rj6, rj2, rj5,0,0 ,1,0)*2+

bk( rj3 ,rj6, rj4, rj5,0,0 ,1,0)*2;

=0.0570022511122276*4+6.244378056112666E-4*2+ 0.03662109375*2=0.3025

因此B原子的基態能級

=-29.1+5.549114061719403-0.3025=-23.85Hartrees=648.72ev

用Hartree-Fork算法得到的值是666.94ev，計算值是Hartree-Fork值的97.26%。

*徐光憲《量子化學》p542

總結

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

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