对同一个构型不同多重度的结构优化和单点能计算

小书童

对原始构型在1  3  5三个自旋多重度下进行了结构优化（OPT+Freq），分别对优化好的结构进行了单点能计算。OPT+Freq用的是 B3PW91 6-31G(d)   单点能计算用的是B2PLYP def2tzvp
能量计算的结果如图所示    B2PLYP +  ZPE =Etot   计算结果表明1自旋多重度的时候结构最稳定，但是观察1对应的结构发现结构扭曲，而另外两个多重度对应的结构没有扭曲，很好。请问这是为什么呢？   输入输出文件如附件，大家可以看看有什么不妥的。

liyuanhe211

排除计算错误的原因，也是有可能的（或者说是一定的）。不同多重态之间成键形式都不一样，乙烯的单重态（基态）是双键，稳定构象是平面的；要是三重态就是个双自由基了，我猜最稳定构象可能是两个CH2平面垂直（没计算过），但肯定和乙烯不一样。这是容易理解的。

小书童

liyuanhe211 发表于 2015-10-17 20:44
排除计算错误的原因，也是有可能的（或者说是一定的）。不同多重态之间成键形式都不一样，乙烯的单重态（基 ...

感谢您的解答，我也希望是我自己的计算错误的原因，但是我自己检查过了，没发现错误，如果您愿意我可以吧我计算的输入和输出文件给您看一下。

sobereva

有机分子体系一般都是单重态最稳定，更高多重度会削弱成键
扭曲就扭曲。也可以尝试其它级别进行优化，并且opt=tight，如果都扭曲说明就是扭曲的。

liyuanhe211

小书童 发表于 2015-10-17 20:57
感谢您的解答，我也希望是我自己的计算错误的原因，但是我自己检查过了，没发现错误，如果您愿意我可以吧 ...

呃没说你计算错了。只是排除一下可能性。如果三种多重态构象一样才不正常。

小书童

sobereva 发表于 2015-10-17 23:27
有机分子体系一般都是单重态最稳定，更高多重度会削弱成键
扭曲就扭曲。也可以尝试其它级别进行优化，并且 ...

sob老师，我计算不同多重度下的结构和单点能是为了找到能量最低，最稳定的结构作为基态，为后面的吸附汞做准备，我用老师的方法试试，计算中没有对ZPE修正，等我修正一次试试结果。谢谢老师的解答。