| 查看: 827 | 回复: 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 ] |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有71人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有2人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
招收博士一名-宁波大学
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
abbott
金虫 (著名写手)
不要用QQ问我东西
- 应助: 16 (小学生)
- 金币: 1071.6
- 散金: 2787
- 红花: 10
- 帖子: 1015
- 在线: 105.2小时
- 虫号: 452267
- 注册: 2007-11-05
- 性别: GG
- 专业: 计算机硬件技术
2楼2009-12-20 17:46:23
3楼2009-12-20 20:59:52
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
4楼2009-12-20 23:12:53
nuaa8222
木虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1655.4
- 帖子: 1405
- 在线: 212.3小时
- 虫号: 559191
- 注册: 2008-05-16
- 专业: 凝聚态物性 II ：电子结构
5楼2009-12-21 09:11:58
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0):我曾经干过4000多基函的opt，嘿嘿 12-22 13:21
小木虫(金币+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
50