计算化学公社

标题: 消除虚频后能量反而升高 [打印本页]

作者
Author: freeEnergy 时间: 2017-12-23 20:56
标题: 消除虚频后能量反而升高
各位老师，同学们好，最近优化结构做振动分析，体系是一个石墨烯的周期性体系（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

作者
Author: sobereva 时间: 2017-12-23 21:39
应当说明具体用的是什么程序、升高了多少

作者
Author: freeEnergy 时间: 2017-12-25 05:04
本帖最后由 freeEnergy 于 2017-12-25 05:06 编辑

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)