关于ONIOM计算之后，QM部分HOMO-LUMO能级的咨询，谢谢指点

小范范1989

背景介绍，想通过ONIOM计算一下中间QM部分的能级，同和是溶剂中单分子的能级进行比较，阐述显示溶剂模型周围分子对中间QM分子的影响。问题如下：
1：通过ONIOM优化之后，读取输出文件
Alpha  occ. eigenvalues --   -0.30438  -0.30259  -0.29995  -0.29746  -0.28833
Alpha  occ. eigenvalues --   -0.27975 -0.26564
Alpha virt. eigenvalues --   -0.11030  -0.09168  -0.07977  -0.07709  -0.06997

所以，这样的话，能够看到HOMO和LUMO的能级。
但是，如果通过优化的S0的结构，做td。会出现如下的输出：
Alpha  occ. eigenvalues --   -0.30485  -0.30376  -0.30102  -0.29856  -0.28830
Alpha  occ. eigenvalues --   -0.28249 -0.26851
Alpha virt. eigenvalues --   -0.11570 -0.09693  -0.08343  -0.08141  -0.07247

问题出现了：通过ONIOM优化得到的结构能级，和基于优化好的结构做td，输出文件的能级不一样。按说都是基于同一个结构，能级不应该是一样的吗？
为了作对比，我计算了溶液中的能级。（PCM模型）

2：溶剂中优化S0得到的能级：
Alpha  occ. eigenvalues --   -0.27727  -0.27662  -0.27554  -0.26639  -0.25854
Alpha  occ. eigenvalues --   -0.25413  -0.23493
Alpha virt. eigenvalues --   -0.07713  -0.05693  -0.04460  -0.04195  -0.03993

给予优化的结构做td，得到的能级：
Alpha  occ. eigenvalues --   -0.25413  -0.23493
Alpha virt. eigenvalues --   -0.07713  -0.05693  -0.04459  -0.04195  -0.03993

可以看到都是一样的。这个应该是没问题吧。

所以，问题如下：
1：我可以通过这种方式来找ONIOM中QM的能级吗？
2：通过ONIOM计算之后，应该读取哪个能级呢？是优化S0的文件，还是做td的文件？还是什么文件？
3：之前咨询sob老师说，这个ONIOM计算完之后，不能通过fchk文件，单独看QM的轨道分布。那如果通过上述方法知道了能级，如何看能级分布呢？背景电荷的方式还是什么方式？
4：谢谢老师指点。

sobereva

描述不是很清楚。
“优化的S0的结构，做td”，结构是ONIOM优化完的整体的结构当中只取出中间分子的结构？还是依然用ONIOM做，只不过high方法又加了个TD？
就直接取ONIOM优化任务最后一步输出当中的QM算的high层计算结果对应的轨道能量就行了。

能级分布不知具体指什么。

小范范1989

sobereva 发表于 2017-10-27 14:09
描述不是很清楚。
“优化的S0的结构，做td”，结构是ONIOM优化完的整体的结构当中只取出中间分子的结构？ ...

sob老师好，我没有描述好。
ONIOM优化完成之后，基于这个，依然用ONIOM做，high加了td关键词。
老师说的这个“就直接取ONIOM优化任务最后一步输出当中的QM算的high层计算结果对应的轨道能量就行了。”是把最后ONIOM优化好的结构，单独把QM取出来，然后做一个单点，来看轨道能级，还有HOMO，LUMO分布吗？
我表述的能积分布，就是看HOMO和LUMO的分布，
谢谢sob老师指点

sobereva

小范范1989 发表于 2017-10-27 16:01
sob老师好，我没有描述好。
ONIOM优化完成之后，基于这个，依然用ONIOM做，high加了td关键词。
老师说 ...

不是重新算一个，就是直接从ONIOM优化任务最后读取QM的轨道能级
你重新算个单点也行，但此时必须也是做ONIOM单点，然后从中读取QM的轨道能。

小范范1989

sobereva 发表于 2017-10-28 01:12
不是重新算一个，就是直接从ONIOM优化任务最后读取QM的轨道能级
你重新算个单点也行，但此时必须也是做O ...

偶，好的，谢谢sob老师指点