计算化学公社

标题: CP2K优化吸附模型固体底部原子冻结失败 [打印本页]

作者
Author: fffff 时间: 2026-2-5 23:35
标题: CP2K优化吸附模型固体底部原子冻结失败
学生使用cif导入Gaussview后保存为gjf，然后用Multiwfn读取gjf生成cp2k输入文件，选择主功能-1 3 优化结构（晶胞固定）然后在主功能9设置原子冻结，冻结原子的序号通过Gaussivew的atom selection功能选择，冻结最下面两层PbI2，观察优化轨迹的xyz文件发现底部原子仍有轻微波动。输入文件及结果见附件。提前感谢各位老师

作者
Author: sobereva 时间: 2026-2-6 04:22
看看文件搞错没有
如果确信没弄错文件，反复定义多次list，每行list只写不超过10个原子序号再试。并且说明当前用的CP2K版本，别用太老的

作者
Author: fffff 时间: 2026-2-6 11:48

sobereva 发表于 2026-2-6 04:22
看看文件搞错没有
如果确信没弄错文件，反复定义多次list，每行list只写不超过10个原子序号再试。并且说明 ...

感谢卢老师，确保文件没搞错。经过测试使用Multiwfn按照z方向重排列序号之后用LIST 1..128的连续序号去冻结就可以正常冻结。如果是类似于下图的非连续序号，则无法正常冻结（需要冻结原子很多的情况下不能正常冻结，不知道冻结少量原子用非连续序号是否可行）
(, 下载次数 Times of downloads: 0)
PS：仅限于CP2K 2025.1 其他版本未做测试
C:\Users\wyf\Desktop\冻结.png

作者
Author: sobereva 时间: 2026-2-6 17:00
北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里opt\Si_surf_recons\opt.inp就是非连续序号的例子，冻结效果没问题

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