计算化学公社

 找回密码 Forget password
 注册 Register
Views: 597|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] orca里面D3BJ-ATM怎么调用,会影响计算速度吗?

[复制链接 Copy URL]

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

D3BJ加上ATM项是不是要自己设置?不是默认的?
优化过程用D3BJ-ATM的计算速度会不会比D3慢?

1万

帖子

0

威望

7398

eV
积分
18160

Level 6 (一方通行)

2#
发表于 Post on 2024-2-6 17:47:38 | 只看该作者 Only view this author
一般不是默认的。
D3BJ-ATM的计算时间是D3BJ的很多倍,但因为D3BJ极其快,所以其实计算时间也不多。一般认为对于DFT,且体系很大(例如几百个原子)的情况,加ATM是值得的。但对于半经验方法往往不值得加
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-2-6 23:29:38 | 只看该作者 Only view this author
wzkchem5 发表于 2024-2-6 17:47
一般不是默认的。
D3BJ-ATM的计算时间是D3BJ的很多倍,但因为D3BJ极其快,所以其实计算时间也不多。一般认 ...

您好,我对ATM在什么时候使用有些困惑,所以提了以下几个问题:

1)ATM在大体系值得加是明显提升精度还是什么其他原因呀?

2)是不是只用在能量计算的时候加ATM就足够了?

3)优化大体系的时候只用D3或者D3BJ,不用ATM,会不会有影响?

4)频率计算的时候是不是可以不用ATM,只用D3或者D3BJ?

1万

帖子

0

威望

7398

eV
积分
18160

Level 6 (一方通行)

4#
发表于 Post on 2024-2-7 00:01:59 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-6 16:29
您好,我对ATM在什么时候使用有些困惑,所以提了以下几个问题:

1)ATM在大体系值得加是明显提升精度 ...

(1)因为ATM是三个原子之间的作用,从整个分子挑出3个原子的方式数和分子大小的三次方成正比,但挑2个原子的方式数和分子大小平方成正比,所以分子越大,ATM项占总色散贡献的比例就越大,当然ATM项的绝对大小也越大
(2)对,一般程序也不支持ATM的解析梯度,所以不方便加ATM优化结构
(3)ATM项即使再大也比D3小得多,一般是D3/D3BJ项的10%左右,所以几乎不会有必须加ATM项才能把结构算准的情形。如果发现D3/D3BJ不够准,建议首先考虑改用D4或含VV10的泛函(wB97M-V等),虽然改善的方面和ATM改善的方面不一样,但是D4、VV10很多程序支持解析梯度,结构优化比较好做。
(4)理论级别和结构优化统一即可
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

199

帖子

2

威望

1526

eV
积分
1765

Level 5 (御坂)

5#
发表于 Post on 2024-2-8 23:10:15 | 只看该作者 Only view this author
wzkchem5 发表于 2024-2-7 00:01
(1)因为ATM是三个原子之间的作用,从整个分子挑出3个原子的方式数和分子大小的三次方成正比,但挑2个原 ...

即使不支持解析梯度用数值梯度,三体项相比SCF还是快不少吧(解析梯度标度也就N^4)

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 13:43 , Processed in 0.183052 second(s), 20 queries , Gzip On.

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