高斯溶剂化计算出错

huage

我最近在做有机分子的溶剂化紫外吸收的计算，这是我的溶剂化吸收的命令%chk=2.chk
%mem=2GB
%nprocshared=12
# b3lyp/6-311g(d)  SCRF=(Solvent=Dichloromethane,Read) Geom=Check Guess=Read

ground

0 1

NonEq=write

--Link1--
%chk=2.chk
%mem=2GB
%nprocshared=12
# b3lyp/6-311g(d)  SCRF=(Solvent=Dichloromethane,externaliteration,Read) TD(NStates=20,Root=1)  Geom=Check Guess=Read

excited

0 1

NonEq=read


在同时提交的问题中有一个出错，信息是Using charges instead of weights in PCMQM.
SC-PCM: The non-equilibrium contribution to the free energy is   -0.022523026 a.u.
SC-PCM: The total energy including the PCM contribution is     -1830.183165912 a.u.
SC-PCM: Erms=  0.42553D-04 Emax=  0.26157D-03 Qconv=  0.10000D-06
SC-PCM: Iterating to achieve self-consistent PCM reaction field, IExtIt=  50.
SC-PCM: Convergence not achieved within  50 iterations ... aborting job.
Error termination via Lnk1e in /opt/g09-D01/g09/l124.exe at Mon Jan 26 05:38:12 2015.
Job cpu time:      35 days 17 hours  9 minutes  4.8 seconds.
File lengths (MBytes):  RWF=   4512 Int=      0 D2E=      0 Chk=    492 Scr=      1
我想请各位老师帮忙看看错误在哪？同时想请问我做的吸收波长与实验测得的相差有将近100nm,请问各位老师又什么建议，谢谢

小范范1989

输出结果是不是写的很清楚？
SC-PCM: Convergence not achieved within  50 iterations ... aborting job.
应该解决这个问题是吧？maxcyc设置？
再就是，我感觉运行这么久，需要这么细致吗？
可以不可以不考虑平衡溶剂非平衡溶剂的影响，按照我们往常的来做，应该有差距，但是不会很大吧？
我也不太明白，只是一家之言，坐等高手指点。

sobereva

溶剂模型计算时出错。虽然是表面电荷迭代不收敛，但具体解决方法不好说，只能多试试，比如不用非平衡溶剂或不用外迭代，基组改一下比如改成6-31G*。
另外强烈建议用#P，多输出些中间信息总是对解决问题有帮助。

huage

小范范1989 发表于 2015-1-27 22:22
输出结果是不是写的很清楚？
SC-PCM: Convergence not achieved within  50 iterations ... aborting job. ...

我也尝试了平衡溶剂与非平衡溶剂，计算的结果两者相差也不大，但是跟实验的数据相差就很大了···谢谢小范范童鞋

huage

sobereva 发表于 2015-1-28 00:23
溶剂模型计算时出错。虽然是表面电荷迭代不收敛，但具体解决方法不好说，只能多试试，比如不用非平衡溶剂或 ...

谢谢sob老师，我再试试加大内存和加大次数试试

hcxytpp@163.com

您好，最近我计算溶剂化能的时候遇上了同样的问题，体系中只有一个结构出现了这样的情况，数据比较急用，如果说是换基组等等去尝试应该不太方便，想请教楼主和Sob老师，有没有好的建议，我的输入是：# ub3lyp/genecp scrf=(smd,solvent=ethanol,read,dovacuum,externaliteration)
radii=uaks 恳请老师指点