 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2537 | 回復(fù): 17 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
donggaomin金蟲 (正式寫手)
|
[求助]
請問VASP5.2中如何設(shè)置INCAR才可以計(jì)算（gamma點(diǎn)的）phonon frequency？
|
||
|
請問VASP5.2中如何設(shè)置INCAR才可以計(jì)算（gamma點(diǎn)的）phonon frequency？ 說明書上只說了設(shè)置IBRION = 6 另外，怎么查看結(jié)果？論壇里提到是動力學(xué)矩陣的對角化，不是很明白。我知道動力學(xué)矩陣對角化后是頻率的平方。 謝謝！ VASP5.2中是不是可以計(jì)算phonon frequency？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 聲子譜 |
» 猜你喜歡
085600，材料與化工321分求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑求院校招收
已經(jīng)有3人回復(fù)
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有8人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有11人回復(fù)
-
0856求調(diào)劑
已經(jīng)有10人回復(fù)
-
266分，求材料冶金能源化工等調(diào)劑
已經(jīng)有9人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有8人回復(fù)
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有16人回復(fù)
-
286求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 急問 現(xiàn)在哪里的計(jì)算中心有vasp5.2.12？
已經(jīng)有6人回復(fù)
- time-on-stream什么意思啊
已經(jīng)有10人回復(fù)
- 全國第二屆理論與計(jì)算化學(xué)軟件Molpro Molcas培訓(xùn)班
已經(jīng)有19人回復(fù)
- 很菜的問題：不設(shè)置參比電極不行么？
已經(jīng)有9人回復(fù)
- 請教有機(jī)分子中O原子和H原子電荷數(shù)怎么計(jì)算，O-H鍵鍵級怎么計(jì)算！謝謝！
已經(jīng)有4人回復(fù)
- 【求助】求VASP計(jì)算ZnO能帶和態(tài)密度的詳細(xì)輸入設(shè)置
已經(jīng)有15人回復(fù)
- 計(jì)劃書，word2007,宏設(shè)置問題
已經(jīng)有5人回復(fù)
- ThinkPad SL410 2842EKC如何恢復(fù)出廠設(shè)置？
已經(jīng)有9人回復(fù)
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時(shí)
- 蟲號: 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
- 散金: 820
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時(shí)
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
補(bǔ)充：這是在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. 這是第二個(gè) 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)在的問題是：是不是第二個(gè)的方式好點(diǎn)？如何查看結(jié)果？ |
2樓2011-05-28 16:27:04
donggaomin
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 565
- 在線: 244.3小時(shí)
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
下面是設(shè)置了IBRION = 6時(shí)的動力學(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強(qiáng)帖: 6
- 應(yīng)助: 91 (初中生)
- 金幣: 1487.6
- 散金: 666
- 紅花: 30
- 沙發(fā): 1
- 帖子: 2224
- 在線: 426.8小時(shí)
- 蟲號: 394087
- 注冊: 2007-06-06
- 專業(yè): 半導(dǎo)體材料
5樓2011-05-29 19:55:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +4 | Micky11223 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考博] 26申博 +3 | 加油沖啊！ 2026-03-26 | 3/150 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +6 | 來好運(yùn)來來來 2026-03-21 | 7/350 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|