版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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計算振動頻率的問題已有5人參與

我想計算一個吸附在硅表面上的水分子的振動頻率。所有的原子都計算的話，我試過，太慢了，由于我只想計算水分子的振動頻率，就把表面硅原子全部固定。方法是在ph.x的輸入文件里加上這樣一項：nat_todo=3，并給出水分子三個原子的原子系數(shù)。能跑起來，但是算完后，在out文件里卻沒有任何振動頻率的結(jié)果。反復(fù)試了幾次都是如此，還以為是bug,換用最新的QE4.2版本算也是同樣的結(jié)果。問題出在哪里？請高人指點迷津...

回復(fù)此樓

» 猜你喜歡

332求調(diào)劑已經(jīng)有5人回復(fù)
085701求調(diào)劑初試286分已經(jīng)有3人回復(fù)
材料工程294求調(diào)劑已經(jīng)有4人回復(fù)
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑已經(jīng)有10人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有5人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有10人回復(fù)
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過已經(jīng)有11人回復(fù)
071000生物學(xué)求調(diào)劑，初試成績343 已經(jīng)有7人回復(fù)
0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有8人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

各位大俠，我可以優(yōu)化結(jié)構(gòu)，但只要一計算頻率就出現(xiàn)下面的問題？已經(jīng)有4人回復(fù)
【交流】振動模式與紅外光譜、拉曼光譜的原理與計算問題的討論已經(jīng)有14人回復(fù)
振動模式，聲子支與對稱性分類--期待熱心的蟲友已經(jīng)有11人回復(fù)
【求助】請教一下關(guān)于Gaussian計算振動頻率的問題已經(jīng)有7人回復(fù)
【求助】gaussian頻率計算問題已經(jīng)有9人回復(fù)
【求助】MS5.0中Dmol3計算頻率問題已經(jīng)有21人回復(fù)
【求助】如何分析ph.x計算振動頻率的結(jié)果已經(jīng)有4人回復(fù)
【求助】用MS中的castep計算聲子振動模頻率，具體怎么算已經(jīng)有10人回復(fù)

1樓 2010-05-31 11:22:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時
蟲號: 1007127
注冊: 2010-04-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
huazhorg(金幣+1):謝謝參與
gavinliu7390(金幣+1):謝謝交流！ 2010-05-31 14:15:52

請貼出你的ph.x的輸入文件。

贊一下(2人)

回復(fù)此樓

2樓2010-05-31 12:52:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

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ù)。
（補充：asr=.true., zue=.true., epsil=.true., trans=.true. 這幾項可以不要，但結(jié)果也是一樣的：沒有聲子。）

[ Last edited by huazhorg on 2010-5-31 at 14:40 ]

贊一下

回復(fù)此樓

3樓2010-05-31 14:38:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強帖: 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

為計算出頻率的原因，可能是你的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 ]

贊一下(1人)

回復(fù)此樓

4樓2010-05-31 16:57:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

等待的心

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 819.7
散金: 620
沙發(fā): 1
帖子: 372
在線: 91.8小時
蟲號: 947336
注冊: 2010-01-23
性別: GG
專業(yè): 物理無機化學(xué)

★
huazhorg(金幣+1):謝謝參與

我想請教下樓主，固定原子能省大概多少時間？
謝謝！
我現(xiàn)在也很困惑，我計算頻率也很花費時間

贊一下(1人)

回復(fù)此樓

waitingisbeautiful!

5樓2010-05-31 18:40:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
qasd(金幣+1):謝謝 2010-05-31 20:42:31

引用回帖:

Originally posted by 等待的心 at 2010-05-31 18:40:17:
我想請教下樓主，固定原子能省大概多少時間？
謝謝！
我現(xiàn)在也很困惑，我計算頻率也很花費時間

如果你的體系很多原子，只需要關(guān)注局部的振動情況（比如吸附體系），可以把其他部分都固定，這樣可以節(jié)省絕大部分時間。如果不固定，N個原子的體系一般需要3N步循環(huán)，如果如果只想讓其中的X個振動，其他的固定，就只需要3X步循環(huán)就夠了。因此固定原子是很有實用意義的。但是現(xiàn)在的問題是，固定原子的問題還沒有解決。

[ Last edited by huazhorg on 2010-5-31 at 18:56 ]

贊一下(1人)

回復(fù)此樓

6樓2010-05-31 18:44:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by valenhou001 at 2010-05-31 16:57:58:
為計算出頻率的原因，可能是你的tr2_ph的收斂精度設(shè)置太高，而未收斂，適當(dāng)降低它。

在做ph.x之前進行的pw.x中，使用多個點。比如下面的：
K_POINTS {automatic}
6 6 1 1 1 1

將ph.x的輸入文件中下面的 ...

我試過了，降低了精度，按照你說的辦法調(diào)整了參數(shù)，也還不行。另外，做SCF運算時，K點只能取在GAMMA點，其他點會出錯的。
嘗試了很多次，現(xiàn)在這個問題好棘手啊……

贊一下

回復(fù)此樓

7樓2010-05-31 18:48:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

如果有人能幫我在pwscf forum 上問問也行，我上不了該論壇。（我注冊了，但是發(fā)現(xiàn)上去很麻煩，不知道怎么上去發(fā)貼）。Thanks all !

贊一下

回復(fù)此樓

8樓2010-05-31 19:34:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強帖: 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 22:09:04
gavinliu7390(金幣+2):謝謝交流！ 2010-06-01 10:48:04

引用回帖:

Originally posted by huazhorg at 2010-05-31 19:34:52:
如果有人能幫我在pwscf forum 上問問也行，我上不了該論壇。（我注冊了，但是發(fā)現(xiàn)上去很麻煩，不知道怎么上去發(fā)貼）。Thanks all !

可以的話，把你pw.x和ph.x的輸入文件發(fā)到valenhou@gmail.com.我再幫你看看。

贊一下(1人)

回復(fù)此樓

9樓2010-05-31 21:33:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時
蟲號: 629899
注冊: 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by valenhou001 at 2010-05-31 21:33:43:

可以的話，把你pw.x和ph.x的輸入文件發(fā)到valenhou@gmail.com.我再幫你看看。

pw.x和ph.x的輸入輸出文件我已經(jīng)發(fā)送到你的郵箱你了，謝謝。

贊一下

回復(fù)此樓

10樓2010-05-31 22:24:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 huazhorg 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來啊
[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 291求調(diào)劑 +6	HanBeiNingZC 2026-03-24	6/300	2026-03-28 07:55 by baoball
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +9	vv迷 2026-03-22	9/450	2026-03-27 15:59 by 不吃魚的貓
[考研] 274求調(diào)劑 +17	顧九笙要謙虛 2026-03-24	23/1150	2026-03-27 15:16 by caszguilin
[考研] 279 分求調(diào)劑 +4	睡個好覺_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風(fēng)里
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 求調(diào)劑 +8	Auroracx 2026-03-22	8/400	2026-03-26 19:55 by 不吃魚的貓
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	10/500	2026-03-26 13:11 by 公瑾逍遙
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 機械學(xué)碩總分317求調(diào)劑�。。。� +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 284求調(diào)劑 +15	Zhao anqi 2026-03-22	15/750	2026-03-25 12:51 by wht0531
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 求調(diào)劑 +6	研研，接電話 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc

24小時熱門版塊排行榜

[交流] 【求助】ph.x計算振動頻率的問題 已有5人參與

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

[交流] 【求助】ph.x計算振動頻率的問題已有5人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: