计算化学公社

标题: 求助,非平衡溶剂化l错误 [打印本页]

作者
Author:
zoey    时间: 2020-1-28 11:30
标题: 求助,非平衡溶剂化l错误
老师你好,我再重复一篇文章中的计算时遇到 “No special actions if energy rises. Inconsistent number of fixed charges for non-equilibrium.
Error termination via Lnk1e in /data/g09/l502.exe at Tue Jan 28 10:40:04 2020.
”问题,计算的结构是一个自由基,文献中提到溶剂重组能的计算“the solvent reorganization energy was calculated using the non〆quilibrium solvent cage implemented in Gaussian09vD01 (scrf=read keyword with NonEq=write/read input option). ”找不到出现上述错误的原因,试了scf=xqc,scf=conver=6等还是一样的错误。输入文件如下:

%chk=Rf-noneq.chk
%nprocs=16
%mem=10GB
# opt=(calcfc) freq wb97xd nosymm scf=xqc SCRF=(PCM,Solvent=Acetonitrile,read) 6-31+g*

v

0 2
6 -0.977181 1.231882 0.016101
6 -0.411449 0.487959 -1.182507
6 1.141668 0.423902 -1.289692
9 -0.863982 -0.782141 -1.138768
9 -0.867293 1.086301 -2.302330
9 1.646622 1.655171 -1.389541
9 1.653433 -0.163883 -0.206222
9 1.498719 -0.275243 -2.366703
9 -0.730559 2.528907 0.041438
9 -0.754524 0.662437 1.186508

NonEq=write





--link1--
%chk=Rf-noneq.chk
%nprocs=16
%mem=10gb
# wb97xd SCRF=(PCM,Solvent=Acetonitrile,read) 6-31+G* guess=read geom=check

v

0 2

NonEq=read




作者
Author:
zoey    时间: 2020-1-28 20:14
已解决,调整结构可以解决
作者
Author:
sobereva    时间: 2020-1-29 00:55
明显应当去掉弥散函数,对当前体系计算毫无意义,还令SCF难收敛得多。跟调整结构没有必然关系

什么时候才需要弥散函数看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

作者
Author:
zoey    时间: 2020-1-29 11:00
sobereva 发表于 2020-1-29 00:55
明显应当去掉弥散函数,对当前体系计算毫无意义,还令SCF难收敛得多。跟调整结构没有必然关系

什么时候 ...

多谢老师,
我还有一个问题不太明白关于单电子转移能量计算的,原文中把溶剂色重组能认为是总的重组能,计算的时候“the solvent reorganization energy was calculated using the non〆quilibrium
solvent cage implemented in Gaussian09vD01 (scrf=read keyword with NonEq=write/read
input option). This feature allows the user to optimize the solvent cage (SC) for the reactant
and use that frozen solvent cage for the energy calculation of the product (Eproduct).”不太理解怎么优化RC与计算PC能量间的关系。文献原文如下图
作者
Author:
sobereva    时间: 2020-1-29 22:55
“色重组能”意义不明
原文中把溶剂色重组能认为是总的重组能


这么使用隐式溶剂模型,没有直接去用量化程序优化solvent cage,只不过通过恰当地利用非平衡溶剂设置体现了等同的效果。相关信息看
乱谈外重组能的计算
http://sobereva.com/333




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3