计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15789|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[CASTEP/Dmol3/MS] 求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办?

[复制链接 Copy URL]

71

帖子

0

威望

356

eV
积分
427

Level 3 能力者

2*2*2超胞,96个原子。采用GGA+U,PBE方法,电脑配置、操作和结果如下所示,已经跑了一天多了,还没有收敛的意思,请问这是这怎么回事?还有优化进行下去的必要吗是因为K点和截断能设置的不够合理吗?有没有什么经验可以预估计算速度啊?麻烦看到的各位给点建议!感谢!感谢!感谢!


GZ{HPCC0`6~~PP%F[F9~18V.png (11.58 KB, 下载次数 Times of downloads: 68)

GZ{HPCC0`6~~PP%F[F9~18V.png

[NP`QBIEJ1IFWJ]N@ISN7DI.png (40.72 KB, 下载次数 Times of downloads: 67)

[NP`QBIEJ1IFWJ]N@ISN7DI.png

44X(AG8NQ3(GZUCL}VKJYWU.png (108.43 KB, 下载次数 Times of downloads: 54)

44X(AG8NQ3(GZUCL}VKJYWU.png

17A]}9R[XS){I5~KR6}U64O.png (16.86 KB, 下载次数 Times of downloads: 67)

17A]}9R[XS){I5~KR6}U64O.png

%)AKN9SHRN]7KTON((L59PN.png (19.25 KB, 下载次数 Times of downloads: 64)

%)AKN9SHRN]7KTON((L59PN.png

]SBZWFK44XIWBAO_)IS5S`V.png (40.97 KB, 下载次数 Times of downloads: 49)

]SBZWFK44XIWBAO_)IS5S`V.png

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2021-4-16 13:11:42 | 只看该作者 Only view this author
优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental,但是要注意自旋极化和设置初始磁矩之类的。也不知道你的是什么金属

71

帖子

0

威望

356

eV
积分
427

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-4-16 14:53:52 | 只看该作者 Only view this author
喵星大佬 发表于 2021-4-16 13:11
优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental,但是要注意自旋极 ...

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

71

帖子

0

威望

356

eV
积分
427

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2021-4-16 14:56:27 | 只看该作者 Only view this author
看到帖子各位前辈,麻烦会的给点建议,我回来仔细想想,感谢!感谢!感谢!

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

5#
发表于 Post on 2021-4-21 17:01:56 | 只看该作者 Only view this author
起码先看看收敛趋势,弄清楚距离收敛大概还有多久,没有收敛希望就早点停了

而且如果研究热至相变,也应该是跑AIMD,除非当前优化的就是相变后的结构
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

45

帖子

0

威望

4453

eV
积分
4498

Level 6 (一方通行)

6#
发表于 Post on 2021-4-21 20:15:44 | 只看该作者 Only view this author
本帖最后由 泡泡媛 于 2021-4-23 13:18 编辑

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

71

帖子

0

威望

356

eV
积分
427

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-4-24 21:21:26 | 只看该作者 Only view this author
泡泡媛 发表于 2021-4-21 20:15
有没有做过收敛性测试,就是k点和截断能的测试,我觉得是你这两个选择不合理,看着能量曲线下降了点,但是 ...

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

71

帖子

0

威望

356

eV
积分
427

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2021-4-24 21:38:30 | 只看该作者 Only view this author
sobereva 发表于 2021-4-21 17:01
起码先看看收敛趋势,弄清楚距离收敛大概还有多久,没有收敛希望就早点停了

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

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

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

9#
发表于 Post on 2021-4-25 00:07:31 | 只看该作者 Only view this author
XWH 发表于 2021-4-24 21:38
您好,我之前模拟这个体系用得是lammps,但是这个没有合适的势函数,对于MS,这种热至相变采用离散式升温 ...

CP2K最好。不仅其DFT是周期性体系程序中最快的之一(对较大体系远快于CASTEP、VASP),而且还可以尝试用里面的GFN1-xTB,能跑的空间和时间尺度比DFT还大得多
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

239

帖子

0

威望

5249

eV
积分
5488

Level 6 (一方通行)

10#
发表于 Post on 2021-4-25 14:48:27 来自手机 | 只看该作者 Only view this author
最简单可靠的方法,查文献,看别人算VO2用什么参数

2

帖子

0

威望

27

eV
积分
29

Level 2 能力者

11#
发表于 Post on 2024-9-13 21:12:39 | 只看该作者 Only view this author
后来怎么解决的呢 友友 我现在在做掺杂了元素的 弹性参数测试 也是一直不收敛 我已经将SCF cycle调到2500,charge跳到了0.02 依旧不收敛

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-19 22:11 , Processed in 0.206708 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list