计算化学公社

标题: 计算1800个原子（Pd/TiO2)体系的结构优化(含4%氧缺陷）时出现报错结构不收敛的问题 [打印本页]

作者
Author: 皮卡丘闪电链 时间: 2025-6-23 15:20
标题: 计算1800个原子（Pd/TiO2)体系的结构优化(含4%氧缺陷）时出现报错结构不收敛的问题
本帖最后由皮卡丘闪电链于 2025-6-23 19:03 编辑

各位老师好，我在进行Pd/TiO2超大体系（含4%氧缺陷）的结构优化。
目前任务是可以跑起来，inp文件如图：
现在不知道这样对不对，想请老师们看一下有没有问题。
目前收敛值是470多跑了94步，距离10-5还有如海般的距离。这样还需要继续跑吗。。。
如果不对的话我也想尽了cp2k能想到的办法。主要体系是1800个原子的超大体系。
还有显示ELPA报错，好像也无法避开。
还有一个主要问题我认为是比较关键的就是表面制造48个氧缺陷，导致表面电荷是不守恒的，模板中对电荷是没有设置的，也可能导致能量不收敛。
请各位老师指点一下。是否该用LAMMPS等其他软件优化。

作者
Author: sobereva 时间: 2025-6-23 19:34
1800原子DFT优化就别想了，除非计算条件超级好。GFN1-xTB或SCC-DFTB优化倒是可以考虑，至少也先用这等便宜级别做个预优化
SCF不收敛问题看下文
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

作者
Author: 皮卡丘闪电链 时间: 2025-6-25 14:30
十分感谢卢老师的解答，目前按照您说的我已经跑完GFN1-xTB的预优化。我还有一个问题想请教您，我想模拟Pd/TiO2在1500K的温度下进行包裹的全过程，这里会用到LASP和PdTiO的NN势，这个程度的预优化可以用来作为初始结构吗，是否还需要DFT的精细优化？

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