计算化学公社

标题: 高斯CCSD(T)T1值在哪里？ [打印本页]

作者
Author: mmjjmmjj 时间: 2018-11-8 08:52
标题: 高斯CCSD(T)T1值在哪里？
请教下用高斯做CCSD(T)单点，结果log文件中并没有找到T1 diag，是不是得在输入行中加特殊关键词？
谢谢！

作者
Author: sobereva 时间: 2018-11-8 09:13
CCSD(T1diag)。就在CCSD能量下面

作者
Author: mmjjmmjj 时间: 2018-11-8 09:18

sobereva 发表于 2018-11-8 09:13
CCSD(T1diag)。就在CCSD能量下面

木有。。。全文搜索也找不到。。。
ccsd(t,maxcycle=200)/gen guess=read geom=allcheck\\Title Card Require
d\\1,4\N,0,0.,-0.000000003,0.94394699\Sc,0,-0.0000000008,1.65646249,-0
.10491298\Sc,0,-1.434538599,-0.8282312502,-0.10491298\Sc,0,1.434538599
7,-0.8282312488,-0.10491298\\Version=ES64L-G09RevD.01\HF=-2333.5652049
\MP2=-2333.9188728\MP3=-2333.9032714\MP4D=-2333.9201781\MP4DQ=-2333.90
76499\PUHF=-2333.567213\PMP2-0=-2333.9205963\PMP3-0=-2333.9045822\MP4S
DQ=-2333.9224005\CCSD=-2333.9401269\CCSD(T)=-2333.9853958\S2=3.816752\
S2-1=3.79747\S2A=3.751793\RMSD=2.426e-09\PG=C03V

作者
Author: sobereva 时间: 2018-11-8 09:39

mmjjmmjj 发表于 2018-11-8 09:18
木有。。。全文搜索也找不到。。。
ccsd(t,maxcycle=200)/gen guess=read geom=allcheck\\Title Card Re ...

你根本就没写t1diag

作者
Author: Scx 时间: 2019-4-11 17:06
楼主，您好！我想请问，您后来找到T1diag的数值在哪个位置了吗？我的是正常结束的，但没找到这个数值在哪里看? 您能帮我看下哪里有问题吗？谢谢您！
Test job not archived.
1\1\GINC-NODE4\SP\RCCSD(T)-FC\Gen\O1Sc1(1+)\\11-Apr-2019\0\
\# CCSD(T,T1DIAG)/gen guess=read geom=allcheck iop(5/13=1)\\Title Card
Required\\1,1\Sc,0,-1.1420153964,-0.74850298,0.\O,0,0.4833327664,-0.7
4850298,0.\\Version=ES64L-G09RevD.01\State=1-SG\HF=-834.4322526\MP2=-8
34.7386732\MP3=-834.6872956\MP4D=-834.7107393\MP4DQ=-834.7079468\MP4SD
Q=-834.7319099\CCSD=-834.7095281\CCSD(T)=-834.7321486\RMSD=1.454e-10\P
G=C*V [C*(O1Sc1)]\\@

作者
Author: sobereva 时间: 2019-4-11 17:14

Scx 发表于 2019-4-11 17:06
楼主，您好！我想请问，您后来找到T1diag的数值在哪个位置了吗？我的是正常结束的，但没找到这个数值在哪里 ...

CCSD迭代之后部分，甭在这里看

作者
Author: Scx 时间: 2019-4-11 17:30

sobereva 发表于 2019-4-11 17:14
CCSD迭代之后部分，甭在这里看

谢谢老师！我找到了
T1 Diagnostic = 0.04097511
Largest amplitude= 4.22D-02
Time for triples= 18.71 seconds.
T4(CCSD)= -0.30176260D-01

作者
Author: lastzealot 时间: 2021-1-9 16:58

sobereva 发表于 2019-4-11 17:14
CCSD迭代之后部分，甭在这里看

老师，请问这个t1diag需要作为关键词输入到命令行吗？
另外，这个值在什么范围是可以认为单参考态结果可信？
谢谢老师

作者
Author: wzkchem5 时间: 2021-1-9 18:51

lastzealot 发表于 2021-1-9 16:58
老师，请问这个t1diag需要作为关键词输入到命令行吗？
另外，这个值在什么范围是可以认为单参考态结果可 ...

需要，要把命令行的CCSD(T)换成CCSD(T,T1diag)
一般认为0.02以内是可信的，不过这个标准不能教条化地去使用，还需要结合其他标准以及化学常识

作者
Author: lastzealot 时间: 2021-1-9 19:29

wzkchem5 发表于 2021-1-9 18:51
需要，要把命令行的CCSD(T)换成CCSD(T,T1diag)
一般认为0.02以内是可信的，不过这个标准不能教条化地去 ...

首先感谢您的细心回答
您所说的其它标准包括哪些？在输出文件里能找到吗？
另外，没有单电子的d或f体系一般可以被认为是单参考体系吗？

作者
Author: sobereva 时间: 2021-1-9 21:07

lastzealot 发表于 2021-1-9 19:29
首先感谢您的细心回答
您所说的其它标准包括哪些？在输出文件里能找到吗？
另外，没有单电子的d或f体系 ...

经验直觉、TAE%、CASSCF最大权重、CISD中HF行列式的贡献、D1指数、B1指数等等。一些文章综合使用不同指数讨论静态相关强弱，如J. Chem. Theory Comput., 11, 5865 (2015)

显然不是那么简单。诸如四个氢摆成一个正方形，让键长距离不远不近，都是一个多参考特征极强的体系

作者
Author: lastzealot 时间: 2021-1-9 21:26

sobereva 发表于 2021-1-9 21:07
经验直觉、TAE%、CASSCF最大权重、CISD中HF行列式的贡献、D1指数、B1指数等等。一些文章综合使用不同指数 ...

谢谢老师
受益匪浅

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