使用cLR-SMD溶剂模型进行激发态优化，请问该从什么地方读取能量？

81521

TDDFT进行激发态构象优化时，应用了cLR-SMD溶剂模型，但不知道该从什么地方读取能量


如果从“Total Energy, E(TD-HF/TD-DFT) =" 中读取，似乎没有考虑校正的作用；
从"Total energy after correction =" 中读取, 最后优化的结构能量又不是最小的（最后一次-967.992574884，之前的有-967.992579734 a.u.)

谢谢大家。

81521

突然想到，有没有可能目前Gaussian 16并不支持校正线性响应（cLR）溶剂模型下的优化，文献好像都是用cLR来算单点能的；

Gaussian只是在线性响应LR的基础上优化了结构，然后再把优化时每步的结构进行校正，得到校正能量

wzkchem5

81521 发表于 2024-4-18 12:21
突然想到，有没有可能目前Gaussian 16并不支持校正线性响应（cLR）溶剂模型下的优化，文献好像都是用cLR来 ...

应该是这样的。支持cLR的程序一般都不支持它的解析梯度（因为难写且计算量大），如果程序没有算数值梯度，那应该就是用的校正之前的能量的梯度来做的计算

81521

wzkchem5 发表于 2024-4-18 19:35
应该是这样的。支持cLR的程序一般都不支持它的解析梯度（因为难写且计算量大），如果程序没有算数值梯度 ...

感谢！