老师，我想请教一下怎么用高斯算H+的能量呢？

zhengry

老师，我想请教一下怎么用高斯算H+的能量呢？在我的机理中有COO2-+H+→COOH-以及COOH-+H+→CO+H2O，这两步反应涉及到质子的转移。
我的液体环境是介电常数80的水溶液，用到的基组是def2-TZVP(Fe,Co),def2-SVP(H,O,N,C)，泛函是B3LYP-D3。
可能是加了水环境或者其他的原因，导致我算的能量特别大，比如我用这个环境算H2的能量是-1.17a.u，也就是-31.69e.V。我算COO2-和COOH-的差值就是-12e.V，所以如果用其他方式算出的H+能量去算的话就不太合理。但是直接算H+能量会因为没有电子报错，所以想问下可以用其他什么方法算H+的能量吗？或者用什么代替H+的能量呢？
附件是我算H2和H+的输入以及输出文件。

dengyh

sobereva 发表于 2024-10-20 22:06
质子真空下单点能为0
稍有量化理论常识就知道这点

好的，谢谢老师的指导。

sobereva

dengyh 发表于 2024-10-20 21:48
你好，问一下这里的提供的是自由能和焓值，质子的单点能应该是多少呢

质子真空下单点能为0
稍有量化理论常识就知道这点

dengyh

Yaocy 发表于 2023-11-2 10:53
298 K, 1 atm下，质子热力学数据：
G(H+)=-26.3101 kJ/mol
H(H+)=6.1394 kJ/mol

你好，问一下这里的提供的是自由能和焓值，质子的单点能应该是多少呢

zhengry

shenzp 发表于 2023-11-2 10:58
建议阅读谈谈该从Gaussian输出文件中的什么地方读电子能量，搞明白什么是电子能量。对稍大点的体系，电子 ...

好的，非常感谢

zhengry

shenzp 发表于 2023-11-2 10:58
建议阅读谈谈该从Gaussian输出文件中的什么地方读电子能量，搞明白什么是电子能量。对稍大点的体系，电子 ...

好的谢谢老师，我认真学习下

shenzp

可能是加了水环境或者其他的原因，导致我算的能量特别大，比如我用这个环境算H2的能量是-1.17a.u，也就是-31.69e.V。

建议阅读谈谈该从Gaussian输出文件中的什么地方读电子能量，搞明白什么是电子能量。对稍大点的体系，电子能量几万个a.u.都很正常，只有差值才有化学意义。

所以想问下可以用其他什么方法算H+的能量吗

博文谈谈隐式溶剂模型下溶解自由能和体系自由能的计算里面有讲如何获得质子的溶解自由能以及计算级别的选择，建议仔细阅读。

另外，你的标题不严谨。能量有电子能量、自由能等，自由能包括电子能量加热力学校正量，需要说清楚是是什么能量。从你的问题来看，应该是指在溶液中的自由能。真空中质子的电子能量为0，有限温度下由于平动会有热力学能，两者之和即为气相下的自由能，再加上溶解自由能就是溶液中的自由能。

Yaocy

298 K, 1 atm下，质子热力学数据：
G(H+)=-26.3101 kJ/mol
H(H+)=6.1394 kJ/mol

0 K下各种热力学数据为0