| 5 | 1/1 | 返回列表 |
| 查看: 2979 | 回復(fù): 15 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
huazhorg鐵桿木蟲 (正式寫手)
|
[交流]
【求助】ph.x計算振動頻率的問題 已有5人參與
|
||
| 我想計算一個吸附在硅表面上的水分子的振動頻率。所有的原子都計算的話,我試過,太慢了,由于我只想計算水分子的振動頻率,就把表面硅原子全部固定。方法是在ph.x的輸入文件里加上這樣一項:nat_todo=3,并給出水分子三個原子的原子系數(shù)。能跑起來,但是算完后,在out文件里卻沒有任何振動頻率的結(jié)果。反復(fù)試了幾次都是如此,還以為是bug,換用最新的QE4.2版本算也是同樣的結(jié)果。問題出在哪里?請高人指點迷津... |
至尊木蟲 (職業(yè)作家)
|
為計算出頻率的原因,可能是你的tr2_ph的收斂精度設(shè)置太高,而未收斂,適當(dāng)降低它。 在做ph.x之前進行的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 ] |
至尊木蟲 (職業(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ù)。 (補充:asr=.true., zue=.true., epsil=.true., trans=.true. 這幾項可以不要,但結(jié)果也是一樣的:沒有聲子。) [ Last edited by huazhorg on 2010-5-31 at 14:40 ] |
金蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +6 | leaflight 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +3 | Koxui 2026-03-28 | 3/150 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +3 | Shern—- 2026-03-28 | 3/150 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +11 | zhubinhao 2026-03-27 | 11/550 |
|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +4 | 26要上岸 2026-03-28 | 4/200 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 5/250 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +6 | 來好運來來來 2026-03-21 | 7/350 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|