计算化学公社

标题: cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正 [打印本页]

作者
Author: Janson_Shi 时间: 2023-3-14 12:53
标题: cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正
cp2k用SCCS模型计算2D COF 收敛，但是 Polarisation energy 为正。[attach]64074[/attach] (, 下载次数 Times of downloads: 11)

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

作者
Author: fineren 时间: 2023-3-14 15:04
RUN_TYPE ENERGY
做过构型优化了吗？

作者
Author: Janson_Shi 时间: 2023-3-14 15:18

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

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

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

(, 下载次数 Times of downloads: 27)
(, 下载次数 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输入文件，超级容易还少走无数弯路！

作者
Author: Janson_Shi 时间: 2023-3-14 16:49

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

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

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