计算化学公社

标题: 求助Mo2C晶胞优化后晶格参数变化过大 [打印本页]

作者
Author: zjf1240291644 时间: 2022-9-12 20:14
标题: 求助Mo2C晶胞优化后晶格参数变化过大
计算对象：Mo2C
task：Geomertry Optimiza
计算参数：GGA PBE
(, 下载次数 Times of downloads: 22)

优化前晶格参数如下
(, 下载次数 Times of downloads: 23)

优化后晶格参数如下
(, 下载次数 Times of downloads: 21)

求助为什么晶格参数会变动这么大

作者
Author: 卡开发发 时间: 2022-9-12 20:24
DMol3优化晶胞的方法比较拉胯，建议慎重使用。我个人平时优化晶格只考虑通过变量扫描的办法，虽然繁琐但确实没啥好办法。

作者
Author: zjf1240291644 时间: 2022-9-12 20:33

卡开发发发表于 2022-9-12 20:24
DMol3优化晶胞的方法比较拉胯，建议慎重使用。我个人平时优化晶格只考虑通过变量扫描的办法，虽然繁琐但确 ...

请问变量扫描是什么方法，可以推荐下相关的资料学习下吗

作者
Author: 卡开发发 时间: 2022-9-12 22:11

zjf1240291644 发表于 2022-9-12 20:33
请问变量扫描是什么方法，可以推荐下相关的资料学习下吗

假定结构是立方格子，那么其实可以简单扫描晶胞长度来找到能量最低的，复杂的情况确实没什么特别好的办法，要么使用单纯形法或者通过差分算stess。我自己一般不会用扫描的方法，而是直接将晶胞长度作为一维变量做一维线搜索，其中效率低下的是二分法，也可以用抛物线逼近，相关内容你可以参考Numerical Receipts。上述实现基本上可能得靠自己写程序来做。客观说，DMol3的晶胞优化功能确实不太行。

作者
Author: zjf1240291644 时间: 2022-9-12 22:35

卡开发发发表于 2022-9-12 22:11
假定结构是立方格子，那么其实可以简单扫描晶胞长度来找到能量最低的，复杂的情况确实没什么特别好的办法 ...

感谢大佬

作者
Author: 脉冲星射线暴 时间: 2022-9-12 23:30
本帖最后由脉冲星射线暴于 2022-9-12 23:33 编辑

试过quality调高一些吗？选的赝势合适不？
不行的话还是用QE的vc-relax计算比较好，可选的赝势更多。至少根据我师兄的话，MS集成的很多模块都过时了，比如CASTEP。

作者
Author: sobereva 时间: 2022-9-15 05:39
推荐用CP2K做，完全免费，速度很快，stress可以解析计算，结合Multiwfn产生输入文件非常方便
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

作者
Author: zjf1240291644 时间: 2022-9-15 10:14

sobereva 发表于 2022-9-15 05:39
推荐用CP2K做，完全免费，速度很快，stress可以解析计算，结合Multiwfn产生输入文件非常方便
使用Multiwfn ...

谢谢推荐

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