激发态计算关键词使用问题

liuyjhx

#p opt CAM-B3LYP/6-31g* scrf=(smd,solvent=dichloromethane)
temperature=333.15 td(root=1,singlet) density
# CAM-B3LYP/6-31g(d) TD(nstates=50) scrf=(solvent=dichloromethane)
老师看到激发态命令有这么两种计算方法，不知道应该用哪一种？

暖空

第一种看命令应该是要做S1态的优化，
温度我没写过，你可以在高斯手册里搜一下关键词看看应该怎写。
除温度外的S1态优化用#p opt CAM-B3LYP/6-31g* scrf=(smd,solvent=dichloromethane) ,最好加上freq一起计算。默认的就是root=1 nstates=3. singlet不用写，也是默认的，scrf是考虑溶剂效应，模型用的SMD。
第二种命令是做50个态的垂直激发，溶剂模型采用的是默认的PCM

暖空

http://sobereva.com/314
可以看一下这个帖子，sob老师写的很全面

三石草祭

暖空 发表于 2019-6-14 20:04
第一种看命令应该是要做S1态的优化，
温度我没写过，你可以在高斯手册里搜一下关键词看看应该怎写。
除温 ...

GAUSSIAN16的话激发态频率计算还不是很耗时，如果其他版本，加freq关键词要耗时多至少一个数量级

暖空

三石草祭 发表于 2019-6-14 20:50
GAUSSIAN16的话激发态频率计算还不是很耗时，如果其他版本，加freq关键词要耗时多至少一个数量级

我做opt的同时一般都会加上freq做频率分析，话说做激发态的优化可以不用计算频率么？

sobereva

写温度根本毫无意义，对优化结果完全没影响，仔细看此文
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423（http://bbs.keinsci.com/thread-10159-1-1.html）

写density莫名其妙，没有特殊理由不要写

只要仔细读了此文，就不可能再有当前问的问题
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

sobereva

暖空 发表于 2019-6-14 21:21
我做opt的同时一般都会加上freq做频率分析，话说做激发态的优化可以不用计算频率么？

完全可以。但是你得有自信确保优化得到的是极小点

暖空

sobereva 发表于 2019-6-14 22:36
完全可以。但是你得有自信确保优化得到的是极小点

谢谢老师的回复，我做激发态优化时（目前就用tddft做过s1的优化）都是直接用的基态构型，或者手动破坏一下对称性再进行opt,td ,freq。
不过在重复一个文献分子的时候确实发现了点儿问题，直接在基态构型下优化的s1和文献并不一样。但是把初猜调成文献上的S1态构型，就能优化出来一样的构型。
请问老师，在优化激发态的时候，是不是需要先扫一下势能面，然后再给初猜做激发态优化呢？

sobereva

暖空 发表于 2019-6-15 14:35
谢谢老师的回复，我做激发态优化时（目前就用tddft做过s1的优化）都是直接用的基态构型，或者手动破坏一 ...

S1势能面可能存在多个极小点，初猜结构不同可能收敛到不同的
如果想找全，可以做扫描

暖空

sobereva 发表于 2019-6-16 00:43
S1势能面可能存在多个极小点，初猜结构不同可能收敛到不同的
如果想找全，可以做扫描

明白了，谢谢老师！