计算绝热能相关问题

jiang012582

                                                                   一，  请问大家在计算绝热能差时，想得到左侧的绝热能差图，采用以下计算方法，首先优化基态结构S0，读取输出文件最后一步的E(SCF Done)能量，记为基态S0极小点能量，2，分别优化S1,T1,T2,T3,T4等激发态，坐标均采用优化后的基态坐标，其中S1的极小点能量读取的是，最后一步输出的S1平衡结构下S1态的能量 Total Energy, E(TD-HF/TD-DFT) =  -517.325970004   ，3，使用优化后的基态坐标，单独优化T1,T2,T3,T4,采用的是TDFT的方法，同时做TDDFT计算，输入文件如下# b3lyp/6-31g(d) opt td(triplets,nstates=10,root=1) scrf=(solvent=aceticacid,pcm) ，T2,T3,T4态改root=2,3,4.电荷自旋多重度均为0，1  ，在读取T1态极小点能量时，我读取的是最后一步输出的Total Energy, E(TD-HF/TD-DFT)，我想问在读取T2,T3,T4态的极小点能量时，是否采用和T1态能量读取相同的办法，T2态也是读取迭代过程最后一步里的 Excited State 2 下方对应的 Total Energy, E(TD-HF/TD-DFT) 。
二，我在看有关重组能计算帖子时看见第二个图，我想问在E1到E2的垂直吸收的过程，计算E2能量，采用优化好的S0结构，输入关键词 # b3lyp/6-31g(d) td(triplets,root=1,nstates=5) 电荷自旋多重度为0，3.  读取Excited State 1 下方对应的 Total Energy, E(TD-HF/TD-DFT) ，E4能量在基于优化好的S0结构进行T1态优化，采用UDFT方法，电荷自旋多重度为0,3计算单点能读 E(SCF Done)。 最后计算E3能量时，采用优化好的T1态结构做TDDFT计算，其中电荷和自旋多重度设置为0，1.这样的方法是否可行？感觉E1到E2采用UDFT，E4到E3采用TDDFT好像不合理。
这点问题真有点糊涂，学到现在也不知道学到啥，谢谢老师们耐心看完。

sobereva

“我想问在读取T2,T3,T4态的极小点能量时，是否采用和T1态能量读取相同的办法”
是

td(triplets)结合0 3的计算没有任何意义

问问题的时候注意排版和条理，要问的问题如果较多的话，应当一条一条写明确，最好能编号，否则看着很累，也不好针对性地回复。

jiang012582

谢谢老师我下次注意，老师我改成下面两种方式您看哪种正确呀   一：1.E1取优化S0结构单点能 2.E2取当前优化好S0结构在，0，3下计算一次单点记为E2   3.基于优化好S0结构0，3下优化T1结构，计算单点记为E3.  4.E4为取T1结构，0，1计算单点记为E4能量。二. 1.E1取优化的S0结构的单点能，2，基于优化好S0结构做TDDFT计算，电荷自旋设置为0，1. 取第一步迭代输出的Excited state1 下面Total energy E(TD-HF/TD-DFT)记为E2.     
3.基于优化好S0结构0，1下加入TD优化T1结构取迭代最后一步输出Excited state1 下面Total energy E(TD-HF/TD-DFT)记为E3     4.取T1结构，做TDDFT计算电荷自旋设置为0，1.取第一步迭代输出的Excited state1 下面Total energy E(TD-HF/TD-DFT)记为E4.
请问老师我这里有没有错误，主要想E1-E2 E3-E4采用相同方法计算，和当前E2与E4取TDDFT或者UDFT下能量。谢谢老师了

sobereva

jiang012582 发表于 2019-5-8 10:27
谢谢老师我下次注意，老师我改成下面两种方式您看哪种正确呀   一：1.E1取优化S0结构单点能 2.E2取当前优化 ...

排版还是很乱，一 和 二 怎么都不换行

只需要算T1的情况， 用 一

jiang012582

忘记分段了谢谢老师