求助，非平衡溶剂化l错误

zoey

老师你好，我再重复一篇文章中的计算时遇到 “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

sobereva

“色重组能”意义不明

原文中把溶剂色重组能认为是总的重组能

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

zoey

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能量间的关系。文献原文如下图

sobereva

明显应当去掉弥散函数，对当前体系计算毫无意义，还令SCF难收敛得多。跟调整结构没有必然关系

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

zoey

已解决，调整结构可以解决