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

xiaoqiu007

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 356.5
散金: 7
帖子: 100
在線: 20.8小時(shí)
蟲號(hào): 707167
注冊(cè): 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ù)此樓

» 猜你喜歡

0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有4人回復(fù)
0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級(jí)證書40+數(shù)一英一287 已經(jīng)有5人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有6人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西南交大，求調(diào)劑已經(jīng)有3人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)
焦慮已經(jīng)有11人回復(fù)
331求調(diào)劑（0703有機(jī)化學(xué) 已經(jīng)有7人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有5人回復(fù)
297求調(diào)劑已經(jīng)有5人回復(fù)

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

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í)變化嗎？已經(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ù)
【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂已經(jīng)有21人回復(fù)

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

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

yinwei_li

木蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 2079.1
散金: 1532
紅花: 10
帖子: 896
在線: 194.2小時(shí)
蟲號(hào): 443243
注冊(cè): 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)度是不能用超軟贗勢(shì)的，很多年不做了，記不清了，你可以檢驗(yàn)下試試！

贊一下(1人)

回復(fù)此樓

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

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

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

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

landaulev

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 39
帖子: 6
在線: 7.1小時(shí)
蟲號(hào): 1259266
注冊(cè): 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的回帖

竇政科

新蟲 (初入文壇)

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 3636656
注冊(cè): 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的回帖

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

[考研] 297求調(diào)劑 +5	戲精丹丹丹 2026-03-17	5/250	2026-03-18 12:46 by Linda Hu
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 收復(fù)試調(diào)劑生 +3	雨后秋荷 2026-03-18	3/150	2026-03-18 11:36 by xdu電子
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[基金申請(qǐng)] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz