求助：隐式溶剂模型计算单点极难收敛的问题

Penson · 发表于 Post on 2021-2-3 22:38:56

使用cp2k7.1，对一个6*6的4层金属锌切面进行真空结构优化后，尝试加隐式溶剂水（DFT/SCCS）进行单点能计算，但是发现调整了两种溶剂化的方法（默认的ANDREUSSI和FATTEBERT-GYGI），都没法收敛，请问还有什么策略可以采用进一步帮助收敛呢？附上我的输入文件和优化后的切面坐标文件。如有大神相助，感激不尽！

wzkchem5 · 发表于 Post on 2021-2-3 23:02:59

第一层用显式水吧
溶剂对金属的极化作用太强，而且也不真实，最外面一层Zn原子肯定是和水配位的，需要显式描述水分子
另外，除非你是在研究刚从真空里放进水的金属Zn，否则表面肯定不是clean surface。刚从空气里放进水的Zn，表面是ZnO；已经在水里停留了一小段时间的Zn，表面是Zn(OH)2

Penson · 发表于 Post on 2021-2-4 00:01:28

wzkchem5 发表于 2021-2-3 23:02
第一层用显式水吧
溶剂对金属的极化作用太强，而且也不真实，最外面一层Zn原子肯定是和水配位的，需要显式 ...

第一层有36个锌原子，应该和显式配位水的数量比例成多少合适呢？感觉会变成好大的体系

wzkchem5 · 发表于 Post on 2021-2-4 11:32:50

Penson 发表于 2021-2-4 00:01
第一层有36个锌原子，应该和显式配位水的数量比例成多少合适呢？感觉会变成好大的体系

比较保险的方法是多加一些水分子，用半经验跑个分子动力学退火，再只保留那些和Zn表面离得最近的水。

Penson · 发表于 Post on 2021-2-4 23:42:15

wzkchem5 发表于 2021-2-4 11:32
比较保险的方法是多加一些水分子，用半经验跑个分子动力学退火，再只保留那些和Zn表面离得最近的水。

了解，十分感谢

时光的邂逅 · 发表于 Post on 2022-2-28 14:14:18

Penson 发表于 2021-2-4 23:42
了解，十分感谢

我也遇到了加溶剂后算单点不收敛的问题，请问您后面怎么处理了呢？

Penson · 发表于 Post on 2022-3-1 10:20:54

时光的邂逅发表于 2022-2-28 14:14
我也遇到了加溶剂后算单点不收敛的问题，请问您后面怎么处理了呢？

暂时也没有很好的解决办法

biogon · 发表于 Post on 2023-4-7 17:55:30

时光的邂逅发表于 2022-2-28 14:14
我也遇到了加溶剂后算单点不收敛的问题，请问您后面怎么处理了呢？

先用1e-6以上收敛限算一次无溶剂的，再加溶剂读原来波函数算，基本就一步收敛

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：隐式溶剂模型计算单点极难收敛的问题