计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办？

求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办？

查看数: 15791 | 评论数: 10 | 收藏 Add to favorites 2

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

XWH

发布时间: 2021-4-16 12:56

正文摘要:

2*2*2超胞，96个原子。采用GGA+U,PBE方法，电脑配置、操作和结果如下所示，已经跑了一天多了，还没有收敛的意思，请问这是这怎么回事？还有优化进行下去的必要吗是因为K点和截断能设置的不够合理吗？有没有什么经验 ...

回复 Reply

Salina 发表于 Post on 2024-9-13 21:12:39

后来怎么解决的呢友友我现在在做掺杂了元素的弹性参数测试也是一直不收敛我已经将SCF cycle调到2500，charge跳到了0.02 依旧不收敛

get-it 发表于 Post on 2021-4-25 14:48:27

最简单可靠的方法，查文献，看别人算VO2用什么参数

sobereva 发表于 Post on 2021-4-25 00:07:31

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

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

XWH 发表于 Post on 2021-4-24 21:38:30

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

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

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

XWH 发表于 Post on 2021-4-24 21:21:26

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

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

泡泡媛 发表于 Post on 2021-4-21 20:15:44

本帖最后由泡泡媛于 2021-4-23 13:18 编辑

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

sobereva 发表于 Post on 2021-4-21 17:01:56

起码先看看收敛趋势，弄清楚距离收敛大概还有多久，没有收敛希望就早点停了

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

XWH 发表于 Post on 2021-4-16 14:56:27

看到帖子各位前辈，麻烦会的给点建议，我回来仔细想想，感谢！感谢！感谢！

XWH 发表于 Post on 2021-4-16 14:53:52

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

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

喵星大佬 发表于 Post on 2021-4-16 13:11:42

优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental，但是要注意自旋极化和设置初始磁矩之类的。也不知道你的是什么金属

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register