24小时热门版块排行榜    

查看: 1029  |  回复: 2

liuxianlv

金虫 (正式写手)

[交流] 【求助】分子A和B的opt+freq结果不同 已有2人参与

以下是分子A和B的opt+freq,1,A优化后是4个YES频率计算后却变为2个YES了;B是优化后是4个YES频率计算后仍为4个YES;2,AB分子优化结束都是-- Stationary point found.而频率计算结束时AB不同,A没有-- Stationary point found.而B有-- Stationary point found.

我的疑问:1、A分子这样的结果可用吗?2、A分子为什么会变成2个YES,(优化时已经算了,频率算时为什么还要再 Optimization completed.
-- Stationary point found.)3、A频率结束没有Stationary point found是不是意味着没找到稳定点。4、AB为什么会有差异?

分子A的opt+freq
%nprocshared=4
Will use up to    4 processors via shared memory.
%rwf=LOC1,2000MB,LOC2,2000MB,LOC3,2000MB,LOC4,2000MB,LOC5,2000MB,LOC6,2000MB,LOC
7, 2000MB,LOC8, 2000MB
%mem=200MW
%chk=D:\Data1\Scratch\ams.chk
Default route:  MaxDisk=40GB
----------------------------------------------------
#p opt freq b3lyp/6-31g(d) geom=connectivity pop=reg
----------------------------------------------------
..........................................................................
省略
...........................................................................
优化计算

Item               Value     Threshold  Converged?
Maximum Force            0.000036     0.000450     YES
RMS     Force            0.000009     0.000300     YES
Maximum Displacement     0.000116     0.001800     YES
RMS     Displacement     0.000026     0.001200     YES
Predicted change in Energy=-8.960091D-07
Optimization completed.
    -- Stationary point found.
..........................................................................
省略
...........................................................................
Leave Link 9999 at Thu Jun 10 03:25:03 2010, MaxMem=  209715200 cpu:       0.0
Job cpu time:  0 days  9 hours 39 minutes 58.0 seconds.
File lengths (MBytes):  RWF=    439 Int=      0 D2E=      0 Chk=     12 Scr=      1
Normal termination of Gaussian 03 at Thu Jun 10 03:25:03 2010.
(Enter C:\G03W\l1.exe)
Link1:  Proceeding to internal job step number  2.
------------------------------------------------------------------
#P Geom=AllCheck Guess=Read SCRF=Check GenChk RB3LYP/6-31G(d) Freq
------------------------------------------------------------------
频率计算
..........................................................................
省略
...........................................................................
Item               Value     Threshold  Converged?
Maximum Force            0.000036     0.000450     YES
RMS     Force            0.000009     0.000300     YES
Maximum Displacement     0.056686     0.001800     NO
RMS     Displacement     0.019196     0.001200     NO
Predicted change in Energy=-2.951346D-06
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
..........................................................................
省略
............................................................................
Job cpu time:  0 days  4 hours 25 minutes 33.0 seconds.
File lengths (MBytes):  RWF=    489 Int=      0 D2E=      0 Chk=     12 Scr=      1
Normal termination of Gaussian 03 at Thu Jun 10 07:50:37 2010.

另一个分子分子B的opt+freq
..........................................................................
省略输入文件
...........................................................................

    Item               Value     Threshold  Converged?
Maximum Force            0.000022     0.000450     YES
RMS     Force            0.000006     0.000300     YES
Maximum Displacement     0.001759     0.001800     YES
RMS     Displacement     0.000337     0.001200     YES
Predicted change in Energy=-3.095954D-08
Optimization completed.
    -- Stationary point found.
..........................................................................
省略
...........................................................................
Job cpu time:  0 days  2 hours 15 minutes 25.0 seconds.
File lengths (MBytes):  RWF=    439 Int=      0 D2E=      0 Chk=      7 Scr=      1
Normal termination of Gaussian 03 at Tue Jun 08 14:37:20 2010.
(Enter C:\G03W\l1.exe)
Link1:  Proceeding to internal job step number  2.
------------------------------------------------------------------
#P Geom=AllCheck Guess=Read SCRF=Check GenChk RB3LYP/6-31G(d) Freq优化结束频率开始
------------------------------------------------------------------
..........................................................................
省略
...........................................................................

Item               Value     Threshold  Converged?
Maximum Force            0.000022     0.000450     YES
RMS     Force            0.000006     0.000300     YES
Maximum Displacement     0.001174     0.001800     YES
RMS     Displacement     0.000298     0.001200     YES
Predicted change in Energy=-2.355722D-08
Optimization completed.
    -- Stationary point found.

..........................................................................
省略
...........................................................................

Job cpu time:  0 days  0 hours 43 minutes  3.0 seconds.
File lengths (MBytes):  RWF=    439 Int=      0 D2E=      0 Chk=      7 Scr=      1
Normal termination of Gaussian 03 at Tue Jun 08 15:20:24 2010.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):谢谢 2010-06-12 10:13:50
New Gaussian users are often surprised to see that the final part frequency calculation output that looks that of a geometry optimization at the beginning of a frequency job:

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Initialization pass.Link 103, which performs geometry optimizations, is executed at the beginning and end of all frequency calculations. This is done so that the quadratic optimization step can be computed using the correct second derivatives.
前途光明,出路难觅!
2楼2010-06-12 08:28:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流!! 2010-06-13 11:02:06
频率算完了再算一步优化时为了验证得到的结构是否真的是最小。因为在优化的过程中,除非使用CalcAll,那么都是用一阶导数去外推二阶导数,得到的二阶导数不一定是正确的。所以程序最后要验证一下。
Anewday,anewhour,anewminute,anewpeople.
3楼2010-06-12 13:17:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuxianlv 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-30 19:26 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:29 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:33 by hZiFeudZyoGR
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见