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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】gaussian03 计算中虚频的处理
汕头大学海洋科学接受调剂
52/1返回列表
查看: 1857  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

majun04

木虫 (著名写手)

majun

[交流] 【求助】gaussian03 计算中虚频的处理

在优化有机物的结构式过程中，发现四个条件都满足了，然后在计算频率的过程中，有虚频出现。一个是N－H的弯曲振动，一个是甲基上C－H键的弯曲振动。虚频都不大，请问怎么消除虚频？
回复此楼

» 猜你喜欢
乐观、执着、勤奋
1楼 2008-05-16 08:42:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingfei101

我确实遇到了个只有第二个是yes的
Berny optimization.
Search for a local minimum.
Step number   1 out of a maximum of   2
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Second derivative matrix not updated -- analytic derivatives used.
     Eigenvalues ---   -0.00015   0.00000   0.00000   0.00135   0.00307
     Eigenvalues ---    0.00380   0.00452   0.00812   0.01054   0.01063
     Eigenvalues ---    0.01316   0.01901   0.01918   0.02749   0.03003
     Eigenvalues ---    0.08853   0.11355   0.12261   0.13353   0.14290
     Eigenvalues ---    0.17766   0.19793   0.27267   0.31526   0.32081
     Eigenvalues ---    0.33146   0.35219   0.36903   0.41685   0.43097
     Eigenvalues ---    0.45464   0.45939   0.45944   0.46027   0.46043
     Eigenvalues ---    0.46099   0.46266   0.46457   0.47521
Eigenvalue   1 out of range, new value =    0.000154 Eigenvector:
                           1
           B1           0.00009
           B2           0.00007
           B3          -0.00408
           B4          -0.00013
           B5          -0.00014
           B6           0.00041
           B7          -0.00089
           B8          -0.00001
           B9          -0.00005
           B10         -0.00002
           B11          0.00043
           B12          0.00000
           B13         -0.00003
           B14          0.00000
           A1           0.00058
           A2          -0.01644
           A3           0.01825
           A4           0.02712
           A5          -0.04557
           A6           0.01676
           A7          -0.00036
           A8           0.00087
           A9          -0.00038
           A10          0.00019
           A11         -0.00052
           A12          0.00002
           A13          0.00067
           D1          -0.57792
           D2           0.17353
           D3           0.20714
           D4           0.19064
           D5           0.74297
           D6           0.01426
           D7           0.01510
           D8           0.01569
           D9          -0.01565
           D10          0.00045
           D11          0.00257
           D12          0.00071
Eigenvalue   2 out of range, new value =    0.000000 Eigenvector:
                           1
           B1          -0.00019
           B2          -0.00015
           B3          -0.00146
           B4           0.00002
           B5           0.00003
           B6          -0.00003
           B7          -0.00500
           B8           0.00005
           B9          -0.00002
           B10         -0.00004
           B11         -0.00758
           B12          0.00004
           B13         -0.00047
           B14         -0.00001
           A1           0.00242
           A2           0.00462
           A3          -0.00123
           A4          -0.00879
           A5           0.01029
           A6           0.00824
           A7          -0.02532
           A8          -0.00377
           A9           0.02687
           A10          0.00203
           A11         -0.00846
           A12          0.00429
           A13          0.01392
           D1           0.44774
           D2           0.46726
           D3           0.46013
           D4           0.46149
           D5          -0.02277
           D6           0.19537
           D7           0.19347
           D8           0.19042
           D9          -0.19950
           D10          0.00267
           D11          0.04681
           D12          0.00982
Quadratic step=1.192D+01 exceeds max=3.000D-01 adjusted using Lamda=-2.966D-05.
Angle between NR and scaled steps=  49.51 degrees.
Angle between quadratic step and forces=  76.36 degrees.
Linear search not attempted -- first point.
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                 (Linear)    (Quad)   (Total)
    B1        1.83918   0.00008   0.00000   0.00023   0.00023   1.83940
    B2        1.84411   0.00033   0.00000   0.00067   0.00067   1.84478
    B3        5.61514   0.00025   0.00000   0.00741   0.00741   5.62255
    B4        1.91915  -0.00002   0.00000   0.00002   0.00002   1.91917
    B5        1.91924  -0.00003   0.00000  -0.00027  -0.00027   1.91897
    B6        1.91934  -0.00006   0.00000  -0.00011  -0.00011   1.91923
    B7        5.48734   0.00005   0.00000  -0.00188  -0.00188   5.48546
    B8        1.91873  -0.00002   0.00000  -0.00034  -0.00034   1.91838
    B9        1.93589   0.00017   0.00000   0.00027   0.00027   1.93616
   B10        1.91849   0.00002   0.00000  -0.00007  -0.00007   1.91842
   B11        5.67358  -0.00016   0.00000  -0.01227  -0.01227   5.66131
   B12        1.91808   0.00004   0.00000   0.00038   0.00038   1.91847
   B13        1.93154   0.00013   0.00000   0.00030   0.00030   1.93184
   B14        1.91883  -0.00006   0.00000  -0.00026  -0.00026   1.91858
    A1        1.87147  -0.00027   0.00000   0.00676   0.00676   1.87822
    A2        0.00044  -0.00097   0.00000  -0.05717  -0.05717  -0.05673
    A3        1.92495   0.00027   0.00000  -0.01593  -0.01593   1.90902
    A4        1.98326  -0.00025   0.00000   0.05957   0.05957   2.04283
    A5        1.92299  -0.00002   0.00000  -0.04540  -0.04540   1.87759
    A6        2.15664   0.00054   0.00000   0.07801   0.07801   2.23465
    A7        2.17168  -0.00002   0.00000   0.01232   0.01232   2.18399
    A8        1.37768  -0.00022   0.00000  -0.00433  -0.00433   1.37335
    A9        2.14772  -0.00021   0.00000  -0.01395  -0.01395   2.13377
   A10        1.12598  -0.00060   0.00000  -0.00189  -0.00189   1.12408
   A11        2.12346  -0.00032   0.00000  -0.04735  -0.04735   2.07610
   A12        0.34481   0.00030   0.00000  -0.00364  -0.00364   0.34117
   A13        1.96959   0.00016   0.00000   0.05402   0.05402   2.02361
    D1        0.32869   0.00019   0.00000  -0.11173  -0.11173   0.21695
    D2        1.17578  -0.00029   0.00000  -0.01063  -0.01063   1.16514
    D3       -2.99893  -0.00041   0.00000   0.01839   0.01839  -2.98054
    D4       -0.89564   0.00052   0.00000   0.02974   0.02974  -0.86590
    D5       -2.22594   0.00020   0.00000   0.06711   0.06711  -2.15883
    D6        1.05762  -0.00029   0.00000   0.08618   0.08618   1.14380
    D7        2.87617   0.00007   0.00000   0.08916   0.08916   2.96533
    D8       -1.58266   0.00006   0.00000   0.08731   0.08731  -1.49535
    D9       -2.86499  -0.00019   0.00000  -0.09798  -0.09798  -2.96296
   D10        2.04016  -0.00004   0.00000  -0.03199  -0.03199   2.00816
   D11        2.93121   0.00004   0.00000   0.12900   0.12900   3.06020
   D12       -1.98387  -0.00009   0.00000  -0.02047  -0.02047  -2.00433
         Item               Value     Threshold  Converged?
Maximum Force            0.000968     0.000450     NO
RMS     Force            0.000282     0.000300     YES
Maximum Displacement     0.128999     0.001800     NO
RMS     Displacement     0.048045     0.001200     NO
Predicted change in Energy=-6.635032D-05
我是用b3lyp/6-311++g**做氨分子的频率时出现的这种情况。请问怎么回事？

[ Last edited by bingfei101 on 2008-5-21 at 21:39 ]
赞一下 回复此楼
11楼2008-05-21 21:33:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

cuihang

铁杆木虫 (正式写手)

VASPer
★ ★ ★
majun04(金币+1,VIP+0):谢谢你
lei0736(金币+2,VIP+0):谢谢
http://muchong.com/bbs/viewthread.php?tid=648127&fpage=6
这篇可能对你有帮助
赞一下(10人) 回复此楼
聪明人知道自己是傻子，傻子以为自己是聪明人
2楼2008-05-16 09:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer
虽然说的是98中的情况，但是大概原理差不多吧
赞一下 回复此楼
聪明人知道自己是傻子，傻子以为自己是聪明人
3楼2008-05-16 09:31:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun
有的时候觉得讲的太泛泛了，对于初学者具体实施起来很困难。
赞一下 回复此楼
乐观、执着、勤奋
4楼2008-05-16 09:33:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +10 汪！？！ 2026-04-09 11/550 2026-04-15 09:17 by fenglj492
[基金申请] 山东省基金2026 +7 jerry681 2026-04-08 10/500 2026-04-15 09:00 by jerry681
[考研] 291 求调剂 +38 化工2026届毕业� 2026-04-09 39/1950 2026-04-15 08:03 by 不我拉绿卡
[考研] 药学305求调剂 +6 玛卡巴卡boom 2026-04-11 6/300 2026-04-14 19:48 by zhouxiaoyu
[考研] 335求调剂 +19 想上岸呀！！ 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 322求调剂，08工科 +4 今天是个小号 2026-04-08 4/200 2026-04-13 00:20 by baobaoye
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 一志愿华中农微生物，288分，三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 求调剂，262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 求调剂，一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
信息提示
关闭
请填处理意见
关闭