楼主第一个castep,没有成功吧,结果给出的是lst max found,而不是TS,所以差别大了 |
zyj19831206 发表于 2015-12-17 00:14 这里疑问点很多。 首先反应能方面做个解释,CASTEP反应能算出来0.02211eV而DMol3反应能算出来-0.691kcal/mol,这两者的误差在2kcal/mol以内。如果实际的反应能真的是这样,可能DFT-GGA的计算也根本无法进行定性,泛函引起的误差可能远大于这个。 我自己做表面比较多,CASTEP我很少用于计算过渡态,不过早上试了一下,发觉一些问题不妨讨论一下: LZ的过程是否是通过LST/QST的方法?TS的图上是否出现CG的过程?如果没出现CG的过程,LZ的计算精度如何?如果出现了CG,按说后续应该有CG的过程,能量应该是参考CG之后的,除此,如果LST/CG之后RMS force没收敛还应该进行QST/CG。 一般来说,CASTEP和DMol3小的误差会有一些,但不会差的很离谱(@jiewei 兄做了些数据的收集,不妨参考http://bbs.keinsci.com/forum.php ... hlight=%C9%AB%C9%A2) |
参与人数Participants 1 | eV +2 | 收起 理由Reason |
---|---|---|
sobereva | + 2 |
minagami 发表于 2015-12-16 23:45 是的,我也是这么理解的,两个模块搜索出问题这么多? |
感觉是参数设置有问题。即使不谈找过渡态,就反应的总自由能变而言,也不可能一个算出来是吸热另一个是放热呀。 |
本帖最后由 卡开发发 于 2015-12-16 22:36 编辑 啥反应,拿出来看看。 还有CASTEP参数经过测试否,截断能和KPOINTS还有SIGMA(Smearing)值大小都会影响force的计算,会造成不准确。另外DMol3的多极展开和Smearing也都有影响。除了这些,模守恒赝势和超软赝势的结果本身也还是有些差异的。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 07:00 , Processed in 0.258975 second(s), 29 queries , Gzip On.