计算化学公社

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

[量化理论] M06-2X加D3校正的问题

[复制链接 Copy URL]

516

帖子

1

威望

5052

eV
积分
5588

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
sob大大你好,刚刚看了你的“简谈量子化学计算中的泛函的选择 ”这篇文章,里面提到对弱相互作用体系比较合适的泛函之一是M06-2X-D3。在M06-2X中已经有了以弱相互作用的拟合,如果再加D3校正会不会过?一般来说是加好还是不加好?有没有相关文献讨论过这个问题?

谢谢!

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
找雷 + 1

查看全部评分 View all ratings

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

2#
发表于 Post on 2018-1-23 11:06:39 | 只看该作者 Only view this author
完全不会过。
D3不应当理解为是用来描述色散作用的,而是对你选择的泛函基础上进行校正以提高精度,只不过是对于色散作用描述极差的诸如B3LYP,D3主要作用才是起到描述色散作用的作用。D3的参数对每种泛函是不同的。M06-2X已经能很好描述色散作用,所以加了D3只会稍微改进一点精度,完全不会double-counting。明尼苏达系列泛函不能用BJ阻尼的D3是因为BJ阻尼下必然导致double-counting。这些问题Grimme都已经充分考虑了。

评分 Rate

参与人数
Participants 3
eV +12 收起 理由
Reason
超限制抱怨 + 2 谢谢
让你变成回忆 + 5 谢谢
Xin + 5 赞!

查看全部评分 View all ratings

北京科音自然科学研究中心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!

516

帖子

1

威望

5052

eV
积分
5588

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2018-1-23 13:18:46 | 只看该作者 Only view this author
sobereva 发表于 2018-1-23 11:06
完全不会过。
D3不应当理解为是用来描述色散作用的,而是对你选择的泛函基础上进行校正以提高精度,只不过 ...

非常感谢!

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

4#
发表于 Post on 2021-1-14 14:22:51 | 只看该作者 Only view this author
本帖最后由 星斗如盘 于 2021-1-14 14:27 编辑
sobereva 发表于 2018-1-23 11:06
完全不会过。
D3不应当理解为是用来描述色散作用的,而是对你选择的泛函基础上进行校正以提高精度,只不过 ...

使用b3lyp-d3(BJ)/6-31g(d) 优化结构,使用b3lyp-d3(BJ)/6-311+G(d,p)计算能量怎么样?和M06/6-311+G(d,p)计算能量哪个更好呢?PS: M06貌似不能加d3(BJ)会报错

188

帖子

2

威望

3202

eV
积分
3430

Level 5 (御坂)

5#
发表于 Post on 2021-1-14 14:38:46 | 只看该作者 Only view this author
星斗如盘 发表于 2021-1-14 14:22
使用b3lyp-d3(BJ)/6-31g(d) 优化结构,使用b3lyp-d3(BJ)/6-311+G(d,p)计算能量怎么样?和M06/6-311+G(d,p ...

你说的是M06-2X?这个只能加D3(也就是D3ZERO)不能加D3BJ. 有机体系算热力学现在用M06-2X好一点。

235

帖子

1

威望

1176

eV
积分
1431

Level 4 (黑子)

6#
发表于 Post on 2021-1-14 15:00:59 | 只看该作者 Only view this author
量化小菜鸡 发表于 2021-1-14 14:38
你说的是M06-2X?这个只能加D3(也就是D3ZERO)不能加D3BJ. 有机体系算热力学现在用M06-2X好一点。

不纯粹是有机体系,有一个过渡金属,所以使用了M06,但是直接写m06/genecp em=gd3bj就报错了

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

7#
发表于 Post on 2021-1-14 17:12:17 | 只看该作者 Only view this author
星斗如盘 发表于 2021-1-14 15:00
不纯粹是有机体系,有一个过渡金属,所以使用了M06,但是直接写m06/genecp em=gd3bj就报错了

M06同理,可以加GD3而不能加GD3BJ
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

8#
发表于 Post on 2021-1-15 07:52:27 | 只看该作者 Only view this author
星斗如盘 发表于 2021-1-14 14:22
使用b3lyp-d3(BJ)/6-31g(d) 优化结构,使用b3lyp-d3(BJ)/6-311+G(d,p)计算能量怎么样?和M06/6-311+G(d,p ...

为什么不能加这里明确说了
DFT-D色散校正的使用
http://sobereva.com/210

如果被研究的问题只涉及有机配体部分,用M06-2X可以,如果牵扯到配位键部分,绝对不能用M06-2X,仔细看
简谈量子化学计算中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!

187

帖子

1

威望

1245

eV
积分
1452

Level 4 (黑子)

9#
发表于 Post on 2022-2-26 01:14:23 | 只看该作者 Only view this author
sobereva 发表于 2018-1-23 11:06
完全不会过。
D3不应当理解为是用来描述色散作用的,而是对你选择的泛函基础上进行校正以提高精度,只不过 ...

老师您好! 请问D3对于M06-2X不会重复计算这个在D3的文献中有写吗? 准备好好看一下然后保存下来,防止有些低水平的审稿人argue使用M06-2X-D3进行计算。谢谢老师!
计算小白 水平很菜
随时可喷 万望指点

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

10#
发表于 Post on 2022-2-26 01:51:00 | 只看该作者 Only view this author
WaterMirror 发表于 2022-2-25 18:14
老师您好! 请问D3对于M06-2X不会重复计算这个在D3的文献中有写吗? 准备好好看一下然后保存下来,防止有 ...

可以引D3的原始文章https://aip.scitation.org/doi/pdf/10.1063/1.3382344和D3BJ的原始文章https://onlinelibrary.wiley.com/doi/epdf/10.1002/jcc.21759
在argue的时候一定要指出,此处给M06-2X加的D3是用的zero damping而不是BJ damping,因此只描述了长程色散作用而没有描述中程色散作用,M06-2X本身则只有中程色散作用而没有长程色散作用,所以虽然两者都描述色散,却没有重复计算任何东西。其他泛函常用的D3BJ是既描述了长程色散作用又描述了中程色散作用,所以D3BJ是不能搭配M06-2X等包含中程色散作用的泛函的。
那些误以为M06-2X不能加D3的人,一般分两类:
(1)只用过D3BJ,没怎么用过基于zero damping的D3,因此误以为D3BJ的结论适用于所有DFT-D方法;
(2)用的不是高斯,而是关键词写法和高斯不同的其他软件。比如在ORCA里面写D3对应于D3BJ,写D3zero才是高斯的D3,所以你如果不把zero写出来,那从纯ORCA用户的视角看,可能就会误认为M06-2X-D3是错误的。

评分 Rate

参与人数
Participants 3
eV +14 收起 理由
Reason
乐平 + 5 我很赞同
joyzhan + 4 谢谢
WaterMirror + 5 谢谢

查看全部评分 View all ratings

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

187

帖子

1

威望

1245

eV
积分
1452

Level 4 (黑子)

11#
发表于 Post on 2022-2-26 12:42:27 | 只看该作者 Only view this author
wzkchem5 发表于 2022-2-26 01:51
可以引D3的原始文章https://aip.scitation.org/doi/pdf/10.1063/1.3382344和D3BJ的原始文章https://onlin ...

谢谢老师详细的解答!
计算小白 水平很菜
随时可喷 万望指点

7

帖子

0

威望

203

eV
积分
210

Level 3 能力者

12#
发表于 Post on 2024-10-29 08:38:45 | 只看该作者 Only view this author
各位老师好,对于D3较正只适用于H~Pu原子,那么对于Pu原子后的Am,Cm原子该怎么加D3较正呢?

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

13#
发表于 Post on 2024-10-29 10:05:30 | 只看该作者 Only view this author
gaojiahao 发表于 2024-10-29 01:38
各位老师好,对于D3较正只适用于H~Pu原子,那么对于Pu原子后的Am,Cm原子该怎么加D3较正呢?

如果没有人拟合过这些元素的D3参数的话,可以不加。
或者也可以用wB97M-V等可以完全考虑长程色散的泛函
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

7

帖子

0

威望

203

eV
积分
210

Level 3 能力者

14#
发表于 Post on 2024-10-30 16:17:11 | 只看该作者 Only view this author
wzkchem5 发表于 2024-10-29 10:05
如果没有人拟合过这些元素的D3参数的话,可以不加。
或者也可以用wB97M-V等可以完全考虑长程色散的泛函

好的,谢谢老师得回答。

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 11:31 , Processed in 0.176378 second(s), 22 queries , Gzip On.

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