|
仔仔细细阅读此文N遍 谈谈量子化学中基组的选择 http://sobereva.com/336 你用比该用的更大的基组审稿人肯定不会批你,因为对结果合理性无害,但真正懂行的审稿人只会这么想:作者水平不行,根本不懂如何划算地选择基组,用这么大基组真是瞎糟蹋计算资源 |
那丢失的 发表于 2016-12-15 09:23 (不了解你们的体系,即使有必要加弥散)基组用6-311+G(d)足以,且3-zeta用B3LYP有些掉价,用M06-2X/ma-TZVP可能高效很多。 |
liyuanhe211 发表于 2016-12-15 09:15 好的,谢谢老师,b3lyp/6-311+g(3df)我们用来看结构,能量用的CBS-QB3,我去查查然后自己做一做计算,有什么问题再请教。 |
那丢失的 发表于 2016-12-15 09:02 自己找个简单反应计算能量差比较一下就知道了、查benchmark也行,不用别人说。DFT根本用不到那么多极化,当前结果本身没问题(这也可能是审稿人不批评的原因)但浪费大量时间、或者说用同等时间本可以获得更好结果,即常说的基组不匹配,。 |
sobereva 发表于 2016-12-14 18:47 结构优化用b3lyp/6-311+g(3df),如果含氢元素则用b3lyp/6-311++g(3df,3pd),算能量用CBS-QB3,现在可以用CCSD(T)了,这都是都是参考之前文章的,也没有审稿人告诉我们说你们不对,我做的都是十原子以下小分子,得到计算结果也与实验结果吻合比较好,所以一直在用,我主攻实验,所以对计算这块没花心思,现在看来还有好长的路要走。NBO的计算结果还是和老板讨论下让老板找他们问问为什么不一样吧,谢谢Sob老师. |
liyuanhe211 发表于 2016-12-14 18:50 看来错的人很多,包括大牛,哈哈,b3lyp/6-311+g(3df)是沿用之前文章的,如果含氢元素则用b3lyp/6-311++g(3df,3pd),也没有审稿人说不对的。 |
那丢失的 发表于 2016-12-14 09:26 误把spin population写成spin density的文献想找有的是,错了就是错了,有文献这么写也是错的。我记得还有IF甚高的文献用高斯算了CCSDT呢。 Σ(spin density)的写法就更奇怪了(虽然大家知道你什么意思),对一个空间分布的三维函数做“求和”。。。 b3lyp/6-311+g(3df)是非常奇怪的组合 |
那丢失的 发表于 2016-12-14 09:26 b3lyp/6-311+g(3df)做优化是不划算的,对于普通有机体系,结果不会比用6-311G**优化更强,耗时还多得多得多。而且由于没考虑氢的极化,牵扯到氢的任务结果会明显更糟 仔细看 谈谈量子化学中基组的选择 http://sobereva.com/336 至于你1L的体系,鉴于体系很小,用def2-TZVP优化耗时不会很高还很理想,比起用pople基组要强。 先把所有计算条件都保证相同了再去比结果,而且注意用的NBO版本是否相同。如果所有地方都相同,你的结果就是和他们的不同,直接管他们要输入输出文件,都是合作者,没道理不给。 不要以为JACS上的计算和说法都对,spin density and NBO atomic charges这种错误说法在哪里都常见,但你绝对不会看到NBO开发者的文章里有NBO atomic charges这种用词,这都是NBO用户自己造出来的词。把原子自旋布居说成spin density在前面的博文专门批了,极为常见的不良用词,体现出作者根本没搞清楚基本概念。 |
那丢失的 发表于 2016-12-13 20:44 我才不信作者用的是UCCSD(T)-F12//cc-pVTZ-F12做的优化 哪有拿这么高级别还没有解析梯度的做优化的 分明是作者没表述清楚,根本就不懂//这个符号的使用习俗 按他那么写成了CCSD(T)优化,DFT算单点,岂不搞笑 而且作者缺乏基本知识,根本没有NBO电荷这个说法,也不能说原子自旋密度,犯了常识性,看 量子化学中的一些常见不良写法和用词 http://sobereva.com/298 另外作者用的是def2-TZVPP你用的是def2-TZVP 几何结构不同、基组不同,都会导致NPA自旋布居不同 |
liyuanhe211 发表于 2016-12-13 20:24 不是啊,UCCSD(T)-F12/cc-pVTZ-F12优化结构并得到单点能,这是一个人做的计算,另外一个人只不过是拿他优化好的这个结构用UBP86/def2-TZVPP做NBO分析 |
那丢失的 发表于 2016-12-13 19:30 我就是说这作者要么缺乏常识不知道“ // ”符号是优化在后单点在前、要么缺乏常识不知道优化用小级别、单点用高级别 |
liyuanhe211 发表于 2016-12-13 19:25 结构优化用的UCCSD(T)-F12/cc-pVTZ-F12,得到的结构再用UBP86/def2-TZVPP做的NBO分析。我也是直接用UCCSD(T)-F12/cc-pVTZ-F12得到的结构做的计算 |
那丢失的 发表于 2016-12-13 18:28 这作者能写出来“UBP86/def2-TZVPP//UCCSD(T)-F12/cc-pVTZ-F12”这种表述就不怎么可信,还拿。。。 以及UCCSD(T)-F12显然不是用高斯做的,因为它这奇怪的表述也不清楚这些分析到底是在前后哪个水平做的,说不定是后者 |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 | 我很赞同 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 10:43 , Processed in 0.206251 second(s), 26 queries , Gzip On.