计算Solvent效应报错"No solvent atoms in DisRep"

Jacken_zzq

用CPCM模型，加入dispersion, repulsion correction时出错。
请教该如何解决，是缺了参数吗？已经设定了数据库中有的solvent为什么还提示"No solvent atom in DisRep"。
谢谢各位老师~

输入文件如下：
# opt() freq m062x/6-311G(2df,p) SCRF(CPCM, solvent=AceticAcid,read)

Title Card Required

0 1
O                  0.00000000    0.00000000    0.11813800
H                  0.00000000    0.75681000   -0.47255200
H                  0.00000000   -0.75681000   -0.47255200

dis
rep
cav

输出文件最后一段：
------------------------------------------------------------------------------
Polarizable Continuum Model (PCM)
=================================
Model                : C-PCM.
Atomic radii         : UFF (Universal Force Field).
Polarization charges : Total charges.
Charge compensation  : None.
Solution method      : On-the-fly selection.
Cavity type          : Scaled VdW (van der Waals Surface) (Alpha=1.100).
Cavity algorithm     : GePol (No added spheres)
                        Default sphere list used, NSphG=    3.
                        Lebedev-Laikov grids with approx.  5.0 points / Ang**2.
                        Smoothing algorithm: Karplus/York (Gamma=1.0000).
                        Polarization charges: spherical gaussians, with
                                              point-specific exponents (IZeta= 3).
                        Self-potential: point-specific (ISelfS= 7).
                        Self-field    : sphere-specific E.n sum rule (ISelfD= 2).
1st derivatives      : Analytical E(r).r(x)/FMM algorithm (CHGder, D1EAlg=3).
                        Cavity 1st derivative terms included.
Solvent              : AceticAcid, Eps=   6.252800 Eps(inf)=   1.882384
------------------------------------------------------------------------------
Spheres list:
ISph  on   Nord     Re0    Alpha      Xe            Ye            Ze
    1   O      1    1.7500  1.100      0.000000      0.000000      0.118138
    2   H      2    1.4430  1.100      0.000000      0.756810     -0.472552
    3   H      3    1.4430  1.100      0.000000     -0.756810     -0.472552
------------------------------------------------------------------------------
Warning! Inconsistent VMol and RSolv for this solvent, using VMol=    0.00 Ang**3.
No solvent atoms in DisRep.
Error termination via Lnk1e in /opt/shared/Gaussian/g09.d01/g09/l301.exe at Sat Sep  2 18:35:29 2017.

kyuu

把后边的
Dis
Rep
Cav
去掉才让算
非要算非极性作用啊

Jacken_zzq

kyuu 发表于 2017-9-6 01:45
把后边的
Dis
Rep

算过不加参数的了，老板想要看看加了非极性作用会不会有影响

kyuu

Jacken_zzq 发表于 2017-9-6 04:01
算过不加参数的了，老板想要看看加了非极性作用会不会有影响

单一默认溶剂，你试试美年达+SMD呗

sobereva

用SMD。此时非极性部分会单独输出，总能量减去这个就是只考虑极性的部分
高斯用户没事绝对不要用CPCM，没任何好处，而且高斯里CPCM对中性体系用的校正参数还是错的。

Jacken_zzq

sobereva 发表于 2017-9-6 05:10
用SMD。此时非极性部分会单独输出，总能量减去这个就是只考虑极性的部分
高斯用户没事绝对不要用CPCM，没 ...

好的 谢谢老师