Dear all,<br>   when I test the ecutwfc for CrO2, I found that the stress at the end of the output file change largely with differenr ecutwfc values, I choose the ultrasoft pseudopotential for Cr and O, the pseudopotential as follow:<br>  Cr  Cr.pbe-sp-van.UPF<br>  O   O.pbe-rrkjus.UPF<br>   what can cause this result and who can tell me?<br>   The input file as below:<br><br><br>. ../../environment/environment_2<br>for ecut in 30 35 40 45 50  <br>do<br> rm -rf $TMP_DIR/*<br><br>cat > cro2.scf.in << EOF<br> &control<br>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='cro2',<br>    tstress = .true.<br>    tprnfor = .true.<br>    pseudo_dir = '$PSEUDO_DIR/',<br>    outdir='$TMP_DIR/'<br> /<br> &system<br>       ibrav=  6, celldm(1) = 8.3544825, celldm(3) =0.65957928, nat= 6, ntyp= 2,<br>    nspin = 2,  starting_magnetization6(1)=0.7,starting_magnetization(2)=0.0,ecutwfc =$ecut,ecutrho=400,<br>    occupations = 'smearing',smearing = 'methfessel-paxton',degauss = 0.03<br><br> /<br> &electrons<br>      mixing_beta = 0.7<br>    conv_thr =  1.0d-8<br> /<br>ATOMIC_SPECIES<br> Cr  52.000  Cr.pbe-sp-van.UPF <br> O   15.999  O.pbe-rrkjus.UPF<br>ATOMIC_POSITIONS {crystal}<br>Cr 0.000000  0.000000  0.000000<br>Cr 0.500000  0.500000  0.500000<br>O  0.301000  0.301000  0.000000<br>O  0.699000  0.699000  0.000000<br>O  0.199000  0.801000  0.500000<br>O  0.801000  0.199000  0.500000<br>K_POINTS {automatic}<br>4 4 6 0 0 0<br><br>EOF<br><br>$PARA_PREFIX /home/pwscf/pwscf_215/bin/pw.x $PARA_POSTFIX < cro2.scf.in >cro2.scf.out_$ecut<br>done<br><br><!-- footer --><br><br><br><br><br><div style="border-bottom:1px solid #999"></div><br>

        <font color="black" style="font-size:14.8px">3G 时 代 来 临 了，坚 决 对 传 统 邮 箱 说 不 ！</font>
        <br>
         <a href="http://www.126.com/" target="_blank" style="font-size:13px;line-height:160%;color:blue">新 一 代 极 速 3G 邮 箱 闪 亮 登 场 ，280 兆 网 盘 免 费 送 ！ 点 击 此 处 注 册
</a>