计算化学公社

标题: 势能曲线扫描发生电子态转变、定位不到TS的问题 [打印本页]

作者
Author: haoguoyuerban 时间: 2026-6-10 12:29
标题: 势能曲线扫描发生电子态转变、定位不到TS的问题
[attach]131933[/attach]大家好，求助一个问题: 使用Gaussian DFT对一个多核金属配合物体系（含有酚自由基）的扫描质子转移坐标的时候，发现扫描过程中有能量断层下降的现象，查看输出文件发现此处发生了金属离子上自旋分布的改变，意味着金属离子被氧化，酚自由基被还原，也就是质子耦合电子转移（PCET），能量断层应该表示电子态（或者势能面）发生了跃迁？但是始终没发现质子转移过程中的过渡态（TS），也就没办法计算这个过程的能垒。请问各位专家，这种情况该如何解决？怎么能定位到这个过程的TS或者计算这个过程的能垒？谢谢！

作者
Author: sobereva 时间: 2026-6-10 12:36
如果不是柔性扫描时几何结构发生了突越，说明电子结构突然变化，大概率是突越前能量高的那个点没收敛到相应结构下的稳定波函数。扫描之前就应当做波函数稳定性测试，必须以稳定的波函数当初猜进行扫描

作者
Author: haoguoyuerban 时间: 2026-6-15 17:16

sobereva 发表于 2026-6-10 12:36
如果不是柔性扫描时几何结构发生了突越，说明电子结构突然变化，大概率是突越前能量高的那个点没收敛到相应 ...

感谢老师回复！！我这个是柔性扫描的,扫描之前也做过波函数稳定性测试(stable=opt)和几何优化的。扫描结果说明金属-配合物体系的电子结构的突然变化和质子转移是强烈耦合在一起的，但是中间没有发现TS，这样合理吗？我发现有一篇文献“Tutorial on computing nonadiabatic proton-coupled electron transfer rate constants”J. Chem. Phys. 163, 091501 (2025）可能设计这个问题的解决办法（但不确定），里面建议使用Constrained DFT 分别固定两种电子态扫描质子转移，但是Gaussian好像没有这个功能，不知老师您有没有好的建议？

作者
Author: sobereva 时间: 2026-6-15 17:25

haoguoyuerban 发表于 2026-6-15 17:16
感谢老师回复！！我这个是柔性扫描的,扫描之前也做过波函数稳定性测试(stable=opt)和几何优化的。扫描结 ...

NWChem、CP2K都支持约束性DFT
谈谈约束性DFT(CDFT)
http://sobereva.com/271（http://bbs.keinsci.com/thread-535-1-1.html）

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