计算化学公社

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

[ORCA] def2-QZVP'是benchmark文献里面采用的基组,但是查找基组的原文献,没有加撇号的基组

[复制链接 Copy URL]

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我打算计算能量的时候采用文献建议的方法:PW6B95/def2-QZVP',但是遇到了一些设置的问题:
1)def2-QZVP'加了撇号,是指的def2-QZVPP还是def2-QZVP呢?我在引用文献里没有查到加撇号的表示方法。

2)PW6B95在orca中是不是需要设置RIJCOSX加速,还是说默认加速方法,不用设置?


5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

2#
发表于 Post on 2024-2-12 17:42:11 | 只看该作者 Only view this author
1 不是标准写法。若无特殊情况,完全没必要效仿
若非要重复文献的数据,看文章里引用的文献

2 毫无必要用PW6B95,ORCA里能用的wB97M-V远比这个强,看下文。ORCA 5开始杂化泛函默认用RIJCOSX
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-2-13 17:42:05 | 只看该作者 Only view this author
sobereva 发表于 2024-2-12 17:42
1 不是标准写法。若无特殊情况,完全没必要效仿
若非要重复文献的数据,看文章里引用的文献

Grimme的benchmark试了wB97M-V和B97M-V等方法,发现B97M-V计算相互作用能的误差最小,略优于wB97M-V。
不过我计算相互作用能也没有实验数据做对比,所以就跟着benchmark文献推荐的方法算。

请问B97M-V是不是自带D3和ATM项?我搜不到相关文献说B97M-V是带了哪些校正项的,好奇怪。

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

4#
发表于 Post on 2024-2-13 19:12:40 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-13 17:42
Grimme的benchmark试了wB97M-V和B97M-V等方法,发现B97M-V计算相互作用能的误差最小,略优于wB97M-V。
...

B97M-V算弱相互作用明显不及wB97M-V,诸多文章都体现了这一点,比如看Mol. Phys., 115, 2315 (2017)图35、36弱相互作用排名,以及J. Chem. Theory Comput. 2018, 14, 5725表3。对于纯泛函SIE问题可能显著的情况,B97M-V更是远不及wB97M-V靠谱。不要简单用“Grimme的benchmark”来描述,Grimme的benchmark文章多了去了

B97M-V带的是VV10色散校正,跟D3是两码事
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-2-13 23:46:21 | 只看该作者 Only view this author
sobereva 发表于 2024-2-13 19:12
B97M-V算弱相互作用明显不及wB97M-V,诸多文章都体现了这一点,比如看Mol. Phys., 115, 2315 (2017)图35 ...

我刚用orca试了计算溶剂化单点能wB97M-V/QZVP, smd, defgrid3, tightscf, noautostart miniprint nopop,给了6700MB每核,16核。跑一个300多个原子的结构,结果是ORCA finished by error termination in SCF。
我设置的内存和核数在计算频率的时候够用,但是当时计算频率用的方法是TPSS-D3/def2-svp。

报错的原因是我基组太大,内存给的太小吗?但是计算能量的时候希望用大基组来避免BSSE,不知道该怎么办才好。

135

帖子

0

威望

1549

eV
积分
1684

Level 5 (御坂)

6#
发表于 Post on 2024-2-14 10:39:10 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-13 23:46
我刚用orca试了计算溶剂化单点能wB97M-V/QZVP, smd, defgrid3, tightscf, noautostart miniprint nopop, ...

你可以先4个或者10个核跑跑看,看是否跑得起来,确实很有可能内存不足。

269

帖子

0

威望

1758

eV
积分
2027

Level 5 (御坂)

7#
发表于 Post on 2024-2-14 10:53:17 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-13 23:46
我刚用orca试了计算溶剂化单点能wB97M-V/QZVP, smd, defgrid3, tightscf, noautostart miniprint nopop, ...

体系太大了,即使你分配再多的内存也不大可能算得动QZVP。
担心BSSE可以用GCP。
见社长这篇文章第三节:http://sobereva.com/214

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

8#
发表于 Post on 2024-2-14 12:02:43 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-13 23:46
我刚用orca试了计算溶剂化单点能wB97M-V/QZVP, smd, defgrid3, tightscf, noautostart miniprint nopop, ...

300原子还QZ就别想了
wB97M-V/def2-TZVP结合gCP倒还差不多,普通双路服务器下一般几个小时内跑完。想更便宜些还可以用def2-TZVP(-f)。如果要更便宜,用r2SCAN-3c(前提是纯泛函对此体系没有明显SIE问题的情况)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2024-2-14 20:35:03 | 只看该作者 Only view this author
renzhogn424 发表于 2024-2-14 10:39
你可以先4个或者10个核跑跑看,看是否跑得起来,确实很有可能内存不足。

试过10核,会报错说内存不够。

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-2-14 20:36:03 | 只看该作者 Only view this author
dnlx 发表于 2024-2-14 10:53
体系太大了,即使你分配再多的内存也不大可能算得动QZVP。
担心BSSE可以用GCP。
见社长这篇文章第三节 ...

真的不能用这么大基组,我换def2-TZVP可以跑了

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-2-14 20:37:21 | 只看该作者 Only view this author
sobereva 发表于 2024-2-14 12:02
300原子还QZ就别想了
wB97M-V/def2-TZVP结合gCP倒还差不多,普通双路服务器下一般几个小时内跑完。想更 ...

那可以只在计算能量的时候加gCP吗?优化和频率计算不加gCP可以吗?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

12#
发表于 Post on 2024-2-15 13:53:54 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-14 20:37
那可以只在计算能量的时候加gCP吗?优化和频率计算不加gCP可以吗?

opt freq加不加gCP没多大影响,加了也没坏处
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2024-2-16 20:33:20 | 只看该作者 Only view this author
sobereva 发表于 2024-2-15 13:53
opt freq加不加gCP没多大影响,加了也没坏处

我找到的几篇同类型的体系都是只在能量计算的时候加gCP。
opt freq不加gCP没有什么影响,那就不加了。

但是想问下是不是opt freq加gCP会明显增加耗时?

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

14#
发表于 Post on 2024-2-16 20:41:53 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-16 20:33
我找到的几篇同类型的体系都是只在能量计算的时候加gCP。
opt freq不加gCP没有什么影响,那就不加了。
...

诸如6-31G*这种档次,带不带gCP对优化出的结构有可能有不可忽视的影响,比如弱相互作用距离相差个<0.2埃,但def2-TZVP这种档次了做opt freq带不带gCP对结果的影响完全可以忽略不计
原理上不增加耗时
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

191

帖子

0

威望

493

eV
积分
684

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2024-2-16 23:12:11 | 只看该作者 Only view this author
sobereva 发表于 2024-2-16 20:41
诸如6-31G*这种档次,带不带gCP对优化出的结构有可能有不可忽视的影响,比如弱相互作用距离相差个

如果优化和频率用def2-svp,那加不加gCP有影响吗?
打算优化和频率用def2-svp,能量用def2-TZVP。
如果优化和频率用def2-TZVP,那能量的基组要增大,QZVP反正跑不动。

本版积分规则 Credits rule

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

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

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