| 查看: 6981 | 回復(fù): 24 | |||||||
[交流]
【求助】VASP頻率分析 已有11人參與
|
|||||||
|
請問高手們,我想得到氫分子的振動(dòng)頻率vibrational frequency(cm-1),設(shè)置了NSW=1, IBRION=5, POTIM=0.015,OUTCAR中得到下面的結(jié)果,請問哪個(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è)頻率對應(yīng)于什么類型的振動(dòng)。第一支頻率就是2個(gè)H原子的相對振動(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ì)算頻率就是能量對位置求二次導(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ù)分子本身的對稱性來做振動(dòng)分析的,跟相對坐標(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方向。 我是猜測你的POSCAR紅,兩個(gè)H原子的x,y坐標(biāo)是相同的,兩個(gè)原子沿著z方向排列。要不你帖出POSCAR看看? 出現(xiàn)這么多頻率是程序自己折騰的,程序試圖給出每一個(gè)(根據(jù)你的設(shè)置)可能的振動(dòng)的頻率。特征頻率是H-H沿鍵長方向的振動(dòng),我猜”碰巧“就是你這里的1f。你的計(jì)算數(shù)值和實(shí)驗(yàn)數(shù)值吻合的很好。。。! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|