xiaoqiu007

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時
蟲號: 707167
注冊: 2009-02-23
專業(yè): 金屬功能材料

[求助] pwscf計算raman光譜報錯

之前的自洽計算能順利完成，但執(zhí)行到ph.x時報錯,輸入文件及報錯信息如下：
輸入文件:
#!/bin/sh

####################################################################

export PW_ROOT=/usr/sharesoft/apps/pwscf/espresso-4.3/bin
export PSEUDO_DIR=../pseudo

export PARA_PREFIX='mpirun  -np 18'
# or export PARA_PREFIX='mpirun' ,export  PARA_POSTFIX= -np 3
export TMP_DIR=../tmp/tmp1

rm -rf $TMP_DIR/*

export name='d194'

export a='0GPa'

cat > $name.scf.in_$a << EOF

&control
calculation = 'scf'
restart_mode='from_scratch',
prefix='$name',
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
tstress=.t.,
tprnfor=.t.
/
&system
ibrav=4,
celldm(1)=6.357227657,
celldm(3)=1.814393151,
  nat= 8, ntyp= 2
ecutwfc =80,
occupations='smearing', smearing='methfessel-paxton', degauss=0.02
/
&electrons
mixing_beta = 0.7
conv_thr =  1.0d-8
/
ATOMIC_SPECIES
H    1.008    h.cpi.UPF
Sc 44.9559 Sc.cpi.UPF
ATOMIC_POSITIONS {crystal}
H 0.3333330000000028  0.6666669999999962  0.5923567327221010
H 0.6666669999999965  0.3333330000000026  0.0923567327221009
H 0.6666669999999972  0.3333330000000038  0.4076432672778991
H 0.3333330000000037  0.6666669999999973  0.9076432672778990
H 0.0000000000000000  0.0000000000000000  0.7500000000000000
H 0.0000000000000000  0.0000000000000000  0.2500000000000000
Sc  0.6666669999999968  0.3333330000000032  0.7500000000000000
Sc  0.3333330000000032  0.6666669999999968  0.2500000000000000
K_POINTS {automatic}
18 18 18  0 0 0
EOF
$PARA_PREFIX $PW_ROOT/pw.x  < $name.scf.in_$a >$name.scf.out_$a

cat > $name.phG.in << EOF
phonons of $name at Gamma
&inputph
tr2_ph=1.0d-12,
prefix='$name',
epsil=.true.,
lraman=.true.,
amass(1)=1.008,
amass(2)=44.955910,
outdir='$TMP_DIR/',
fildyn='$name.dynG',
/
0.0 0.0 0.0
EOF
$PARA_PREFIX $PW_ROOT/ph.x  < $name.phG.in >$name.phG.out

報錯信息：

Program PHONON v.4.3    starts on 28Jun2011 at  2:51:23

   This program is part of the open-source Quantum ESPRESSO suite
   for quantum simulation of materials; please cite
      "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
      URL http://www.quantum-espresso.org",
   in publications or presentations arising from this work. More details at
   http://www.quantum-espresso.org/ ... ng_Quantum-ESPRESSO

   Parallel version (MPI), running on 18 processors
   R & G space division:  proc/pool = 18

   Ultrasoft (Vanderbilt) Pseudopotentials

Info: using nr1, nr2, nr3 values from input

Info: using nr1s, nr2s, nr3s values from input
   file h.cpi.UPF: wavefunction(s)  1 3 renormalized
   file Sc.cpi.UPF: wavefunction(s)  3 renormalized

Stick Mesh
----------
nst = 889,  nstw = 253, nsts = 889
            n.st n.stw n.sts n.g n.gw n.gs
min       49    14    49 2165    338 2165
max       50    15    50 2170    347 2170
      889    253    889 39005 6131 39005

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from phq_readin : error #       1
   no elec. field with metals
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from phq_readin : error #       1
   no elec. field with metals
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

回復(fù)此樓

» 猜你喜歡

312求調(diào)劑已經(jīng)有7人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有16人回復(fù)
280求調(diào)劑已經(jīng)有3人回復(fù)
266求調(diào)劑已經(jīng)有7人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有5人回復(fù)
環(huán)境工程調(diào)劑已經(jīng)有8人回復(fù)
302求調(diào)劑已經(jīng)有9人回復(fù)
297求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有4人回復(fù)
焦慮已經(jīng)有10人回復(fù)

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

PWscf計算激發(fā)態(tài)能量已經(jīng)有13人回復(fù)
PWscf計算homo lumo或者valence conduction band，問題已解決已經(jīng)有14人回復(fù)
PWSCF計算PDOS時，選哪列數(shù)據(jù)？已經(jīng)有15人回復(fù)
pwscf算費米面為什么計算很慢已經(jīng)有5人回復(fù)
PWSCF可以計算半導(dǎo)體材料費米能級變化嗎？已經(jīng)有19人回復(fù)
【求助】pwscf計算得到的能帶圖很亂已經(jīng)有16人回復(fù)
【求助】pwscf計算中的warning要緊嗎？如何消除呢？已經(jīng)有9人回復(fù)
【其他】Pwscf不能并行計算已經(jīng)有4人回復(fù)
【求助】關(guān)于PWSCF介電常數(shù)和聲子譜的計算已經(jīng)有7人回復(fù)
【求助】PWSCF計算介電極化時為何出現(xiàn)如下提示？已經(jīng)有19人回復(fù)
【求助】用PWSCF中的LDA+U做反鐵磁絕緣體系的弛豫計算已經(jīng)有23人回復(fù)
【求助】請教PWSCF晶體建模和自洽計算，不收斂已經(jīng)有21人回復(fù)

1樓 2011-06-28 03:25:11

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

landaulev

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 39
帖子: 6
在線: 7.1小時
蟲號: 1259266
注冊: 2011-04-08
性別: GG
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

★
uuv2010(金幣+1): 歡迎參加交流 2011-09-11 16:23:52

這里計算含金屬體系的介電性質(zhì)：
epsil=.true.,
lraman=.true.
說明文件里有解釋：
Do not set epsil to .true. if you have a metallic system or q/=0: the code will complain and stop.

贊一下(1人)

回復(fù)此樓

2樓2011-09-10 22:05:17

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

yinwei_li

木蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2079.1
散金: 1532
紅花: 10
帖子: 896
在線: 194.2小時
蟲號: 443243
注冊: 2007-10-27
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★ ★
zzy870720z(金幣+2): 謝謝指教 2011-09-21 13:09:54

引用回帖:

1樓: Originally posted by xiaoqiu007 at 2011-06-28 03:25:11:
之前的自洽計算能順利完成，但執(zhí)行到ph.x時報錯,輸入文件及報錯信息如下：
輸入文件:
#!/bin/sh

####################################################################

export PW_ROOT=/usr/sharesoft/ ...

第一：金屬材料是無法計算拉曼強度的，如果真的想算，可以去掉occupations='smearing', smearing='methfessel-paxton', degauss=0.02近似的計算試試。
第二：似乎記得算拉曼強度是不能用超軟贗勢的，很多年不做了，記不清了，你可以檢驗下試試！

贊一下(1人)

回復(fù)此樓

與其臨淵羨魚，不如退而結(jié)網(wǎng)！

3樓2011-09-21 10:59:43

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

竇政科

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 48
紅花: 1
帖子: 22
在線: 148.1小時
蟲號: 5065517
注冊: 2016-09-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

我也遇到了同樣的問題，請問樓主你的問題解決了嗎，如果解決了，方便說一下嗎，謝謝，急啊

贊一下

回復(fù)此樓

龍翔天宇

4樓2017-02-22 11:55:53

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

W.Z.Zhou

捐助貴賓 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 527.8
帖子: 71
在線: 58.8小時
蟲號: 3636656
注冊: 2015-01-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

樓主解決這個問題了嗎？同樣遇到這個問題，新人求教！

贊一下

回復(fù)此樓

幸福也罷，苦難也罷，都來自外界，堅毅則來自內(nèi)心，來自一個人自身的努力。

5樓2018-05-03 09:34:20

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

湘粵cxy

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 118.3
帖子: 36
在線: 13.5小時
蟲號: 4800129
注冊: 2016-06-27
性別: MM
專業(yè): 基礎(chǔ)物理學(xué)

金屬是不能計算拉曼的，所以相應(yīng)scf.in文件中不可以有smearing, degauss這樣的字樣

贊一下

回復(fù)此樓

6樓2021-12-12 23:55:23

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

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

[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321