计算化学公社

标题: 激发态振动分析疑问，先opt，再读取chk进行振动分析，出现NO [打印本页]

作者
Author: 倪宇晴 时间: 2018-8-31 10:18
标题: 激发态振动分析疑问，先opt，再读取chk进行振动分析，出现NO
本帖最后由倪宇晴于 2018-8-31 18:32 编辑

请教各位老师，我进行激发态振动分析（gaussian 16）的时候遇到了这样一个问题：
1. 先opt，再读取chk进行振动分析出现NO；
关键词是：# opt td=(root=1,triplet)/6-311g(d,p) m062x density
# freq td=(root=1,triplet)/6-311g(d,p) m062x density guess=read geom=check

2.利用关键词# opt freq td=(root=1,triplet)/6-311g(d,p) m062x density进行计算，最后的结果收敛了。

PS:看2的计算结果的时候发现其振动分析的关键词是 #N Geom=AllCheck Guess=TCheck SCRF=Check GenChk RM062X TD(FC,Read)/6-311G(d,p) Freq。在手册的TD下面只找到了read，不知道这个FC是什么意思。FC是Frozen Core么？

想请教一下各位老师，1.这两个方法的区别在什么地方？是哪一步出现了问题？
2.优化激发态时,如果想获取激发态的波函数是不是同样需要加density。

作者
Author: sobereva 时间: 2018-8-31 18:21
“振动分析不收敛”意义不明，SCF不收敛还是CPHF不收敛需要说清楚
仔细看输出文件弄清楚原因，至少把情况表达清楚

opt freq任务自动生成的freq那一步的关键词你不用管

尽量别把不确定能否用在一起的关键词用在一起，否则可能导致耗时莫名其妙增加或者结果诡异或者报错。获得激发态波函数，在优化和振动分析完了，单独再做一个任务，心里放心得多

作者
Author: 倪宇晴 时间: 2018-8-31 18:31

sobereva 发表于 2018-8-31 18:21
“振动分析不收敛”意义不明，SCF不收敛还是CPHF不收敛需要说清楚
仔细看输出文件弄清楚原因，至少把情况 ...

谢谢老师，不收敛是做振动分析是出现了NO的情况。Maximum Force，Maximum Displacement等四个全都是NO。

作者
Author: sobereva 时间: 2018-8-31 18:55

倪宇晴发表于 2018-8-31 18:31
谢谢老师，不收敛是做振动分析是出现了NO的情况。Maximum Force，Maximum Displacement等四个全都是NO。

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

作者
Author: 倪宇晴 时间: 2018-9-2 16:16

sobereva 发表于 2018-8-31 18:55
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

好的，谢谢老师。我先看看帖子。

作者
Author: 倪宇晴 时间: 2018-9-3 11:30
本帖最后由倪宇晴于 2018-9-3 17:10 编辑

sobereva 发表于 2018-8-31 18:55
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

老师，我发现先用opt td=(root=1,triplet)/6-311g(d,p) m062x，之后用 freq td=(root=1,triplet)/6-311g(d,p) m062x guess=read geom=check计算的时候会出现NO，用 freq td=(root=1,triplet,read)/6-311g(d,p) m062x geom=check计算的时候结果就和opt freq一致，没有出现问题。这是不是因为guess=read读取的是chk文件中基态的波函数作为初猜进行振动分析，而TD=read读取的是chk文件中激发态的波函数进行振动分析。

作者
Author: sobereva 时间: 2018-9-4 03:47

倪宇晴发表于 2018-9-3 11:30
老师，我发现先用opt td=(root=1,triplet)/6-311g(d,p) m062x，之后用 freq td=(root=1,triplet)/6-311g( ...

程序做TDDFT的时候，先对当前结构做基态DFT计算产生基态轨道，然后再做电子激发计算，这个过程中要求解组态系数，两个过程都需要迭代求解。guess=read读的是对第一步而言的，TD=read读的是对第二步而言的

出现NO不代表有问题，而且即便opt freq直接联用，照样可能freq之后有NO。

之后单独做freq的时候，guess=read TD=read一起用，有助于减少耗时，而且能容易保证波函数承接优化过程最后一步的波函数

作者
Author: 倪宇晴 时间: 2018-9-4 09:06

sobereva 发表于 2018-9-4 03:47
程序做TDDFT的时候，先对当前结构做基态DFT计算产生基态轨道，然后再做电子激发计算，这个过程中要求解组 ...

谢谢老师！

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