计算化学公社

标题: 求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办？ [打印本页]

作者
Author: XWH 时间: 2021-4-16 12:56
标题: 求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办？
2*2*2超胞，96个原子。采用GGA+U,PBE方法，电脑配置、操作和结果如下所示，已经跑了一天多了，还没有收敛的意思，请问这是这怎么回事？还有优化进行下去的必要吗是因为K点和截断能设置的不够合理吗？有没有什么经验可以预估计算速度啊？麻烦看到的各位给点建议！感谢！感谢！感谢！

作者
Author: 喵星大佬 时间: 2021-4-16 13:11
优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental，但是要注意自旋极化和设置初始磁矩之类的。也不知道你的是什么金属

作者
Author: XWH 时间: 2021-4-16 14:53

喵星大佬发表于 2021-4-16 13:11
优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental，但是要注意自旋极 ...

谢谢，我模拟的是VO2热至相变，这种金属氧化物341k会发生相变，由绝缘性变为金属性，这样也不需要勾选metal吗？我之前试过660ev,k点1，不勾选metal的，可是结构优化以后本该是有能带间隙的，显示没有。文献中给出的模拟相关信息细节比较少，太多不懂

作者
Author: XWH 时间: 2021-4-16 14:56
看到帖子各位前辈，麻烦会的给点建议，我回来仔细想想，感谢！感谢！感谢！

作者
Author: sobereva 时间: 2021-4-21 17:01
起码先看看收敛趋势，弄清楚距离收敛大概还有多久，没有收敛希望就早点停了

而且如果研究热至相变，也应该是跑AIMD，除非当前优化的就是相变后的结构

作者
Author: 泡泡媛 时间: 2021-4-21 20:15
本帖最后由泡泡媛于 2021-4-23 13:18 编辑

有没有做过收敛性测试，就是k点和截断能的测试，我觉得是你这两个选择不合理，看着能量曲线下降了点，但是其他的标准都震荡的很厉害

作者
Author: XWH 时间: 2021-4-24 21:21

泡泡媛发表于 2021-4-21 20:15
有没有做过收敛性测试，就是k点和截断能的测试，我觉得是你这两个选择不合理，看着能量曲线下降了点，但是 ...

谢谢，我再去看看收敛性测试

作者
Author: XWH 时间: 2021-4-24 21:38

sobereva 发表于 2021-4-21 17:01
起码先看看收敛趋势，弄清楚距离收敛大概还有多久，没有收敛希望就早点停了

而且如果研究热至相变，也应 ...

您好，我之前模拟这个体系用得是lammps，但是这个没有合适的势函数，对于MS,这种热至相变采用离散式升温，模拟300K,400K,500K,600K下弛豫，前一个温度下跑完的模型是后一个温度下的开始，以这种模拟方式观察相变可行吗？我最后要模拟的是镀膜后材料之间的热应力变化（需要算各层材料的弹性模量，泊松比，热膨胀系数这些），用MS的话不知道这个想法可以吗？还是要用vasp比较好？我目前只接触过lammps\MS\VASP.可以给些建议吗？

作者
Author: sobereva 时间: 2021-4-25 00:07

XWH 发表于 2021-4-24 21:38
您好，我之前模拟这个体系用得是lammps，但是这个没有合适的势函数，对于MS,这种热至相变采用离散式升温 ...

CP2K最好。不仅其DFT是周期性体系程序中最快的之一（对较大体系远快于CASTEP、VASP），而且还可以尝试用里面的GFN1-xTB，能跑的空间和时间尺度比DFT还大得多

作者
Author: get-it 时间: 2021-4-25 14:48
最简单可靠的方法，查文献，看别人算VO2用什么参数

作者
Author: Salina 时间: 2024-9-13 21:12
后来怎么解决的呢友友我现在在做掺杂了元素的弹性参数测试也是一直不收敛我已经将SCF cycle调到2500，charge跳到了0.02 依旧不收敛

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