 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2980 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
huazhorg鐵桿木蟲 (正式寫手)
|
[交流]
【求助】ph.x計(jì)算振動頻率的問題 已有5人參與
|
||
| 我想計(jì)算一個吸附在硅表面上的水分子的振動頻率。所有的原子都計(jì)算的話，我試過，太慢了，由于我只想計(jì)算水分子的振動頻率，就把表面硅原子全部固定。方法是在ph.x的輸入文件里加上這樣一項(xiàng)：nat_todo=3，并給出水分子三個原子的原子系數(shù)。能跑起來，但是算完后，在out文件里卻沒有任何振動頻率的結(jié)果。反復(fù)試了幾次都是如此，還以為是bug,換用最新的QE4.2版本算也是同樣的結(jié)果。問題出在哪里？請高人指點(diǎn)迷津... |
回復(fù)此樓» 猜你喜歡
一志愿西電085401數(shù)一英一299求調(diào)劑 六級521
已經(jīng)有3人回復(fù)
-
生物考研337分求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有14人回復(fù)
-
285求調(diào)劑
已經(jīng)有10人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有29人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有5人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國科學(xué)院大學(xué)265求調(diào)劑
已經(jīng)有7人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 各位大俠，我可以優(yōu)化結(jié)構(gòu)，但只要一計(jì)算頻率就出現(xiàn)下面的問題？
已經(jīng)有4人回復(fù)
- 【交流】振動模式與紅外光譜、拉曼光譜的原理與計(jì)算問題的討論
已經(jīng)有14人回復(fù)
- 振動模式，聲子支與對稱性分類--期待熱心的蟲友
已經(jīng)有11人回復(fù)
- 【求助】請教一下關(guān)于Gaussian計(jì)算振動頻率的問題
已經(jīng)有7人回復(fù)
- 【求助】gaussian頻率計(jì)算問題
已經(jīng)有9人回復(fù)
- 【求助】MS5.0中Dmol3計(jì)算頻率問題
已經(jīng)有21人回復(fù)
- 【求助】如何分析ph.x計(jì)算振動頻率的結(jié)果
已經(jīng)有4人回復(fù)
- 【求助】用MS中的castep計(jì)算聲子振動模頻率，具體怎么算
已經(jīng)有10人回復(fù)
huazhorg
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 21 (小學(xué)生)
- 金幣: 4953.8
- 散金: 5
- 紅花: 2
- 帖子: 381
- 在線: 363.2小時
- 蟲號: 629899
- 注冊: 2008-10-18
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2010-05-31 18:48:03
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-05-31 12:52:42
huazhorg
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 21 (小學(xué)生)
- 金幣: 4953.8
- 散金: 5
- 紅花: 2
- 帖子: 381
- 在線: 363.2小時
- 蟲號: 629899
- 注冊: 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 最下面的三個數(shù)，即是水的三個原子在超胞中的系號.（即1-24號原子都是硅原子) 其中輸出文件最后的部分內(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小時
- 蟲號: 1007127
- 注冊: 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中，使用多個點(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)整是針對4.1.2版本的。試試看。 [ Last edited by valenhou001 on 2010-5-31 at 17:00 ] |
4樓2010-05-31 16:57:58
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290分調(diào)劑求助 +4 | 吉祥止止陳 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +8 | 材料學(xué)257求調(diào)劑 2026-03-28 | 9/450 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
|
[考研] 合肥區(qū)域性重點(diǎn)一本招收調(diào)劑 +4 | 6266jl 2026-03-30 | 4/200 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +8 | hyl3153942 2026-03-29 | 8/400 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|