小木虫

当前位置：首页 > 量子化学 >【求助】opt+freq之后能量为什么不是最低？

【求助】opt+freq之后能量为什么不是最低？

作者 hbnuwangshuling: 来源: 小木虫 600 12 举报帖子

我在计算并三苯分子时使用 opt freq关键字没有出现虚频，得到能量是E(RB+HF-LYP) = -539.546299015，我对优化后的结构稍作了一下调整，竟然能量（E(RB+HF-LYP) = -539.546359942）比我优化得到的能量还低，这是为什么？优化之后的结构和我做调整之后的结构放到了附件中！

下载提醒：APP中打开可直接下载，点击下载

structure.doc

今日热帖

精华评论

hbnuwangshuling

我重新上传了一下结构：
% MEM=2GB
% NPROCLIND=1
% NPROCSHARED=4

# B3LYP/6-31G** nosymm

0opt

0 1
C             3.660044 -0.713030 0.000015
C             2.479245 -1.406933 -0.000045
C             1.223658 -0.722595 -0.000028
C             1.223660 0.722596 -0.000019
C             2.479248 1.406933 0.000024
C             3.660046 0.713026 0.000058
C             0.000000 -1.403623 -0.000022
C             0.000001 1.403625 -0.000026
C          -1.223658 0.722598 -0.000022
C          -1.223659 -0.722595 -0.000005
C          -2.479246 -1.406933 0.000025
H          -2.476926 -2.493823 0.000064
C          -3.660045 -0.713029 0.000052
C          -3.660045 0.713027 0.000005
C          -2.479247 1.406933 -0.000028
H          -0.000006 -2.491343 -0.000043
H             4.606168 -1.246226 0.000041
H             2.476923 -2.493823 -0.000093
H             2.476928 2.493822 0.000025
H             4.606165 1.246230 0.000113
H          -0.000005 2.491345 -0.000045
H          -4.606166 -1.246230 0.000086
H          -4.606166 1.246228 0.000020
H          -2.476931 2.493823 -0.000071

调整后的结构：

% MEM=2GB
% NPROCLIND=1
% NPROCSHARED=4

# B3LYP/6-31G** nosymm

0opt tiaozheng

0 1
C       3.660096 -0.712985  0.000020
C       2.479166 -1.406940 -0.000050
C       1.223657 -0.722698 -0.000030
C       1.223657  0.722698 -0.000020
C       2.479166  1.406940  0.000020
C       3.660106  0.712985  0.000060
C       0.000000 -1.403568 -0.000020
C       0.000000  1.403578 -0.000030
C       -1.223657  0.722698 -0.000020
C       -1.223657 -0.722698 -0.000010
C       -2.479166 -1.406940  0.000030
H       -2.476989 -2.493834  0.000060
C       -3.660106 -0.712985  0.000050
C       -3.660106  0.712985  0.000010
C       -2.479166  1.406940 -0.000030
H       -0.000010 -2.491284 -0.000040
H       4.606205 -1.246230  0.000040
H       2.476979 -2.493834 -0.000090
H       2.476989  2.493834  0.000030
H       4.606195  1.246230  0.000110
H       -0.000010  2.491294 -0.000050
H       -4.606205 -1.246230  0.000090
H       -4.606205  1.246230  0.000020
H       -2.476989  2.493834 -0.000070，
mengfc

初试构型肯定有影响，高斯没有那么智能，你加上OPT+FREQ就找到的一定是基态！其实他找到的是一个稳定结构！
jeanie163

可以看看分子势能面。初始结构对得到的结构有很大的影响。得到的只是一个范围内相对能量最低的结构。
hbnuwangshuling

我在调整后的结构基础上继续优化，结果得到能量是SCF Done: E(RB+HF-LYP) = -539.546298975，比我在调整的结构上只求单点能（ SCF Done: E(RB+HF-LYP) = -539.546359942 ）高，那岂不是越优化，越能量升高了？即使是局域最小，也不应该吧？而且比我最初没有调整时的结构优化的能量还高呢。这是怎么一回事？
aaacjj

# B3LYP/6-31G** nosymm,你这个是计算优化的？好像算的是电子密度吧
erylingjet

引用回帖:
Originally posted by aaacjj at 2009-11-2 09:52:
# B3LYP/6-31G** nosymm,你这个是计算优化的？好像算的是电子密度吧

单点吧