<div dir="auto">Dear users,<div dir="auto">I'm trying to perform a convergence test for CaSiO3 tetragonal phase. Following is my input :</div><div dir="auto">#!/bin/sh</div><div dir="auto"><br></div><div dir="auto">for aa in 5.6131</div><div dir="auto">do</div><div dir="auto">for cc in 7.181 7.182 7.183 7.184 7.185</div><div dir="auto">do</div><div dir="auto">for k in 8</div><div dir="auto">do</div><div dir="auto">for ecut in 30</div><div dir="auto">do</div><div dir="auto">cat > ${aa}_${cc}_${k}_${ecut}.in << EOF</div><div dir="auto"> &control</div><div dir="auto">    calculation = 'scf',</div><div dir="auto">    prefix = '${aa}_${cc}_${k}_${ecut}'</div><div dir="auto">    tprnfor=.true.</div><div dir="auto">    tstress=.true.</div><div dir="auto">    outdir = '/home/user/PNV/CaSiO3_tetragonal/'</div><div dir="auto">    pseudo_dir = '/home/user/PNV/CaSiO3_tetragonal/pseudo/'</div><div dir="auto"> /</div><div dir="auto"> &system</div><div dir="auto">    ibrav =  7,</div><div dir="auto">    A=$aa</div><div dir="auto">    C=$cc</div><div dir="auto">    nat =  10,</div><div dir="auto">    ntyp = 3,</div><div dir="auto">    ecutwfc = $ecut,</div><div dir="auto">/    </div><div dir="auto">&electrons</div><div dir="auto">    diagonalization='david'</div><div dir="auto">    mixing_mode = 'plain'</div><div dir="auto">    mixing_beta = 0.7</div><div dir="auto">    conv_thr =  1.0d-12</div><div dir="auto">/</div><div dir="auto">ATOMIC_SPECIES</div><div dir="auto"> </div><div dir="auto">Ca 40.078  Ca.pz-spn-kjpaw_psl.1.0.0.UPF</div><div dir="auto">Si 28.0855 Si.pz-n-kjpaw_psl.0.1.UPF</div><div dir="auto"> O 15.999   O.pz-n-kjpaw_psl.0.1.UPF</div><div dir="auto"><br></div><div dir="auto">ATOMIC_POSITIONS {crystal}</div><div dir="auto">Ca 0.250 0.750 0.500</div><div dir="auto">Ca 0.750 0.250 0.500</div><div dir="auto">Si 0.000 0.000 0.000</div><div dir="auto">Si 0.500 0.500 0.000</div><div dir="auto">O  0.222 0.722 0.9441</div><div dir="auto">O  0.250 0.250 0.000</div><div dir="auto">O  0.278 0.222 0.500</div><div dir="auto">O  0.722 0.778 0.500</div><div dir="auto">O  0.750 0.750 0.000</div><div dir="auto">O  0.778 0.278 0.0559</div><div dir="auto"><br></div><div dir="auto">K_POINTS (automatic)</div><div dir="auto">  $k $k $k 1 1 1</div><div dir="auto"><br></div><div dir="auto">EOF</div><div dir="auto">mpirun -np 40 -machinefile x2 /apps/codes/qe/6.5/bin/pw.x < ${aa}_${cc}_${k}_${ecut}.in >  ${aa}_${cc}_${k}_${ecut}scf.out</div><div dir="auto">done</div><div dir="auto">done</div><div dir="auto">done</div><div dir="auto">done</div><div dir="auto"><br></div><div dir="auto">To perform the test, I fixed c=7.182 angstrom and varied 'a'. The value of 'a' for which I attained zero pressure was 5.6131 angstrom, which is quite different from other theoretical finding which is a=4.9812 angstrom (<a href="https://sci-hub.se/https://doi.org/10.1107/s0108768106035762">https://sci-hub.se/https://doi.org/10.1107/s0108768106035762</a>). The error in 'a' is about 10-12%, which is quite large.</div><div dir="auto"><br></div><div dir="auto">Can I know why I am getting this much discrepancy in lattice constant? And how can I get the correct value of 'a'?</div></div>