24小时热门版块排行榜    

查看: 702  |  回复: 5
当前主题已经存档。

nuaa8222

木虫 (著名写手)

[交流] 【求助】高斯计算电荷分布,出现问题(急)

大家好,我在用高斯计算电荷分布时,输入文件如下:

%chk=title
%mem=2GB
#p b3lyp/3-21g charge scf=tight pop=chelpg

title

0 1
分子描述

背景电荷


输出文件中,出现了下列信息:
Convergence on energy, delta-E=      NaN
SCF Done:  E(RB+HF-LYP) =                 NaN A.U. after    1 cycles
             Convg  =    0.4391D+01             -V/T =     NaN
             S**2   =   0.0000
KE= 9.517492208094D+03 PE=-1.709445150761D+05 EE=                NaN
Leave Link  502 at Sun Dec 20 12:19:45 2009, MaxMem=  268435456 cpu:    3193.2

能量变成NaN  咋回事呀 查了一些资料 找不出原因  虽然最后Normal termination of Gaussian 03     但对计算的准确性 不敢确定  请教高手  什么情况会导致能量变成NaN   这样计算得到的结果可靠吗?谢谢大家乐

[ Last edited by nuaa8222 on 2009-12-21 at 15:59 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abbott

金虫 (著名写手)

不要用QQ问我东西


小木虫(金币+0.5):给个红包,谢谢回帖交流
NaN这个貌似是程序中的....用来表达, 当前值已经不在预先设定的范内.....
Chemistry[]==[]Chem[]is[]try!!!
2楼2009-12-20 17:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hunanxiang

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
不知道你的分子结构是否有问题,以至于优化没有完成
3楼2009-12-20 20:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+3,VIP+0):感谢回答! 12-21 09:29
NaN是not a number的意思,多数是计算中出现了除0的错误。
在此错误提示之前,应该出现了一大堆坐标或者矩阵元已经是NaN了,如果你在#后面写了p的话,可以明确知道是哪个步骤出来的错误,就可以看看出错前是怎么跑的,来估计出错的原因。你只给出最后这几行出错信息,没有监控的过程,大家很难猜的。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
4楼2009-12-20 23:12:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nuaa8222

木虫 (著名写手)

多谢各位的解答,我贴出了详细的输出文件,麻烦大家再帮我看看
5楼2009-12-21 09:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3,VIP+0):我曾经干过4000多基函的opt,嘿嘿 12-22 13:21
嗯,你的体系很生猛啊。2200多个基函。好在只是单点算。
你的输出文件在301就出现了“Nuclear-Charge attraction  =  NaN”
应该是根本还没跑起来就死翘了。
背景电荷当然不能放在原子自身身上了。原子自身的核和电子都是带电的,
你在原子所在的同一位置放置电荷,是什么意思呢?是要把原子核挤飞了,还是要把
电子挤到核里面去?背景电荷我的理解是应该放在分子所处的环境中——根据环境
确定背景电荷的位置。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
6楼2009-12-21 11:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nuaa8222 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见