计算化学公社

标题: 激发态计算关键词使用问题 [打印本页]

作者
Author:
liuyjhx    时间: 2019-6-14 18:54
标题: 激发态计算关键词使用问题
#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)
老师看到激发态命令有这么两种计算方法,不知道应该用哪一种?


作者
Author:
暖空    时间: 2019-6-14 20:04
第一种看命令应该是要做S1态的优化,
温度我没写过,你可以在高斯手册里搜一下关键词看看应该怎写。
除温度外的S1态优化用#p opt CAM-B3LYP/6-31g* scrf=(smd,solvent=dichloromethane) ,最好加上freq一起计算。默认的就是root=1 nstates=3. singlet不用写,也是默认的,scrf是考虑溶剂效应,模型用的SMD。
第二种命令是做50个态的垂直激发,溶剂模型采用的是默认的PCM
作者
Author:
暖空    时间: 2019-6-14 20:07
http://sobereva.com/314
可以看一下这个帖子,sob老师写的很全面
作者
Author:
三石草祭    时间: 2019-6-14 20:50
暖空 发表于 2019-6-14 20:04
第一种看命令应该是要做S1态的优化,
温度我没写过,你可以在高斯手册里搜一下关键词看看应该怎写。
除温 ...

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

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

作者
Author:
sobereva    时间: 2019-6-14 22:35
写温度根本毫无意义,对优化结果完全没影响,仔细看此文
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423http://bbs.keinsci.com/thread-10159-1-1.html

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

只要仔细读了此文,就不可能再有当前问的问题
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314http://bbs.keinsci.com/thread-2413-1-1.html
作者
Author:
sobereva    时间: 2019-6-14 22:36
暖空 发表于 2019-6-14 21:21
我做opt的同时一般都会加上freq做频率分析,话说做激发态的优化可以不用计算频率么?

完全可以。但是你得有自信确保优化得到的是极小点
作者
Author:
暖空    时间: 2019-6-15 14:35
sobereva 发表于 2019-6-14 22:36
完全可以。但是你得有自信确保优化得到的是极小点

谢谢老师的回复,我做激发态优化时(目前就用tddft做过s1的优化)都是直接用的基态构型,或者手动破坏一下对称性再进行opt,td ,freq。
不过在重复一个文献分子的时候确实发现了点儿问题,直接在基态构型下优化的s1和文献并不一样。但是把初猜调成文献上的S1态构型,就能优化出来一样的构型。
请问老师,在优化激发态的时候,是不是需要先扫一下势能面,然后再给初猜做激发态优化呢?
作者
Author:
sobereva    时间: 2019-6-16 00:43
暖空 发表于 2019-6-15 14:35
谢谢老师的回复,我做激发态优化时(目前就用tddft做过s1的优化)都是直接用的基态构型,或者手动破坏一 ...

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

作者
Author:
暖空    时间: 2019-6-16 14:57
sobereva 发表于 2019-6-16 00:43
S1势能面可能存在多个极小点,初猜结构不同可能收敛到不同的
如果想找全,可以做扫描

明白了,谢谢老师!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3