<div dir="ltr"><div>Dear all,<br><br></div><div>It appears that there is a problem restarting a vc-relax calculation. The problem happens when convergence is not achieved in the first scf cycle and this calculation is restarted. The vc-relax calculation will proceed. However, it will yield this error at the end of a calculation<br><br><pre> A final scf calculation at the relaxed structure.
 The G-vectors are recalculated for the final unit cell
 Results may differ from those at the preceding step.
....
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine read_rho_xml (1):
     dimensions do not match
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...<br><br></pre><pre>I have tried using wf_collect = .true. I also used the same number of processors for the first calculation and the restart. <br>I have checked the output directory and  I have also used the same machine for the first calculation and restart.<br>Increasing the electron_maxstep = 200 allowed for the calculation to complete the calculation without the need for a restart and also without any errors.<br></pre><pre><br></pre><pre>Sincerely,<br><br></pre><pre>Dae Kwang Jun<br></pre></div></div>