agent99 发表于 2018-1-6 02:38 谢谢啦,我试试 |
sobereva 发表于 2018-1-6 06:49 谢谢sob老师,我研究一下。 |
|
B3LYP对这种体系根本没有极小点,即便你说的优化结果发现和实验差不多,那也只是数值巧合,因为受力已经非常小了就当成优化收敛了,而不是因为B3LYP下极小点就在那里 算这个体系要拿来发表起码CCSD(T) 就算用DFT,D3校正在你的版本里不支持,起码也得用诸如wB97XD、M06-2X这样能描述色散作用的 算弱相互作用基本常识知识看 乱谈DFT-D http://sobereva.com/83 这篇文章专门是理论研究稀有气体二聚体的,你一定要看看:DOI 10.1002/qua.25358 |
|
这种体系的能量极小点应该很浅,你把曲线在极小点附近放大试试。 另外,B3LYP不加色散校正的话对这种很弱的范德华作用很可能无法描述。高斯版本旧没有色散校正的话,M06-2X总有吧?那个至少对弱相互作用描述好一些。 |
yflchx 发表于 2018-1-5 23:25 你好,能给点建议,算一下么,我看文献上加入了矫正,貌似还不错 |
|
本帖最后由 yflchx 于 2018-1-5 23:28 编辑 像这种色散主导的体系,B3LYP无能为力。 看曲线这么平坦,优化出个结果也算正常。 也许有人试图重新优化B3LYP的参数来做非共价相互作用(原理上是说不通的),问题会有两个: 1. 适合氢键体系的系数肯定不适合色散主导的相互作用体系。 2. 针对色散主导的体系优化参数,最终的平均误差大概会在0.5-0.7 kcal/mol左右,然后再也下不去了。 如果做类似的事,写文章投稿会遇到很大的阻力。即便是发出来也很可能被comment。 |
|
本帖最后由 liyuanhe211 于 2018-1-5 21:41 编辑 你觉得图哪儿有问题?以及不能仅凭键长差不多就说级别没问题 |
|
横坐标从2.9 或者3.0开始就好 纵坐标改成相对能量 单位改成KJ/mol |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-23 08:21 , Processed in 0.170190 second(s), 25 queries , Gzip On.