计算化学公社

标题: 请教MS中castep优化构型的问题 [打印本页]

作者
Author: 菩城 时间: 2016-7-31 21:12
标题: 请教MS中castep优化构型的问题
1：castep中提供好几种计算方法，请问有相关文献说明什么方法适合算什么样体系吗？

2：我算的是类似金刚石结构的体系，是需要先收敛测试，再结构优化吗？

3：收敛测试是不是说选择计算能量，然后一步步调节K-point和E-cut至能量变化很小？

4：优化构型是否一定要优化晶格参数？

请懂的朋友解答一下，感谢！

作者
Author: 卡开发发 时间: 2016-8-1 10:46
1、有，泛函的选择暂时可以看Truhlar的一些测试方面的文献，简单点来说可以这样选择：
（1）PBE，各方面较为均衡（或者说啥都能算但算啥都不是最好），最为常用，PW91的结果与之接近但稍有差异；
（2）RPBE较为适合计算表面体系，不过高估晶格常数是常有的事情，不是很推荐使用；
（3）PBEsol和Wu-Cohen比较适合处理晶格常数，PBEsol使用更广泛，个人觉得WC的文献中测试例子太少；但对于某些情况这两种泛函不见得比PBE来得更好，如金属体系可能不如LSDA；
（4）杂化泛函的使用大概就是求得为了一个好的带隙，HSE06是个好的选择，不过非得使用模守恒赝势加上杂化泛函本身计算慢，即便不用All bands/EDFT，也慢得相当惊人；
额外的修正：
（1）DFT-D的修正这个就不说了；
（2）DFT+U的修正，用于修正d、f电子自作用改善带隙，但对于d、f全满、半满系统没有改善，强行使用可能会得到非物理的结果，计算量比杂化泛函小；U值的确定严格来说需要线性响应方法确定，CASTEP暂时不支持。

2、当然要先测试再优化，否则截断能、k点等怎么确定呢？

3、对，但个人建议除了能量最好还要参考力和压力的收敛，尤其是变胞优化；如果涉及到性质计算，可能还得考虑所感兴趣的性质的收敛情况；

4、晶格参数是晶体构型的一个指标，对于固体的计算显然应该考虑；有些情况没必要优化结构（如计算一个分子，得把分子放到一个大的格子中，或者切出的表面等），这种情况不需要优化，但某些1D、2D系统只是难于直接优化晶格，这得想想其他的办法。

作者
Author: 菩城 时间: 2016-8-1 11:27

卡开发发发表于 2016-8-1 10:46
1、有，泛函的选择暂时可以看Truhlar的一些测试方面的文献，简单点来说可以这样选择：
（1）PBE，各方面较 ...

赞！很详细！在网上也经常看到您的回复，很受益，感谢！

作者
Author: lwj 时间: 2016-10-12 14:14

卡开发发发表于 2016-8-1 10:46
1、有，泛函的选择暂时可以看Truhlar的一些测试方面的文献，简单点来说可以这样选择：
（1）PBE，各方面较 ...

请问老师，测试K点和截断能选用精度需要很高还是很低就行了？

作者
Author: 卡开发发 时间: 2016-10-12 19:00

lwj 发表于 2016-10-12 14:14
请问老师，测试K点和截断能选用精度需要很高还是很低就行了？

测试稍高点没关系，低了的话对某些测试结果可能不可靠。

作者
Author: lwj 时间: 2016-10-12 21:37

卡开发发发表于 2016-10-12 19:00
测试稍高点没关系，低了的话对某些测试结果可能不可靠。

谢谢老师。

作者
Author: alonewolfyang 时间: 2019-6-6 21:15
老师，针对用MS建模，我只是用来建模，构型优化和其他计算均用VASP，还需要以上繁琐的测试吗？

作者
Author: 如遇古剑 时间: 2020-7-18 17:16

卡开发发发表于 2016-8-1 10:46
1、有，泛函的选择暂时可以看Truhlar的一些测试方面的文献，简单点来说可以这样选择：
（1）PBE，各方面较 ...

请问老师收敛测试是怎么做的呢？

作者
Author: 卡开发发 时间: 2020-7-19 02:15

如遇古剑发表于 2020-7-18 17:16
请问老师收敛测试是怎么做的呢？

现有大部分人的标准采用的是50eV动能截断的间隔增加使得能量变化小于1meV/atom，当然这不是唯一的衡量方式。

作者
Author: 如遇古剑 时间: 2020-7-19 09:16

卡开发发发表于 2020-7-19 02:15
现有大部分人的标准采用的是50eV动能截断的间隔增加使得能量变化小于1meV/atom，当然这不是唯一的衡量方 ...

谢谢卡开发发老师的指点！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)