| 查看: 6988 | 回復(fù): 24 | |||||||
[交流]
【求助】VASP頻率分析 已有11人參與
|
|||||||
|
請(qǐng)問高手們,我想得到氫分子的振動(dòng)頻率vibrational frequency(cm-1),設(shè)置了NSW=1, IBRION=5, POTIM=0.015,OUTCAR中得到下面的結(jié)果,請(qǐng)問哪個(gè)才是我所需要的氫分子的頻率? 1 f = 127.990789 THz 804.189843 2PiTHz 4269.313033 cm-1 529.327544 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.000036 -0.000036 0.707238 0.000000 0.000000 0.750000 0.000283 -0.000011 -0.706976 2 f = 4.690886 THz 29.473704 2PiTHz 156.471101 cm-1 19.399951 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.748316 -0.000024 0.000172 0.000000 0.000000 0.750000 -0.663342 0.000014 -0.000131 3 f = 4.396660 THz 27.625032 2PiTHz 146.656802 cm-1 18.183133 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.000022 -0.759554 -0.000017 0.000000 0.000000 0.750000 0.000015 0.650444 0.000012 4 f/i= 0.023611 THz 0.148351 2PiTHz 0.787570 cm-1 0.097646 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.004764 -0.000534 -0.706959 0.000000 0.000000 0.750000 0.005330 -0.000629 -0.707218 5 f/i= 0.283258 THz 1.779760 2PiTHz 9.448456 cm-1 1.171460 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.663323 0.001788 0.004920 0.000000 0.000000 0.750000 0.748294 0.002094 0.005187 6 f/i= 0.340255 THz 2.137885 2PiTHz 11.349685 cm-1 1.407182 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.001819 0.650441 -0.000575 0.000000 0.000000 0.750000 -0.002059 0.759551 -0.000620 Eigenvectors after division by SQRT(mass) Eigenvectors and eigenvalues of the dynamical matrix ---------------------------------------------------- 1 f = 127.990789 THz 804.189843 2PiTHz 4269.313033 cm-1 529.327544 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.000036 -0.000036 0.707238 0.000000 0.000000 0.750000 0.000283 -0.000011 -0.706976 2 f = 4.690886 THz 29.473704 2PiTHz 156.471101 cm-1 19.399951 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.748316 -0.000024 0.000172 0.000000 0.000000 0.750000 -0.663342 0.000014 -0.000131 3 f = 4.396660 THz 27.625032 2PiTHz 146.656802 cm-1 18.183133 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.000022 -0.759554 -0.000017 0.000000 0.000000 0.750000 0.000015 0.650444 0.000012 4 f/i= 0.023611 THz 0.148351 2PiTHz 0.787570 cm-1 0.097646 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.004764 -0.000534 -0.706959 0.000000 0.000000 0.750000 0.005330 -0.000629 -0.707218 5 f/i= 0.283258 THz 1.779760 2PiTHz 9.448456 cm-1 1.171460 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 0.663323 0.001788 0.004920 0.000000 0.000000 0.750000 0.748294 0.002094 0.005187 6 f/i= 0.340255 THz 2.137885 2PiTHz 11.349685 cm-1 1.407182 meV X Y Z dx dy dz 0.000000 0.000000 0.000000 -0.001819 0.650441 -0.000575 0.000000 0.000000 0.750000 -0.002059 0.759551 -0.000620 Finite differences POTIM= 1.500000000000000E-002 |
gassian | vasp 學(xué)習(xí)貼 | 有關(guān)零點(diǎn)能計(jì)算的一些帖子 | 零點(diǎn)能 |
金蟲 (著名寫手)
金蟲 (著名寫手)
|
頻率下面給出了本征矢量,分析本征矢量可以知道這個(gè)頻率對(duì)應(yīng)于什么類型的振動(dòng)。第一支頻率就是2個(gè)H原子的相對(duì)振動(dòng),即你所要的振動(dòng)頻率。因?yàn)橐粋(gè)分子的振動(dòng)模式數(shù)是3N,N為原子個(gè)數(shù),所以這里有6支頻率。后面的5支頻率是旋轉(zhuǎn)和剛性平移的頻率,為0。 你計(jì)算的值和實(shí)驗(yàn)值差別較大,可能是如下原因: 最可能的原因出在POTIM。 POTIM=0.015太小了,這樣計(jì)算的位移后構(gòu)型的能量和位移前構(gòu)型的能量可能沒差別哦,可以設(shè)為POTIM=0.4。你可以查查資料,有限差分法計(jì)算頻率就是能量對(duì)位置求二次導(dǎo)數(shù),你POTIM太小,VASP計(jì)算的位移和沒位移的能量沒明顯的差別。理想的POTIM為0.4-1.0。 計(jì)算頻率時(shí)H原子間距是取的優(yōu)化的還是實(shí)驗(yàn)的? 你是采用中心差分法嗎?即NFREE=2嗎? 最好把INCAR貼出來。 [ Last edited by stractor on 2011-6-4 at 02:05 ] |
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
做頻率計(jì)算的時(shí)候,H2可以沿著軸放置。表面吸附分子方向一般不是這么”規(guī)則“的,是不是要求矢量和呢? ————————————————————————————————————— 自己算了一個(gè)例子,實(shí)際上無論分子如何放置,頻率計(jì)算是根據(jù)分子本身的對(duì)稱性來做振動(dòng)分析的,跟相對(duì)坐標(biāo)無關(guān)。 [ Last edited by cenwanglai on 2011-2-23 at 14:27 ] |
|
我按照你說的做了設(shè)置,計(jì)算的結(jié)果是 1 f = 127.818139 THz 803.105052 2PiTHz 4263.554054 cm-1 528.613521 meV 2 f = 3.840658 THz 24.131568 2PiTHz 128.110567 cm-1 15.883692 meV 3 f = 3.610350 THz 22.684500 2PiTHz 120.428321 cm-1 14.931214 meV 4 f/i= 0.003690 THz 0.023185 2PiTHz 0.123083 cm-1 0.015260 meV 5 f/i= 0.584949 THz 3.675344 2PiTHz 19.511806 cm-1 2.419156 meV 6 f/i= 0.654984 THz 4.115385 2PiTHz 21.847909 cm-1 2.708796 meV 1 f = 127.818139 THz 803.105052 2PiTHz 4263.554054 cm-1 528.613521 meV 2 f = 3.840658 THz 24.131568 2PiTHz 128.110567 cm-1 15.883692 meV 3 f = 3.610350 THz 22.684500 2PiTHz 120.428321 cm-1 14.931214 meV 4 f/i= 0.003690 THz 0.023185 2PiTHz 0.123083 cm-1 0.015260 meV 5 f/i= 0.584949 THz 3.675344 2PiTHz 19.511806 cm-1 2.419156 meV 6 f/i= 0.654984 THz 4.115385 2PiTHz 21.847909 cm-1 2.708796 meV |
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
呵呵,可能是我沒有說清楚我的意思。 不是1f就代表z方向。 我是猜測(cè)你的POSCAR紅,兩個(gè)H原子的x,y坐標(biāo)是相同的,兩個(gè)原子沿著z方向排列。要不你帖出POSCAR看看? 出現(xiàn)這么多頻率是程序自己折騰的,程序試圖給出每一個(gè)(根據(jù)你的設(shè)置)可能的振動(dòng)的頻率。特征頻率是H-H沿鍵長(zhǎng)方向的振動(dòng),我猜”碰巧“就是你這里的1f。你的計(jì)算數(shù)值和實(shí)驗(yàn)數(shù)值吻合的很好啊。。。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +20 | 不知道叫什么! 2026-03-15 | 22/1100 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +3 | NIFFFfff 2026-03-20 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|