<div dir="ltr"><div><div>When I calculate the elastic constants for supercells of my system using thermo_pw.x. I found the compatibility problem of ibrav. My system is tetragonal <br>Please how can I solve the issue and continue my calculations? <br></div><div>Thank you in advance for any assistance <br></div><div>The input / out files of calculations are shown as follows<br></div><div><br></div><div>&INPUT_THERMO</div>  what='mur_lc_elastic_constants',<br>/<br>&CONTROL<br>  Â  Â  Â  Â  Â  Â  Â  Â calculation = 'scf' ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â pseudo_dir = './',<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â outdir='temp/'<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  prefix = 'psto50' ,<br> /<br> &SYSTEM<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ibrav = 6,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â celldm(1) = 14.74,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â celldm(3) = 1.06,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â nat = 20,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ntyp = 4,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ecutwfc = 40.0 ,<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  nbnd = 62,<br> /<br> &ELECTRONS<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  conv_thr = 1.0d-8 ,<br>  Â  Â  Â  Â  Â  Â  Â  Â mixing_mode = 'plain' ,<br>  Â  Â  Â  Â  Â  Â  Â  Â mixing_beta = 0.7 ,<br>  Â  Â  Â  Â  Â  Â diagonalization = 'davidson' ,<br> /<br>ATOMIC_SPECIES<br>  Â Pb Â 207.20000 Â Pb.pbe-mt_fhi.UPF<br>  Â Sn Â 118.70000 Â Sn.pbe-mt_fhi.UPF<br>  Â Ti Â  47.86700 Â Ti.pbe-mt_fhi.UPF<br>  Â  O Â  15.99940 Â O.pbe-mt_fhi.UPF<br>ATOMIC_POSITIONS (angstrom)<br>Sn Â  Â  Â  0.000000000 Â  0.000000000 Â  0.000000000<br>Ti Â  Â  Â  1.995794146 Â  2.008568413 Â  2.009579978<br>O Â  Â  Â  Â 1.984411263 Â  2.008567974 Â -0.000000000<br>O Â  Â  Â  Â 0.000000000 Â  2.008568082 Â  2.009579978<br>O Â  Â  Â  Â 1.985281957 Â  0.000000000 Â  2.009579978<br>Sn Â  Â  Â  0.000000000 Â  4.017139736 Â  0.000000000<br>Ti Â  Â  Â  1.995794146 Â  6.025711058 Â  2.009579978<br>O Â  Â  Â  Â 1.984411263 Â  6.025711497 Â -0.000000000<br>O Â  Â  Â  Â 0.000000000 Â  6.025711389 Â  2.009579978<br>O Â  Â  Â  Â 1.985283246 Â  4.017139736 Â  2.009579978<br>Pb Â  Â  Â  4.016267350 Â  0.000000000 Â  0.000000000<br>Ti Â  Â  Â  6.036740554 Â  2.008568413 Â  2.009579978<br>O Â  Â  Â  Â 6.048123437 Â  2.008567974 Â -0.000000000<br>O Â  Â  Â  Â 4.016267350 Â  2.008566335 Â  2.009579978<br>O Â  Â  Â  Â 6.047252743 Â  0.000000000 Â  2.009579978<br>Pb Â  Â  Â  4.016267350 Â  4.017139736 Â  0.000000000<br>Ti Â  Â  Â  6.036740554 Â  6.025711058 Â  2.009579978<br>O Â  Â  Â  Â 6.048123437 Â  6.025711497 Â -0.000000000<br>O Â  Â  Â  Â 4.016267350 Â  6.025713136 Â  2.009579978<br>O Â  Â  Â  Â 6.047251454 Â  4.017139736 Â  2.009579978<br>K_POINTS automatic<br> 4  4  4 Â 1  1  1</div><div><br></div><div><b>Output</b></div><div><br></div><div> Reading input from _temporary_1<br>  Â  Â  Â  Â  Â  Â  Â file Pb.pbe-mt_fhi.UPF: wavefunction(s) Â 5f renormalized<br>  Â  Â  Â  Â  Â  Â  Â file Sn.pbe-mt_fhi.UPF: wavefunction(s) Â 5d 4f renormalized<br>  Â  Â  Â  Â  Â  Â  Â file Ti.pbe-mt_fhi.UPF: wavefunction(s) Â 4f renormalized<br>  Â  Â  Â  Â  Â  Â  Â file O.pbe-mt_fhi.UPF: wavefunction(s) Â 4f renormalized<br><br>  Â Info: using nr1, nr2, nr3 values from input<br><br>  Â Info: using nr1, nr2, nr3 values from input<br><br>  Â  Â C_1 (1) is incompatible with the tetragonal Bravais lattice<br>  Â  Â It is compatible with the<br>  Â  Â triclinic Bravais lattice; ibrav= Â  14<br>  Â  Â You might want to change the Bravais lattice or to<br>  Â  Â understand why the symmetries are wrong before continuing<br>  Â  Â The point group or the Laue class are not used to reduce the number of<br>  Â  Â computed tensor components</div><div>Info: using nr1, nr2, nr3 values from input<br>  Â  Â Computing the frozen ions elastic constants at the minimum volume<br>  Â  Â FFT mesh: ( Â  60, Â  60, Â  64 )<br><br>  Â  Â Bravais lattice:<br><br>  Â  Â ibrav= Â 6: tetragonal<br>  Â  Â Cell parameters:<br><br>  Â  Â alat= Â 14.740000 a.u., c/a= Â  1.060000<br><br><br>  Â  Â Starting primitive lattice vectors:<br>  Â  Â crystal axes: (cart. coord. in units of alat)<br><br>  Â  Â  Â  Â  Â  Â  Â a(1) = ( Â  1.000000 Â  0.000000 Â  0.000000 )<br>  Â  Â  Â  Â  Â  Â  Â a(2) = ( Â  0.000000 Â  1.000000 Â  0.000000 )<br>  Â  Â  Â  Â  Â  Â  Â a(3) = ( Â  0.000000 Â  0.000000 Â  1.060000 )<br><br>  Â  Â Starting reciprocal lattice vectors:<br>  Â  Â reciprocal axes: (cart. coord. in units 2 pi/alat)<br><br>  Â  Â  Â  Â  Â  Â  Â b(1) = ( Â 1.000000 Â 0.000000 Â 0.000000 )<br>  Â  Â  Â  Â  Â  Â  Â b(2) = ( Â 0.000000 Â 1.000000 Â 0.000000 )<br>  Â  Â  Â  Â  Â  Â  Â b(3) = ( Â 0.000000 Â 0.000000 Â 0.943396 )<br>  Â Â  ......</div><div><br></div><div><br></div><div>Shiferaw G.</div><div>PhD student</div><div>Department of physics<br></div></div>