24小时热门版块排行榜    

查看: 230  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +14 梁梁校校 2026-03-19 16/800 2026-03-26 18:13 by 乐呵呵的追梦人
[考研] 化学工程085602 305分求调剂 +16 RichLi_ 2026-03-25 16/800 2026-03-26 16:18 by bruceku07
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 311求调剂 +4 勇敢的小吴 2026-03-20 4/200 2026-03-25 18:12 by xcjcqu
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 各位老师您好:本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
信息提示
请填处理意见