<div dir="ltr"><div>Dear users,</div><div>I'm trying to do scf calculation of teragonal <span class="gmail-value">I4/mcm</span> CaSiO3. Following is my input. The parameters A and C as well as atomic positions are taken from materials project <a href="https://materialsproject.org/materials/mp-3387/">https://materialsproject.org/materials/mp-3387/</a> . However, energy doesn't converge after 100 iterations. <br></div><div>Â </div><div>&control<br>Â Â calculation = 'scf',<br>Â Â prefix = '5.115_5.115_2_40'<br>Â Â tprnfor=.true.<br>Â Â tstress=.true.<br>Â Â outdir = '/home/nkb/QE/qe6.5/PW/examples/example01/CaSiO3_tetragonal/'<br>Â Â pseudo_dir = '/home/nkb/QE/qe6.5/pseudo/'<br>Â /<br>Â &system<br>Â Â ibrav = Â 7, <br>Â Â A=5.115<br>Â Â C=5.115 <br>Â Â nat = Â 10, <br>Â Â ntyp = 3,<br>Â Â ecutwfc = 40,<br>/ Â Â <br>&electrons<br>Â Â diagonalization='david'<br>Â Â mixing_mode = 'plain'<br>Â Â mixing_beta = 0.7<br>Â Â conv_thr = Â 1.0d-12<br>Â /<br><br>ATOMIC_SPECIES<br>Â <br>Ca 40.078 Â Ca.pbe-spn-kjpaw_psl.1.0.0.UPF<br>Si 28.0855 Si.pbe-n-kjpaw_psl.1.0.0.UPF<br>Â O 15.999 Â O.pbe-n-kjpaw_psl.1.0.0.UPF<br><br>ATOMIC_POSITIONS {crystal}<br>Ca 0.250 0.750 0.500<br>Ca 0.750 0.250 0.500<br>Si 0.000 0.000 0.000<br>Si 0.500 0.500 0.000<br>O Â 0.222 0.722 0.9441<br>O Â 0.250 0.250 0.000<br>O Â 0.278 0.222 0.500<br>O Â 0.722 0.778 0.500<br>O Â 0.750 0.750 0.000<br>O Â 0.778 0.278 0.0559<br><br>K_POINTS (automatic)<br>Â 2 2 2 1 1 1</div><div><br></div><div>Can anyone let me know what's wrong with my input. Any kind of help is appreciated.</div><div>Thanks.<br></div></div>