zhangby 发表于 2026-2-10 13:34
卢老师好，
1、关于泛函的选择，我是在计算SMD模型下该络合物的溶解自由能，按照您在“谈谈隐式溶剂模型 ...

如果你的最终目的是算分子在溶剂下的自由能，那种算溶解自由能的做法对当前情况（非普通有机分子）未必理想，步骤还更多。直接按下文的过程算溶剂下的自由能就完了，直接带着溶剂模型做opt freq和高精度单点
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

至于你的第2个问题，一个可能是溶剂效应影响了能量本来就相近的不同自旋态之间的相对能量，导致真空和隐式溶剂模型下最稳定波函数对应的自旋多重度不同。也有可能只是stable=opt巧合优化到了不同波函数（注：不代表stable提示稳定就真的是能量最低波函数）。

你做的是波函数稳定性测试+发现不稳定时自动试图优化出稳定的波函数，而不是波函数稳定性测试这一项。当前说明程序发现另外的自旋态能量更低
配体、计算级别等因素都会影响配合物不同自旋态之间相对能量的高低。别把配合物里的Eu(III)简单视为孤立的Eu(III)原子。

M05-2X算镧系配合物是糟糕透顶的选择
认真看下文
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

重新编辑帖子，好好写标题，别用汉语拼音