當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

xiaoqiu007

新蟲 (小有名氣)

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

[求助] pwscf計(jì)算raman光譜報(bào)錯(cuò)

之前的自洽計(jì)算能順利完成，但執(zhí)行到ph.x時(shí)報(bào)錯(cuò),輸入文件及報(bào)錯(cuò)信息如下：
輸入文件:
#!/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

報(bào)錯(cuò)信息：

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ù)此樓

» 猜你喜歡

265求調(diào)劑已經(jīng)有8人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有21人回復(fù)
295材料求調(diào)劑，一志愿武漢理工085601專碩已經(jīng)有4人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有9人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？已經(jīng)有7人回復(fù)
材料學(xué)碩297已過四六級求調(diào)劑推薦已經(jīng)有6人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有7人回復(fù)
本人考085602 化學(xué)工程專碩已經(jīng)有20人回復(fù)

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

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

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

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

W.Z.Zhou

捐助貴賓 (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

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

查看全部 6 個(gè)回答

landaulev

銅蟲 (初入文壇)

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

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

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

這里計(jì)算含金屬體系的介電性質(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小時(shí)
蟲號: 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:
之前的自洽計(jì)算能順利完成，但執(zhí)行到ph.x時(shí)報(bào)錯(cuò),輸入文件及報(bào)錯(cuò)信息如下：
輸入文件:
#!/bin/sh

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

export PW_ROOT=/usr/sharesoft/ ...

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

贊一下(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小時(shí)
蟲號: 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的回帖

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

[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+3	babangida 2026-03-19	7/350	2026-03-20 14:13 by hejicker
[考研] 279分求調(diào)劑一志愿211 +9	chaojifeixia 2026-03-19	10/500	2026-03-20 12:29 by lature00
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	14/700	2026-03-20 09:30 by kkcoco25
[考研] 304求調(diào)劑 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike