计算化学公社

标题: 求助：使用sobMECP程序搜索极小能量交叉点的疑惑 [打印本页]

作者
Author: myyang 时间: 2019-6-19 11:20
标题: 求助：使用sobMECP程序搜索极小能量交叉点的疑惑
我的计算的体系是一个双自由基，想用sobMECP程序sobMECP程序来搜索S=0 和S=1的极小能量交叉点。在做Input_Header_A/B文件是存在疑惑：

Input_Header_A:

%mem=6GB
%nproc=24
%chk=singlet.chk
#n UM062X/6-31G** force guess(read)

First State

0 1

Input_Header_B

%mem=6GB
%nproc=24
%chk=triplet.chk
#n UM062X/6-31G** force guess(read)

Second State

0 3

输入文件中 UM062X是否正确，看Sob老师的例子中都没用U？
guess(read) 需要变成 guess(read, mix)吗? 当我没加mix时，计算得到的MECP点我单独拿出来计算（#p UM062x/6-31g** guess=(mix,always)）发现S=0 和S=1 时的能量不相等

谢谢大家

作者
Author: Accelerator 时间: 2019-6-19 12:52
先自行做一次单点得到chk并改成相应的文件名之后计算就可以读取了

作者
Author: 小范范1989 时间: 2019-6-19 15:23
Ｓ０和Ｓ１势能面交叉点，不都是Ｃａｓｓｃｆ做吗？
ｓｏｂ－－ＭＥＣＰ合适吗？

作者
Author: zjxitcc 时间: 2019-6-19 15:30

小范范1989 发表于 2019-6-19 15:23
Ｓ０和Ｓ１势能面交叉点，不都是Ｃａｓｓｃｆ做吗？
ｓｏｂ－－ＭＥＣＰ合适吗？

题主的是S=0与S=1，不是S0, S1.

作者
Author: myyang 时间: 2019-6-19 17:42

Accelerator 发表于 2019-6-19 12:52
先自行做一次单点得到chk并改成相应的文件名之后计算就可以读取了

谢谢
我这样测试了， Step=1 的时候，就有4个YES(除能量没有收敛)，后续优化100步也未全部收敛，感觉奇怪，遗憾

作者
Author: sobereva 时间: 2019-6-19 21:36
对单重态、自旋多重度大于1的情况分别默认用R和U
自己写U除非是算对称破缺单重态，否则几乎总是多余的。我之前给的sobMECP例子里没有一个是算双自由基体系的

guess(read, mix)没有任何意义
算之间手动按照此文得到期望的对称破缺态的chk：
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
sobMECP跑的时候照常用guess=read就完了

作者
Author: 小范范1989 时间: 2019-6-20 08:23

zjxitcc 发表于 2019-6-19 15:30
题主的是S=0与S=1，不是S0, S1.

偶，抱歉，思维定势，看错了。

作者
Author: myyang 时间: 2019-6-20 10:47

sobereva 发表于 2019-6-19 21:36
对单重态、自旋多重度大于1的情况分别默认用R和U
自己写U除非是算对称破缺单重态，否则几乎总是多余的。我 ...

谢谢Sob 老师的指导！

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