24小时热门版块排行榜    

查看: 1855  |  回复: 6

taipoboy

新虫 (初入文坛)

[交流] 求國內同胞幫幫港澳兄弟解決Quantum Espresso 的運行問題 已有2人参与

求幫忙:當我試驗運行一些以下文件,發現運行不了,主要問題顯示:

     Error in routine pw_readfile (1):
     error opening xml data file

我的文件:
#!/bin/sh

###############################################################################
##
##  HIGH VERBOSITY EXAMPLE
##
###############################################################################

# run from directory where this script is
cd `echo $0 | sed 's/\(.*\)\/.*/\1/'` # extract pathname
EXAMPLE_DIR=`pwd`

# check whether echo has the -e option
if test "`echo -e`" = "-e" ; then ECHO=echo ; else ECHO="echo -e" ; fi

$ECHO
$ECHO "$EXAMPLE_DIR : starting"
$ECHO
$ECHO "This example tests pw.x with USPS and DFT+U in the noncollinear, spin-orbit case."
$ECHO "It calculates the band structure of ferromagnetic LiZnSb."
$ECHO

# set the needed environment variables
. ../../../environment_variables

# required executables and pseudopotentials
BIN_LIST="pw.x"
PSEUDO_LIST="Li.pw-mt_fhi, Zn.pw-mt_fhi, Sb.pw-mt_fhi.UPF"

$ECHO
$ECHO "  executables directory: $BIN_DIR"
$ECHO "  pseudo directory:      $PSEUDO_DIR"
$ECHO "  temporary directory:   $TMP_DIR"
$ECHO "  checking that needed directories and files exist...\c"

# check for directories
for DIR in "$BIN_DIR" "$PSEUDO_DIR" ; do
    if test ! -d $DIR ; then
        $ECHO
        $ECHO "ERROR: $DIR not existent or not a directory"
        $ECHO "Aborting"
        exit 1
    fi
done
for DIR in "$TMP_DIR" "$EXAMPLE_DIR/results" ; do
    if test ! -d $DIR ; then
        mkdir $DIR
    fi
done
cd $EXAMPLE_DIR/results

# check for executables
for FILE in $BIN_LIST ; do
    if test ! -x $BIN_DIR/$FILE ; then
        $ECHO
        $ECHO "ERROR: $BIN_DIR/$FILE not existent or not executable"
        $ECHO "Aborting"
        exit 1
    fi
done

# check for pseudopotentials
for FILE in $PSEUDO_LIST ; do
    if test ! -r $PSEUDO_DIR/$FILE ; then
       $ECHO
       $ECHO "Downloading $FILE to $PSEUDO_DIR...\c"
            $WGET $PSEUDO_DIR/$FILE $NETWORK_PSEUDO/$FILE 2> /dev/null
    fi
    if test $? != 0; then
        $ECHO
        $ECHO "ERROR: $PSEUDO_DIR/$FILE not existent or not readable"
        $ECHO "Aborting"
        exit 1
    fi
done
$ECHO " done"

# how to run executables
PW_COMMAND="$PARA_PREFIX $BIN_DIR/pw.x $PARA_POSTFIX"
$ECHO
$ECHO "  running pw.x     as: $PW_COMMAND"
$ECHO

# self-consistent calculation for LiZnSb with fully relativistic US-PP
cat > LiZnSb.band_pbe.in << EOF
&control
    calculation = 'bands'
    prefix='LiZnSb',
    pseudo_dir = '$PSEUDO_DIR/',
    outdir='$TMP_DIR/'
/
&system   
    ibrav     = 4,
    nat       = 6,
    celldm(1) = 4.428346,
    celldm(3) = 1.616938,
    ntyp      = 3,
    ecutwfc   = 18.D0,
    ecutrho   = 180.0,
    noncolin=.true.
    lspinorb=.true.
    starting_magnetization(1)=0.5,
    angle1(1)=90.0
    angle2(1)=0.0
    occupations='smearing',
    smearing='mp',
    degauss=0.04,
    lda_plus_u=.true.
    lda_plus_u_kind=1
    Hubbard_U(1)=2.2
    Hubbard_J(1,1)=1.75
    Hubbard_J(2,1)=0.0
/
&electrons
    conv_thr =  1.0d-12,
/
ATOMIC_SPECIES
Li    6.94100  Li.pbe-mt_fhi.UPF
Zn   65.38000  Zn.pbe-mt_fhi.UPF
Sb  121.76000  Sb.pbe-hgh.UPF
ATOMIC_POSITIONS {bohr}
Li      0.999990977    0.999983308    0.175222936   
Li      0.000007824    0.000015649    0.675222656   
Zn      0.333351242    0.666702484    0.394028851   
Zn      0.666647559    0.333296472    0.894028572   
Sb      0.333362979    0.666725957    0.011748006   
Sb      0.666635822    0.333273000    0.511747727  
K_POINTS tpiba_b
3
1.0 0.0 0.0 40
0.0 0.0 0.0 40
0.0 0.0 1.0 1
EOF
$ECHO "  running the band calculation for LiZnSb with USPP and DFT+U and spin-orbit...\c"
$PW_COMMAND < LiZnSb.band_pbe.in > LiZnSb.band_pbe.out
check_failure $?
$ECHO " done"

# clean TMP_DIR
$ECHO "  cleaning $TMP_DIR...\c"
#rm -rf $TMP_DIR/LiZnSb*
$ECHO " done"


$ECHO
$ECHO "$EXAMPLE_DIR: done"

Sample Text
Program PWSCF v.5.2.0 starts on 11Jul2015 at  5:50:32

     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/quote

     Serial version
     Waiting for input...
     Reading input from standard input

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3

     Atomic positions and unit cell read from directory:
     /home/kayi/Desktop/Program/espresso-5.2.0/tempdir/LiZnSb.save/

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine pw_readfile (1):
     error opening xml data file
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...


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

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

条件错误?
不要说话
2楼2015-07-11 21:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这个脚本需要预先计算好的scf输出文件才能计算,单单只有这个脚本是不能计算的。
3楼2015-07-11 21:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taipoboy

新虫 (初入文坛)

那應該怎麼樣做~~是不是先把SCF算好,才可以算下一步Bands
4楼2015-07-12 09:53:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taipoboy

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by KalaShayminS at 2015-07-11 21:55:47
这个脚本需要预先计算好的scf输出文件才能计算,单单只有这个脚本是不能计算的。

那應該怎麼樣做~~是不是先把SCF算好,才可以算下一步Bands
5楼2015-07-12 10:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

引用回帖:
5楼: Originally posted by taipoboy at 2015-07-12 10:34:07
那應該怎麼樣做~~是不是先把SCF算好,才可以算下一步Bands...

是的
6楼2015-07-12 19:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taipoboy

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by KalaShayminS at 2015-07-12 19:06:21
是的...

謝謝~~真的解決了。
7楼2015-07-15 16:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 taipoboy 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 6/300 2026-08-28 01:46 by alongwaytogo
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见