<div dir="ltr"><div><div><div>Hello all<br></div>Am running the phonon calculation for Ti3O4N. I have successfully run the scf calculation using the file below. However, i did not manage to obtain force constants on each discrete q-point from 17 q points.  Does anyone has a similar issue or has an idea on how to solve it?<br></div>My inputs:<br>&control<br>     calculation='scf',<br>     prefix='Ti3O4N',<br>     pseudo_dir = './'<br>     outdir='./tmp/'<br>/<br>&system    <br>    ibrav= 9, <br>    celldm(1)= 7.2871, <br>    celldm(2)= 2.4839, <br>    celldm(3)= 2.6290, <br>    nat=  16, <br>    ntyp= 3, <br>    ecutwfc = 40.0, <br>    ecutrho = 320.0,<br>    occupations = 'smearing'<br>    smearing = 'mv'<br>    degauss = 0.005<br>    nbnd = 97,<br>/<br>&ELECTRONS<br>     diagonalization = 'david',<br>     mixing_beta = 0.3,<br>     conv_thr = 1.0d-10<br>/<br>&IONS<br>     ion_dynamics='bfgs',<br>/<br>&CELL<br>     press=0.d0,<br>     wmass=0.005,<br>     cell_dynamics='bfgs',<br> /<br>ATOMIC_SPECIES<br> Ti  47.90    Ti.pbe-sp-van_ak.UPF<br> O   15.9994   O.pbe-van_ak.UPF<br> N   14.007    N.pbe-van_ak.UPF<br><br>ATOMIC_POSITIONS {alat}<br>Ti       0.000000000   1.596317050   1.505773265<br>Ti       0.000000000  -1.596317050   2.820290719<br>Ti       0.000000000   1.596317050  -0.191255810<br>Ti       0.000000000  -1.596317050  -1.505773265<br>Ti       0.000000000   1.736319689   0.657258727<br>Ti       0.000000000  -1.736319689   1.971776182<br>O        0.000000000   0.454425827   1.477033616<br>O        0.000000000  -0.454425827   2.791551070<br>O        0.000000000   0.454425827  -0.162516161<br>O        0.000000000  -0.454425827  -1.477033616<br>O        0.000000000   1.360457568   0.998830060<br>O        0.000000000  -1.360457568   2.313347515<br>O        0.000000000   1.360457568   0.315687394<br>O        0.000000000  -1.360457568  -0.998830060<br>N        0.000000000   1.815338206   1.971776182<br>N        0.000000000  -1.815338206   0.657258727<br>K_POINTS {automatic}<br>5 2 2 1 1 1<br><br><br>Phonon of Ti3O4N<br>&inputph<br>  ldisp=.true.,<br>  nq1=4, nq2=4, nq3=3<br>  amass(1)=47.90,<br>  amass(2)=15.9994,<br>  amass(3)=14.007 <br>  prefix='Ti3O4N',<br>  outdir='./tmp/'<br>  fildyn='Ti3O4N.dyn',<br>  tr2_ph=1.0d-14,<br>  epsil=.false.<br>  recover=.true.<br> /<br><br></div>Best regards,<br><div style="font-size:12.8px">Perpetua Muchiri</div><div style="font-size:12.8px">Computational Material Science Group,</div><div style="font-size:12.8px">University of Eldoret,</div><div style="font-size:12.8px">Eldoret, Kenya</div></div>