计算化学公社

标题: TD优化第一激发态结构不收敛，疑似出现势能面交叉 [打印本页]

作者
Author: 木子叶 时间: 2023-10-18 11:12
标题: TD优化第一激发态结构不收敛，疑似出现势能面交叉
本帖最后由木子叶于 2023-10-21 22:14 编辑

%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的结构无法优化了吗？

作者
Author: sobereva 时间: 2023-10-19 06:06
问题交代清楚，什么叫前一个、后一个

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

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

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

作者
Author: 木子叶 时间: 2023-10-19 19:22
本帖最后由木子叶于 2023-10-19 19:29 编辑

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

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

作者
Author: tiandikuoyuan 时间: 2023-10-19 20:06

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

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

作者
Author: 木子叶 时间: 2023-10-19 21:07

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

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

作者
Author: 木子叶 时间: 2023-10-21 15:25

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关键词？

作者
Author: zjxitcc 时间: 2023-10-21 20:36

木子叶发表于 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和坐标底下的连接关系毫无用处，建议删去。

作者
Author: 木子叶 时间: 2023-10-21 22:17

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

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

作者
Author: zjxitcc 时间: 2023-10-21 22:58

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

本帖6楼的意思

作者
Author: 木子叶 时间: 2023-10-21 23:03

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

谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)