应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助 高斯计算过程中一直不收敛 scf(maxcyc=500,xqc)
51/1返回列表
查看: 3201  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

忧伤的小猪

银虫 (小有名气)

[求助] 求助 高斯计算过程中一直不收敛 scf(maxcyc=500,xqc)

高斯计算过程中一直不收敛scf(maxcyc=500,xqc) l508错误

[ Last edited by 忧伤的小猪 on 2013-4-1 at 10:31 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-04-01 10:30:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家
gaussian默认的初始轨道占据对于镧系、锕系原子计算来说通常是错的，很可能是难收敛的激发态，或者即使收敛也是高能量的激发态。解决办法：
0，前期工作：查资料，确定基态的多重度和轨道占据情况。
1，给原子加、减足够的电荷，使(n)f轨道占据数为0或14，(n+1)d轨道占据为0，(n+2)s占据为0或2，让原子变成容易收敛的闭壳层离子。保存chk文件，并且做备份。
2，计算中性原子，用guess=read读入1产生的初始占据，参照文献中的轨道占据，必要的话用guess=(read,alter)交换轨道。
3，遇到能量振荡，远未达到收敛标准的情况，先恢复原来备份的chk文件，然后交换不同编号的简并轨道，重复2。

例如，计算NpF6的原子化能。在B3LYP级别，使用gaussian默认产生的占据，得到的结果是780 kcal/mol，实验值为680 kcal/mol。人工调整占据之后，为706 kcal/mol。再加上BSSE、旋轨耦合、零点能校正，得到结果703 kcal/mol。
如果换成精度更高的UCCSD(T)，默认占据得到840 kcal/mol，人工调整占据并加上各种校正后，结果为684 kcal/mol，已经非常好了。
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
9楼2013-11-15 07:48:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

yongleli

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2013-04-01 20:14:40
自洽场计算出错。
你用什么算的，如果是DFT建议增加格点密度。

还有请重新看看输入的结构是否合理,
中间每一步的能量（#p的话就会在输出里）
要打出来做个曲线看看，是不是正在收敛，
还是在震荡。
赞一下(1人) 回复此楼
2楼2013-04-01 10:43:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忧伤的小猪

银虫 (小有名气)
引用回帖:
2楼: Originally posted by yongleli at 2013-04-01 10:43:50
自洽场计算出错。
你用什么算的，如果是DFT建议增加格点密度。

还有请重新看看输入的结构是否合理,
中间每一步的能量（#p的话就会在输出里）
要打出来做个曲线看看，是不是正在收敛，
还是在震荡。

计算单个原子，出这样的错
赞一下 回复此楼
3楼2013-04-02 08:49:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujialing

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
是不是你选择的基组不太好？我一般都会换一个大点的基组，或者采用scf=qc命令。
赞一下 回复此楼
越努力，越幸运！
4楼2013-04-02 09:57:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭