5個(gè)電子軌道的徑向部分是，

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電子的波函數(shù)，一個(gè)是α位一個(gè)是β位，對(duì)應(yīng)兩種自旋。

計(jì)算相互作用能用公式

比如計(jì)算k15，兩個(gè)電子的量子數(shù)分別是（1，0，0）（2，1，0）

計(jì)算k=0-4的bk值

因?yàn)橹挥衎1不為0，因此只計(jì)算G1

因此

用同樣的辦法計(jì)算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原子的基態(tài)能級(jí)

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

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

*徐光憲《量子化學(xué)》p542

總結(jié)

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

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。