24小时热门版块排行榜    

查看: 1853  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 6/300 2026-08-28 01:46 by alongwaytogo
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-28 01:23 by xxniao123
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 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
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
信息提示
请填处理意见