冰释之川 发表于 2017-3-19 15:07 嗯嗯，明白了！谢谢老师！ |
本帖最后由 冰释之川 于 2017-3-19 15:08 编辑 万卷书万里路 发表于 2017-3-19 14:35 一般情况下，溶剂效应对分子的几何结构影响较小（参见 http://sobereva.com/327 1.3节<隐式溶剂模型对体系性质的影响>）；如果你有实验值对照，那么最好考虑溶剂效应。 优化几何结构的过程和后面的泛函调控没有必然联系。泛函调控是针对你优化好的几何结构下做相应泛函的优化，此时你的分子几何结构已经确定下来了。 |
ggdh 发表于 2017-3-14 00:23 老师您好！ 这几天看了帖子后又计算产生了疑问， 第一条您说，不准，因为没加溶剂没有调控，我想问，如果是先用CAM-B3LYP泛函优化激发态结构，然后再用调控的泛函在加溶剂效应的情况下计算激发态单点能，那在前一步优化的过程（CAM-B3LYP）中我还需要加溶剂效应吗？ 也就是说，优化结构的过程溶剂对最后的结构有没有影响，优化结构的过程中是否加调控又有没有影响？ 谢谢了！ |
赵云跳槽 发表于 2017-3-14 16:39 好的，谢谢！我再仔细看看 |
本帖最后由 赵云跳槽 于 2017-3-14 16:42 编辑 万卷书万里路 发表于 2017-3-14 14:28 S1的绝热能=E(TD-HF/TD-KS)+HF能量 td opt或tda opt计算 T1的绝热能=E(TD-HF/TD-KS)+HF能量 td(triplets) opt或tda(triplets) opt计算，或 opt使用0 3 对于HF/DFT计算 HF能量=SCF Done的能量 垂直指的是不优化构型，你图中的E_S1和E_T1都是优化后的，都是绝热的 还是认真看sob的博文吧 Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法 http://sobereva.com/314 优化长程校正泛函w参数的简便工具optDFTw http://sobereva.com/346 |
ggdh 发表于 2017-3-14 00:23 谢谢老师的解答！ 我还有两个问题， 1.“不写TD， 0 1的情况下是基态， 0 3的情况下是T1” 那我想，求delta ST应该是S1和T1的能量用同样的方法水平下计算出来、意义对应的上的才对吧， 我单纯只在0 3 下优化结构得到的SCF Done能量（也就是文件末尾的HF能量）就是ET1，对于求绝热delta_ST，对应要用的S1态能量应该是# opt td=singlets 的E(TD-HF/TD-KS)数值吗？还是别的？ 2.要是我还想更精确点，“用调控了参数的泛函再计算T1态的单点能”，只想计算T1的单点能，我拿 0 3 下优化过的结构再： # LC-wPBE/tzvp IOP（···）IOP（···） 0 3 我这样写对吗？ 第一种方法，用非调控的泛函（比如B3LYP）优化结构之后再用调控了的来计算单点能，和第二种方法：直接用调控参数的泛函优化激发态结构并直接读取能量，这两种哪一种更能准确一点？ 3.和这句“如果是0 3算的T1,那么scf done 是ET1，如果是0 1+ td/tda算的T1，那么scf done 是B点能量。” 这里的“如果是0 1+ td/tda算的T1，”具体指哪种： （1）# td=triplets 0 1 （输入文件的结构是优化了的 0 3 结构） （2）# td=triplets 0 1 （输入文件的结构是优化了的 0 1 结构，基态s0结构） （3）# opt td=triplets 0 1 （输入文件的结构是未经优化的初始结构） 还是这三个都不是， 还想知道，这三种计算各自的意义，能不能在势能面上找到各自输出文件里的E(TD-HF/TD-KS)的位置？ 感觉我问的也有点问题，但是也说不清楚，麻烦老师能指导！ |
赵云跳槽 发表于 2017-3-14 11:12 您好！ 我不能理解这个：“S1/T1的绝热能=E(TD-HF/TD-KS)+HF能量” 这个绝热能的意思是绝热发射（激发）能吗？ 如果是的话，是不是应该是减号？HF能量等于SCF Done的能量是不是。 还有这句：“如果是垂直能，那就直接等于两者E(TD-HF/TD-KS)能量差，因为此时的S1/T1对应的基态结构是完全相同的。” 两者E(TD-HF/TD-KS)的时候不是上面那个图里的ES1和ET1吗，怎么着对应的基态结构就一样了？还请解释一下，谢谢啦！ |
本帖最后由 ggdh 于 2017-3-14 11:31 编辑 冰释之川 发表于 2017-3-14 06:53 w跟结构有关，也跟轨道的形状有关。一般在优化激发态过程中。轨道的基本性质不会变（共轭PI轨道），那么w变化就不会太大。统一用基态结构的就好。 |
wyf 发表于 2017-3-14 10:34 delta_ST绝热能=S1的绝热能-T1的绝热能 S1/T1的绝热能=E(TD-HF/TD-KS)+HF能量 如果是垂直能，那就直接等于两者E(TD-HF/TD-KS)能量差，因为此时的S1/T1对应的基态结构是完全相同的 |
冰释之川 发表于 2017-3-14 06:53 文献的做法是：基态分子进行这一系列的w参数调控 |
冰释之川 发表于 2017-3-14 10:48 非常感谢老师的解答，这边还有一个问题，那我们在计算单-三线态的驱动能△G时，也是用到的这个绝热能差么？还是用到的两个HF能量相减得到的能量呢？因为之前听到一个老师讲的是，这个驱动能等于二者零点矫正能相减，如果我没有做频率计算的话，那这个能量是不是就是两个HF相减得到的能量呢？这个驱动能跟绝热能差有什么区别呢？ |
wyf 发表于 2017-3-14 10:34 如果都用TDDFT算，则绝热ΔEst 为两者E(TD-HF/TD-KS)能量差，并不是HF能量差（你还要考虑激发能） |
| 这边还有一个问题想请教一下sob老师，如果我都用td的方法去计算S1态和T1态的能量，那我的绝热能应该是两个total energy相减，还是两个HF能量相减呢？ |
本帖最后由 冰释之川 于 2017-3-14 07:44 编辑 ggdh 发表于 2017-3-14 00:23 钟老师的回答我很赞同。另外咨询一下，关于用调控泛函算绝热ΔEst，由于涉及到S1和T1的平衡结构，那么问题来了，调控是分别对这两个结构进行优化w值呢？还是大家都取基态结构下（或者是S1结构下？）的最优w值来算E（S1）和E（T1）呢？ |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-20 21:47 , Processed in 0.175053 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites