基态为开壳层三重态，如何计算更高的三重激发态

cdw32156

我想对一个基态为开壳层三重态的三角烯分子，垂直激发到更高的三重激发态，计算吸收光谱，请问是这样计算吗：
第一步：opt ub3lyp/6-311+g(d,p) nosymm pop=full iop(9/40=4)
0 3
第二步：采取上一步优化好的结构 ub3lyp/6-311+g(d,p) td(Nstates=12,root=1) nosymm pop=full iop(9/40=4)
0 3

zjxitcc

zxmre20 发表于 2024-11-28 06:28
老师您好，我是在macos上使用homebrew编译的mokit， V1.2.5-with-option 版本
在运行您的例子时（三重态 ...

版本过低，请至少使用MOKIT-v1.2.6rc28。当然，用现在最新的v1.2.6 Release版更好。GitLab上可联网安装MOKIT，下载Linux预编译版MOKIT，下载MOKIT的Windows可执行程序等等多种选择，任选其一。

zxmre20

zjxitcc 发表于 2024-4-19 13:30
建议使用GAMESS的MRSF-TDDFT，从三重态RODFT出发计算T0, T1, T2, ... （几乎自旋纯态，肉眼不可察觉自旋污 ...

老师您好，我是在macos上使用homebrew编译的mokit， V1.2.5-with-option 版本
在运行您的例子时（三重态O2）第二步
$ fch2inp O2.fch -mrsf 便出现错误：
ERROR in subroutine fch2inp: wrong command line arguments!
Example 1 (R(O)HF, UHF, CAS): fch2inp a.fch
Example 2 (GVB)             : fch2inp a.fch -gvb [npair]
Example 3 (ROGVB)           : fch2inp a.fch -gvb [npair] -open [nopen]

请问这是版本导致的吗？ 谢谢老师

天天121

zjxitcc 发表于 2024-5-12 22:53
把inp文件中这几处改成如下样式
RUNTYP=OPTIMIZE
$TDDFT NSTATE=5 MULT=3 IROOT=2 $END

老师这个MRSF-TDDFT能做NTO吗

猫头鹰

zjxitcc 发表于 2024-5-12 22:53
把inp文件中这几处改成如下样式
RUNTYP=OPTIMIZE
$TDDFT NSTATE=5 MULT=3 IROOT=2 $END

老师好，想请教个问题：通过GAMESS的输出文件可以看激发态(T1)的电子分布以及轨道贡献吗

猫头鹰

zjxitcc 发表于 2024-5-12 22:53
把inp文件中这几处改成如下样式
RUNTYP=OPTIMIZE
$TDDFT NSTATE=5 MULT=3 IROOT=2 $END

谢谢老师的耐心指导，我尝试一下

zjxitcc

猫头鹰 发表于 2024-5-12 22:25
老师您说的正是我想表达的意思。

把inp文件中这几处改成如下样式
RUNTYP=OPTIMIZE
$TDDFT NSTATE=5 MULT=3 IROOT=2 $END
清理下上次计算的残留文件。把这个inp文件提交计算。优化过程中，往往会发生态顺序的变化，如果你的STATE 2变成了其他态，需要杀掉进程，提取最后一帧结构，重新做ROBHandHLYP计算、生成inp文件做结构优化。

猫头鹰

zjxitcc 发表于 2024-5-12 22:20
看不懂。能否换用更严谨的语言表示一下。我举个例子：假设有个体系基态是三重态，用MRSF-TDDFT算出来基态 ...

老师您说的正是我想表达的意思。

zjxitcc

猫头鹰 发表于 2024-5-12 21:56
老师好，换了一个小体系之后，这个问题解决了。想问一下如何得到最低能量的三重态激发态的几何构型吗呢？

看不懂。能否换用更严谨的语言表示一下。我举个例子：假设有个体系基态是三重态，用MRSF-TDDFT算出来基态也是三重态，它当然称为T0，我现在想得到T1态的极小点结构。

猫头鹰

zjxitcc 发表于 2024-5-9 23:16
假设你一直用的是BHandHLYP，且GAMESS SCF在5圈内收敛，随后在MRSF对角化步骤报错，那你可以尝试在$SYSTE ...

老师好，换了一个小体系之后，这个问题解决了。想问一下如何得到最低能量的三重态激发态的几何构型吗呢？

zjxitcc

猫头鹰 发表于 2024-5-9 23:12
明白了，多谢老师指点，比较傻，不知道BHANDHLYPH和BHHLYP是一回事

假设你一直用的是BHandHLYP，且GAMESS SCF在5圈内收敛，随后在MRSF对角化步骤报错，那你可以尝试在$SYSTEM中添加MEMDDI=2000。

猫头鹰

zjxitcc 发表于 2024-5-9 23:07
fch2inp a.inp -mrsf会在inp文件里指定用BHHLYP，不管gjf和fch文件里是啥泛函。这里需要第三次强调，SF计 ...

明白了，多谢老师指点，比较傻，不知道BHANDHLYPH和BHHLYP是一回事

zjxitcc

猫头鹰 发表于 2024-5-9 22:58
我gaussian用B3LYP的时候，转化为inp后，显示的是BHHLYP，以为这个是需要自己改成B3LYP。按mokit自动改的 ...

fch2inp a.inp -mrsf会在inp文件里指定用BHHLYP，不管gjf和fch文件里是啥泛函。这里需要第三次强调，SF计算最常用的泛函是BHandHLYP（这是gjf文件中的写法），不建议用别的。如果你想用别的泛函（不推荐），需要自己修改inp文件里的泛函名称。

猫头鹰

zjxitcc 发表于 2024-5-9 22:53
fch2inp a.fch -mrsf产生的inp文件中写的是DFTTYP=BHHLYP，这就是正确合理的，用起来也很方便，为啥你要 ...

我gaussian用B3LYP的时候，转化为inp后，显示的是BHHLYP，以为这个是需要自己改成B3LYP。按mokit自动改的话，那用B3LYP的时候在inp文件中也是BHHLYP吗？

zjxitcc

猫头鹰 发表于 2024-5-9 22:51
抱歉，老师，在原帖中已修改为图片。

fch2inp a.fch -mrsf产生的inp文件中写的是DFTTYP=BHHLYP，这就是正确合理的，用起来也很方便，为啥你要把它手动改成错的DFTTYP=BHANDHLYP我不理解。明明是很方便的运行一下，偏要给它改成报错的

“建议使用SF这类方法最常用泛函BHandHLYP”这句话的意思是你在Gaussian的gjf文件中要用BHandHLYP泛函。而GAMESS的inp文件中所有信息都是自动写好的。