单点能计算入门:水分子示例
下载 Gaussian 输入文件,复算固定构型水分子的电子能。
准备固定构型
使用下方可下载的 water.gjf,在 Gaussian 16 Rev. C.02 中运行。输入坐标单位为 Å,电荷为 0,自旋多重度为 1。采用气相 B3LYP/def2-SVP,并启用严格 SCF 收敛条件。
运行计算
在可使用 Gaussian 的环境执行 g16 < water.gjf > water.log。在 CEMP 单点能页面可查看实际提供的输入模式和模板。如果所选流程在求能量前优化构型,运行的就是另一种计算。
读取结果
该输入的最后一条 SCF Done 记录为 −76.3581208721 Hartree,SCF 迭代 9 次。还应确认日志正常结束。此值是固定构型的电子能,未加入零点能和热校正。
示例输入或输出
g16 < water.gjf > water.log
SCF Done: E(RB3LYP) = -76.3581208721 A.U. after 9 cycles