计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 用MRCC软件对一些耦合簇方法的测试

用MRCC软件对一些耦合簇方法的测试

查看数: 18855 | 评论数: 26 | 收藏 Add to favorites 4

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

Accelerator

发布时间: 2017-7-22 09:08

正文摘要:

MRCC是由匈牙利布达佩斯大学开发的主攻高精度量化计算的软件包，其亮点在于可用于任意阶耦合簇（包括其各种微扰近似）和组态相互作用，并有对应的解析梯度。笔者将对该软件的耦合簇功能进行试用。限于笔者的电脑 ...

回复 Reply

ldatea 发表于 Post on 2020-8-19 22:02:17

本帖最后由 ldatea 于 2020-8-19 22:28 编辑

zjxitcc 发表于 2020-8-19 10:00
不好意思我想问得详细点儿，是哪部分没重复出楼主的计算结果？某个电子能量？某个相对能量？打开表格没看 ...

我用的计算方法和基组没有严格对应，现在直接给出不能重复的地方。先前的表格把基组从aug-cc-pvtz改成jun-cc-pvqz，把冻芯去掉，对于水和氢氧根也没啥影响（单点能差398kcal/mol附近）。这里冻芯对反应能的结果大一点。（为了一眼就能看懂，没贴额外信息）

1.楼主的图中，CCSD（T）和CCSDT(Q)相差无几
2.第 4 张条形图中M06-2X/aug-cc-pvtz//b3lyp/6-31g*与CCSDT(Q)/aug-cc-pvtz//M06-2X/def2-tzvp相差12.1kcal/mol

3.我算出来的M06-2X/aug-cc-pvtz//b3lyp/6-31g*与CCSD(T)/aug-cc-pvtz//M06-2X/def2-tzvp相差1.1kcal/mol
这三条里无论如何，至少有一条是错的。

喵星大佬 发表于 Post on 2020-8-19 15:35:31

ldatea 发表于 2020-8-18 22:58
不知道为什么，我用G16无法重复你的结果。我使用的能量参考值时CCSD(T)/jun-cc-pvqz下的能量，和CCSDT(Q)差 ...

这很容易理解，因为直线往任意方向弯的受力是平衡的，所以用受力优化怎么优化都是直线，这也是几何优化之后有虚频，消虚频的时候为什么要手动破坏对称性。对称的结构一定是能量极值点，但未必是极小值点。

zjxitcc 发表于 Post on 2020-8-19 10:00:38

ldatea 发表于 2020-8-18 22:58
不知道为什么，我用G16无法重复你的结果。我使用的能量参考值时CCSD(T)/jun-cc-pvqz下的能量，和CCSDT(Q)差 ...

不好意思我想问得详细点儿，是哪部分没重复出楼主的计算结果？某个电子能量？某个相对能量？打开表格没看出所说的问题在哪儿

ldatea 发表于 Post on 2020-8-18 22:58:13

本帖最后由 ldatea 于 2020-8-18 23:52 编辑

不知道为什么，我用G16无法重复你的结果。我使用的能量参考值时CCSD(T)/jun-cc-pvqz下的能量，和CCSDT(Q)差不了多少。
我的结论是:
此体系中，
1.  只要结构是接近的，不同的几何优化结果对最终能量的影响很小，大概就是0点几kcal/mol
2.  算阴离子体系能量不加弥散函数结果会很烂
这两条通常都是成立的

如果你的结果是对的，就说明  小基组几何优化+大基组单点在这个体系会严重不符。虽然不是说一定不可能，但在这么一个很正常的体系里出这种事会非常非常怪。

另外还发现一个很神奇的事，如果把水摆成直线几何优化（键拉长的情况下，键压缩的情况我还没试过），即便用nosymm，（没算Hessian时，没尝试过算Hessian的情况）在M06-2X/def2-tzvp下总是会优化出一个直线的构型。

能量差.xlsx

11.97 KB, 阅读权限: 10, 下载次数 Times of downloads: 2

lidanhui2 发表于 Post on 2018-2-9 16:37:04

您好，能发一下MRCC的安装软件吗？940469319@qq.com

Warm_Cloud 发表于 Post on 2017-7-24 22:31:18

MRCC慢在了积分变换以及排序上，它的CCSD(T)的算法本身更快。高斯的积分变化出奇的快，CCSD(T)的话，molpro更稳健一些。

scf 发表于 Post on 2017-7-24 01:02:46

beefly 发表于 2017-7-22 09:15
CCSDT和CCSDT(Q)可以用cfour 2.0做。人写的代码比计算机产生的代码效率高

展开说说？automatical implementation的挑战在哪里？

我本是个娃娃 发表于 Post on 2017-7-23 08:15:35

sobereva 发表于 2017-7-22 22:33
你直接下载预编译好的包，把路径加入到$PATH里直接就能用

感谢大博士

sobereva 发表于 Post on 2017-7-22 23:10:28

beefly 发表于 2017-7-22 22:34
可能是版本问题，需要用2014年的2.0beta版

是我搞错了，之前机子里有cfour v1版，明明编译了v2 beta，结果不知怎么回事跑的是v1。删了v1版现在CCSDTQ能正常计算了。谢谢提示。

beefly 发表于 Post on 2017-7-22 22:34:21

可能是版本问题，需要用2014年的2.0beta版

sobereva 发表于 Post on 2017-7-22 22:33:53

我本是个娃娃发表于 2017-7-22 09:27
MRCC能做极化率和超极化率计算吗？

另外能把安装的过程写一点吗？

你直接下载预编译好的包，把路径加入到$PATH里直接就能用

sobereva 发表于 Post on 2017-7-22 22:12:04

beefly 发表于 2017-7-22 17:32
需要设定CC模块：CC_PROGRAM=NCC
现在只能算闭壳层单点能
http://slater.chemie.uni-mainz.de/cfour/in ...

设了NCC后提示Option NCC for keyword CC_PROGRAM not known. Default used.，结果自动用的还是VCC，最终CCSD和CCSDTQ结果相同，输出文件附上了
官网上看貌似NCC在公开版本里没有：http://slater.chemie.uni-mainz.de/cfour/index.php?n=Main.CCModules

CCSD

36.77 KB, 下载次数 Times of downloads: 5

CCSDTQ

36.69 KB, 下载次数 Times of downloads: 5

beefly 发表于 Post on 2017-7-22 17:32:22

sobereva 发表于 2017-7-22 10:26
你说的2.0是非公开版本么？
用公开版本，写CALC=CCSDT(Q)时候发现结果和CCSD一样。貌似公开版本最高就 ...

需要设定CC模块：CC_PROGRAM=NCC
现在只能算闭壳层单点能
http://slater.chemie.uni-mainz.d ... hp?n=CFOURV2.00beta

jiangning198511 发表于 Post on 2017-7-22 15:42:55

sobereva 发表于 2017-7-22 10:26
你说的2.0是非公开版本么？
用公开版本，写CALC=CCSDT(Q)时候发现结果和CCSD一样。貌似公开版本最高就 ...

CFOUR貌似在今年秋天之前发布新版本

A new release of CFOUR should happen no later than autumn 2017

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register