计算化学公社

标题: CASTEP计算吸附时的对称性问题 [打印本页]

作者
Author:
tessal    时间: 2018-11-2 10:28
标题: CASTEP计算吸附时的对称性问题
新手一枚,请教下各位,我用castep计算吸附能,先结构优化再单点能计算,结构优化的时候总是会提醒是否需要提高对称性以加快速度,我试过提高跟不提高的情况,计算出来的能量差别很大,请问大家一般需不需要提高对称性计算呢?什么情况下需要提高对称性呢?谢谢大家

作者
Author:
卡开发发    时间: 2018-11-3 18:59
1、关于能量相差很大,你得先确定一下你转化成primitive cell之后原子数目和未转化或未使用对称性之前是否一致,如果原子数不一样的话,两者应该是差原子数目比例的整数倍。
2、如果这样进行换算能量差别仍然不合理,可能需要提高计算精度,这可能是数值误差导致的。
3、如果(能够确定)体系本身确实有对称性,当然主张采用对称性来进行加速,但如果体系没有对称性而被强制要求对称会导致一些问题。比如规则的晶体的bulk确实有好的对称性,优化晶格结构的时候就应当考虑对称性,但要注意磁序对空间对称性的破坏,MS对这种情况没办法识别(我不清楚新的版本如何)。
作者
Author:
tessal    时间: 2018-11-6 09:45
卡开发发 发表于 2018-11-3 18:59
1、关于能量相差很大,你得先确定一下你转化成primitive cell之后原子数目和未转化或未使用对称性之前是否 ...

谢谢卡开发发,我看了下,相差的能量确实是原子数目差的整数倍。还想请教下,那如果我是计算氮气在氧化银(111)面的吸附能,这种情况下需要使用对称性吗?我看到CASTEP自带的例子中,CO在Pd(110)面的吸附,由于CO的加入,使得Pd的对称性改变,所以要加入对称性,我不明白这个对称性加入的意义是什么?是加快计算吗?
作者
Author:
卡开发发    时间: 2018-11-6 14:44
tessal 发表于 2018-11-6 09:45
谢谢卡开发发,我看了下,相差的能量确实是原子数目差的整数倍。还想请教下,那如果我是计算氮气在氧化银 ...

加快计算的同时也造成了自由度的约束,但这样的约束到底是否有真实的物理意义就得凭你自己的判断了。
作者
Author:
tessal    时间: 2018-11-7 17:27
卡开发发 发表于 2018-11-6 14:44
加快计算的同时也造成了自由度的约束,但这样的约束到底是否有真实的物理意义就得凭你自己的判断了。

好的,万分感谢




欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3