计算化学公社

标题: 吸附能计算中sccs溶剂参数设置问题 [打印本页]

作者
Author: madhatter 时间: 2022-3-12 18:53
标题: 吸附能计算中sccs溶剂参数设置问题
本帖最后由 madhatter 于 2022-3-12 21:04 编辑

我做离子液体中某种分子在锌表面的吸附能计算，用的cp2k，为了找到标准参数，重复文献中锌表面水的吸附能，我在真空下对锌002晶面和水分子优化。然后计算能量

文献给的是-0.299eV（top），我真空下得到的是0.1eV（正值），设置sccs，用Multiwfn初始设置，介电常数是78.36，其他都是默认，pbe算的能量，吸附能的结果是-0.06eV，我想问下这个差距大吗？，其次关于sccs中的设置都是什么，例如ALPHA BETA GAMMA这些，在算我的溶剂体系时，我直接用量化计算中的溶剂介电常数和表面张力参数就可以吗？其他的就按照默认的？
&SCCS

   ALPHA [N*m^-1] 0.0
   BETA [kbar] 0.0
   GAMMA [mN/m] 0.0
   DIELECTRIC_CONSTANT 78 #Water
   EPS_SCCS 1E-6 #Default. Requested accuracy for the SCCS iteration cycle
   EPS_SCF 0.5 #Default. SCCS iteration is activated only if SCF iteration is converged to this threshold
   MAX_ITER 100 #Default. Maximum number of SCCS iteration steps
   DERIVATIVE_METHOD FFT #Default. Method for calculation of numerical derivatives. Can also be CD3, CD5, CD7
   &ANDREUSSI
      RHO_MAX 0.0035 #Default
      RHO_MIN 0.0001 #Default
   &END ANDREUSSI
&END SCCS

---------------------------------------------------------------
补充刚刚有算了一下  对分子不用sccs，其他加sccs，结果时0.36eV，所以吸附能计算过程中是否需要都加隐式溶剂

我对分子是直接用算总能量的文件把锌的坐标删除，直接算的，我用Multiwfn生成文件算分子不加sccs，结果是0.36eV,所以有人知道具体的计算过程吗？这两种不同的方法能量不一样吗？

作者
Author: sobereva 时间: 2022-3-13 07:39
真空下算出来吸附能为正这明显不对。先把这个算对再说SCCS的事。

ALPHA BETA GAMMA这些参数去看SCCS原文便知含义，影响的是溶剂模型的非极性部分。原文并没有对不同溶剂直接给出相应的值，CP2K里也没预置各种溶剂的参数，只能根据其物理意义自己根据溶剂属性去设。它们对结果的影响都远不如决定溶剂模型极性部分表现的EPS大。

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