计算化学公社

标题: 计算Solvent效应报错"No solvent atoms in DisRep" [打印本页]

作者
Author: Jacken_zzq 时间: 2017-9-6 00:58
标题: 计算Solvent效应报错"No solvent atoms in DisRep"
用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.

作者
Author: kyuu 时间: 2017-9-6 01:45
本帖最后由 kyuu 于 2017-9-6 02:49 编辑

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

作者
Author: Jacken_zzq 时间: 2017-9-6 04:01

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

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

作者
Author: kyuu 时间: 2017-9-6 04:08

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

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

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

作者
Author: Jacken_zzq 时间: 2017-9-6 08:29

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

好的谢谢老师

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