计算化学公社

标题: 求助CP2K优化结构中原子析出的问题 [打印本页]

作者
Author: 萧瑟 时间: 2025-12-11 10:59
标题: 求助CP2K优化结构中原子析出的问题
本帖最后由萧瑟于 2025-12-11 11:01 编辑

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

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

作者
Author: sobereva 时间: 2025-12-11 17:32
图片描述具体，说明哪个是你说的“脱出的现象”的图
第一个图的结构没问题

作者
Author: 萧瑟 时间: 2025-12-12 11:25

sobereva 发表于 2025-12-11 17:32
图片描述具体，说明哪个是你说的“脱出的现象”的图
第一个图的结构没问题

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

作者
Author: sobereva 时间: 2025-12-13 04:02

萧瑟发表于 2025-12-12 11:25
好的社长，第一个是优化前的初态是相对正常的，第二张是优化完成后的末态单个Cd原子就往上远离晶格结 ...

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

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

(, 下载次数 Times of downloads: 0)

作者
Author: 萧瑟 时间: 2025-12-15 15:24

sobereva 发表于 2025-12-13 04:02
监控优化过程，看SCF是否有异常

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)