| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2982 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
huazhorg鐵桿木蟲 (正式寫手)
|
[交流]
【求助】ph.x計算振動頻率的問題 已有5人參與
|
||
| 我想計算一個吸附在硅表面上的水分子的振動頻率。所有的原子都計算的話,我試過,太慢了,由于我只想計算水分子的振動頻率,就把表面硅原子全部固定。方法是在ph.x的輸入文件里加上這樣一項:nat_todo=3,并給出水分子三個原子的原子系數(shù)。能跑起來,但是算完后,在out文件里卻沒有任何振動頻率的結(jié)果。反復(fù)試了幾次都是如此,還以為是bug,換用最新的QE4.2版本算也是同樣的結(jié)果。問題出在哪里?請高人指點迷津... |
金蟲 (正式寫手)

至尊木蟲 (職業(yè)作家)
鐵桿木蟲 (正式寫手)
|
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. 這幾項可以不要,但結(jié)果也是一樣的:沒有聲子。) [ Last edited by huazhorg on 2010-5-31 at 14:40 ] |
至尊木蟲 (職業(yè)作家)
|
為計算出頻率的原因,可能是你的tr2_ph的收斂精度設(shè)置太高,而未收斂,適當(dāng)降低它。 在做ph.x之前進(jìn)行的pw.x中,使用多個點。比如下面的: 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 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +28 | 洛神哥哥 2026-03-31 | 28/1400 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +21 | 新時代材料 2026-03-27 | 23/1150 |
|
|
[考研] 304求調(diào)劑 +12 | 素年祭語 2026-03-31 | 15/750 |
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-04-01 | 3/150 |
|
|
[考研] 英一數(shù)二生物信息學(xué)287分,本科生物科學(xué),求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +5 | 大饞小子 2026-03-27 | 7/350 |
|
|
[考研] 321求調(diào)劑 +8 | y-yh 2026-04-01 | 8/400 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-30 | 7/350 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[考研] 322求調(diào)劑 +8 | 三水sss 2026-04-01 | 8/400 |
|
|
[考研] 282求調(diào)劑 +6 | 呼吸都是減肥 2026-04-01 | 6/300 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +7 | 鐵憨憨123425 2026-03-31 | 7/350 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +11 | dongfang59 2026-03-27 | 11/550 |
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|