hi everyone, <br>     When I do raman calculation of Li2O, I found the "second order response" process is not convergent. The log is below:<br><br>      iter #   1   <a href="http://av.it">av.it</a>.:   8.1      thresh= 
0.100E-01 alpha_mix =  0.700 |ddv_scf|^2 =  0.208E-03<br>      iter #   2   <a href="http://av.it">av.it</a>.:  13.1     thresh= 0.144E-02 alpha_mix =  0.700 |ddv_scf|^2 =  0.905E-01<br>      iter #   3   <a href="http://av.it">
av.it</a>.:  11.0     thresh= 0.100E-01 alpha_mix =  0.700 |ddv_scf|^2 =  0.228E-03<br>      iter #   4   <a href="http://av.it">av.it</a>.:   8.8      thresh= 0.151E-02 alpha_mix =  0.700 |ddv_scf|^2 =  0.178E-03<br>      iter #   5   
<a href="http://av.it">av.it</a>.:  10.9      thresh= 0.133E-02 alpha_mix =  0.700 |ddv_scf|^2 =  0.118E-08      <br>      iter #   6   <a href="http://av.it">av.it</a>.:  12.6      thresh= 0.343E-05 alpha_mix =  0.700 |ddv_scf|^2 =  
0.155E-07<br>      iter #   7   <a href="http://av.it">av.it</a>.:   9.6      thresh= 0.124E-04 alpha_mix =  0.700 |ddv_scf|^2 =  0.279E-08<br>      iter #   8   <a href="http://av.it">av.it</a>.:  11.9      thresh= 0.529E-05
 alpha_mix =  0.700 |ddv_scf|^2 =  0.191E-08<br>      iter #   9   <a href="http://av.it">av.it</a>.:  11.9      thresh= 0.437E-05 alpha_mix =  0.700 |ddv_scf|^2 =  0.402E-07<br>      iter #  10   <a href="http://av.it">av.it
</a>.:  10.9      thresh= 0.201E-04 alpha_mix =  0.700 |ddv_scf|^2 =  0.437E-08<br>.............................(omitted)<br>      iter #  42   <a href="http://av.it">av.it</a>.:  10.7      thresh= 0.855E-07 alpha_mix =  
0.700 |ddv_scf|^2 =  0.201E-11<br>      iter #  43   <a href="http://av.it">av.it</a>.:  11.9      thresh= 0.142E-06 alpha_mix =  0.700 |ddv_scf|^2 =  0.181E-12<br>      iter #  44   <a href="http://av.it">av.it</a>.:  13.6
      thresh= 0.425E-07 alpha_mix =  0.700 |ddv_scf|^2 =  0.546E-10<br>      iter #  45   <a href="http://av.it">av.it</a>.:  10.2      thresh= 0.739E-06 alpha_mix =  0.700 |ddv_scf|^2 =  0.190E-11<br>      iter #  46   <a href="http://av.it">
av.it</a>.:  11.0      thresh= 0.138E-06 alpha_mix =  0.700 |ddv_scf|^2 =  0.111E-11<br><br>and the ph input file as:<br>&inputph<br>   tr2_ph=1.0d-14,<br>   prefix='li2O',<br>   amass(1)= 6.941,<br>   amass(2)=
15.999,<br>   fildyn='li2O.dynG',<br>   epsil=.true.,<br>   trans=.true.,<br>   lraman=.true.,<br>   outdir='/home/qgr/_work_/li2O/raman/temp',<br>/<br>0 0 0 <br>    And in the scf input file, conv_thr =  
1.0d-12 with my own Li and O PPs generated using fhi98PP code.<br>    How can I solve the problem?<br>    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>