楼主第一个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）

minagami 发表于 2015-12-16 23:45
感觉是参数设置有问题。即使不谈找过渡态，就反应的总自由能变而言，也不可能一个算出来是吸热另一个是放热 ...

是的，我也是这么理解的，两个模块搜索出问题这么多？

感觉是参数设置有问题。即使不谈找过渡态，就反应的总自由能变而言，也不可能一个算出来是吸热另一个是放热呀。

啥反应，拿出来看看。

还有CASTEP参数经过测试否，截断能和KPOINTS还有SIGMA（Smearing）值大小都会影响force的计算，会造成不准确。另外DMol3的多极展开和Smearing也都有影响。除了这些，模守恒赝势和超软赝势的结果本身也还是有些差异的。