计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Tinker
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 请问M06-L泛函包含色散矫正吗

[复制链接 Copy URL]

72

帖子

1

威望

561

eV
积分
653

Level 4 (黑子)

16#
发表于 Post on 2021-9-9 14:24:10 | 只看该作者 Only view this author
借楼问一下老师:M06-2X已经包括了良好的色散描述作用,请问如果优化计算弱相互作用时单纯使用MO6-2X,以及加上色散矫正gd3后的计算结果相比会有多大差距呢?相比较下哪种方式个更可信呢?

7

帖子

0

威望

25

eV
积分
32

Level 2 能力者

17#
发表于 Post on 2021-9-10 09:57:13 | 只看该作者 Only view this author
liustone 发表于 2021-9-9 14:24
借楼问一下老师:M06-2X已经包括了良好的色散描述作用,请问如果优化计算弱相互作用时单纯使用MO6-2X,以及 ...

最近发表的一篇相关的论文,可供参考一下.
J. Chem. Theory Comput. 2021, 17, 7, 3967–3973

简单说的结论就是M06-L, M06-2X不适合加D3.
如果有审稿人问到M06-L和M06-2X加D3的问题,可以引用这篇文章.

150

帖子

0

威望

3921

eV
积分
4071

Level 6 (一方通行)

18#
发表于 Post on 2021-9-10 10:09:07 | 只看该作者 Only view this author
wudihua 发表于 2021-9-10 09:57
最近发表的一篇相关的论文,可供参考一下.
J. Chem. Theory Comput. 2021, 17, 7, 3967–3973

但这篇文章里面似乎也没泛函什么事儿啊

7

帖子

0

威望

25

eV
积分
32

Level 2 能力者

19#
发表于 Post on 2021-9-10 10:18:15 | 只看该作者 Only view this author
RAL 发表于 2021-9-10 10:09
但这篇文章里面似乎也没泛函什么事儿啊

文章里有比较M06-L和M06-2X 加D3和不加D3对结果的影响.

150

帖子

0

威望

3921

eV
积分
4071

Level 6 (一方通行)

20#
发表于 Post on 2021-9-10 10:39:29 | 只看该作者 Only view this author
本帖最后由 RAL 于 2021-9-10 10:43 编辑
wudihua 发表于 2021-9-10 10:18
文章里有比较M06-L和M06-2X 加D3和不加D3对结果的影响.我复制粘贴引用居然搜到的是另一篇文献
我复制粘贴引用居然搜到了另一篇文献

2021-09-10_103811.jpg (190.75 KB, 下载次数 Times of downloads: 72)

2021-09-10_103811.jpg

453

帖子

0

威望

8172

eV
积分
8625

Level 6 (一方通行)

21#
发表于 Post on 2021-9-10 11:24:11 | 只看该作者 Only view this author
wudihua 发表于 2021-9-10 09:57
最近发表的一篇相关的论文,可供参考一下.
J. Chem. Theory Comput. 2021, 17, 7, 3967–3973

M系列的泛函在参数拟合的时候考虑了非共价相互作用。我的理解是通过参数化让排斥小了一点,所以就表现出了吸引。

在这种情况下加D3、D4是属于引入复杂,效果不会好。

7

帖子

0

威望

25

eV
积分
32

Level 2 能力者

22#
发表于 Post on 2021-9-10 11:56:54 | 只看该作者 Only view this author
yflchx 发表于 2021-9-10 11:24
M系列的泛函在参数拟合的时候考虑了非共价相互作用。我的理解是通过参数化让排斥小了一点,所以就表现出 ...

meta-GGA可以处理部分色散作用,所以能够对非共价体系给出比较好的结果也很正常.

538

帖子

1

威望

5767

eV
积分
6325

Level 6 (一方通行)

23#
发表于 Post on 2021-9-10 15:41:35 | 只看该作者 Only view this author
RAL 发表于 2021-9-10 10:39
我复制粘贴引用居然搜到了另一篇文献

因为一般不写期号,19楼写了期号

72

帖子

1

威望

561

eV
积分
653

Level 4 (黑子)

24#
发表于 Post on 2021-9-10 19:18:41 | 只看该作者 Only view this author
wudihua 发表于 2021-9-10 09:57
最近发表的一篇相关的论文,可供参考一下.
J. Chem. Theory Comput. 2021, 17, 7, 3967–3973

好的,谢谢版主

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

25#
发表于 Post on 2021-9-11 04:39:18 | 只看该作者 Only view this author
wudihua 发表于 2021-9-10 11:56
meta-GGA可以处理部分色散作用,所以能够对非共价体系给出比较好的结果也很正常.

“meta-GGA可以处理部分色散作用”这么表述不太准确
必须看具体的泛函,诸如TPSS虽然也是meta-GGA,但描述色散作用是非常烂的。下文有个TPSS优化DNA片段的例子,可见多么糟糕
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557http://bbs.keinsci.com/thread-17899-1-1.html
26L里的图也体现出TPSS连个极小点都没有,这甚至还不如GGA级别的PBE,起码还勉强有个极小点

但可以说如果专门在设计、参数化泛函的时候考虑对弱相互作用的描述的话,meta-GGA比GGA或许有可能能做到更好、潜力更大,毕竟其自由度更大、引入的信息更多。

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

26#
发表于 Post on 2021-9-11 04:43:34 | 只看该作者 Only view this author
liustone 发表于 2021-9-9 14:24
借楼问一下老师:M06-2X已经包括了良好的色散描述作用,请问如果优化计算弱相互作用时单纯使用MO6-2X,以及 ...

加不加D3对M06-2X优化弱相互作用体系的结果不会有什么差别,对于大多数常研究的体系大小范畴,计算相互作用能也不会有多大差别。但计算大面积pi-pi作用我建议加D3

出自GMTKN55文章






相关内容我在此文也说过了
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413http://bbs.keinsci.com/thread-9772-1-1.html


另外关于前面有人提到的JCTC文章,我之前粗略看过,个人认为,对于严格讨论D3的影响,如果用相对小一些的体系,使得带弥散函数的情况(或者QZ级别的基组)也能算得动,再进行对比更接近于一般研究者的实际研究场景。而且作为参考标准的DLPNO-CCSD(T)如果基组不上到QZ的话(或只用不带弥散的cc-pVTZ),难免可能由于误差巧合抵消原因不能如实反应廉价方法本身的精度差异,更何况后HF对基组的档次要求比普通泛函明显高一个档。另外测试时最好能加入目前已知算弱相互作用最好的非双杂化泛函wB97M-V。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

27#
发表于 Post on 2021-9-11 04:58:26 | 只看该作者 Only view this author

我是量化版的版主


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

72

帖子

1

威望

561

eV
积分
653

Level 4 (黑子)

28#
发表于 Post on 2021-9-12 10:17:37 | 只看该作者 Only view this author
sobereva 发表于 2021-9-11 04:58
我是量化版的版主

谢谢版主老师,我学到了很多,你前面提到的帖子我都收藏并认真阅读过了
借楼请教一个问题:就是我前面提到的溶剂-甲硫醇体系结构优化+振动分析,整个体系并不复杂,20-30个原子的体系,按照这两天的学习,我运行的命令是#p opt=calcall freq m062x/6-311+g(d) em=gd3 int=ultrafine,在进行计算时我分别尝试甲硫醇巯基对着溶剂的N O S原子构建分子体系,可是总遇到这个错误导致计算停止:
Cartesian Forces:  Max     0.000212907 RMS     0.000042805
Linear angle in Bend.
Error termination via Lnk1e in /gdata/lun7/software/G09E01/g09/l716.exe at Sun Sep 12 04:41:26 2021.
Job cpu time:       9 days  1 hours 17 minutes 39.9 seconds.
File lengths (MBytes):  RWF=    304 Int=      0 D2E=      0 Chk=      7 Scr=      2
我尝试过使用opt=cartesian指令运行,计算时间会变得很长,而且有的体系还是会报出这个错误,请问该怎么办?

1万

帖子

0

威望

9871

eV
积分
22113

Level 6 (一方通行)

29#
发表于 Post on 2021-9-12 16:17:48 | 只看该作者 Only view this author
liustone 发表于 2021-9-12 03:17
谢谢版主老师,我学到了很多,你前面提到的帖子我都收藏并认真阅读过了
借楼请教一个问题:就是我前面提 ...

以优化的最后一个结构为初始结构,继续优化即可。
这个问题是所有在内坐标下优化的算法的共性,就是如果初始结构里有一个键角不是180度,但优化过程中变成了180度,就会报错,因为会导致除以0等问题。因为程序无法预知哪个键角会在优化过程中变成180度,所以无法提前避免这个问题,只能是你发现这个报错以后手动继续优化
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

72

帖子

1

威望

561

eV
积分
653

Level 4 (黑子)

30#
发表于 Post on 2021-9-12 22:04:50 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-12 16:17
以优化的最后一个结构为初始结构,继续优化即可。
这个问题是所有在内坐标下优化的算法的共性,就是如果 ...

好嘞,我以最后的文件保存为新文件继续优化,优化命令需要更换嘛?是不是最从opt=calcall换成cartesian啊?

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 02:12 , Processed in 0.297349 second(s), 22 queries , Gzip On.

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