计算化学公社

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

[CP2K] 关于晶胞优化能量收敛但是应力张量不收敛的问题求助

[复制链接 Copy URL]

21

帖子

0

威望

236

eV
积分
257

Level 3 能力者

各位老师好,本人在使用CP2K做施主-受主掺杂钛酸钡体系的晶胞优化时,出现了能量收敛但达到迭代步数之前应力张量始终达不到收敛精度(100bar)的问题。选用的初始结构是geo_opt之后的。请问此时利用cell_opt优化得到的晶胞结构足够可信吗?此外,用这个优化后的结构进行NPT系综下的BOMD模拟始终报错Diagolisation Failed,请问应该如何解决呢?非常感谢各位老师!附件为输入文件,输出文件由于体积过大我只截取关键部分截图。非常感谢各位老师!!

能量的收敛情况.png (46.89 KB, 下载次数 Times of downloads: 8)

能量随迭代步数的收敛情况

能量随迭代步数的收敛情况

初始结构.png (439.45 KB, 下载次数 Times of downloads: 7)

初始结构示意图(锂铝掺杂的钛酸钡)

初始结构示意图(锂铝掺杂的钛酸钡)

达到迭代步数后最后一步的收敛情况.png (124.12 KB, 下载次数 Times of downloads: 7)

达到迭代步数极限前最后一步的收敛信息

达到迭代步数极限前最后一步的收敛信息

1.inp

2.49 KB, 下载次数 Times of downloads: 5

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2021-12-3 01:43:49 | 只看该作者 Only view this author
把晶胞尺寸、压力关系绘制成曲线看看
也可以考虑对晶胞尺寸做个扫描,看压力怎么变化的,取个可以接受的值
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 15:05 , Processed in 0.255440 second(s), 23 queries , Gzip On.

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