计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2402|回复 Reply: 39

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

[复制链接 Copy URL]

133

帖子

2

威望

803

eV
积分
976

Level 4 (黑子)

发表于 Post on 2022-7-13 23:29:21 | 显示全部楼层 Show all |阅读模式 Reading model
由于DLPNO-CCSD(T)没有解析梯度,常规的优化方式没法用,ORCA又不能利用对称性在内坐标下优化。于是我自己写了一个调用ORCA的小程序实现D9h结构的18碳环的结构优化。由于DLPNO-CCSD(T)结合像样的基组在我的电脑上算个单点都很花时间,我换成了PWPB95-D4/def2-TZVP这个级别(也没有解析梯度)。优化算法是NEWUOA。
18C.f90 (51.76 KB, 下载次数 Times of downloads: 6)

评分 Rate

参与人数
Participants 4
eV +18 收起 理由
Reason
biogon + 5
sobereva + 5
Freeman + 5 赞!
wzkchem5 + 3

查看全部评分 View all ratings

6691

帖子

0

威望

4012

eV
积分
10703

Level 6 (一方通行)

发表于 Post on 2022-7-14 01:34:54 | 显示全部楼层 Show all
近两年ORCA可能会实现DLPNO-CCSD(T)解析梯度,我们这边有一个博士生正在写。在此之前会实现canonical CCSD(T)解析梯度。
此外有兴趣的话可以看看NEWUOA和ORCA内置的基于数值梯度的结构优化算法哪个快。对于18碳环这样的高对称性分子肯定前者快,但是可以测一下非对称性分子,假如也是NEWUOA快的话,可以考虑写一个基于NEWUOA 的通用结构优化脚本,上传到这里https://gitlab.gwdg.de/orca-helpers/orca-helpers/-/tree/master/(这里有很多普通orca用户写的脚本),这样其他人也可以用

评分 Rate

参与人数
Participants 2
eV +10 收起 理由
Reason
北大-陶豫 + 5 赞!
Freeman + 5 赞!

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

4万

帖子

99

威望

4万

eV
积分
89946

管理员

公社社长+计算化学玩家

发表于 Post on 2022-7-14 06:29:37 | 显示全部楼层 Show all
顺带一提,一些范围分离双杂化泛函可以正确优化C18,如SOS1-QIDH和RSX-QIDH
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

604

帖子

1

威望

4899

eV
积分
5523

Level 6 (一方通行)

发表于 Post on 2022-7-14 09:56:12 | 显示全部楼层 Show all
本帖最后由 hebrewsnabla 于 2022-7-14 10:06 编辑

对于DLPNO系列方法,包括DLPNO-MP2和DLPNO-CC,数值梯度可能有比较大的问题,参见 https://orcaforum.kofo.mpg.de/vi ... 7610&p=32831#p32798 ,原因是数值扰动时如果轨道domain改变,势能面就不够光滑。而解析梯度(如果存在的话)不会有这个问题。

虽然NEWUOA不是基于数值梯度的,但是搞不好也有这个问题。

256

帖子

0

威望

964

eV
积分
1220

Level 4 (黑子)

发表于 Post on 2022-7-14 10:34:52 | 显示全部楼层 Show all
wzkchem5 发表于 2022-7-14 01:34
近两年ORCA可能会实现DLPNO-CCSD(T)解析梯度,我们这边有一个博士生正在写。在此之前会实现canonical CCSD( ...

这个要写两年吗

2609

帖子

1

威望

5816

eV
积分
8445

Level 6 (一方通行)

MOKIT开发者

发表于 Post on 2022-7-14 10:52:24 | 显示全部楼层 Show all

两年能写出来就已经很好了,这巨复杂。比如MOLCAS的CASPT2解析导数做了N年也没做出来
自动做多参考态计算的程序MOKIT

256

帖子

0

威望

964

eV
积分
1220

Level 4 (黑子)

发表于 Post on 2022-7-14 10:54:16 | 显示全部楼层 Show all
zjxitcc 发表于 2022-7-14 10:52
两年能写出来就已经很好了,这巨复杂。比如MOLCAS的CASPT2解析导数做了N年也没做出来

DLPNO-CCSD(T)相比CCSD(T)的解析导数是更容易还是复杂

2609

帖子

1

威望

5816

eV
积分
8445

Level 6 (一方通行)

MOKIT开发者

发表于 Post on 2022-7-14 11:00:21 | 显示全部楼层 Show all
mfdsrax2 发表于 2022-7-14 10:54
DLPNO-CCSD(T)相比CCSD(T)的解析导数是更容易还是复杂

即使是正则CCSD(T)解析导数也很复杂。而局域相关方法的解析导数复杂得多,即使写出来了,还要考虑解析导数是不是也实现了线性标度的问题。
自动做多参考态计算的程序MOKIT

256

帖子

0

威望

964

eV
积分
1220

Level 4 (黑子)

发表于 Post on 2022-7-14 11:07:01 | 显示全部楼层 Show all
zjxitcc 发表于 2022-7-14 11:00
即使是正则CCSD(T)解析导数也很复杂。而局域相关方法的解析导数复杂得多,即使写出来了,还要考虑解析导 ...

那为什么不先尝试CCSD(T)的解析导数,而是花两年做更难的DLPNO-CCSD(T)导数。CCSD(T)的解析导数对于计算化学的意义也很大

2609

帖子

1

威望

5816

eV
积分
8445

Level 6 (一方通行)

MOKIT开发者

发表于 Post on 2022-7-14 11:11:02 | 显示全部楼层 Show all
本帖最后由 zjxitcc 于 2022-7-14 11:13 编辑
mfdsrax2 发表于 2022-7-14 11:07
那为什么不先尝试CCSD(T)的解析导数,而是花两年做更难的DLPNO-CCSD(T)导数。CCSD(T)的解析导数对于计算 ...

2L里有说“在此之前会实现canonical CCSD(T)解析梯度”,这就是正则CCSD(T),或者说传统CCSD(T)。传统CCSD(T)解析导数已经有好几个量化程序做出来了,ORCA实现它只能是增添一个小功能(因为具有可替代性,且传统CCSD(T)标度是N^7,原理上就没法算大体系),作为自己写代码时的验证。重要的是DLPNO-CCSD(T)解析导数。
自动做多参考态计算的程序MOKIT

6691

帖子

0

威望

4012

eV
积分
10703

Level 6 (一方通行)

发表于 Post on 2022-7-14 17:00:16 | 显示全部楼层 Show all
mfdsrax2 发表于 2022-7-14 03:54
DLPNO-CCSD(T)相比CCSD(T)的解析导数是更容易还是复杂

基本上任何近似方法的代码都比相应的非近似方法的代码复杂。至少当你只需要做对而不需要做快的时候是如此的。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

214

帖子

0

威望

4838

eV
积分
5052

Level 6 (一方通行)

跳跳猪

发表于 Post on 2022-7-15 09:10:26 | 显示全部楼层 Show all
zjxitcc 发表于 2022-7-14 10:52
两年能写出来就已经很好了,这巨复杂。比如MOLCAS的CASPT2解析导数做了N年也没做出来

现在的OpenMOLCAS有CASPT2的解析梯度了呀
流年似水,浮生如梦。

1158

帖子

1

威望

2801

eV
积分
3979

Level 5 (御坂)

发表于 Post on 2022-7-15 09:39:24 | 显示全部楼层 Show all
Mikasa 发表于 2022-7-15 09:10
现在的OpenMOLCAS有CASPT2的解析梯度了呀

我看22jul08的手册还说没有啊

1158

帖子

1

威望

2801

eV
积分
3979

Level 5 (御坂)

发表于 Post on 2022-7-15 09:40:51 | 显示全部楼层 Show all
wzkchem5 发表于 2022-7-14 17:00
基本上任何近似方法的代码都比相应的非近似方法的代码复杂。至少当你只需要做对而不需要做快的时候是如此 ...

orca的NEVPT2梯度现在有没有人在做,以前论坛上说打算做了

214

帖子

0

威望

4838

eV
积分
5052

Level 6 (一方通行)

跳跳猪

发表于 Post on 2022-7-15 09:43:37 | 显示全部楼层 Show all
biogon 发表于 2022-7-15 09:39
我看22jul08的手册还说没有啊

可能只是手册没更新,我们已经用上了~
流年似水,浮生如梦。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2023-2-6 04:51 , Processed in 0.641075 second(s), 26 queries .

快速回复 返回顶部 返回列表 Return to list