 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2534 | 回復: 17 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
donggaomin金蟲 (正式寫手)
|
[求助]
請問VASP5.2中如何設置INCAR才可以計算（gamma點的）phonon frequency？
|
||
|
請問VASP5.2中如何設置INCAR才可以計算（gamma點的）phonon frequency？ 說明書上只說了設置IBRION = 6 另外，怎么查看結果？論壇里提到是動力學矩陣的對角化，不是很明白。我知道動力學矩陣對角化后是頻率的平方。 謝謝！ VASP5.2中是不是可以計算phonon frequency？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 聲子譜 |
» 猜你喜歡
總分293求調劑
已經有3人回復
-
生物技術與工程
已經有3人回復
-
299求調劑
已經有9人回復
-
材料學碩，求調劑
已經有8人回復
-
284求調劑
已經有14人回復
-
292求調劑
已經有4人回復
-
287求調劑
已經有11人回復
-
086003食品工程求調劑
已經有6人回復
-
295求調劑
已經有8人回復
-
考研化學308分求調劑
已經有11人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 急問 現(xiàn)在哪里的計算中心有vasp5.2.12？
已經有6人回復
- time-on-stream什么意思啊
已經有10人回復
- 全國第二屆理論與計算化學軟件Molpro Molcas培訓班
已經有19人回復
- 很菜的問題：不設置參比電極不行么？
已經有9人回復
- 請教有機分子中O原子和H原子電荷數(shù)怎么計算，O-H鍵鍵級怎么計算！謝謝！
已經有4人回復
- 【求助】求VASP計算ZnO能帶和態(tài)密度的詳細輸入設置
已經有15人回復
- 計劃書，word2007,宏設置問題
已經有5人回復
- ThinkPad SL410 2842EKC如何恢復出廠設置？
已經有9人回復
stractor
金蟲 (著名寫手)
- 1ST強帖: 6
- 應助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導體材料
【答案】應助回帖
★ ★ ★
donggaomin(金幣+5): 謝謝stractor!學習了！ 2011-05-29 21:30:19
zzy870720z(金幣+3): 謝謝指教 2011-05-29 23:32:50
donggaomin(金幣+5): 謝謝stractor!學習了！ 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不能太大。 太小的話位移后構型的能量和沒位移的構型的能量沒差別，也會導致誤差。 還有一個誤差，有些體系會遇到，就是VASP中計算頻率是簡諧近似，忽略了能量對位移泰勒展開的高階項。 我測試過H2,CO分子，合理選擇位移的話，計算值和實驗值差別在10cm-1以內。 |
9樓2011-05-29 21:18:00
donggaomin
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
補充：這是在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)在的問題是：是不是第二個的方式好點？如何查看結果？ |
2樓2011-05-28 16:27:04
donggaomin
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2011-05-28 20:42:54
donggaomin
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 819
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
下面是設置了IBRION = 6時的動力學矩陣本征值，按我的理解是本征值就是頻率的平方，但... 請問怎么看下面的結果呢？ 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物技術與工程 +3 | 1294608413 2026-03-25 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 294求調劑材料與化工專碩 +16 | 陌の森林 2026-03-18 | 16/800 |
|
|
|
[考研] 303求調劑 +5 | 元夕元 2026-03-20 | 6/300 |
|
|
|
[考研] 化學調劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
|
[考研] 300求調劑，材料科學英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
|
[考研] 327求調劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 384求調劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +5 | 葵梓衛(wèi)隊 2026-03-18 | 7/350 |
|
|
|
[考研]
 一志愿東華大學化學070300，求調劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 311求調劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 303求調劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
|
[考研] 354求調劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 321求調劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] A區(qū)線材料學調劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 0856調劑，是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|