计算化学公社

标题: 关于S1与S2态的交叉 [打印本页]

作者
Author: 天天121 时间: 2022-4-9 08:46
标题: 关于S1与S2态的交叉
读过Sob老师写的关于搜索不同多重度之间交叉点的博文，现在遇到一个问题，就是我在优化S1（nstate=3,root=1）和S2（nstate=4,root=2）时发现S1的能量比S2更高。查看输出文件发现，S1态优化的过程中并没有发现交叉，但是S2的优化过程中却发现随着优化步数的增加S1和S2发生了交叉。我私以为这是由于圆锥交差点多维度的问题导致S1态优化的时候没有经过那个点，但却是在S2态优化的时候经过了。我想找到那个圆锥交叉点请问有什么好的方法吗谢谢各位了！

作者
Author: kay 时间: 2022-4-9 09:21
gv中读取不准，具体看sob老师帖子http://bbs.keinsci.com/thread-13450-1-1.html

作者
Author: 喵星大佬 时间: 2022-4-9 09:24
本帖最后由喵星大佬于 2022-4-9 09:25 编辑

看跃迁的轨道，交叉并没有导致跟踪错，优化并不需要经过CI点，只要经过交叉超曲面的任意一个点就行
用TDDFT可以找交叉点，用sobMECP+瓜丝，orca都可以

作者
Author: 天天121 时间: 2022-4-9 09:26

kay 发表于 2022-4-9 09:21
gv中读取不准，具体看sob老师帖子http://bbs.keinsci.com/thread-13450-1-1.html

老师您好，我觉得这不是读不读的准的问题老师(Sob老师这个博文我也看过，还是S1态更高)，确实是出现S1和S2的交叉点，但是在S1的优化过程中却并未发现，在S2的优化过程中却发现了

作者
Author: 天天121 时间: 2022-4-9 09:30

喵星大佬发表于 2022-4-9 09:24
看跃迁的轨道，交叉并没有导致跟踪错，优化并不需要经过CI点，只要经过交叉超曲面的任意一个点就行
用TDDF ...

好的谢谢老师那为什么我在优化S1的时候没有发现有交叉，只是在S2优化过程中才发现呢？还有就是平时大多数时候我只需要优化S1态，我现在有点不放心仅凭S1态输出文件是否能判断出有没有交叉了。因为我确实不常用到S2态。

作者
Author: 天天121 时间: 2022-4-9 09:34

喵星大佬发表于 2022-4-9 09:24
看跃迁的轨道，交叉并没有导致跟踪错，优化并不需要经过CI点，只要经过交叉超曲面的任意一个点就行
用TDDF ...

还有就是SobMECP可以用于相同多重度交叉点的搜索吗

作者
Author: 喵星大佬 时间: 2022-4-9 09:39

天天121 发表于 2022-4-9 09:30
好的谢谢老师那为什么我在优化S1的时候没有发现有交叉，只是在S2优化过程中才发现呢？还有就是平时大多 ...

S1/S2是相对于基态结构来说的，在S1极小点没经过交叉超曲面，而从基态极小点去S2极小点要经过那个交叉超曲面。换句话说S1和S2极小点结构都不一样，激发能排序不一样不是很正常吗。。。。

作者
Author: 天天121 时间: 2022-4-9 10:05

喵星大佬发表于 2022-4-9 09:39
S1/S2是相对于基态结构来说的，在S1极小点没经过交叉超曲面，而从基态极小点去S2极小点要经过那个交叉超 ...

但是老师从S1态的输出文件来看 SI态就始终是S1态而S2态给出的结果却是原本的S2态最后变成了最低绝热态也就是卡沙规则里说的那个最低绝热态荧光从哪个结构发射呢到底

作者
Author: 天天121 时间: 2022-4-9 10:47
整理了一下如下图所示，那么到底从哪个态发射荧光呢

作者
Author: kay 时间: 2022-4-9 11:10

天天121 发表于 2022-4-9 10:47
整理了一下如下图所示，那么到底从哪个态发射荧光呢

我个人理解，你优化S1走的绝热路径，优化S2走的透热路径，可以根据SDNTO跟踪一下路径判断吧。但无论走什么路径，应该是绝热的S1态发光。你图中S1的绝热能量为什么比S2的要高？你输出文件中不是低么？交叉点应该往前画，才应该符合你的这个情况吧。

作者
Author: Freeman 时间: 2022-4-9 11:52
楼主看过轨道没有？有没有可能只是152和153的顺序变了？如果是的话，就不涉及交叉；否则用g16+sobmecp或者orca都可以找交叉点，见orca手册189页。

作者
Author: 天天121 时间: 2022-4-9 12:49

kay 发表于 2022-4-9 11:10
我个人理解，你优化S1走的绝热路径，优化S2走的透热路径，可以根据SDNTO跟踪一下路径判断吧。但无论走什 ...

老师您好图片1和2给出的只是优化某一个步数的能量最后能量就是S1比较高因为从S1结构优化的输出文件中并未发现有任何交叉只是在S2的优化文件中看见了交叉所以交叉点推测在后面

作者
Author: 天天121 时间: 2022-4-9 13:02

Freeman 发表于 2022-4-9 11:52
楼主看过轨道没有？有没有可能只是152和153的顺序变了？如果是的话，就不涉及交叉；否则用g16+sobmecp或者o ...

老师您好我看了一下轨道 S1态对不管是从OPT-S1.fch看还是从OPT-S2.fch看都对应的是一个态，S2态也一样，所以应该是出现了交叉点我一会用sobMECP加高斯试一下

作者
Author: Freeman 时间: 2022-4-9 13:09

天天121 发表于 2022-4-9 13:02
老师您好我看了一下轨道 S1态对不管是从OPT-S1.fch看还是从OPT-S2.fch看都对应的是一个态，S2态也一样， ...

如果两个态都用tddft算的话，我还是建议用orca，因为只要一个输入文件就够了，而且比高*快一点。

作者
Author: 天天121 时间: 2022-4-9 13:36

Freeman 发表于 2022-4-9 13:09
如果两个态都用tddft算的话，我还是建议用orca，因为只要一个输入文件就够了，而且比高*快一点。

好的老师十分感谢

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