计算化学公社

标题: 有关《外电场对18碳环调控作用》中能量计算的咨询 [打印本页]

作者
Author: 恒星建造指南 时间: 2021-11-11 13:43
标题: 有关《外电场对18碳环调控作用》中能量计算的咨询
各位老师，我在阅读卢老师的博文“一篇文章深入揭示外电场对18碳环的超强调控作用”时，有关能量计算部分有一个地方不太理解。卢老师在在文章的Support information中给出了能量的计算方法。对于其中一个能量的计算方法我不太明白，计算方法见图一。
我粗浅的理解是这样的：假如要研究的分子为A，那么首先将A在外加电场下进行几何优化，得到优化后的结构，将优化后的结构进行外加电场下的单点能计算，取最后一步收敛的波函数作为初猜，再次进行单点能计算，取第一次自洽场迭代输出的能量作为最终结果。
我的疑问如下：
1、对于红字第一部分的操作，单点能计算好像没有类似于优化的收敛过程，因此在“取最后一步收敛的波函数作为初猜”这一点上我不太理解。
2、对于红字第二部分的操作，在哪里能看到第一次自洽场迭代输出的能量？
谢谢各位老师。

作者
Author: sobereva 时间: 2021-11-11 17:17
你的红字写得不对

算你截图里这个量分为三步：
(1)：先在电场下优化结构，得到R_F
(2)：在(1)的结构下，在没有外电场的情况做一次单点计算，最后存到chk里的收敛的波函数就是Ψ0(R_F)
(3)：在(1)的结构下，设置加外电场的关键词（此时对应带电场的哈密顿H_F），通过guess=read读(2)的chk，即用(2)收敛的波函数当做初猜波函数，第一轮的SCF能量便是截图里这个量（计算时候用#P，读第一次输出的E=后面的值）

作者
Author: 恒星建造指南 时间: 2021-11-11 20:02

sobereva 发表于 2021-11-11 17:17
你的红字写得不对

算你截图里这个量分为三步：

谢谢sob老师

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