计算化学公社

标题: Gaussian计算初始构型得电子到末态，算出的能量与预期不符，而且过渡态没有负的虚频 [打印本页]

作者
Author: Shuhan_Lu 时间: 2022-2-28 19:42
标题: Gaussian计算初始构型得电子到末态，算出的能量与预期不符，而且过渡态没有负的虚频
本帖最后由 Shuhan_Lu 于 2022-2-28 20:28 编辑

想重复文献（ http://dx.doi.org/10.1002/cssc.201601356)里氟代碳酸乙烯酯（FEC)的单电子/双电子还原过程。比如构型47得电子到构型49，文献里49的能量比47高。而我算的47能量=-441.767071Hartree,49能量=-441.840309，明显是49能量远低于47。其它有些构型也有类似情况（其他构型的能量见图片，能量都取自输出文件里的Sum of electronic and thermal Free Energies)。47的构型是从网上找的FEC优化后得到的(#p opt freq b3lyp/6-311+g(d,p) scrf=(solvent=generic,read,pcm) empiricaldispersion=gd3bj)。不知道是47的构型没优化好还是加电子的方法不对(我是通过在gjf文件里改电荷的数值实现的)。另外几个过渡态虚频算出来没有负值，跟文献也不一样。47、49、TS11的输入、输出文件在附件里。Gaussian初学者，面临毕业，谁能帮忙解答下，万分感激！

作者
Author: sobereva 时间: 2022-3-2 03:22
要和文献比就用和文献精确相同的计算级别和设置
如果精确相同了，而且也确信没其它不符的地方（几何结构、程序版本等），直接发邮件问作者

作者
Author: Shuhan_Lu 时间: 2022-3-17 18:37

sobereva 发表于 2022-3-2 03:22
要和文献比就用和文献精确相同的计算级别和设置
如果精确相同了，而且也确信没其它不符的地方（几何结构、 ...

谢谢老师回复，现在在用相同精度试试，

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