计算化学公社

标题: 自旋极化单重态的输入问题 [打印本页]

作者
Author: 星斗如盘 时间: 2020-9-22 16:10
标题: 自旋极化单重态的输入问题
本帖最后由星斗如盘于 2020-9-22 19:21 编辑

看过sob的帖子谈谈片段组合波函数与自旋极化单重态http://sobereva.com/82，有以下两个问题：（1）guess和guess=(mix,always)有什么区别

（2）先计算三线态的单点能然后读取波函数作为初猜，寻找自旋极化单重态的过渡态和guess=mix有什么区别？那一个方法更好？

作者
Author: zjxitcc 时间: 2020-9-22 16:26
（1）always的解释官方就有http://gaussian.com/guess/
简单来说就是（就算要用，也）仅用于结构优化中，与算单点无关。

（2）guess=mix没有统一的定义，不同软件可以采取不同的做法，意思是同一个意思，具体做法可以不一样。
没有哪个更好之说。一般而言谁能量低，谁波函数稳定，谁就更好。

作者
Author: sobereva 时间: 2020-9-22 18:03
1 guess只是一个控制初猜的关键词，必须接具体选项才有意义，仔细看手册相应部分
guess(mix,always)在你提到的文章里，以及这里有解释
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

2 “先计算三线态的单点能然后读取波函数作为初猜来计算自旋极化单重态”语义不明，这种做法根本没法算单重态

作者
Author: 星斗如盘 时间: 2020-9-22 19:09

zjxitcc 发表于 2020-9-22 16:26
（1）always的解释官方就有http://gaussian.com/guess/
简单来说就是（就算要用，也）仅用于结构优化中， ...

在结构优化的中，如果只使用guess=mix不加always，在优化过程中有可能会从自旋极化单重态变成普通的单重态吗？

作者
Author: 星斗如盘 时间: 2020-9-22 19:20

sobereva 发表于 2020-9-22 18:03
1 guess只是一个控制初猜的关键词，必须接具体选项才有意义，仔细看手册相应部分
guess(mix,always)在你提 ...

先计算三线态的单点能然后读取波函数作为初猜寻找自旋极化单重态的过渡态

作者
Author: zjxitcc 时间: 2020-9-22 19:37

星斗如盘发表于 2020-9-22 19:09
在结构优化的中，如果只使用guess=mix不加always，在优化过程中有可能会从自旋极化单重态变成普通的单重 ...

若从共价键拉长的位置做结构优化，跑向距离缩短的平衡位置，这是一个从自旋极化单重态 -> 闭壳层单重态的过程，这种过程无论你加或不加mix或always，缩短一会儿都会收敛到闭壳层单重态。

作者
Author: zjxitcc 时间: 2020-9-22 19:38

星斗如盘发表于 2020-9-22 19:20
先计算三线态的单点能然后读取波函数作为初猜寻找自旋极化单重态的过渡态

这是一种有效的手段，但并不总是奏效。还要看波函数稳定性等其他因素。

作者
Author: sobereva 时间: 2020-9-22 22:04

星斗如盘发表于 2020-9-22 19:09
在结构优化的中，如果只使用guess=mix不加always，在优化过程中有可能会从自旋极化单重态变成普通的单重 ...

除非跑到的相应结构本来稳定的波函数就是闭壳层单重态，否则一般不会

作者
Author: sobereva 时间: 2020-9-22 22:05

星斗如盘发表于 2020-9-22 19:20
先计算三线态的单点能然后读取波函数作为初猜寻找自旋极化单重态的过渡态

你读取三重态波函数当初猜，会导致你之后的任务还是在算三重态，怎么可能得到自旋极化单重态，除非你像ORCA那样做spin flip，但Gaussian不支持

作者
Author: 星斗如盘 时间: 2020-9-23 16:05

sobereva 发表于 2020-9-22 22:05
你读取三重态波函数当初猜，会导致你之后的任务还是在算三重态，怎么可能得到自旋极化单重态，除非你像OR ...

读取三线态波函数作为初猜，但是新的任务自旋多重度设为1

作者
Author: sobereva 时间: 2020-9-24 12:06

星斗如盘发表于 2020-9-23 16:05
读取三线态波函数作为初猜，但是新的任务自旋多重度设为1

这没有意义。比如你用C4H8双自由基测试一下就发现行不通

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