 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2977 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
huazhorg鐵桿木蟲 (正式寫手)
|
[交流]
【求助】ph.x計(jì)算振動(dòng)頻率的問題 已有5人參與
|
||
| 我想計(jì)算一個(gè)吸附在硅表面上的水分子的振動(dòng)頻率。所有的原子都計(jì)算的話，我試過，太慢了，由于我只想計(jì)算水分子的振動(dòng)頻率，就把表面硅原子全部固定。方法是在ph.x的輸入文件里加上這樣一項(xiàng)：nat_todo=3，并給出水分子三個(gè)原子的原子系數(shù)。能跑起來(lái)，但是算完后，在out文件里卻沒有任何振動(dòng)頻率的結(jié)果。反復(fù)試了幾次都是如此，還以為是bug,換用最新的QE4.2版本算也是同樣的結(jié)果。問題出在哪里？請(qǐng)高人指點(diǎn)迷津... |
回復(fù)此樓» 猜你喜歡
086502化學(xué)工程342求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有5人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有7人回復(fù)
-
328求調(diào)劑
已經(jīng)有7人回復(fù)
-
340求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有15人回復(fù)
-
291求調(diào)劑
已經(jīng)有19人回復(fù)
-
331環(huán)境科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 各位大俠，我可以優(yōu)化結(jié)構(gòu)，但只要一計(jì)算頻率就出現(xiàn)下面的問題？
已經(jīng)有4人回復(fù)
- 【交流】振動(dòng)模式與紅外光譜、拉曼光譜的原理與計(jì)算問題的討論
已經(jīng)有14人回復(fù)
- 振動(dòng)模式，聲子支與對(duì)稱性分類--期待熱心的蟲友
已經(jīng)有11人回復(fù)
- 【求助】請(qǐng)教一下關(guān)于Gaussian計(jì)算振動(dòng)頻率的問題
已經(jīng)有7人回復(fù)
- 【求助】gaussian頻率計(jì)算問題
已經(jīng)有9人回復(fù)
- 【求助】MS5.0中Dmol3計(jì)算頻率問題
已經(jīng)有21人回復(fù)
- 【求助】如何分析ph.x計(jì)算振動(dòng)頻率的結(jié)果
已經(jīng)有4人回復(fù)
- 【求助】用MS中的castep計(jì)算聲子振動(dòng)模頻率，具體怎么算
已經(jīng)有10人回復(fù)
huazhorg
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 21 (小學(xué)生)
- 金幣: 4953.8
- 散金: 5
- 紅花: 2
- 帖子: 381
- 在線: 363.2小時(shí)
- 蟲號(hào): 629899
- 注冊(cè): 2008-10-18
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
ph.x 的輸入是這樣的： vibrations of H2O &inputph nat_todo=3, tr2_ph=4.0d-17, prefix='H2O', outdir='～/tmp', amass(1)=28.d0, amass(2)=16.d0, amass(3)=1.d0, fildyn='H2O.dyn', asr=.true., zue=.true., epsil=.true., trans=.true., / 0.0 0.0 0.0 25 26 27 最下面的三個(gè)數(shù)，即是水的三個(gè)原子在超胞中的系號(hào).（即1-24號(hào)原子都是硅原子) 其中輸出文件最后的部分內(nèi)容為： ………………………………………… iter # 10 total cpu time : 54288.4 secs av.it.: 21.0 thresh= 0.256E-07 alpha_mix = 0.700 |ddv_scf|^2 = 0.362E-13 iter # 11 total cpu time : 54509.6 secs av.it.: 21.0 thresh= 0.190E-07 alpha_mix = 0.700 |ddv_scf|^2 = 0.429E-14 End of self-consistent calculation Convergence has been achieved Stopping because representation 1 is not done PHONON : 15h 8m CPU time, 15h11m wall time INITIALIZATION: phq_setup : 1.27s CPU phq_init : 6207.15s CPU phq_init : 6207.15s CPU ……………………………………………………………… 可以看出，循環(huán)結(jié)束后，按照常規(guī)，本應(yīng)該給出聲子結(jié)果的，卻沒有，而是直接結(jié)束了，沒有給出任何聲子的數(shù)據(jù)。 （補(bǔ)充：asr=.true., zue=.true., epsil=.true., trans=.true. 這幾項(xiàng)可以不要，但結(jié)果也是一樣的：沒有聲子。） [ Last edited by huazhorg on 2010-5-31 at 14:40 ] |
3樓2010-05-31 14:38:00
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-05-31 12:52:42
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★
huazhorg(金幣+2): 2010-05-31 18:40:35
qasd(金幣+2):謝謝 2010-05-31 20:42:24
huazhorg(金幣+2): 2010-05-31 18:40:35
qasd(金幣+2):謝謝 2010-05-31 20:42:24
|
為計(jì)算出頻率的原因，可能是你的tr2_ph的收斂精度設(shè)置太高，而未收斂，適當(dāng)降低它。 在做ph.x之前進(jìn)行的pw.x中，使用多個(gè)點(diǎn)。比如下面的： K_POINTS {automatic} 6 6 1 1 1 1 將ph.x的輸入文件中下面的： vibrations of H2O &inputph tr2_ph=4.0d-14, prefix='H2O', outdir='～/tmp', amass(1)=28.d0, amass(2)=16.d0, amass(3)=1.d0, asr=.true., epsil=.false., trans=.true., lraman=.false., elop=.false., fildyn='H2O.dyn', fildrho='H2O.drho', nat_todo=3, / 0.0 0.0 0.0 25 26 27 上面的調(diào)整是針對(duì)4.1.2版本的。試試看。 [ Last edited by valenhou001 on 2010-5-31 at 17:00 ] |
4樓2010-05-31 16:57:58
等待的心
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 819.7
- 散金: 620
- 沙發(fā): 1
- 帖子: 372
- 在線: 91.8小時(shí)
- 蟲號(hào): 947336
- 注冊(cè): 2010-01-23
- 性別: GG
- 專業(yè): 物理無(wú)機(jī)化學(xué)
5樓2010-05-31 18:40:17
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +8 | 新時(shí)代材料 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0856，材料與化工321分求調(diào)劑 +8 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+16
|
孅華 2026-03-22 | 16/800 |
|
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +5 | Sxy112 2026-03-21 | 7/350 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|