计算化学公社

标题: 高斯软件计算溶液中H+能量出错 [打印本页]

作者
Author: xingyuan 时间: 2018-6-4 20:43
标题: 高斯软件计算溶液中H+能量出错
各位老师好，我想求如下图中在溶液中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

麻烦各位老师啦

作者
Author: xingyuan 时间: 2018-6-4 20:44
[attach]14183[/attach]

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

都是常识性问题

作者
Author: xingyuan 时间: 2018-6-5 09:19

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

作者
Author: njfuzjs 时间: 2018-6-5 09:37

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

为什么还用HF...

作者
Author: xingyuan 时间: 2018-6-5 09:42

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

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

作者
Author: sobereva 时间: 2018-6-5 10:02

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

作者
Author: sobereva 时间: 2018-6-5 10:03

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

是

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