Views: 15327|回复 Reply: 10
[CASTEP/Dmol3/MS]
求助过渡组金属氧化物CASTEP几何优化速度太慢怎么办？
[复制链接 Copy URL]
XWH
XWH
当前离线
UID26350
在线时间124 小时
注册时间2021-3-9
最后登录2023-4-28
好友0
71
帖子0
威望356
eV
Level 3 能力者, 积分 credits 427, 距离下一级还需 73 积分 credits
积分427
Level 3 能力者
发消息 Send PM
电梯直达 Go to
楼主
发表于 Post on 2021-4-16 12:56:45
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
2*2*2超胞，96个原子。采用GGA+U,PBE方法，电脑配置、操作和结果如下所示，已经跑了一天多了，还没有收敛的意思，请问这是这怎么回事？还有优化进行下去的必要吗是因为K点和截断能设置的不够合理吗？有没有什么经验可以预估计算速度啊？麻烦看到的各位给点建议！感谢！感谢！感谢！
GZ{HPCC0`6~~PP%F[F9~18V.png
(11.58 KB, 下载次数 Times of downloads: 67)
下载附件 Download
2021-4-16 12:56 上传 Uploaded
[NP`QBIEJ1IFWJ]N@ISN7DI.png
(40.72 KB, 下载次数 Times of downloads: 65)
下载附件 Download
2021-4-16 12:54 上传 Uploaded
44X(AG8NQ3(GZUCL}VKJYWU.png
(108.43 KB, 下载次数 Times of downloads: 50)
下载附件 Download
2021-4-16 12:54 上传 Uploaded
17A]}9R[XS){I5~KR6}U64O.png
(16.86 KB, 下载次数 Times of downloads: 66)
下载附件 Download
2021-4-16 12:53 上传 Uploaded
%)AKN9SHRN]7KTON((L59PN.png
(19.25 KB, 下载次数 Times of downloads: 61)
下载附件 Download
2021-4-16 12:53 上传 Uploaded
]SBZWFK44XIWBAO_)IS5S`V.png
(40.97 KB, 下载次数 Times of downloads: 47)
下载附件 Download
2021-4-16 12:53 上传 Uploaded
收藏 Add to favorites2
回复 Reply
举报 Report
喵星大佬
喵星大佬
当前离线
UID19855
在线时间2275 小时
注册时间2020-5-14
最后登录2025-8-15
好友10
1657
帖子5
威望4562
eV
积分6319
Level 6 (一方通行)
喵星人
发消息 Send PM
2#
发表于 Post on 2021-4-16 13:11:42
|
只看该作者 Only view this author
优化任务一定要把截断能和SCF收敛限设置到fine以上。过渡金属氧化物一般不需要开Mental，但是要注意自旋极化和设置初始磁矩之类的。也不知道你的是什么金属
回复 Reply
举报 Report
XWH
XWH
当前离线
UID26350
在线时间124 小时
注册时间2021-3-9
最后登录2023-4-28
好友0
71
帖子0
威望356
eV
Level 3 能力者, 积分 credits 427, 距离下一级还需 73 积分 credits
积分427
Level 3 能力者
发消息 Send PM
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的，可是结构优化以后本该是有能带间隙的，显示没有。文献中给出的模拟相关信息细节比较少，太多不懂
回复 Reply
举报 Report
XWH
XWH
当前离线
UID26350
在线时间124 小时
注册时间2021-3-9
最后登录2023-4-28
好友0
71
帖子0
威望356
eV
Level 3 能力者, 积分 credits 427, 距离下一级还需 73 积分 credits
积分427
Level 3 能力者
发消息 Send PM
4#
楼主 Author|
发表于 Post on 2021-4-16 14:56:27
|
只看该作者 Only view this author
看到帖子各位前辈，麻烦会的给点建议，我回来仔细想想，感谢！感谢！感谢！
回复 Reply
举报 Report
sobereva
sobereva
当前离线
UID1
在线时间6122 小时
注册时间2014-10-7
最后登录2025-8-17
好友67
6万
帖子99
威望5万
eV
积分120181
管理员
公社社长
发消息 Send PM
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
回复 Reply
举报 Report
泡泡媛
泡泡媛
当前离线
UID12989
在线时间2127 小时
注册时间2019-4-2
最后登录2025-8-17
好友0
45
帖子0
威望4309
eV
积分4354
Level 6 (一方通行)
发消息 Send PM
6#
发表于 Post on 2021-4-21 20:15:44
|
只看该作者 Only view this author
本帖最后由 泡泡媛 于 2021-4-23 13:18 编辑
有没有做过收敛性测试，就是k点和截断能的测试，我觉得是你这两个选择不合理，看着能量曲线下降了点，但是其他的标准都震荡的很厉害
回复 Reply
举报 Report
XWH
XWH
当前离线
UID26350
在线时间124 小时
注册时间2021-3-9
最后登录2023-4-28
好友0
71
帖子0
威望356
eV
Level 3 能力者, 积分 credits 427, 距离下一级还需 73 积分 credits
积分427
Level 3 能力者
发消息 Send PM
7#
楼主 Author|
发表于 Post on 2021-4-24 21:21:26
|
只看该作者 Only view this author
泡泡媛 发表于 2021-4-21 20:15
有没有做过收敛性测试，就是k点和截断能的测试，我觉得是你这两个选择不合理，看着能量曲线下降了点，但是 ...
谢谢，我再去看看收敛性测试
回复 Reply
举报 Report
XWH
XWH
当前离线
UID26350
在线时间124 小时
注册时间2021-3-9
最后登录2023-4-28
好友0
71
帖子0
威望356
eV
Level 3 能力者, 积分 credits 427, 距离下一级还需 73 积分 credits
积分427
Level 3 能力者
发消息 Send PM
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.可以给些建议吗？
回复 Reply
举报 Report
sobereva
sobereva
当前离线
UID1
在线时间6122 小时
注册时间2014-10-7
最后登录2025-8-17
好友67
6万
帖子99
威望5万
eV
积分120181
管理员
公社社长
发消息 Send PM
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
回复 Reply
举报 Report
get-it
get-it
当前离线
UID2981
在线时间1621 小时
注册时间2016-4-27
最后登录2025-8-15
好友2
236
帖子0
威望5069
eV
积分5305
Level 6 (一方通行)
发消息 Send PM
10#
发表于 Post on 2021-4-25 14:48:27
来自手机
|
只看该作者 Only view this author
最简单可靠的方法，查文献，看别人算VO2用什么参数
回复 Reply
举报 Report
Salina
Salina
当前离线
UID68133
在线时间0 小时
注册时间2024-9-6
最后登录2024-9-27
好友0
2
帖子0
威望27
eV
Level 2 能力者, 积分 credits 29, 距离下一级还需 121 积分 credits
积分29
Level 2 能力者
发消息 Send PM
11#
发表于 Post on 2024-9-13 21:12:39
|
只看该作者 Only view this author
后来怎么解决的呢 友友 我现在在做掺杂了元素的 弹性参数测试 也是一直不收敛 我已经将SCF cycle调到2500，charge跳到了0.02 依旧不收敛
回复 Reply
举报 Report