 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2533 | 回復(fù): 17 | |||
donggaomin金蟲 (正式寫手)
|
[求助]
請問VASP5.2中如何設(shè)置INCAR才可以計算（gamma點的）phonon frequency？
|
|
請問VASP5.2中如何設(shè)置INCAR才可以計算（gamma點的）phonon frequency？ 說明書上只說了設(shè)置IBRION = 6 另外，怎么查看結(jié)果？論壇里提到是動力學(xué)矩陣的對角化，不是很明白。我知道動力學(xué)矩陣對角化后是頻率的平方。 謝謝！ VASP5.2中是不是可以計算phonon frequency？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 聲子譜 |
» 猜你喜歡
材料學(xué)碩333求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有7人回復(fù)
-
收08調(diào)劑生
已經(jīng)有6人回復(fù)
-
材料專碩找調(diào)劑
已經(jīng)有5人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有5人回復(fù)
-
化工專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
274求調(diào)劑
已經(jīng)有5人回復(fù)
-
調(diào)劑
已經(jīng)有3人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有15人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有16人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 急問 現(xiàn)在哪里的計算中心有vasp5.2.12？
已經(jīng)有6人回復(fù)
- time-on-stream什么意思啊
已經(jīng)有10人回復(fù)
- 全國第二屆理論與計算化學(xué)軟件Molpro Molcas培訓(xùn)班
已經(jīng)有19人回復(fù)
- 很菜的問題：不設(shè)置參比電極不行么？
已經(jīng)有9人回復(fù)
- 請教有機分子中O原子和H原子電荷數(shù)怎么計算，O-H鍵鍵級怎么計算！謝謝！
已經(jīng)有4人回復(fù)
- 【求助】求VASP計算ZnO能帶和態(tài)密度的詳細輸入設(shè)置
已經(jīng)有15人回復(fù)
- 計劃書，word2007,宏設(shè)置問題
已經(jīng)有5人回復(fù)
- ThinkPad SL410 2842EKC如何恢復(fù)出廠設(shè)置？
已經(jīng)有9人回復(fù)
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
補充：這是在VASP論壇中找到的： you can calculate the gamma-point phonons with vasp.4.6, using a frozen phonon approach please use IBRION=5 in that case, and set POTIM = 0.015 (as the harmonic limit is assumed) in case you do not need all modes, please use the Selective dymanics key in POSCAR and just allow the vibrations you are intrested in (by setting F/T ) in order to save CPU time. it is not possible to continue a broken phonon job (keeping the info) it has to be restarted from the beginning. relativistic effects can be simulated by using spin-orbit splitting (please have a look at the LSORBIT tag, and the additional parameters you have to set. 這是第二個 To do \Gamma point phonon calculate Firstly, fully relax your structure and save WAVECAR Secondly, you can start a phonon calculation by set INCAR as PREC= Accurate ADDGRID = .TRUE. ENCUT= 500.00 eV ISTART=1 ICHARG=0 IBRION= 6 NFREE=4 POTIM=0.015 !displacement of atoms in Ang. IALGO= 48 LREAL=Auto !needed for good convergence NSIM = 4 LCHARG = .False. LPLANE=.TRUE. ~ 現(xiàn)在的問題是：是不是第二個的方式好點？如何查看結(jié)果？ |
2樓2011-05-28 16:27:04
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2011-05-28 20:42:54
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
下面是設(shè)置了IBRION = 6時的動力學(xué)矩陣本征值，按我的理解是本征值就是頻率的平方，但... 請問怎么看下面的結(jié)果呢？ Eigenvectors and eigenvalues of the dynamical matrix ---------------------------------------------------- 1 f = 15.937321 THz 100.137142 2PiTHz 531.611795 cm-1 65.911486 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.000519 -0.000028 0.081397 1.609613 0.929267 2.583248 -0.000519 -0.000028 -0.081397 -0.000113 1.858561 1.976701 0.001581 0.000057 -0.702404 1.609613 0.929267 4.559948 0.001581 0.000057 0.702404 2 f = 13.134427 THz 82.526040 2PiTHz 438.117317 cm-1 54.319644 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.092238 -0.254007 0.000249 1.609613 0.929267 2.583248 0.092238 0.254007 0.000249 -0.000113 1.858561 1.976701 -0.222415 -0.614413 -0.000695 1.609613 0.929267 4.559948 0.222415 0.614413 -0.000695 3 f = 13.116478 THz 82.413261 2PiTHz 437.518595 cm-1 54.245412 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.254559 -0.092022 0.000439 1.609613 0.929267 2.583248 -0.254559 0.092022 0.000439 -0.000113 1.858561 1.976701 0.614185 -0.222503 -0.000963 1.609613 0.929267 4.559948 -0.614185 0.222503 -0.000963 4 f = 12.231166 THz 76.850683 2PiTHz 407.987772 cm-1 50.584055 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.128882 -0.287868 -0.000033 1.609613 0.929267 2.583248 -0.128882 -0.287868 0.000033 -0.000113 1.858561 1.976701 0.258437 0.577694 0.000731 1.609613 0.929267 4.559948 0.258437 0.577694 -0.000731 5 f = 12.215895 THz 76.754733 2PiTHz 407.478389 cm-1 50.520900 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.288034 0.128850 -0.000113 1.609613 0.929267 2.583248 -0.288034 0.128850 0.000113 -0.000113 1.858561 1.976701 0.577610 -0.258454 0.001474 1.609613 0.929267 4.559948 0.577610 -0.258454 -0.001474 6 f = 10.927284 THz 68.658148 2PiTHz 364.494939 cm-1 45.191629 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.000189 0.000636 0.313532 1.609613 0.929267 2.583248 0.000189 -0.000636 0.313532 -0.000113 1.858561 1.976701 -0.000728 0.000873 -0.633795 1.609613 0.929267 4.559948 0.000728 -0.000873 -0.633795 7 f = 7.595291 THz 47.722620 2PiTHz 253.351625 cm-1 31.411609 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.000310 0.000022 0.702406 1.609613 0.929267 2.583248 0.000310 0.000022 -0.702406 -0.000113 1.858561 1.976701 0.000072 -0.000016 0.081397 1.609613 0.929267 4.559948 0.000072 -0.000016 -0.081397 8 f = 2.947281 THz 18.518311 2PiTHz 98.310698 cm-1 12.188977 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.129670 0.640455 -0.000783 1.609613 0.929267 2.583248 -0.129670 -0.640455 -0.000783 -0.000113 1.858561 1.976701 -0.053713 -0.264796 -0.000086 1.609613 0.929267 4.559948 0.053713 0.264796 -0.000086 9 f = 2.934466 THz 18.437793 2PiTHz 97.883244 cm-1 12.135980 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.640217 -0.129724 0.000119 1.609613 0.929267 2.583248 -0.640217 0.129724 0.000119 -0.000113 1.858561 1.976701 -0.265373 0.053582 0.000117 1.609613 0.929267 4.559948 0.265373 -0.053582 0.000117 10 f/i= 0.003591 THz 0.022565 2PiTHz 0.119792 cm-1 0.014852 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.000007 -0.000665 -0.633795 1.609613 0.929267 2.583248 -0.000007 0.000665 -0.633795 -0.000113 1.858561 1.976701 -0.000005 0.000373 -0.313533 1.609613 0.929267 4.559948 0.000005 -0.000373 -0.313533 11 f/i= 0.036145 THz 0.227107 2PiTHz 1.205674 cm-1 0.149485 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.063064 0.629725 -0.000046 1.609613 0.929267 2.583248 0.063064 0.629725 0.000046 -0.000113 1.858561 1.976701 0.031392 0.313821 0.000019 1.609613 0.929267 4.559948 0.031392 0.313821 -0.000019 12 f/i= 0.039877 THz 0.250556 2PiTHz 1.330159 cm-1 0.164919 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 0.629641 -0.063053 -0.000333 1.609613 0.929267 2.583248 0.629641 -0.063053 0.000333 -0.000113 1.858561 1.976701 0.313989 -0.031415 0.000203 1.609613 0.929267 4.559948 0.313989 -0.031415 -0.000203 |
4樓2011-05-29 17:04:52
stractor
金蟲 (著名寫手)
- 1ST強帖: 6
- 應(yīng)助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導(dǎo)體材料
5樓2011-05-29 19:55:42
stractor
金蟲 (著名寫手)
- 1ST強帖: 6
- 應(yīng)助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導(dǎo)體材料
【答案】應(yīng)助回帖
★ ★
donggaomin(金幣+5): 謝謝！請問1f,2f 3f 4f 5f ... f/i 等分別表示什么？如何判斷振動類型？ 2011-05-29 21:21:12
zzy870720z(金幣+2): 謝謝指教 2011-05-29 23:32:18
donggaomin(金幣+5): 謝謝！請問1f,2f 3f 4f 5f ... f/i 等分別表示什么？如何判斷振動類型？ 2011-05-29 21:21:12
zzy870720z(金幣+2): 謝謝指教 2011-05-29 23:32:18
|
舉例 1 f = 15.937321 THz 100.137142 2PiTHz 531.611795 cm-1 65.911486 meV X Y Z dx dy dz -0.000113 1.858561 0.000000 -0.000519 -0.000028 0.081397 1.609613 0.929267 2.583248 -0.000519 -0.000028 -0.081397 -0.000113 1.858561 1.976701 0.001581 0.000057 -0.702404 1.609613 0.929267 4.559948 0.001581 0.000057 0.702404 第一行就是頻率，不同的單位給出。 后面X Y Z是原子坐標 在后面的 dx dy dz是動力學(xué)矩陣的本征矢量，判斷振動類型的時候用到。 |
6樓2011-05-29 20:00:30
dxcharlary
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +83 - 1ST強帖: 2
- 應(yīng)助: 251 (大學(xué)生)
- 金幣: 10420.1
- 散金: 155
- 紅花: 65
- 帖子: 3137
- 在線: 757.6小時
- 蟲號: 936718
- 注冊: 2010-01-04
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
- 管轄: 第一性原理
7樓2011-05-29 20:21:14
dxcharlary
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +83 - 1ST強帖: 2
- 應(yīng)助: 251 (大學(xué)生)
- 金幣: 10420.1
- 散金: 155
- 紅花: 65
- 帖子: 3137
- 在線: 757.6小時
- 蟲號: 936718
- 注冊: 2010-01-04
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
- 管轄: 第一性原理
8樓2011-05-29 20:23:02
stractor
金蟲 (著名寫手)
- 1ST強帖: 6
- 應(yīng)助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導(dǎo)體材料
【答案】應(yīng)助回帖
★ ★ ★
donggaomin(金幣+5): 謝謝stractor!學(xué)習(xí)了！ 2011-05-29 21:30:19
zzy870720z(金幣+3): 謝謝指教 2011-05-29 23:32:50
donggaomin(金幣+5): 謝謝stractor!學(xué)習(xí)了！ 2011-05-29 21:30:19
zzy870720z(金幣+3): 謝謝指教 2011-05-29 23:32:50
|
是的，有限差分的可信度取決于原子的位移 因為有限差分計算Hessen矩陣是 Hij=d2E/dxidxj=[E(ri0+dxi,rj0+dxj)+E(ri0-dxi,rj0-dxj)-2E(ri0,rj0)]/dxidxj 上面的d是微分，rjo就是j原子的平衡位置，E是能量。 上式反映出位移xi, xj不能太大。 太小的話位移后構(gòu)型的能量和沒位移的構(gòu)型的能量沒差別，也會導(dǎo)致誤差。 還有一個誤差，有些體系會遇到，就是VASP中計算頻率是簡諧近似，忽略了能量對位移泰勒展開的高階項。 我測試過H2,CO分子，合理選擇位移的話，計算值和實驗值差別在10cm-1以內(nèi)。 |
9樓2011-05-29 21:18:00
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
10樓2011-05-29 21:41:07
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 donggaomin 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 招08考數(shù)學(xué) +7 | laoshidan 2026-03-20 | 16/800 |
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+3
|
糖葫蘆888ll 2026-03-22 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +7 | A-哆啦Z夢 2026-03-23 | 12/600 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|