一个用DLPNO-CCSD(T)优化18碳环的结构的思路

zjxitcc · 发表于 Post on 2022-8-10 23:09:43

ionexchangeC 发表于 2022-8-10 22:39
用DLPNO-CCSD(T)-M06-2X/cc-pVTZ（M06-2X为参考态的DLPNO-CCSD(T)）按照这个方法优化，最后的结果是长碳碳 ...

为啥不用基于自旋极化UHF的DLPNO-UCCSD(T)。。

sobereva · 发表于 Post on 2022-8-11 07:09:55

zjxitcc 发表于 2022-8-10 23:09
为啥不用基于自旋极化UHF的DLPNO-UCCSD(T)。。

DLPNO-UCCSD(T)比DLPNO-CCSD(T)昂贵甚巨
而且UHF错误地高估了这个体系的对称破缺特征，用U还不如用R
另外U和R我实际用CCSD优化的结果基本没差别

ionexchangeC · 发表于 Post on 2022-8-15 11:48:37

又用CASSCF和NEVPT2优化了18碳环，取的活性空间是（8，8）。
CASSCF：长键1.356埃，短键1.199埃
NEVPT2：长键等于短键，1.289埃

CASSCF.txt (2.34 KB, 下载次数 Times of downloads: 4)

FIC-NEVPT2.txt (2.58 KB, 下载次数 Times of downloads: 7)

wjc404 · 发表于 Post on 2022-10-25 11:44:51

本帖最后由 wjc404 于 2022-10-25 15:34 编辑

是否可以用CCSD(T)-F12(RI)/cc-pVDZ-F12试一下优化(看情况启用DLPNO；这个级别的精度相当于不带F12时的TZ基组结果)？

另外想请教一下，现在有程序支持“DFT的基态上做CC方法”的解析梯度吗？

ionexchangeC · 发表于 Post on 2022-11-4 23:15:38

wjc404 发表于 2022-10-25 11:44
是否可以用CCSD(T)-F12(RI)/cc-pVDZ-F12试一下优化(看情况启用DLPNO；这个级别的精度相当于不带F12时的TZ基 ...

按理来说，你这个计算级别就算加了DLPNO也比普通地用DLPNO-CCSD(T)/cc-pVTZ贵吧

wjc404 · 发表于 Post on 2022-11-7 23:47:51

本帖最后由 wjc404 于 2022-11-8 04:10 编辑

ionexchangeC 发表于 2022-8-15 11:48
又用CASSCF和NEVPT2优化了18碳环，取的活性空间是（8，8）。
CASSCF：长键1.356埃，短键1.199埃
NEVPT2： ...

看了一下（R1=1.22 Å，R2=1.34 Å时） RHF/def2TZVP的轨道hessian的（共3个）本征值小于零对应的向量中大于0.1的元素，发现牵涉最多的轨道为最高的8个占有轨道（均为π）和最低的8个未占有π轨道。它们在每个向量中都有较高贡献，且互相之间没有绝对的主次。猜测在轨道优化时仅考虑（8，8）的活性空间可能偏小了。但再往大就比较昂贵了

。

ionexchangeC · 发表于 Post on 2022-11-8 11:38:55

cyclo[18]carbon.molden.input (7.48 MB, 下载次数 Times of downloads: 3)

ionexchangeC · 发表于 Post on 2022-11-8 11:39:20

wjc404 发表于 2022-11-7 23:47
看了一下（R1=1.22 Å，R2=1.34 Å时） RHF/def2TZVP的轨道hessian的（共3个）本征值小于零对应 ...

已上传CASSCF波函数，欢迎讨论

wjc404 · 发表于 Post on 2022-12-5 10:16:24

本帖最后由 wjc404 于 2022-12-5 10:24 编辑

ionexchangeC 发表于 2022-11-8 11:39
已上传CASSCF波函数，欢迎讨论

谢谢。我用ORCA跑cas(8,8)出来的波函数和这个差不多，能够收敛，但是后续的MCRPA计算发现casscf波函数不稳定（以及nevpt2计算中警告Koopmans矩阵存在负的本征值（结果未上传））。
计算结果: https://pan.baidu.com/s/1mL_uUeRGpgR_oknlpXv9rg 提取码: ari5

Kalinite · 发表于 Post on 2022-12-19 11:01:54

本帖最后由 Kalinite 于 2022-12-19 11:04 编辑

Mikasa 发表于 2022-7-15 09:43
可能只是手册没更新，我们已经用上了~

您好，请问应该通过什么方式实现呢，我在输入里直接调ALASKA，输出文件中显示仍然是numerical gradient。

Numerical gradient, root 1
---------------------------------------------
X Y Z
---------------------------------------------
O1 -0.046735 -0.066128 -0.000000
H1 0.086732 -0.011718 -0.000000
H2 -0.039997 0.077846 0.000000
---------------------------------------------
--- Stop Module: numerical_gradient at Mon Dec 19 10:04:31 2022 /rc=_RC_ALL_IS_WELL_ ---
*** files: xmldump
saved to directory /export/home/chli/xjw/test/molcas/caspt2grad
--- Module numerical_gradient spent 17 seconds ---

复制代码

ionexchangeC · 发表于 Post on 2023-12-12 19:35:08

2023.12.12更新：
用RHF为参考态，DLPNO-CCSD(T)/cc-pVTZ级别优化18碳环，结果是长键1.359埃，短键1.232埃

chrinide · 发表于 Post on 2023-12-14 23:20:33

ionexchangeC 发表于 2023-12-12 19:35
2023.12.12更新：
用RHF为参考态，DLPNO-CCSD(T)/cc-pVTZ级别优化18碳环，结果是长键1.359埃，短键1.232埃

这是算了多长时间了？不过结果可以用来引用的

ionexchangeC · 发表于 Post on 2023-12-14 23:23:46

chrinide 发表于 2023-12-14 23:20
这是算了多长时间了？不过结果可以用来引用的

我是为了测试超算做的这个任务，而我现在甚至忘记了测试账号的密码

Eclair · 发表于 Post on 2024-2-10 00:14:52

mfdsrax2 发表于 2022-7-14 10:34
这个要写两年吗

你觉得慢可以自己来写，ccsd（T）的解析梯度很早以前就推出来了，

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

[算法与编程] 一个用DLPNO-CCSD(T)优化18碳环的结构的思路