TD优化第一激发态结构不收敛，疑似出现势能面交叉

木子叶

%nprocshared=28
%mem=25GB
%chk=20231012-2ohcz-alkali-is-S.chk
# opt=(tight,noeigen) freq td=(singlets,nstates=6,root=1)/6-31g(d)
scrf=(iefpcm,solvent=methanol) geom=connectivity em=gd3 int=superfine
m062x

Title Card Required

-1 1

（高斯09版本，超算计算），请问只能尝试优化S2的结构吗？S1的结构无法优化了吗？

sobereva

问题交代清楚，什么叫前一个、后一个

没事甭用tight，明显更难达到收敛限

输出文件压缩后上传，否则别人根本不知道优化轨迹、结构什么样

tiandikuoyuan

你这不是势能面交叉，这是结构优化振荡，可以参考http://sobereva.com/164。

木子叶

tiandikuoyuan 发表于 2023-10-19 12:12
你这不是势能面交叉，这是结构优化振荡，可以参考http://sobereva.com/164。

好的，谢谢，我按这篇博文尝试了4次不同方式都没有成功，我再试试

tiandikuoyuan

木子叶 发表于 2023-10-19 19:22
好的，谢谢，我按这篇博文尝试了4次不同方式都没有成功，我再试试

你这个分子很小，直接尝试用opt=calcall每一步计算Hessian矩阵

木子叶

tiandikuoyuan 发表于 2023-10-19 20:06
你这个分子很小，直接尝试用opt=calcall每一步计算Hessian矩阵

td下calcall使用无法使用，我手动添加calcall关键词，计算报错

木子叶

tiandikuoyuan 发表于 2023-10-19 20:06
你这个分子很小，直接尝试用opt=calcall每一步计算Hessian矩阵

%nprocshared=28
%mem=25GB
%chk=20231012-2ohcz-alkali-is-S.chk
# opt=(tight,noeigen,recalc=1) freq
td=(singlets,nstates=6,root=1)/6-31g(d) scrf=(iefpcm,solvent=methanol)
geom=connectivity em=gd3 int=ultrafine m062x

Title Card Required

-1 1
尝试加recalc=1，也无法计算，是不是td不能用calcall关键词？

zjxitcc

木子叶 发表于 2023-10-21 15:25
%nprocshared=28
%mem=25GB
%chk=20231012-2ohcz-alkali-is-S.chk

应直接在6L中进行编辑，无需再开一层。另外：
（1）当你说报错时一定要顺带说出报错原因（如果你不懂，就截图展示输出文件最后15行，让别人帮你看报错原因），你这种情况还要说出Gaussian版本，不然的话回答者只能像算命一样给你瞎猜。
（2）opt=tight并不能帮助收敛，对你当前问题毫无用处。
（3）td()括号里singlets,root=1是默认的，不用写。scrf()里iefpcm是默认的，不用写。对于g16来说int=ultrafine是默认的，不用写。geom=connectivity和坐标底下的连接关系毫无用处，建议删去。

木子叶

zjxitcc 发表于 2023-10-21 20:36
应直接在6L中进行编辑，无需再开一层。另外：
（1）当你说报错时一定要顺带说出报错原因（如果你不懂， ...

不好意思，请问“直接在6L中进行编辑，无需再开一层”这句话什么意思？6L是什么软件吗？我百度也没有找到

zjxitcc

木子叶 发表于 2023-10-21 22:17
不好意思，请问“直接在6L中进行编辑，无需再开一层”这句话什么意思？6L是什么软件吗？我百度也没有找到

本帖6楼的意思

木子叶

zjxitcc 发表于 2023-10-21 22:58
本帖6楼的意思

谢谢