结构优化出现虚频+势能扫描

whtu · 发表于 Post on 2016-5-13 17:21:50

首先我对图片中1-2-3-4二面角进行势能扫描，输入文件为OcSO4SCAN1.gjf，我从扫描得到的势能面曲线得到我初始的输入结构（即第一个扫描点）的能量最低（ -1014.52386920 a.u.），于是我选择这个结构在BVP86/TZVP水平重新进行结构优化和频率分析，输入文件为OcSO4num1.gjf，结果发现存在虚频，且能量变高（-1014.52379093 a.u.），不知道为什么？按势能扫描的结果，选用初始结构进行结构优化和势能扫面应该会得到一个能量相差不大的稳定构型。另外，我同样选用第一个扫面点，方法和基组换成：b3lyp/6-311++g(d,p)进行结构优化和频率分析，输入文件为OcSO4B3LYP.gjf，结果可以得到稳定的构型，没有虚频，但是能量变大（-1014.42338430 a.u.），这里可以得到稳定的构型，请问这是为什么？
量化知识实在匮乏，希望能得到大家的帮助，谢谢。

alwens · 发表于 Post on 2016-5-13 18:51:09

不同基组，不能比较的

whtu · 发表于 Post on 2016-5-13 19:45:54

alwens 发表于 2016-5-13 18:51
不同基组，不能比较的

嗯，但是为什么根据势能扫描的结果选的结构进行优化和频率分析却得不到稳定的构型，而且得到的能量还偏高。这里的势能扫描的基组与优化的基组是一样的。

liyuanhe211 · 发表于 Post on 2016-5-13 21:18:44

前两个能量仅差0.2 kJ/mol，分析这个差异在当前方法精度下毫无意义，完全可以认为是一样的。
将不同水平的b3lyp/6-311++g(d,p)与BVP86/TZVP得到的能量相比较也毫无意义。
你的结构在一般的温度下都会是一大堆低能构象的玻尔兹曼分布，用任何一个单一构想表述都可能会造成对性质分析的显著误差。

whtu · 发表于 Post on 2016-5-15 09:26:19

liyuanhe211 发表于 2016-5-13 21:18
前两个能量仅差0.2 kJ/mol，分析这个差异在当前方法精度下毫无意义，完全可以认为是一样的。
将不同水平的 ...

其实我主要的目的不是比较能量的差异，我不明白的是根据势能扫描的结果，我选用第一个扫描点进行结构优化理论上应该能得到一个稳定的构型，可是实际上计算却出现了虚频？

sobereva · 发表于 Post on 2016-5-15 10:19:59

whtu 发表于 2016-5-15 09:26
其实我主要的目的不是比较能量的差异，我不明白的是根据势能扫描的结果，我选用第一个扫描点进行结构优化 ...

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

另外，优化这个体系，用BVP86不会有任何好处。明显不如用B3LYP。

whtu · 发表于 Post on 2016-5-16 08:30:41

sobereva 发表于 2016-5-15 10:19
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

谢谢！

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

[新手求助] 结构优化出现虚频+势能扫描

评分 Rate

浏览过的版块