cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正

Janson_Shi · 发表于 Post on 2023-3-14 12:53:05

cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正。

相关的inp, out, cif 文件附上，麻烦看下参数哪里出了问题？谢谢

Janson_Shi · 发表于 Post on 2023-3-14 16:49:24

sobereva 发表于 2023-3-14 16:35
Multiwfn载入你的cif文件，然后输入
cp2k
t.inp

谢谢群主大佬解答,主要这里涉及到好几个模型,尤其是COF+H2O的模型,所以之前都是手动调参数使COF+H2O SCCS 收敛。结果返回来用相同的参数计算2D COF 就出现上述问题。我先统一下参数，如果对其他模型有相关问题再请教。

sobereva · 发表于 Post on 2023-3-14 16:35:00

Multiwfn载入你的cif文件，然后输入
cp2k
t.inp
7 //开启SCCS
0
得到的输入文件对应PBE/DZVP-MOLOPT-SR-GTH在SCCS描述的水环境下计算单点
产生t.inp和计算完的结果在附件里，结果完全合理。

t.inp (8.44 KB, 下载次数 Times of downloads: 27)

t.out (40.21 KB, 下载次数 Times of downloads: 13)

SCF最后一轮：
  SCCS| Iteration cycle converged in 27 steps
  SCCS|
  SCCS| Polarisation charge                                     -0.005985169559
  SCCS| Hartree energy of solute and solvent [a.u.]          392.757709746063
  SCCS| Hartree energy of the solute only [a.u.]             392.781239084831
  SCCS| Polarisation energy [a.u.]                         -0.023529338768
  SCCS|                      [kcal/mol]                               -14.765
  SCCS| Dispersion free energy [a.u.]                         -0.131855531402
  SCCS|                      [kcal/mol]                               -82.741
  SCCS| Repulsion free energy  [a.u.]                            0.107670747519
  SCCS|                      [kcal/mol]                               67.564
21 Broy./Diag. 0.40E+00 10.1    0.00000346    -410.0624906437  1.01E-04

始终强烈建议用Multiwfn产生CP2K输入文件，超级容易还少走无数弯路！

Janson_Shi · 发表于 Post on 2023-3-14 15:18:24

fineren 发表于 2023-3-14 15:04
RUN_TYPE ENERGY
做过构型优化了吗？

你的意思用SCCS模型来跑结构优化吗？

fineren · 发表于 Post on 2023-3-14 15:04:19

RUN_TYPE ENERGY
做过构型优化了吗？

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

[CP2K] cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正

浏览过的版块