求助cp2k 结构优化后频率分析有很多虚频

elpa · 发表于 Post on 2022-5-18 00:11:41

计算流程：初始结构从文献（[size=7.223pt]DOI: 10.1039/d1ma00330e）获取并进行相应修改（替换一部分基团），使用含384个原子的原胞做CELL_OPT，再做GEO_OPT，然后VIBRATIONAL_ANALYSIS，发现有100个虚频。

CELL_OPT及GEO_OPT精度：
   MAX_FORCE  8.0000000E-005
   MAX_DR    8.00000000E-004
   RMS_DR    5.0000000E-004
   RMS_FORCE  5.00000000E-005

其他参数详见输入文件，都在附件。

现在比较困惑的是，虚频产生的原因是参数设置精度不够，还是泛函赝势等问题，还是体系本身的问题。

cellopt.inp (2.45 KB, 下载次数 Times of downloads: 13)

celloptinput.xyz (24.8 KB, 下载次数 Times of downloads: 3)

cp2kfreq.inp (2.24 KB, 下载次数 Times of downloads: 44)

freq.dat (9.34 KB, 下载次数 Times of downloads: 26)

freqinput.xyz (24.8 KB, 下载次数 Times of downloads: 8)

geoopt.inp (2.15 KB, 下载次数 Times of downloads: 9)

geooptinput.xyz (24.42 KB, 下载次数 Times of downloads: 2)

elpa · 发表于 Post on 2022-5-19 15:16:53

开发者在groups里写了一段话

我看懂了但又没完全懂。

啊这就离谱 · 发表于 Post on 2023-1-2 16:44:25

遇到了相似的问题，请问楼主最后解决了么

elpa · 发表于 Post on 2023-1-6 10:11:19

啊这就离谱发表于 2023-1-2 16:44
遇到了相似的问题，请问楼主最后解决了么

没解决，不管怎么提高精度最后还是有虚频。不过如果是算很少几个原子而不是所有原子的频率，提高精度有希望解决。我最后换vasp，发现会好很多。唉离谱

啊这就离谱 · 发表于 Post on 2023-1-9 00:35:22

elpa 发表于 2023-1-6 10:11
没解决，不管怎么提高精度最后还是有虚频。不过如果是算很少几个原子而不是所有原子的频率，提高精度有希 ...

我也是还没解决，试过了Google group里面所有回复的方法都解决不了这个问题，现在换成簇模型了

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

[CP2K] 求助cp2k 结构优化后频率分析有很多虚频

浏览过的版块