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

zjxitcc · 发表于 Post on 2024-5-9 13:23:41

本帖最后由 zjxitcc 于 2024-5-9 22:40 编辑

猫头鹰发表于 2024-5-9 11:10
实验室的服务器是126G的内存。有两个问题：
1、在输入文件中已经修改MWORDS， ...

（1）继续加大内存。
（2）不建议使用B3LYP，建议使用SF这类方法最常用泛函BHandHLYP。这句话的意思是你在Gaussian和GAMESS中都要用BHandHLYP，随后fch2inp -mrsf会自动在inp文件里写好这个泛函，不用改。你不能在Gaussian中用B3LYP，然后在GAMESS中用BHandHLYP，这样传过去的轨道效果很差。
（3）若体系为中性，基组强烈不建议使用6-311+g(d,p)，建议改用def2TZVP。由于体系较大，可以重要的原子用def2TZVP，不重要的原子用6-31G(d,p)。自定义基组在gjf文件中很容易操作。基组设置要合理，不要偷懒或乱设置。
（4）注意监测计算进度，比如RODFT几圈收敛，时间是花在SCF上，还是花在Davidson迭代对角化上，要有更为明确的概念，而非“计算时收敛了5个小时然后进入MRSF部分出错”。模糊用语，无法对症下药，没有人能帮忙。

猫头鹰 · 发表于 Post on 2024-5-9 22:32:40

本帖最后由猫头鹰于 2024-5-9 22:48 编辑

zjxitcc 发表于 2024-5-9 13：23
（1）继续加大内存。
（2）不建议使用B3LYP，建议使用SF这类方法最常用泛函BHandHLYP。
（3）若体系为 ...

在了解完MWORDS的涵义之后，我用了44个核，MWORDS=3000，然后我用exetyp=check之后是通过，然而在计算时在同一位置又出错了。

另外，我用BHANDHLYP【基组是混合基组，ROBHANDHLYP/GEN......】用gaussian算完之后转化，计算时出错，以下是inp文件和log文件。

zjxitcc · 发表于 Post on 2024-5-9 22:41:41

猫头鹰发表于 2024-5-9 22:32
在了解完MWORDS的涵义之后，我用了44个核，MWORDS=3000，然后我用exetyp=check之后是通过，然而在计算时 ...

看不出来。至少要有fch文件、inp文件或gms文件，最起码也要有报错信息的截图（不要复制粘贴文字，空格没对齐，会乱）

猫头鹰 · 发表于 Post on 2024-5-9 22:51:22

zjxitcc 发表于 2024-5-9 22:41
看不出来。至少要有fch文件、inp文件或gms文件，最起码也要有报错信息的截图（不要复制粘贴文字，空格没 ...

抱歉，老师，在原帖中已修改为图片。

zjxitcc · 发表于 Post on 2024-5-9 22:53:36

本帖最后由 zjxitcc 于 2024-5-9 22:55 编辑

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

fch2inp a.fch -mrsf产生的inp文件中写的是DFTTYP=BHHLYP，这就是正确合理的，用起来也很方便，为啥你要把它手动改成错的DFTTYP=BHANDHLYP

我不理解。明明是很方便的运行一下，偏要给它改成报错的

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

猫头鹰 · 发表于 Post on 2024-5-9 22:58:57

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

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

zjxitcc · 发表于 Post on 2024-5-9 23:07:53

本帖最后由 zjxitcc 于 2024-5-9 23:10 编辑

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

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

猫头鹰 · 发表于 Post on 2024-5-9 23:12:23

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

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

zjxitcc · 发表于 Post on 2024-5-9 23:16:04

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

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

猫头鹰 · 发表于 Post on 2024-5-12 21:56:33

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

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

zjxitcc · 发表于 Post on 2024-5-12 22:20:41

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

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

猫头鹰 · 发表于 Post on 2024-5-12 22:25:34

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

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

zjxitcc · 发表于 Post on 2024-5-12 22:53:28

本帖最后由 zjxitcc 于 2024-5-13 10:34 编辑

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

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

猫头鹰 · 发表于 Post on 2024-5-12 23:06:02

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

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

猫头鹰 · 发表于 Post on 2024-5-16 16:10:24

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 基态为开壳层三重态，如何计算更高的三重激发态