计算化学公社

标题: 求助：隐式溶剂模型计算单点极难收敛的问题 [打印本页]

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

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

作者
Author: Penson 时间: 2021-2-4 00:01

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

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

作者
Author: wzkchem5 时间: 2021-2-4 11:32

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

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

作者
Author: Penson 时间: 2021-2-4 23:42

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

了解，十分感谢

作者
Author: 时光的邂逅 时间: 2022-2-28 14:14

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

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

作者
Author: Penson 时间: 2022-3-1 10:20

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

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

作者
Author: biogon 时间: 2023-4-7 17:55

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

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

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