计算化学公社

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

[综合交流] 坚持使用B3LYP算有机体系能量的人的下场

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

跳转到指定楼层 Go to specific reply
#
前阵子接到一个J Mol Struct的审稿,算的是3+2环加成,作者用的B3LYP/6-311+G**(优化和单点能都是)。有实验有计算,内容挺丰富,发JMS怎么都够了。虽然B3LYP算有机体系能量早就是老黄历了,早该淘汰了,可有些人顽固不化,但在我来看还不足以成为需要重算所有数据的理由,所以我的审稿意见也没提这点。没想到修回的稿子里,第二个审稿人提出:
- It is well-know that any calculation using b3lyp functional must be  done including empirical dispersion correction. All the calculations  should be repeated coonsidering dispersion correction

作者很不幸,碰到个较真的审稿人。作者修回的稿子里的回复是:
All calculations using B3LYP functional were repeated consider empirical dispersion correction B3LYP-D3. As well as calculations were performed using the M06-2X and M11 functionals. Corrected in accordance with the Reviewer’s suggestions (see in the text, page 6 and in Supplementary Materials).
多惨,数据重算了一遍。(用M11这性能平庸的泛函计算是纯粹多余的,可见作者对泛函的发展现状明显缺乏关注)

我在下文里曾明确表态什么泛函适合算有机体系:
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html
文中明确说了,用M06-2X算有机体系能量是值得推荐的(还有更好的wB97M-V)。当初这作者要是用M06-2X,哪用得着如今重算?这年头还用B3LYP算有机体系能量问题已经严重底气不足了,虽然当前计算加D3不是有绝对必要性,但加了也比直接用干巴巴的B3LYP强、也不至于招致这个审稿人的这顿批。我在此文里也强调了加D3可以改进反应能的计算:
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413http://bbs.keinsci.com/thread-9772-1-1.html
不过,这个审稿人过度重视色散校正、甚至不是非得加色散校正的情况还要让作者加了重算,这点我是不赞同的。

我之前在此文中鼓励使用def或def2基组代替尺寸相近的Pople基组:
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
这审稿人还有个意见
- Pople basis  sets are not the best. A solid justification of their use and comparison with other levels of theory must be included
作者的回复是
In the text, calculations have been added at various levels of theory (see in the text, page 6 and in Supplementary Materials – Tables S3 and S4).
可见,当初要是用def或def2,特别是计算能量时候用def2-TZVP,哪有这麻烦?

网上答疑时,对方用B3LYP/6-311++G**这个垃圾组合(精度又差、无意义地乱加弥散导致耗时大为增加)做有机体系计算的情况见过好几次,我总是苦口婆心地说别没事瞎加弥散(尤其是氢的)、算有机体系能量B3LYP过时了、建议改用def/def2,可对方还总是不以为然(理由往往是在之前的文章中见多了,或者是实验室祖传的,等等),似乎我这些建议都是多管闲事似的,这些人以后注定要吃大亏。像这种体系(几十个原子),用wB97XD或B3LYP-D3(BJ)结合6-311G**或def-TZVP优化,然后M06-2X或wB97M-V或双杂化泛函结合def2-TZVP算单点,妥妥的,耗时不很高精度又满足一般研究需要,哪个审稿人也挑不出毛病(真碰见有找茬的可以轻易怼回去)。

注:注意我只是说B3LYP算有机体系能量如今已经难以发表了,而优化是另一码事,见《谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的》(http://sobereva.com/557)。单点和几何优化完全可以用不同的泛函,见《量子化学研究中切换泛函应当注意的问题》(http://sobereva.com/415)、《浅谈为什么优化和振动分析不需要用大基组》(http://sobereva.com/387

评分 Rate

参与人数
Participants 43
eV +156 收起 理由
Reason
SharkYYX2025 + 2 233333
18257553618 + 3 赞!
cyh时代 + 4 精品内容
Emilia923 + 3 我很赞同
人间划水艺术家 + 5 GJ!
phur + 1 233333
hahaha333 + 3 谢谢
longqiyang + 4 好物!
啾啾啾0424 + 3 精品内容
石头砸鱼 + 2 谢谢
davi + 4 谢谢
caiqianyu + 4 正解
gengle + 5 赞!
Denny2716 + 4 谢谢
忘不掉的回忆 + 5 赞!
zhaoxiaoni + 5
chessrein + 4 谢谢
xiaoxu + 4 谢谢分享
玉米猫 + 3 我很赞同
drosin + 2

查看全部评分 View all ratings

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

76

帖子

0

威望

622

eV
积分
698

Level 4 (黑子)

41#
发表于 Post on 2025-11-23 13:05:21 | 只看该作者 Only view this author

好滴,谢谢老师

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

40#
 楼主 Author| 发表于 Post on 2025-11-23 00:56:04 | 只看该作者 Only view this author
cyh时代 发表于 2025-11-22 22:36
方便问一下卢老师,我最近用高斯计算了一下mobius碳纳米带,84个原子也计算了同类型的碳纳米管,用的是wb97 ...

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

76

帖子

0

威望

622

eV
积分
698

Level 4 (黑子)

39#
发表于 Post on 2025-11-22 22:36:01 | 只看该作者 Only view this author
方便问一下卢老师,我最近用高斯计算了一下mobius碳纳米带,84个原子也计算了同类型的碳纳米管,用的是wb97xd/def2-svp优化是否合理,本来想着用def-tzvp优化,我这电脑没跑动,就用了这个,http://sobereva.com/387看你这帖子说了,可以用,因为这个帖子上面都是推荐的是3zeta基组,不知道是否可行

71

帖子

0

威望

689

eV
积分
760

Level 4 (黑子)

38#
发表于 Post on 2023-1-5 23:31:58 | 只看该作者 Only view this author
收藏了

310

帖子

3

威望

6404

eV
积分
6774

Level 6 (一方通行)

37#
发表于 Post on 2020-7-27 13:56:14 | 只看该作者 Only view this author
coolrainbow 发表于 2020-7-27 11:10
ORCA算Pople基组时,那种SP基组的积分计算没有优化过,相当于分别算SX和PX,所以效率很低。

是啊,张老师最近在忙啥?
欢迎使用量子化学软件Amesp

249

帖子

13

威望

3593

eV
积分
4102

Level 6 (一方通行)

36#
发表于 Post on 2020-7-27 11:11:52 | 只看该作者 Only view this author
唉,有的人是另一个极端,拿个100原子的体系二话不说MP2/ccpvtz就几何优化,白白浪费一堆机时还算不出来

249

帖子

13

威望

3593

eV
积分
4102

Level 6 (一方通行)

35#
发表于 Post on 2020-7-27 11:10:17 | 只看该作者 Only view this author
Warm_Cloud 发表于 2019-4-12 18:32
pople基组确实该退休了,很不解的是用ORCA还用pople基组,耗时比def2更多(主要是6-31g**和def2-SVP)而且精 ...

ORCA算Pople基组时,那种SP基组的积分计算没有优化过,相当于分别算SX和PX,所以效率很低。

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

34#
 楼主 Author| 发表于 Post on 2020-7-26 20:45:09 | 只看该作者 Only view this author
TANGYIHAN 发表于 2020-7-26 19:26
社长,有机体系优化过渡态,使用M06-2X有没有必要加D3呢

可加可不加,加了没坏处但也没显著改进
北京科音自然科学研究中心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

15

帖子

0

威望

707

eV
积分
722

Level 4 (黑子)

33#
发表于 Post on 2020-7-26 19:26:19 | 只看该作者 Only view this author
社长,有机体系优化过渡态,使用M06-2X有没有必要加D3呢

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

32#
 楼主 Author| 发表于 Post on 2020-6-11 23:52:16 | 只看该作者 Only view this author
leo1992 发表于 2020-6-11 22:45
社长,计算有机体系的热力学数据,1.用wB97XD/6-311g**(用B3LYP一直有大的虚频),结构优化振动分析获得 ...

1 合理
2 不用
北京科音自然科学研究中心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

18

帖子

0

威望

648

eV
积分
666

Level 4 (黑子)

31#
发表于 Post on 2020-6-11 22:45:28 | 只看该作者 Only view this author
sobereva 发表于 2020-2-3 04:51
分清楚M06和M06-2X
不要随意把B3LYP简写,仔细看
计算化学中的一些常见不良写法和用词

社长,计算有机体系的热力学数据,1.用wB97XD/6-311g**(用B3LYP一直有大的虚频),结构优化振动分析获得热力学校正量,再用M06-2X/Def2TZVPP算单点,是否合理?
2.M06-2X/Def2TZVPP计算单点时是否需要加em=D3?谢谢老师。

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

30#
 楼主 Author| 发表于 Post on 2020-2-3 04:51:39 | 只看该作者 Only view this author
loovfnd 发表于 2020-2-1 11:11
感谢社长推荐这个帖子。有个疑问,如果用B3LYP优化结构,用M06-2X计算能量(单点能)呢。若这个可行用零点 ...

分清楚M06和M06-2X
不要随意把B3LYP简写,仔细看
计算化学中的一些常见不良写法和用词
http://sobereva.com/298http://bbs.keinsci.com/thread-1358-1-1.html

此问题在此文说了
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415http://bbs.keinsci.com/thread-9878-1-1.html
北京科音自然科学研究中心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

98

帖子

0

威望

2481

eV
积分
2579

Level 5 (御坂)

29#
发表于 Post on 2020-2-1 11:11:56 | 只看该作者 Only view this author
感谢社长推荐这个帖子。有个疑问,如果用B3LYP优化结构,用M06-2X计算能量(单点能)呢。若这个可行用零点校正获得自由能可以用B3的零点加上M06的单点吗?

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

28#
 楼主 Author| 发表于 Post on 2019-5-23 00:39:16 | 只看该作者 Only view this author
胡桃木 发表于 2019-5-22 16:22
请问sob老师,用b3lyp做优化,然后用m062x做激发态计算,这样可以被接受吗

如果你是指B3LYP优化基态,没问题
北京科音自然科学研究中心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

4

帖子

0

威望

158

eV
积分
162

Level 3 能力者

27#
发表于 Post on 2019-5-22 16:22:27 | 只看该作者 Only view this author

请问sob老师,用b3lyp做优化,然后用m062x做激发态计算,这样可以被接受吗

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

GMT+8, 2026-1-25 03:01 , Processed in 0.186589 second(s), 22 queries , Gzip On.

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