如何用高斯做CCSD（T）单点能计算？

lei234 · 发表于 Post on 2019-5-2 03:22:29

最近有个体系，不算大，不到60个原子，元素有C,N,O,H，结合了几个水分子。用M06-2X/6-31++G(d,p)优化读的能量，审稿人让做一下CC计算，（ I highlyrecommend to run CC calculations (as SP on the DFT structures) to have muchmore definitive conclusions.），还提到相同的体系别人已经用CCSD(T)计算过了。我打算用CCSD（T）做一下单点计算，对比一下能量。请问关键字应该怎么写？还有：
1.基组怎么选？我看别人用的是aug-cc-pVDZ
2.CCSD和CCSD(T)的区别？T是指激发三重态吗？我的体系都是基态，没有激发态。
3.内存应该设置多大？
如果计算单点能的话，直接CCSD(T)/aug-cc-pVDZ 就可以吗？

期待大神回复~~~

Lluo · 发表于 Post on 2023-10-25 11:21:51

zjxitcc 发表于 2023-9-21 17:42
CCSD(T)是先算完CCSD，在CCSD波函数基础上，以微扰的方式考虑三激发。t1和t2振幅在CCSD中是迭代求解的（ ...

好的明白了，谢谢老师

zjxitcc · 发表于 Post on 2023-9-21 17:42:41

本帖最后由 zjxitcc 于 2023-9-21 17:43 编辑

Lluo 发表于 2023-9-21 17:36
所以楼主知道CCSD和CCSD(T)的区别了吗？最近我看文献也疑惑这个，在gaussian的计算方法选择里有CCSD，但是 ...

CCSD(T)是先算完CCSD，在CCSD波函数基础上，以微扰的方式考虑三激发。t1和t2振幅在CCSD中是迭代求解的（给一组初猜，迭代至收敛），t3振幅在(T)中是一次求解、无迭代的。
如何设置：就是不要把事情想复杂，直接写比如CCSD(T)/cc-pVTZ。

Lluo · 发表于 Post on 2023-9-21 17:36:48

所以楼主知道CCSD和CCSD(T)的区别了吗？最近我看文献也疑惑这个，在gaussian的计算方法选择里有CCSD，但是没明白这个（T）的含义，和如何设置？

MiaoZuo · 发表于 Post on 2023-2-27 11:21:36

mfdsrax2 发表于 2023-2-13 20:23
48核2TB内存，RWF文件最大2.8TB，放在SSD里，算了34天

原文提及接近三十个结构都用这个级别的精度计算，而且是为了验证实验化学机理的，我感觉数据可能不太靠谱

wzkchem5 · 发表于 Post on 2023-2-13 21:49:34

ogrebit 发表于 2023-2-13 14:20
谢谢，内存单条128G好贵，看来玩不起ccsd

可以用ORCA做DLPNO-CCSD(T)计算，普通16核机子算100个原子cc-pVTZ都可以，和普通CCSD(T)比误差1kcal/mol左右

mfdsrax2 · 发表于 Post on 2023-2-13 21:49:12

ogrebit 发表于 2023-2-13 21:20
谢谢，内存单条128G好贵，看来玩不起ccsd

1TB内存应该也能算，单条64的还可以接受

ogrebit · 发表于 Post on 2023-2-13 21:20:23

mfdsrax2 发表于 2023-2-13 20:23
48核2TB内存，RWF文件最大2.8TB，放在SSD里，算了34天

谢谢，内存单条128G好贵，看来玩不起ccsd

mfdsrax2 · 发表于 Post on 2023-2-13 20:23:31

本帖最后由 mfdsrax2 于 2023-2-13 20:25 编辑

ogrebit 发表于 2023-2-13 19:58
请问你算三十个原子，用了多少核，硬件是多大内存，硬盘是多大？需要pcie3.0以上的nvme的ssd吗？

48核2TB内存，RWF文件最大2.8TB，放在SSD里，算了34天

ogrebit · 发表于 Post on 2023-2-13 19:58:24

mfdsrax2 发表于 2023-2-13 19:10
SSCD(T)是什么方法？CCSD(T)的话我算过三十多个原子的，还含氯，也能正常算，大概要算一周左右

请问你算三十个原子，用了多少核，硬件是多大内存，硬盘是多大？需要pcie3.0以上的nvme的ssd吗？

mfdsrax2 · 发表于 Post on 2023-2-13 19:10:35

MiaoZuo 发表于 2022-11-30 19:57
我最近审稿看到一个课题组G16用SSCD(T)/ccpVTZ级别下以及SMD溶解模型计算单点能，IM和TS加起来二十几个结构 ...

SSCD(T)是什么方法？CCSD(T)的话我算过三十多个原子的，还含氯，也能正常算，大概要算一周左右

MiaoZuo · 发表于 Post on 2022-11-30 20:15:08

wjc404 发表于 2022-11-30 20:07
可以看一下平均每个任务的基函数数量和价层电子总数，以及对称性的情况。光从原子数看有时候不一定能判断 ...

好的，谢谢指导~

wjc404 · 发表于 Post on 2022-11-30 20:07:57

本帖最后由 wjc404 于 2022-11-30 20:10 编辑

MiaoZuo 发表于 2022-11-30 19:57
我最近审稿看到一个课题组G16用SSCD(T)/ccpVTZ级别下以及SMD溶解模型计算单点能，IM和TS加起来二十几个结构 ...

可以看一下平均每个任务的基函数数量和价层电子总数，以及对称性的情况。光从原子数看有时候不一定能判断计算开销。

MiaoZuo · 发表于 Post on 2022-11-30 19:57:13

我最近审稿看到一个课题组G16用SSCD(T)/ccpVTZ级别下以及SMD溶解模型计算单点能，IM和TS加起来二十几个结构，每个结构都在25原子左右，我觉得很扯。结果后来还是发表在其他期刊了，还两篇类似的。社长感觉这数据是否有问题

bd1225an · 发表于 Post on 2020-3-25 17:06:26

sobereva 发表于 2020-3-25 15:43
这是实际计算经验，没有什么文献

好的
谢谢老师

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

[Gaussian/gview] 如何用高斯做CCSD（T）单点能计算？

评分 Rate