计算化学公社

标题: 对同一个构型不同多重度的结构优化和单点能计算 [打印本页]

作者
Author: 小书童 时间: 2015-10-17 19:47
标题: 对同一个构型不同多重度的结构优化和单点能计算
本帖最后由小书童于 2015-10-17 22:11 编辑

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

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

作者
Author: 小书童 时间: 2015-10-17 20:57

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

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

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

作者
Author: liyuanhe211 时间: 2015-10-17 23:28

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

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

作者
Author: 小书童 时间: 2015-10-18 08:56

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

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

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