应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】opt+freq之后能量为什么不是最低？
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 1165  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

hbnuwangshuling

铜虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-11-01 16:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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
赞一下 回复此楼
2楼2009-11-01 17:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 11-1 20:12
初试构型肯定有影响，高斯没有那么智能，你加上OPT+FREQ就找到的一定是基态！其实他找到的是一个稳定结构！
赞一下(11人) 回复此楼
3楼2009-11-01 17:15:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★
lei0736(金币+1,VIP+0):谢谢 11-1 20:12
很简单：上次找到的是局域最小点
这次找到的……可能是局域最小点也可能是全局最小点
赞一下(10人) 回复此楼
4楼2009-11-01 17:20:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jeanie163

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 11-2 08:18
可以看看分子势能面。初始结构对得到的结构有很大的影响。得到的只是一个范围内相对能量最低的结构。
赞一下(11人) 回复此楼
5楼2009-11-01 20:21:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbnuwangshuling

铜虫 (小有名气)

为什么优化之后能量仍然高呢？

我在调整后的结构基础上继续优化，结果得到能量是SCF Done:  E(RB+HF-LYP) =  -539.546298975，比我在调整的结构上只求单点能（ SCF Done:  E(RB+HF-LYP) =  -539.546359942 ）高，那岂不是越优化，越能量升高了？即使是局域最小，也不应该吧？而且比我最初没有调整时的结构优化的能量还高呢。这是怎么一回事？
赞一下 回复此楼
6楼2009-11-02 09:05:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
erylingjet(金币+1,VIP+0):感谢交流 11-2 09:58
# B3LYP/6-31G** nosymm,你这个是计算优化的？好像算的是电子密度吧
赞一下(11人) 回复此楼
7楼2009-11-02 09:52:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by aaacjj at 2009-11-2 09:52:
# B3LYP/6-31G** nosymm,你这个是计算优化的？好像算的是电子密度吧

单点吧
回复此楼
8楼2009-11-02 09:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by hbnuwangshuling at 2009-11-2 09:05:
我在调整后的结构基础上继续优化，结果得到能量是SCF Done:  E(RB+HF-LYP) =  -539.546298975，比我在调整的结构上只求单点能（ SCF Done:  E(RB+HF-LYP) =  -539.546359942 ）高，那岂不是越优化，越能量升高了？ ...

很简单啊。上次你的构型落在一个局域最小点中，这次落在另一个局域最小点中了呗
赞一下(1人) 回复此楼
9楼2009-11-02 10:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaacjj

铜虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
搞错了，呵呵，不好意思撒
赞一下(1人) 回复此楼
10楼2009-11-02 10:11:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hbnuwangshuling 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323分（计算机视觉和大模型项目）能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 085601，一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 电子信息调剂交叉学科有推荐吗 +6 jhtfeybgj 2026-04-01 8/400 2026-04-04 07:52 by 1753564080
[考研] 求调剂不挑专业 +3 xrh030412 2026-04-01 3/150 2026-04-03 14:40 by 氮气气气
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 建环，能源，土木老师路过看一看！！！ +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 085600，320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 材料工程322分 +8 哈哈哈吼吼吼哈 2026-04-01 8/400 2026-04-02 11:53 by 3041
[考研] 272求调剂，接受跨专业调剂！ +4 闲鱼卢 2026-03-31 4/200 2026-04-02 11:18 by guyan1000
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 生物学296求调剂 +10 汤圆包 2026-03-29 14/700 2026-04-01 10:44 by 求调剂zz
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 一志愿中海洋材料357 +4 麦恩莉. 2026-03-30 4/200 2026-03-31 14:35 by 记事本2026
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
信息提示
关闭
请填处理意见
关闭