求助：关于使用CP2K配合Multiwfn和Sobtop生成非标准氨基酸力场参数的疑问

Leafia

各位老师们好, 最近我正在尝试为GROMACS生成非标准氨基酸残基的拓扑文件，参考了Sobtop官方FAQ5 (http://sobereva.com/soft/Sobtop/#download)的流程，我打算使用CP2K完成前期的量子化学计算。但在实操中遇到了一些细节疑惑，特来请教：

1. cp2k是否可以完全代替Gaussian和ORCA进行几何优化和振动分析呢，因为我看到http://sobereva.com/651这里Sobereva老师有提到可以用cp2k的 .molden文件进行处理后输入进Multiwfn?

2. 我已经使用CP2K进行了几何优化和振动分析，但在振动分析中开启.Molden 输出后，得到了很多文件。请问用于拟合RESP电荷的波函数，我该选择用哪个文件呢？我是否应该在几何优化完成后，单独拿优化好的最优坐标，跑一次 RUN_TYPE=ENERGY 的单点计算，来提取唯一正确的.Molden 输入给Multiwfn？

感谢各位前辈百忙之中答疑！

牧生

http://sobereva.com/soft/Sobtop/#download  中说的很清楚，用量子化学程序来计算。

常用的CP2K是第一性原理程序，Gaussian和ORCA是量子化学程序。