计算Be原子基态能级

用类氢轨道计算Be原子基态能级，Be原子基态4个核外电子的轨道分别是1s和2s，1s轨道的半径是1， 2s轨道的半径是4。轨道分别是

计算能级

因为E1=E2，E3=E4，J13=J14=J23=J24,K13=K14=K23=K24

所以原式可以化简

分别展开

计算库仑积分,因为都是s轨道可以用公式

由此计算J12

同样计算J13=0.24，J34=0.15

计算交换积分

因此Be原子基态能级为=2*(-8-0.875)+2.5+0.24*4+0.15-0.0456*4=-14.317Hartrees，Hartree-Fock算法的值是-14.57，是Hartree-Fock算法的98.2%。

import sympy
import math
from sympy import symbols, cancela = sympy.Symbol('a')
e = sympy.Symbol('e')
m = sympy.Symbol('m')
h = sympy.Symbol('h')
l = sympy.Symbol('l')
lp = sympy.Symbol('lp')
r = sympy.Symbol('r')
EE = sympy.Symbol('EE')
R = sympy.Symbol('R')
r1 = sympy.Symbol('r1')
r2 = sympy.Symbol('r2')
r3 = sympy.Symbol('r3')c1 = sympy.Symbol('c1')
c2 = sympy.Symbol('c2')
c3 = sympy.Symbol('c3')μ = sympy.Symbol('μ')
v = sympy.Symbol('v')α = sympy.Symbol('α')
β = sympy.Symbol('β')x = sympy.Symbol('x')
y = sympy.Symbol('y')
z = sympy.Symbol('z')θ1= sympy.Symbol('θ1')
θ2= sympy.Symbol('θ2')
Φ1= sympy.Symbol('Φ1')
Φ2= sympy.Symbol('Φ2')θ= sympy.Symbol('θ')
Ψ= sympy.Symbol('Ψ')
Φ= sympy.Symbol('Φ')
pi=sympy.Symbol('pi')
E=sympy.Symbol('E')
I=sympy.Symbol('I')
sin=sympy.Symbol('sin')
cos=sympy.Symbol('cos')
diff=sympy.Symbol('diff')
integrate=sympy.Symbol('integrate')pi=sympy.pi
E=sympy.E
sin=sympy.sin
cos=sympy.cos
diff=sympy.diff
integrate=sympy.integrate#Be原子基态能级def hin( fx1 ,fx2, z ):fx = fx1z=z# 拉普拉斯算符f1 = (1 / (r * r)) * diff((r * r * diff(fx, r)), r)f2 = (1 / (r * r * sin(θ))) * diff((sin(θ) * diff(fx, θ)), θ)f3 = (1 / (r * r * sin(θ) * sin(θ))) * diff(fx, Φ, Φ)f8 = fx2*(-1 / 2) * (f1 + f2 + f3)# print   (   f1 )# print   (   f2 )# print   (   f3 )# print    ( f8 )# 球坐标积分  动能f9 = (integrate((integrate(integrate(f8 * r * r * sin(θ), (r, 0, float('inf'))), (θ, 0, pi))), (Φ, 0, 2 * pi)))# print(f9)f10 = fx2 * (-z / r) * fx# 势能f11 = (integrate((integrate(integrate(f10 * r * r * sin(θ), (r, 0, float('inf'))), (θ, 0, pi))), (Φ, 0, 2 * pi)))# print(f11)print("H", f9 + f11)return f9 + f11def jin (fr1 ,fr2 ):f21 = fr1 * fr2 * (1 / r1) * fr1 * fr2 * r1 * r1 * r2 * r2f22 = fr1 * fr2 * (1 / r2) * fr1 * fr2 * r1 * r1 * r2 * r2f23 = (integrate(f21, (r2, 0, r1)))f24 = (integrate(f22, (r2, r1, float('inf'))))f25 = (16 * pi ** (2) * integrate(f24 + f23, (r1, 0, float('inf'))))# print("f23",f23)# print("f24",f24)print("J", f25)return f25def kin (fr1,fr2,fr3,fr4 ):# 交换积分#fr1 = (z) ** (1.5) * sympy.exp(-z * r1) * pi ** (-0.5)#fr2 = (z) ** (1.5) * sympy.exp(-z * r2) * pi ** (-0.5)#fr3 = (z) ** (1.5) * sympy.exp(-z * r2) * pi ** (-0.5)#fr4 = (z) ** (1.5) * sympy.exp(-z * r1) * pi ** (-0.5)f21 = fr1 * fr2 * (1 / r1) * fr3 * fr4 * r1 * r1 * r2 * r2f22 = fr1 * fr2 * (1 / r2) * fr3 * fr4 * r1 * r1 * r2 * r2f23 = (integrate(f21, (r2, 0, r1)))f24 = (integrate(f22, (r2, r1, float('inf'))))f36 = (16 * pi ** (2) * integrate(f24 + f23, (r1, 0, float('inf'))))# print("f23",f23)# print("f24",f24)print("K", f36)return f36def sab (fr1 ,fr2 ):# print("f23",f23)# print("f24",f24)f9 = (integrate((integrate(integrate( fr1*fr2 * r * r * sin(θ), (r, 0, float('inf'))), (θ, 0, pi))), (Φ, 0, 2 * pi)))print("S", f9)return f9z=4a0=1
a1=4fx1=(z/a0)**(1.5)*2*sympy.exp(-z*r/a0  )*(4*pi)**(-0.5)
fx2=(z/( 2*a1))**(1.5)*(2-z*r/a1)*sympy.exp(-z*r/(2*a1)  )*(4*pi)**(-0.5)fj1 = (z/a0)**(1.5)*2*sympy.exp(-z*r1/a0  )*(4*pi)**(-0.5)
fj2 = (z/a0)**(1.5)*2*sympy.exp(-z*r2/a0  )*(4*pi)**(-0.5)fj3=(z/( 2*a1))**(1.5)*(2-z*r1/a1)*sympy.exp(-z*r1/(2*a1)  )*(4*pi)**(-0.5)
fj4=(z/( 2*a1))**(1.5)*(2-z*r2/a1)*sympy.exp(-z*r2/(2*a1)  )*(4*pi)**(-0.5)fh=hin(fx1, fx1, z)*2+ hin(fx2, fx2, z)*2fj=jin(fj1,fj2)+jin(fj1,fj4)*4+jin(fj3,fj4)fk=kin( fj1 , fj4 , fj2 , fj3 )*4print(   fh+fj-fk , (fh+fj-fk)/14.57   )

计算Be原子基态能级相关推荐

计算硼原子基态能级B---动能和势能
设B原子核外5个电子的轨道分别是这5个电子的主,角,磁量子数分别是(1,0,0), (1,0,0), (2,0,0), (2,0,0), (2,1,0). 硼原子的能量为其中E1=E2,E3=E4 ...
计算硼原子的基态能级的java程序
这个程序由2个java程序,1个python程序共3部分构成,具体包括 public class eneB:用于计算硼原子的动能,势能,库仑能和相互作用能 public class overlap3: ...
计算碳原子系综的能级C
C原子核外6个电子的排布是轨道的径向部分是 String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )"; String rj2 =&q ...
计算硼原子的基态能级B---交换能
5个电子轨道的径向部分是, String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )"; String rj2 ="(z/a0 ...
计算硼原子的基态能级B---库仑排斥能
5个电子的轨道的径向部分是, String rj1 ="(z/a0)**(1.5)*2*sympy.exp(-z*r1/a0 )"; String rj2 ="(z/a ...
计算氦原子的基态能级
把He核外的两个电子看作有库仑斥力的类H核外电子,由H的径向函数和球谐函数得到He核外电子的波函数得到总的波函数于是能级为其中哈密顿量电子1的动能+势能+电子2的动能+势能+两个电子的库仑排 ...
轨道半径对氦原子基态能级的影响He
设He原子的两个核外电子的轨道分别是 He原子的基态能级为其中a0是半径,让a0由0.2增加到2.2步长是0.01,看看a0的变化对E有什么影响得到的数据画成图当a0=1.19时,E取最小值,E ...
计算类氦离子基态能级z=1-103从氢到铹
让原子核外有两个电子,让核电荷数z=1-103从H到Lr. 两个1s电子波函数为能级为分成三部分,两个电子的动能和势能,加上两个电子之间的排斥能将第一部分打开第二部分与第一部分操作方法相同打 ...
python ovito模块计算某一类原子的MSD均方位移
此文章为用python ovito模块计算某一类原子的MSD ovito.modifiers.CalculateDisplacementsModifier 基于ovito.pipline.Refer ...

计算Be原子基态能级

计算Be原子基态能级相关推荐

最新文章

热门文章