Gaussian计算初始构型得电子到末态，算出的能量与预期不符，而且过渡态没有负的虚频

Shuhan_Lu · 发表于 Post on 2022-2-28 19:42:45

本帖最后由 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初学者，面临毕业，谁能帮忙解答下，万分感激！

sobereva · 发表于 Post on 2022-3-2 03:22:58

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

Shuhan_Lu · 发表于 Post on 2022-3-17 18:37:20

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

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

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

[Gaussian/gview] Gaussian计算初始构型得电子到末态，算出的能量与预期不符，而且过渡态没有负的虚频

浏览过的版块