高斯软件计算溶液中H+能量出错

xingyuan

各位老师好，我想求如下图中在溶液中O-H之间的键能，即EEACOO--H。按照平时做法，对EEACOOH计算结构优化和能量，EEACOO-计算结构优化和能量，对H+计算结构优化和能量。但是奇怪的是计算H+结构优化时出错，显示利用了UAO模型，但是H上面没有键所以出错。但是我就想计算H+单独的能量。如何计算？
这个是输入文件：
%chk=h.chk
#p opt freq ub3lyp/6-31+g(d,p) scrf=(solvent=dmso) geom=connectivity

Title Card Required

1 1
H                 -0.97701148    0.61576354    0.00000000

1

输出文件后几行：
United Atom Topological Model (UA0  parameters set).
UA0: Hydrogen    1 is unbound. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.
Error termination via Lnk1e in d:\gauss03\l301.exe at Mon Jun 04 17:04:34 2018.
Job cpu time:  0 days  0 hours  0 minutes  1.0 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      1 Scr=      1

麻烦各位老师啦

sobereva

xingyuan 发表于 2018-6-5 09:42
那我在计算另一个结构EEACOO-时，也要算自由基，然后电荷设为0，自旋多重度为2？麻烦老师了

是

sobereva

xingyuan 发表于 2018-6-5 09:19
好的，多谢老师。我在网上看到有人这样算的您看对吗？麻烦了。
%chk=h.chk
# HF/6-311++G** POP=FULL I ...

又是严重的以讹传讹
别信其它论坛的讨论，90%以上都是错的
仔细阅读
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

xingyuan

sobereva 发表于 2018-6-5 06:34
你算键能，应该算的是氢自由基，怎么算的是H+？
G03默认用的是联合原子模型，构造PCM孔洞时候不显式地考虑 ...

那我在计算另一个结构EEACOO-时，也要算自由基，然后电荷设为0，自旋多重度为2？麻烦老师了

njfuzjs

xingyuan 发表于 2018-6-5 09:19
好的，多谢老师。我在网上看到有人这样算的您看对吗？麻烦了。
%chk=h.chk
# HF/6-311++G** POP=FULL I ...

为什么还用HF...

xingyuan

sobereva 发表于 2018-6-5 06:34
你算键能，应该算的是氢自由基，怎么算的是H+？
G03默认用的是联合原子模型，构造PCM孔洞时候不显式地考虑 ...

好的，多谢老师。我在网上看到有人这样算的您看对吗？麻烦了。
%chk=h.chk
# HF/6-311++G** POP=FULL IOP(6/7=3)   

H  

0,2
1   0.0 0.0 0.0

sobereva

你算键能，应该算的是氢自由基，怎么算的是H+？
G03默认用的是联合原子模型，构造PCM孔洞时候不显式地考虑氢，因此也不能直接这么算
而且就一个原子/离子，显然不可能做opt

都是常识性问题

xingyuan