<div dir="ltr">Dear Ari, thanks for you reply and kind attention. By checking my input files, I managed to get the phonon code running under the <span style="font-size:13px"> 'tr2_ph' </span> threshold value of 1.0D-12. So, the code successfully ends the electric fields calculation, dielectric constant and effective charges (dF/dE) stages. Nevertheless, I still have an issue: the second order response computation does not converge; for example, after 11 iterations, calculations the |ddv_scf|^2 term  gives as high as 1E+06.<div><div><br></div><div>Here is my <a href="http://pw.in" target="_blank">pw.in</a> and <a href="http://ph.in" target="_blank">ph.in</a> files, so you or any in the pwscf forum can give me some light of advice...</div><div><br></div><div>-------------</div><div><a href="http://pw.in" target="_blank">pw.in</a>:</div><div><br></div><div><div> &CONTROL</div><div>                       title = 'MAPbI3-no-H-NC-PZ' ,</div><div>                 calculation = 'scf' ,</div><div>                      outdir = '/home/quantum/Descargas/mapbi-real/' ,</div><div>                  pseudo_dir = '/home/quantum/Descargas/espresso-5.2.0-GPU/pseudo_PZ_NC/' ,</div><div>                      prefix = 'MAPbI3-real' ,</div><div> /</div><div> &SYSTEM</div><div>                       ibrav = 6,</div><div>                   celldm(1) = 16.722185273,</div><div>                   celldm(3) = 1.428636004,</div><div>                         nat = 24,</div><div>                        ntyp = 4,</div><div>                     ecutwfc = 70 ,</div><div>                     ecutrho = 280 ,</div><div>                    noncolin = .false. ,</div><div> /</div><div> &ELECTRONS</div><div>            electron_maxstep = 1000,</div><div>                 startingpot = 'atomic' ,</div><div>                 startingwfc = 'atomic+random' ,</div><div>                 mixing_beta = 0.7 ,</div><div> /</div><div>ATOMIC_SPECIES</div><div>   Pb  207.20000  Pb.rel-pz-n-nc.UPF </div><div>    I  126.90400  I.rel-pz-n-nc.UPF </div><div>    C   12.01070  C.rel-pz-nc.UPF </div><div>    N   14.00670  N.rel-pz-nc.UPF </div><div>ATOMIC_POSITIONS alat </div><div>    C      0.500000000    0.000000000    0.352000000    </div><div>    C      0.000000000    0.500000000    0.852000000    </div><div>    C      0.000000000    0.500000000    0.352000000    </div><div>    C      0.500000000    0.000000000    0.852000000    </div><div>    I      0.000000000    0.000000000    0.247200000    </div><div>    I      0.500000000    0.500000000    0.747200000    </div><div>    I      0.000000000    0.000000000    0.747200000    </div><div>    I      0.500000000    0.500000000    0.247200000    </div><div>    I      0.214200000    0.714200000    0.004600000    </div><div>    I      0.714200000    0.214200000    0.504600000    </div><div>    I      0.785800000    0.285800000    0.004600000    </div><div>    I      0.285800000    0.785800000    0.504600000    </div><div>    I      0.285800000    0.214200000    0.004600000    </div><div>    I      0.785800000    0.714200000    0.504600000    </div><div>    I      0.714200000    0.785800000    0.004600000    </div><div>    I      0.214200000    0.285800000    0.504600000    </div><div>    N      0.500000000    0.000000000    0.242000000    </div><div>    N      0.000000000    0.500000000    0.742000000    </div><div>    N      0.000000000    0.500000000    0.242000000    </div><div>    N      0.500000000    0.000000000    0.742000000    </div><div>   Pb      0.000000000    0.000000000    0.000000000    </div><div>   Pb      0.500000000    0.500000000    0.500000000    </div><div>   Pb      0.000000000    0.000000000    0.500000000    </div><div>   Pb      0.500000000    0.500000000    0.000000000    </div><div>K_POINTS automatic </div><div>  4 4 4   0 0 0 </div></div><div>----------END OF FILE--------</div><div><br></div><div><a href="http://ph.in">ph.in</a>:</div><div><br></div><div><div>'MAPbI3-no-H-NC-PZ'        </div><div> &INPUTPH</div><div>                      outdir = '/home/quantum/Descargas/mapbi-real/' ,</div><div>                      prefix = 'MAPbI3-real' ,</div><div>                      fildyn = 'dmat.MAPbI3' ,</div><div>                       trans = .true.,</div><div>                       epsil = .true.,</div><div>                      lraman = .true.,</div><div>                     recover = .false. ,</div><div>                       qplot = .false. ,</div><div>                    amass(1) = 207.20000,</div><div>                    amass(2) = 126.90400,</div><div>                    amass(3) = 12.01070,</div><div>                    amass(4) = 14.00670,</div><div>                 max_seconds = 700000,</div><div> /</div><div>0 0 0 </div></div><div><br></div><div>------END OF FILE------</div><div><br></div><div><br></div><div><br></div><div>Best regards, </div><div><br></div><div>Josué Clavijo</div><div><br></div></div></div>