如果计算THF环上的各个键的键能？

Dreamer1 · 发表于 Post on 前天 21:55

老师好，我正在计算THF溶剂分子上面的每个键的键能。我现在已经对整个分子结构进行了 #P B3LYP/def2TZVP em=GD3BJ opt freq 几何优化和振动分析并且能得到整体的电子能量。但是环上的键比如C-C键断开后，还是一个完整连续的分子，并没有分成两个片段，那么这该怎么去几何优化和计算能量呢？电荷和自旋多重度该如何设置呢？

sobereva · 发表于 Post on yesterday 04:14

扭转合适的二面角，使得在相应构型下断键的两个位点之间充分远离，然后视为双自由基体系优化和算能量。注意需要当做自旋极化单重态计算：
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

但小环体系，键能本身就没有严格定义，和环张力有显著耦合
对于仅仅需要算C-C键能的情况，最省事也很合理的做法是根据Multiwfn计算的拉普拉斯键级估计，看下文拉普拉斯键级部分：
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471（http://bbs.keinsci.com/thread-12488-1-1.html）

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 如果计算THF环上的各个键的键能？