计算化学公社

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

[CP2K] 求助CP2K优化结构中原子析出的问题

[复制链接 Copy URL]

77

帖子

0

威望

999

eV
积分
1076

Level 4 (黑子)

本帖最后由 萧瑟 于 2025-12-11 11:01 编辑

如题,我在计算以CdS为基础的一些结构模型时,偶尔会出现表面的Cd原子脱出的现象,之前还出现过靠中间的Cd原子优化完往上飘起来的情况一个孤立原子在真空中。计算使用的inp文件是用multiwfn生成的,参数变动的地方不多,现在附上相关文件,想请教各位老师,出现这种情况的原因有哪位老师知道是因为什么吗?目前都是Cd元素出现这种问题,在inp文件中需要调整哪些关键词和参数可以解决这种问题?



ps:之前的算例找不到了   昨天又出现这种类似的情况,想来请教下各位老师

屏幕截图 2025-12-11 105446.png (131.71 KB, 下载次数 Times of downloads: 3)

初

屏幕截图 2025-12-11 105457.png (106.27 KB, 下载次数 Times of downloads: 4)

末

ZCS111-half-HCOOH.cif

17.87 KB, 下载次数 Times of downloads: 1

ZCS111-half-HCOOH.inp

12.04 KB, 下载次数 Times of downloads: 2

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2025-12-11 17:32:26 | 只看该作者 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

77

帖子

0

威望

999

eV
积分
1076

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-12-12 11:25:06 | 只看该作者 Only view this author
sobereva 发表于 2025-12-11 17:32
图片描述具体,说明哪个是你说的“脱出的现象”的图
第一个图的结构没问题

好的社长,第一个是优化前的初态是相对正常的 ,第二张是优化完成后的末态   单个Cd原子就往上远离晶格结构,跑出去的Cd与最近的S原子距离超过3.5A了  应该没有化学成键作用了。

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

4#
发表于 Post on 2025-12-13 04:02:55 | 只看该作者 Only view this author
萧瑟 发表于 2025-12-12 11:25
好的社长,第一个是优化前的初态是相对正常的 ,第二张是优化完成后的末态   单个Cd原子就往上远离晶格结 ...

监控优化过程,看SCF是否有异常

另外,羧基结合到表面上,质子应该是解离掉的状态,不会连着吸附的氧。例如在TiO2上吸附:



北京科音自然科学研究中心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

77

帖子

0

威望

999

eV
积分
1076

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-12-15 15:24:48 | 只看该作者 Only view this author
sobereva 发表于 2025-12-13 04:02
监控优化过程,看SCF是否有异常

另外,羧基结合到表面上,质子应该是解离掉的状态,不会连着吸附的氧 ...

好的谢谢社长, 这个结构是中间态, 不过您说得我以后也会多注意

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 04:58 , Processed in 0.169703 second(s), 23 queries , Gzip On.

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