<div>Dear PWSCF users,</div>
<div>      I have a question of how to efficiently restart the calculation of LSDA. It is a little bit complicated.</div>
<div>      The system I study is always metallic.</div>
<div>      When I relax a system using LDA, usually I can first relax it at a &#39;high temperature&#39; (a wider broading width for smearing) and a coarse k grid. Then I use the relaxed atomic positions and the charge density (charge-density.dat) to restart the relaxation at a &#39;lower temperature&#39; and a more dense k grid. </div>

<div>      However if I run a spin-polarized calculation (LSDA), I find that a system could be paramagnetic at &quot;high temperature&quot; and ferromagnetic at &#39;low temperature&#39;. However if I relax a big system using very small broading width and a very dense k grid, it is too computationally intensive. But if I first relax it using a wider broading width and a coarse k grid, the system will relax to a paramagnetic state. Then if I restart it, the system will always be in the paramagnetic state.  </div>

<div>      So is there a clear trick so that I can restart the calculation using the previous charge density but make a new initial guess for the magnetization rather than reading from the file &quot;spin-polarization.dat&#39;?</div>

<div>      Thank you very much.</div>
<div> </div>
<div>Hanghui Chen</div>
<div>Department of Physics</div>
<div>Yale University</div>