计算化学公社

标题: orca计算自旋极化单重态如何表示UCCSD（T） [打印本页]

作者
Author: BinWang 时间: 2022-11-16 16:01
标题: orca计算自旋极化单重态如何表示UCCSD（T）
本帖最后由 BinWang 于 2022-11-16 16:09 编辑

各位老师：
   大家好！我计算单重态双自由基，用高斯先产生一个UHF的初猜，并且用opt=stable得到稳定的波函数，然后根据sob老师的帖子http://sobereva.com/517将高斯的输出波函数喂给orca。我发现orca里面不允许输入UCCSD（T）。请问写成UHF CCSD（T）是一个意思吗？
   同样的问题，如果我想用orca里面的双杂化泛函，比如高斯里面的ub2blyp在orca里面是uks b2plyp ?
   我已经用高斯里面的uccsd（t）算过的体系，想用orca里面的uccsd（t）计算，直接调用高斯UCCSD（T）的波函数作为初猜。该怎么输入？如果直接输入UHF CCSD（T）感觉浪费了，直接写CCSD（T）则会优化到闭壳层。

作者
Author: wzkchem5 时间: 2022-11-16 17:45
对，UHF CCSD(T)就是UCCSD(T)的意思，UKS B2PLYP就是UB2PLYP
据我所知应该没有办法把高斯的CCSD(T)波函数传给orca，只能传HF部分。
此外，不要把高斯产生初猜传给orca当作标准计算流程，大多数情况下用orca直接找稳定的波函数更快。orca又不是没有稳定性分析功能

作者
Author: BinWang 时间: 2022-11-16 18:17

wzkchem5 发表于 2022-11-16 17:45
对，UHF CCSD(T)就是UCCSD(T)的意思，UKS B2PLYP就是UB2PLYP
据我所知应该没有办法把高斯的CCSD(T)波函数 ...

非常感谢老师的回复

作者
Author: qzxchem 时间: 2023-6-20 20:51
您好，想请问一下，UCCSSD(T) 计算OSS的输入应该怎么写呢？方便贴一下嘛？看到一些帖子上写加上flipspin [翻转的原子序号]

作者
Author: wzkchem5 时间: 2023-6-20 21:18

qzxchem 发表于 2023-6-20 13:51
您好，想请问一下，UCCSSD(T) 计算OSS的输入应该怎么写呢？方便贴一下嘛？看到一些帖子上写加上flipspin [ ...

在orca手册里搜flipspin就知道怎么写了。
建议先做flipspin的HF计算，计算完成且检查发现波函数确实正确后，再读这个波函数做UCCSD(T)

作者
Author: qzxchem 时间: 2023-6-20 21:32

wzkchem5 发表于 2023-6-20 21:18
在orca手册里搜flipspin就知道怎么写了。
建议先做flipspin的HF计算，计算完成且检查发现波函数确实正确 ...

谢谢您！老师，我去做一下！

作者
Author: qzxchem 时间: 2023-6-21 17:34
本帖最后由 qzxchem 于 2023-6-22 16:39 编辑

wzkchem5 发表于 2023-6-20 21:18
在orca手册里搜flipspin就知道怎么写了。
建议先做flipspin的HF计算，计算完成且检查发现波函数确实正确 ...

老师您好，根据您给我的指点我先做了flipspin的HF
正常结束计算后，在对称性破缺部分显示如下
请问老师，是否可以从这里的S**2判断已经收敛到了OSS正确的波函数了呢？或者说，我该如何判断已经收敛到了正确的波函数了呢？谢谢老师！

作者
Author: qzxchem 时间: 2023-6-21 21:24

wzkchem5 发表于 2022-11-16 17:45
对，UHF CCSD(T)就是UCCSD(T)的意思，UKS B2PLYP就是UB2PLYP
据我所知应该没有办法把高斯的CCSD(T)波函数 ...

老师您好，还想请问一下，使用基组外推的方法CCSD(T)/CBS (cc-pVTZ->QZ extrapolation)，orca在做动外推功能是否能保证对每个基组的计算都是对对称破缺态的呀？谢谢您！

作者
Author: wzkchem5 时间: 2023-6-21 21:58

qzxchem 发表于 2023-6-21 10:34
老师您好，根据您给我的指点我先做了flipspin的HF，orca版本5.0.4输入如下
正常结束计算后，在对称性破 ...

如果已知这个体系只可能有一种自旋破缺的波函数，那么可以。否则不行，必须结合自旋布居、自旋密度等判断，因为有的体系有不止一个自旋破缺波函数，其中只有一个是正确的

作者
Author: wzkchem5 时间: 2023-6-21 22:01

qzxchem 发表于 2023-6-21 14:24
老师您好，还想请问一下，使用基组外推的方法CCSD(T)/CBS (cc-pVTZ->QZ extrapolation)，orca在做动外推 ...

应该是大概率可以但不能严格保证。首先其中的TZ计算应该读你提前准备好的对称性破缺态波函数作为初猜；其次QZ计算会自动读取TZ波函数作为初猜。读对称性破缺波函数可以使得收敛的波函数有较大概率是对称性破缺，但当初猜的基组和实际计算基组不符时，这个概率不到100%。所以还是需要检查自旋布居等结果才能完全确定

作者
Author: qzxchem 时间: 2023-6-21 22:32

wzkchem5 发表于 2023-6-21 22:01
应该是大概率可以但不能严格保证。首先其中的TZ计算应该读你提前准备好的对称性破缺态波函数作为初猜；其 ...

感谢老师细致的解答，还有三个问题想要请教您。
1⃣️我的体系超级小，所以我直接在QZ计算波函数可以嘛？
2⃣️在orca的输出中怎么检查自旋布居分布呢？我在orca的输出中没有查到有spin population呢；
3⃣️我是现在高斯中用UB3LYP stable=opt方法得到稳定波函数的结构，通过自旋布局确认了结构是OSS的,S**2为0.7，想用orca CCSD(T)方法计算高精度的单点能量，其算出的S**2=1.1，老师这两个应该选用哪个数值呢，非常感谢老师！

作者
Author: wzkchem5 时间: 2023-6-22 00:17

qzxchem 发表于 2023-6-21 15:32
感谢老师细致的解答，还有三个问题想要请教您。
1⃣️我的体系超级小，所以我直接在QZ计算波 ...

（1）ORCA应该不支持TZ、QZ分别读取两个不同的初猜。假设你5Z算不动的话，用TZ准备初猜、然后手动或写脚本检查QZ有没有收敛对，我觉得是最优方法了。此外也可以考虑做F12计算，也能得到CBS结果但无需外推。
（2）你的输入文件不能写nopop。在此前提下，在输出文件里搜大写SPIN
（3）两者各有各的意义，不能说一定选用哪个，得看你报<S**2>值的目的（注意不能漏写"<"和">"，单写S**2是错误的写法）

作者
Author: qzxchem 时间: 2023-6-22 16:42

wzkchem5 发表于 2023-6-22 00:17
（1）ORCA应该不支持TZ、QZ分别读取两个不同的初猜。假设你5Z算不动的话，用TZ准备初猜、然后手动或写脚 ...

非常感谢老师耐心细致的讲解！谢谢您啦！

作者
Author: qzxchem 时间: 2023-6-22 19:25
本帖最后由 qzxchem 于 2023-10-9 10:55 编辑

wzkchem5 发表于 2023-6-22 00:17
（1）ORCA应该不支持TZ、QZ分别读取两个不同的初猜。假设你5Z算不动的话，用TZ准备初猜、然后手动或写脚 ...

老师您好，感谢您的指点～
1.我现在先在UHF方法下，使用flipspin，检查输出中的SPIN，得到了正确的波函数结果；(看MULLIKEN ATOMIC CHARGES AND SPIN POPULATIONS 部分的SPIN POPULATIONS是否分布正确)
2.使用UCCSDT方法，TZ基组读UHF下的波函数计算，在这一步中的输入如下：
! UHF CCSD(T) cc-pVTZ tightSCF MORead
%moinp "hf.gbw"
%maxcore xxx
%pal nprocs xxx end
*xyz 0 1 #(OSS开壳层)

老师，我想请问的问题是
1⃣️在2.的输入中是不是不需要再使用flipspin了呀？只需要读取UHF计算下的波函数即可～
2⃣️老师，看到手册中会加%SCF GuessMode xxxx end 的关键词，这部分对%SCF 的设置是不是如果没有特殊的计算需要就不需要加呀？谢谢您！

作者
Author: wzkchem5 时间: 2023-6-22 19:59

qzxchem 发表于 2023-6-22 12:25
老师您好，感谢您的指点～
1.我现在先在UHF方法下，使用flipspin，检查输出中的SPIN，得到了正确的波函 ...

（1）对，因为flipspin的作用在于准备自旋破缺的初猜，如果初猜本来就是自旋破缺的，普通SCF算法就足以让这个波函数完全收敛，不需要也不应该再flip一次
（2）一般不用写

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