sobereva 发表于 2022-10-2 17:18
如果你是指m062x/6-31G*结合SMD算溶解自由能，可以。如果你是指拿它算电子能量部分，太low

打扰了sob老师了因为我用的是orca 所以不会写noneq=xxx 所以直接自定义的smd orca中优化可以得到溶剂cpcm的文件，但是如何把这部分信息像高斯一样放到单点计算那里去我就不会了，所以我就在所有过程里都加了隐式溶剂（smd自定义的溶剂）

我目前先算的外重组能(a) 优化中性状态结构 （cpcm）
(b) 在(a)的结构下算中性状态的能量E1 （smd）
(c) 在(a)的结构下算阳离子状态的能量E2 （smd）
(d) 优化阳离子状态结构 （cpcm）
(e) 在(d)的结构下算中性状态的能量E3 （smd）
(f) 在(d)的结构下算阳离子状态的能量E4 （smd）
重组能为
λ(I)=|E3-E1|
λ(II)=|E2-E4|
λh= λ(I)+λ(II)
我优化结构都是用b3lyp def2-TZVP+隐式溶剂cpcm（具有离子特性），算能量的时候是M062X 6-31G*结合smd得到了单点能，直接相加减算λh
因为加了smd所以我用的这个基组，我不太清楚这里是电子能量还是溶解自由能去相加减求λh，目的是用λh去算还ΔG = 0.25 (ΔrG + λ)2/λ，所以这里是用电子能量应该用更高精度的泛函加基组？有推荐吗？

如果你是指m062x/6-31G*结合SMD算溶解自由能，可以。如果你是指拿它算电子能量部分，太low