24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 2621  |  回复: 5
【悬赏金币】回答本帖问题,作者xiaoqiu007将赠送您 10 个金币

xiaoqiu007

新虫 (小有名气)

[求助] pwscf计算raman光谱报错

之前的自洽计算能顺利完成,但执行到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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

landaulev

铜虫 (初入文坛)

【答案】应助回帖


uuv2010(金币+1): 欢迎参加交流 2011-09-11 16:23:52
这里计算含金属体系的介电性质:
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.
2B
2楼2011-09-10 22:05:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinwei_li

木虫 (正式写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 谢谢指教 2011-09-21 13:09:54
引用回帖:
1楼: Originally posted by xiaoqiu007 at 2011-06-28 03:25:11:
之前的自洽计算能顺利完成,但执行到ph.x时报错,输入文件及报错信息如下:
输入文件:
#!/bin/sh

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

export PW_ROOT=/usr/sharesoft/ ...

第一:金属材料是无法计算拉曼强度的,如果真的想算,可以去掉occupations='smearing', smearing='methfessel-paxton', degauss=0.02近似的计算试试。
第二:似乎记得算拉曼强度是不能用超软赝势的,很多年不做了,记不清了,你可以检验下试试!
与其临渊羡鱼,不如退而结网!
3楼2011-09-21 10:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

窦政科

新虫 (初入文坛)

我也遇到了同样的问题,请问楼主你的问题解决了吗,如果解决了,方便说一下吗,谢谢,急啊
龙翔天宇
4楼2017-02-22 11:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

W.Z.Zhou

捐助贵宾 (小有名气)

楼主解决这个问题了吗?同样遇到这个问题,新人求教!
幸福也罢,苦难也罢,都来自外界,坚毅则来自内心,来自一个人自身的努力。
5楼2018-05-03 09:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

湘粤cxy

新虫 (初入文坛)

金属是不能计算拉曼的,所以相应scf.in文件中不可以有smearing, degauss这样的字样
6楼2021-12-12 23:55:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoqiu007 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 工科学硕305求调剂 +3 屋顶橙子味~ 2024-04-10 4/200 2024-04-16 23:29 by waxx333
[基金申请] 基金和生小孩 +27 Ausy 2024-04-15 27/1350 2024-04-16 23:10 by 博csc分享
[考研] 土木工程281求调剂 +4 乔兮木 2024-04-13 4/200 2024-04-16 21:40 by zjl渐行渐远
[基金申请] 颜宁:基础研究应该“立地”,即从需求端作为问题的出发点,立项背景就是面向真实需求 +3 zju2000 2024-04-16 4/200 2024-04-16 21:01 by sjtu2012
[论文投稿] The total environment格式 25+4 lyyyyds0 2024-04-13 8/400 2024-04-16 20:02 by 浅梦觉得
[找工作] 杭州国企和浙江高校如何选择? +10 restart2024 2024-04-15 10/500 2024-04-16 16:01 by 剧情需要203
[基金申请] 做什么工作,能够在不social,仅仅靠学习,就获得50w经费支持? +9 ECOblance 2024-04-10 12/600 2024-04-16 10:46 by icefireice
[教师之家] 美国教授开除中国研究生:我就不该录取你! +4 苏东坡二世 2024-04-15 4/200 2024-04-16 07:59 by lizhiweiouc
[考研] 调剂 +8 人间十三月@ 2024-04-10 10/500 2024-04-15 22:01 by zzzhonghl
[考研] 求调剂化工学硕276 +12 星星陨落 2024-04-13 12/600 2024-04-15 17:36 by 奇悠悠
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[晶体] 单晶仪在50kv 30mA ,电压电流没降,能开门上样吗? +3 15665390996 2024-04-10 3/150 2024-04-15 07:36 by gyliu
[考研] 290 生物与医药求调剂 +7 GPX4 2024-04-11 9/450 2024-04-13 09:40 by lincunhui
[考研] 280求调剂 +3 黑皮冰棒 2024-04-12 6/300 2024-04-12 16:18 by 黑皮冰棒
[考研] 267求调剂 +4 Wang可 2024-04-11 7/350 2024-04-11 22:43 by Wang可
[考研] 考研调剂 +5 夏末的晨曦12 2024-04-10 6/300 2024-04-11 18:51 by 17814699553
[考研] 生物与医药调剂! +9 夏天快过去了吗 2024-04-10 10/500 2024-04-11 17:29 by 王有田有学上
[考研] 329化学调剂 +8 瞧瞧乔儿 2024-04-10 9/450 2024-04-11 09:46 by yqWu武
[考研] 329求调剂 +9 瞧瞧乔儿 2024-04-10 9/450 2024-04-11 00:10 by xuxiaodai
[考研] 生物与医药304求调剂 +6 王有田有学上 2024-04-10 9/450 2024-04-10 15:41 by hunanzang
信息提示
请填处理意见