24小时热门版块排行榜    

查看: 1686  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

majun04

木虫 (著名写手)

majun

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

在优化有机物的结构式过程中,发现四个条件都满足了,然后在计算频率的过程中,有虚频出现。一个是N-H的弯曲振动,一个是甲基上C-H键的弯曲振动。虚频都不大,请问怎么消除虚频?
回复此楼
乐观、执着、勤奋
已阅   回复此楼   关注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):谢谢
聪明人知道自己是傻子,傻子以为自己是聪明人
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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见