24小时热门版块排行榜    

查看: 220  |  回复: 0

CQ_fly

新虫 (小有名气)

[求助] 关于核吸引力的积分

用Hartree Fock 近似,现在写到分子积分的那一部分,用的是STO-3G基组,写到核吸引能的时候发现一个问题:

根据多中心的核吸引能计算公式:(A|-Zc /r1c |B)= S * F0[(alpha+beta)|Rp-Rc|^2],其中S很好计算出来,关键问题是后面的F0函数的值。这个函数:F0(t)=1/2(pi/t)^(1/2)*erf(t^(1/2))。其中的erf是error function,那么我的问题是:

如果积分两边都是同一个核,RA=RB,同时我还想计算A核对电子的吸引作用,即Rc=RA,做这个积分就会得到F0函数内的因变量是0(Rp=0,Rc=0),如此一来,整个积分的值就变成0了(因为erf(0)=0)。但是这肯定是不对的。不知道问题出在哪里,金币不多,菜鸟一枚,希望大家能帮我解决一下疑惑。谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 CQ_fly 的主题更新
信息提示
请填处理意见