hi everyone, <br>&nbsp; &nbsp;&nbsp; When I do raman calculation of Li2O, I found the &quot;second order response&quot; process is not convergent. The log is below:<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 1&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp;&nbsp; 8.1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 
0.100E-01 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.208E-03<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 2&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 13.1&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.144E-02 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.905E-01<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 3&nbsp;&nbsp; <a href="http://av.it">
av.it</a>.:&nbsp; 11.0&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.100E-01 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.228E-03<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 4&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp;&nbsp; 8.8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.151E-02 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.178E-03<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 5&nbsp;&nbsp; 
<a href="http://av.it">av.it</a>.:&nbsp; 10.9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.133E-02 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.118E-08&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 6&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 12.6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.343E-05 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 
0.155E-07<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 7&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp;&nbsp; 9.6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.124E-04 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.279E-08<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 8&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 11.9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.529E-05
 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.191E-08<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp;&nbsp; 9&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 11.9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.437E-05 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.402E-07<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp; 10&nbsp;&nbsp; <a href="http://av.it">av.it
</a>.:&nbsp; 10.9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.201E-04 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.437E-08<br>.............................(omitted)<br>&nbsp; &nbsp; &nbsp; iter #&nbsp; 42&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 10.7 &nbsp; &nbsp;&nbsp; thresh= 0.855E-07 alpha_mix =&nbsp; 
0.700 |ddv_scf|^2 =&nbsp; 0.201E-11<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp; 43&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 11.9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.142E-06 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.181E-12<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp; 44&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 13.6
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.425E-07 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.546E-10<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp; 45&nbsp;&nbsp; <a href="http://av.it">av.it</a>.:&nbsp; 10.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.739E-06 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.190E-11<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iter #&nbsp; 46&nbsp;&nbsp; <a href="http://av.it">
av.it</a>.:&nbsp; 11.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; thresh= 0.138E-06 alpha_mix =&nbsp; 0.700 |ddv_scf|^2 =&nbsp; 0.111E-11<br><br>and the ph input file as:<br>&amp;inputph<br>&nbsp;&nbsp; tr2_ph=1.0d-14,<br>&nbsp;&nbsp; prefix=&#39;li2O&#39;,<br>&nbsp;&nbsp; amass(1)= 6.941,<br>&nbsp;&nbsp; amass(2)=
15.999,<br>&nbsp;&nbsp; fildyn=&#39;li2O.dynG&#39;,<br>&nbsp;&nbsp; epsil=.true.,<br>&nbsp;&nbsp; trans=.true.,<br>&nbsp;&nbsp; lraman=.true.,<br>&nbsp;&nbsp; outdir=&#39;/home/qgr/_work_/li2O/raman/temp&#39;,<br>/<br>0 0 0 <br>&nbsp;&nbsp;&nbsp; And in the scf input file, conv_thr =&nbsp; 
1.0d-12 with my own Li and O PPs generated using fhi98PP code.<br>&nbsp;&nbsp;&nbsp; How can I solve the problem?<br>&nbsp; &nbsp; What is more wried, this situation never happened days before when using the exactly same input file and PPs. That is , I cannot repeat the result again with all the same input file though I check all the input file carefully. The scf result are the same, but the 
ph.x gave different convergence process. Besides, I got different raman frequencies.<br><br>Thanks in advance.<br><br>Bests, Qian<br><br>