计算化学公社

标题: UDFT和UTDDFT计算开壳三线态和激发态的问题 [打印本页]

作者
Author: sqpc 时间: 2018-10-28 21:59
标题: UDFT和UTDDFT计算开壳三线态和激发态的问题
各位老师同学们好！看了sob 老师的文章，但仍不知道gaussian计算开壳三线态的关键词是opt ub3lyp/6-31G** guess(mix,always) 0  3 吗？
进行吸收光谱的计算是 td(nstates=10) ub3lyp/6-31G** guess(mix,always) 0  1 吗?
进行T1到Tn的计算是td(triplets,nstates=10) ub3lyp/6-31G** guess(mix,always) 0  1吗？
还有一个问题就是计算单线态和三线态的双自由基能级差需要消除自旋污染吗？还有他们的绝对能量需要消除自选污染吗？
麻烦老师们指点一下。谢谢谢谢。

作者
Author: sobereva 时间: 2018-10-28 22:43
“开壳三线态”用词不当，三重态本来就是开壳层

如果你是指参考态是开壳层单重态，要算它的三重态激发态，你的关键词不对。只有参考态是闭壳层时triplet才有意义，这里明说了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

算吸收光谱那么些没错

DFT计算不需要管自旋污染

作者
Author: sqpc 时间: 2018-10-28 23:00

sobereva 发表于 2018-10-28 22:43
“开壳三线态”用词不当，三重态本来就是开壳层

如果你是指参考态是开壳层单重态，要算它的三重态激发态 ...

谢谢老师的回复，1.有的双自由基体系它的基态是三线态的双自由基，请问这种应该用什么关键词算。2.还有从T1到Tn的计算是td(triplets,nstates=10) ub3lyp/6-31G** guess(mix,always) 0 1这么写对吗？麻烦老师您了，还有3.DFT为什么不需要管自旋污染呀?新手小白，恳请老师多多指点

作者
Author: sobereva 时间: 2018-10-28 23:26

sqpc 发表于 2018-10-28 23:00
谢谢老师的回复，1.有的双自由基体系它的基态是三线态的双自由基，请问这种应该用什么关键词算。2.还有从 ...

1 三重态一般不叫它双自由基。写0 3 TD(nstates=x)

2 我已经说了，参考态是开壳层时用triplets毫无意义，写了也白写

3 DFT引入轨道波函数的目的和HF完全不同，用KS-DFT波函数是否很满足S^2算符本征函数这点来衡量结果既没理论意义，也没什么实际意义

作者
Author: sqpc 时间: 2018-10-29 15:04

sobereva 发表于 2018-10-28 23:26
1 三重态一般不叫它双自由基。写0 3 TD(nstates=x)

2 我已经说了，参考态是开壳层时用triplets毫无意 ...

谢谢老师指点

作者
Author: sqpc 时间: 2018-10-29 15:14

sobereva 发表于 2018-10-28 23:26
1 三重态一般不叫它双自由基。写0 3 TD(nstates=x)

2 我已经说了，参考态是开壳层时用triplets毫无意 ...

老师，对于您的第一个回复，还是不是十分明确只用0 3td不用guess关键词吗？还有此时不需要加opt吗？因为我想得到最低三线态能量。如果我只想得到三线态的最低能量用uDFT关键词opt ub3lyp/6-31G** guess(mix,always) 0 3不行吗？谢谢老师您了

作者
Author: sobereva 时间: 2018-10-29 16:30

sqpc 发表于 2018-10-29 15:14
老师，对于您的第一个回复，还是不是十分明确只用0 3td不用guess关键词吗？还有此时不需要加opt吗？因为 ...

guess(mix,always)是结合0 1算开壳层单重态的，给0 3用根本毫无意义
仔细看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

0 1还是0 3，取决于基态是什么自旋多重度

作者
Author: sqpc 时间: 2018-10-29 16:42

sobereva 发表于 2018-10-29 16:30
guess(mix,always)是结合0 1算开壳层单重态的，给0 3用根本毫无意义
仔细看
谈谈片段组合波函数与自旋 ...

已经仔细看啦，就是说对于基态是双自由基体系他们的三线态的计算就和闭壳层一样吧，只不过都用UDFT和UTDDFT,而且不用加triplets，只改成0 3 就行。老师我这么理解对吗?

作者
Author: liyuanhe211 时间: 2018-10-29 17:17

sqpc 发表于 2018-10-29 16:42
已经仔细看啦，就是说对于基态是双自由基体系他们的三线态的计算就和闭壳层一样吧，只不过都用UDFT和UTDD ...

说阅读，“仔细看了”明显没读懂、没理解；说表达，写了三句话两句是病句、剩下一句读着贼费劲。与其抽象的说不明白，还不如把你要算的具体体系、和要算的物理量发出来，把你想写的关键词写全，然后让人修改就完了。

作者
Author: sqpc 时间: 2018-10-29 20:57

liyuanhe211 发表于 2018-10-29 17:17
说阅读，“仔细看了”明显没读懂、没理解；说表达，写了三句话两句是病句、剩下一句读着贼费劲。与其抽象 ...

我在发帖的时候就已经写了啊，语文不好，抱歉污了您的眼

作者
Author: sqpc 时间: 2018-10-29 21:03

sobereva 发表于 2018-10-29 16:30
guess(mix,always)是结合0 1算开壳层单重态的，给0 3用根本毫无意义
仔细看
谈谈片段组合波函数与自旋 ...

老师那我要计算双自由基体系的最低三线态的关键词是应该opt ub3lyp/6-31G(d) 0 3 是吧
计算双自由基体系的T1 到Tn的激发能关键词是td (nstates=n) ub3lyp/6-31G（d）0 3?

作者
Author: sobereva 时间: 2018-10-29 21:32

sqpc 发表于 2018-10-29 21:03
老师那我要计算双自由基体系的最低三线态的关键词是应该opt ub3lyp/6-31G(d) 0 3 是吧
计算双自由基体系 ...

1 是。没必要自己特意写U
2 显然不是。基态是单重态，显然不能写0 3，而应当guess=mix TD UB3LYP/6-31G* 0 1

作者
Author: sqpc 时间: 2018-10-29 22:12

sobereva 发表于 2018-10-29 21:32
1 是。没必要自己特意写U
2 显然不是。基态是单重态，显然不能写0 3，而应当guess=mix TD UB3LYP/6-31G* ...

谢谢老师！

作者
Author: sqpc 时间: 2018-10-29 22:17

sobereva 发表于 2018-10-29 21:32
1 是。没必要自己特意写U
2 显然不是。基态是单重态，显然不能写0 3，而应当guess=mix TD UB3LYP/6-31G* ...

老师我要算的是T1到Tn的激发能呀，怎么感觉您写的是S0到SN的呢？真的麻烦老师了，基础比较薄弱，太感谢老师您了

作者
Author: sobereva 时间: 2018-10-29 23:45

sqpc 发表于 2018-10-29 22:17
老师我要算的是T1到Tn的激发能呀，怎么感觉您写的是S0到SN的呢？真的麻烦老师了，基础比较薄弱，太感谢老 ...

参考态是开壳层的时候你没法指定算什么自旋多重度的激发态，只能根据算出来的激发态的<S**2>判断是几重态，把其中能算得上是三重态的都取出来

作者
Author: sqpc 时间: 2018-10-30 17:25
好的，谢谢老师！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)