sai77 发表于 2025-6-4 07:39
Tools-Point Group里调整

好的，谢谢老师的帮助

TUT 发表于 2025-6-3 20:38
其实我尝试过用mo62x/def2TZVP算单点能，只是结果还是背离预期的a的能量比b更低，所以就没有继续用了。请 ...

Tools-Point Group里调整

sai77 发表于 2025-5-31 20:51
为什么单点能的基组跟优化保持一样的？至少应该提升一到两个档次吧
另外结构a应该是C2v对称的，你给的结果 ...

其实我尝试过用mo62x/def2TZVP算单点能，只是结果还是背离预期的a的能量比b更低，所以就没有继续用了。请问老师调整a的对称性应该怎么操作呢？

sai77 发表于 2025-5-31 20:51
为什么单点能的基组跟优化保持一样的？至少应该提升一到两个档次吧
另外结构a应该是C2v对称的，你给的结果 ...

单点能的级别是要高点，我同意。
对称性方面，Gaussian默认的对称性的阈值比较严，差一丁点都不行，我的经验是Gaussian经常判低了对称性，除非更改Tolerance。

为什么单点能的基组跟优化保持一样的？至少应该提升一到两个档次吧
另外结构a应该是C2v对称的，你给的结果只优化到了C2，还可以进一步调整对称性继续优化

chands 发表于 2025-5-31 12:15
要计算自由能差值，因为自由能差值（等温等压下是Gibbs自由能）与平衡系数直接相关。
我按照你gjf的级别 ...

好的，那我再算一下试试，感谢老师的帮助

TUT 发表于 2025-5-31 11:53
原来还要计算自由能，谢谢老师指导！附件里的gif是真空优化出的结构，请问老师是在m062x-D3/6-31G+(d,p) ...

要计算自由能差值，因为自由能差值（等温等压下是Gibbs自由能）与平衡系数直接相关。
我按照你gjf的级别算的，没动过。
除非溶质分子与溶剂分子有相当强的作用（比如形成氢键），真空优化的结构和加溶剂模型优化的结构应该差不多。你用真空优化也没啥问题。考虑到你的分子有对称性，优化和振动分析时加上对称性会节省时间。

chands 发表于 2025-5-30 20:44
我用你附件里的结构、级别和溶剂模型直接计算单点能
A=-1196.620869 Hartree
B=-1196.621324 Hartree

原来还要计算自由能，谢谢老师指导！附件里的gif是真空优化出的结构，请问老师是在m062x-D3/6-31G+(d,p)下计算的能量吗？不经过溶剂模型做结构优化直接用溶剂模型计算单点能的结果也是可用的吗？

我用你附件里的结构、级别和溶剂模型直接计算单点能
A=-1196.620869 Hartree
B=-1196.621324 Hartree
B低了11 kJ/mol 左右吧。不过严格来说要比较自由能。