消除虚频后能量反而升高

freeEnergy

各位老师，同学们好，最近优化结构做振动分析，体系是一个石墨烯的周期性体系（31个原子），中间掺杂了一个金属原子，优化和振动的各项参数都是一致的，优化结构后发现存在一个虚频，虚频的方向是金属原子向面外运动，于是调整结构将原子稍微移出表面一点重新优化后，虚频消失了，但是奇怪的是体系能量反而升高了，这是因为什么？麻烦sob老师看看！
下面是有一个虚频的具体：
    1   -0.004770       -24.5       -3.04      -0.735（虚频）
    2   -0.000066        -0.3       -0.04      -0.010
    3    0.000033         0.2        0.02       0.005
    4    0.000349         1.8        0.22       0.054
    5    0.013967        71.8        8.90       2.152
    6    0.023446       120.5       14.94       3.613
    7    0.026470       136.1       16.87       4.079
    8    0.027784       142.8       17.71       4.282
    9    0.028294       145.4       18.03       4.360
   10    0.037013       190.3       23.59       5.704
   11    0.056013       287.9       35.70       8.632
   12    0.059216       304.4       37.74       9.126
   13    0.059424       305.5       37.87       9.158
   14    0.060029       308.6       38.26       9.251
   15    0.066746       343.1       42.54      10.286
   16    0.068483       352.0       43.65      10.554
   17    0.071285       366.4       45.43      10.986
   18    0.077413       397.9       49.34      11.930
   19    0.077896       400.4       49.65      12.004
   20    0.078699       404.6       50.16      12.128
   21    0.079386       408.1       50.60      12.234
   22    0.086174       443.0       54.92      13.280
   23    0.087066       447.6       55.49      13.418
   24    0.089862       461.9       57.27      13.848
   25    0.100695       517.6       64.18      15.518
   26    0.102824       528.6       65.53      15.846
   27    0.105173       540.6       67.03      16.208
   28    0.108385       557.2       69.08      16.703
   29    0.116748       600.1       74.41      17.992
   30    0.118140       607.3       75.30      18.206
   31    0.120451       619.2       76.77      18.562
   32    0.120989       621.9       77.11      18.645
   33    0.121754       625.9       77.60      18.763
   34    0.122142       627.9       77.85      18.823
   35    0.123465       634.7       78.69      19.027
   36    0.126029       647.8       80.32      19.422
   37    0.129089       663.6       82.27      19.894
   38    0.133383       685.7       85.01      20.555
   39    0.136146       699.9       86.77      20.981
   40    0.141070       725.2       89.91      21.740
   41    0.142917       734.7       91.09      22.025
   42    0.144598       743.3       92.16      22.284
   43    0.147330       757.3       93.90      22.705
   44    0.149588       769.0       95.34      23.053
   45    0.151603       779.3       96.62      23.363
   46    0.151881       780.7       96.80      23.406
   47    0.152440       783.6       97.16      23.492
   48    0.153423       788.7       97.78      23.644
   49    0.155776       800.8       99.28      24.006
   50    0.157367       808.9      100.30      24.252
   51    0.157781       811.1      100.56      24.315
   52    0.159627       820.6      101.74      24.600
   53    0.160373       824.4      102.21      24.715
   54    0.169281       870.2      107.89      26.088
   55    0.179975       925.2      114.71      27.736
   56    0.182739       939.4      116.47      28.162
   57    0.193667       995.5      123.43      29.846
   58    0.199300      1024.5      127.02      30.714
   59    0.204033      1048.8      130.04      31.443
   60    0.212478      1092.2      135.42      32.744
   61    0.217261      1116.8      138.47      33.482
   62    0.233225      1198.9      148.64      35.942
   63    0.242800      1248.1      154.75      37.417
   64    0.244335      1256.0      155.73      37.654
   65    0.245136      1260.1      156.24      37.777
   66    0.248066      1275.2      158.10      38.229
   67    0.250200      1286.2      159.46      38.558
   68    0.252409      1297.5      160.87      38.898
   69    0.252702      1299.0      161.06      38.943
   70    0.255312      1312.4      162.72      39.346
   71    0.256660      1319.4      163.58      39.553
   72    0.257109      1321.7      163.87      39.623
   73    0.258178      1327.2      164.55      39.787
   74    0.261430      1343.9      166.62      40.288
   75    0.263101      1352.5      167.69      40.546
   76    0.265911      1366.9      169.48      40.979
   77    0.268544      1380.4      171.15      41.385
   78    0.271405      1395.2      172.98      41.826
   79    0.272279      1399.6      173.54      41.960
   80    0.272940      1403.0      173.96      42.062
   81    0.273835      1407.6      174.53      42.200
   82    0.274823      1412.7      175.16      42.352
   83    0.276990      1423.9      176.54      42.686
   84    0.278871      1433.5      177.74      42.976
   85    0.286448      1472.5      182.57      44.144
   86    0.289377      1487.5      184.43      44.595
   87    0.289456      1487.9      184.48      44.607
   88    0.295428      1518.6      188.29      45.528
   89    0.298001      1531.9      189.93      45.924
   90    0.298539      1534.6      190.27      46.007
   91    0.305818      1572.1      194.91      47.129
   92    0.307803      1582.3      196.18      47.435
   93    0.311243      1599.9      198.37      47.965

sobereva

应当说明具体用的是什么程序、升高了多少

freeEnergy

sobereva 发表于 2017-12-23 21:39
应当说明具体用的是什么程序、升高了多少

用的程序是DMol3，在里面直接做的振动分析，升高了大概1ev，并且在没有虚频的那个结构中，我稍微将那个原子拉近表面，结果优化后原子会回到面内同时能量下降了，做振动分析发现就有虚频了，并且振动模式与最开始的一致，还是那个原子的振动，这样是不是已经可以确定原子不在面内是稳定的结构，但是能量为什么会升高呢？