<div>Dear users,</div><div> </div><div>I am performing MD simulations of a material at T=300K. I've noticed that if I use the following options:</div><div> </div><div><div> pot_extrapolation = 'second_order'</div><div>  pot_extrapolation = 'second_order'</div><div> </div><div>SCF steps are quite unstable (but converge anyway), starting from first MD step:</div><div> </div><div><div>scf convergence threshold =      1.0E-10</div><div>     estimated scf accuracy    <       0.00000080 Ry</div><div>     estimated scf accuracy    <       0.00000044 Ry</div><div>     estimated scf accuracy    <       0.00000016 Ry</div><div>     estimated scf accuracy    <          9.5E-09 Ry</div><div>     estimated scf accuracy    <          2.1E-09 Ry</div><div>     estimated scf accuracy    <          2.0E-10 Ry</div><div>     estimated scf accuracy    <          1.0E-10 Ry</div><div>     estimated scf accuracy    <          2.0E-11 Ry</div><div>     estimated scf accuracy    <       0.00982726 Ry</div><div>     estimated scf accuracy    <       0.00335630 Ry</div><div>     estimated scf accuracy    <       0.00036212 Ry</div><div>     estimated scf accuracy    <       0.00002729 Ry</div><div>     estimated scf accuracy    <       0.00001055 Ry</div><div>     estimated scf accuracy    <       0.00001287 Ry</div><div>     estimated scf accuracy    <       0.00000066 Ry</div><div>     estimated scf accuracy    <       0.00000097 Ry</div><div>     estimated scf accuracy    <       0.00000007 Ry</div><div>     estimated scf accuracy    <       0.00000007 Ry</div><div>     estimated scf accuracy    <       0.00000002 Ry</div><div>     estimated scf accuracy    <          8.7E-09 Ry</div><div>     estimated scf accuracy    <          2.2E-09 Ry</div><div>     estimated scf accuracy    <          1.1E-10 Ry</div><div>     estimated scf accuracy    <          9.8E-12 Ry</div><div>     estimated scf accuracy    <     124.55557104 Ry</div><div>     estimated scf accuracy    <     566.20308384 Ry</div><div>     estimated scf accuracy    <      24.52678095 Ry</div><div>     estimated scf accuracy    <      11.49168590 Ry</div><div>     estimated scf accuracy    <       4.74770768 Ry</div><div>     estimated scf accuracy    <       2.84584072 Ry</div><div>     estimated scf accuracy    <       1.82872297 Ry</div><div>     estimated scf accuracy    <       1.13633834 Ry</div><div>     estimated scf accuracy    <       0.85163312 Ry</div><div>     estimated scf accuracy    <       0.75225475 Ry</div><div>     estimated scf accuracy    <       0.23630576 Ry</div><div>     estimated scf accuracy    <       0.06657040 Ry</div><div>     estimated scf accuracy    <       0.12356423 Ry</div><div>     estimated scf accuracy    <       0.10745706 Ry</div><div>     estimated scf accuracy    <       0.07377183 Ry</div><div>     estimated scf accuracy    <       0.03204377 Ry</div><div>     estimated scf accuracy    <       0.01108520 Ry</div><div>     estimated scf accuracy    <       0.00237517 Ry</div><div>     estimated scf accuracy    <       0.00051729 Ry</div><div>     estimated scf accuracy    <       0.00010551 Ry</div><div>     estimated scf accuracy    <       0.00001862 Ry</div><div>     estimated scf accuracy    <       0.00000216 Ry</div><div>     estimated scf accuracy    <       0.00000081 Ry</div><div>     estimated scf accuracy    <       0.00000112 Ry</div><div>     estimated scf accuracy    <       0.00000038 Ry</div><div>     estimated scf accuracy    <       0.00000017 Ry</div><div>     estimated scf accuracy    <       0.00000002 Ry</div><div>     estimated scf accuracy    <          2.1E-09 Ry</div><div>     estimated scf accuracy    <          5.2E-09 Ry</div><div>     estimated scf accuracy    <          7.1E-10 Ry</div><div>     estimated scf accuracy    <          2.1E-11 Ry</div><div>     estimated scf accuracy    <      82.01884485 Ry</div><div>     estimated scf accuracy    <     541.29488792 Ry</div><div>     estimated scf accuracy    <      30.68489908 Ry</div><div>     estimated scf accuracy    <       9.32622744 Ry</div><div>     estimated scf accuracy    <       5.74410066 Ry</div><div>     estimated scf accuracy    <       2.78523601 Ry</div><div>...</div><div> </div><div>This is a how SCF accuracy looks without those steps:</div><div> </div><div><div>scf convergence threshold =      1.0E-10</div><div>     estimated scf accuracy    <       0.00000012 Ry</div><div>     estimated scf accuracy    <       0.00000003 Ry</div><div>     estimated scf accuracy    <          4.4E-09 Ry</div><div>     estimated scf accuracy    <          3.9E-10 Ry</div><div>     estimated scf accuracy    <          6.1E-11 Ry</div><div>     estimated scf accuracy    <       0.00832406 Ry</div><div>     estimated scf accuracy    <       0.00236207 Ry</div><div>     estimated scf accuracy    <       0.00032964 Ry</div><div>     estimated scf accuracy    <       0.00002279 Ry</div><div>     estimated scf accuracy    <       0.00000535 Ry</div><div>     estimated scf accuracy    <       0.00000077 Ry</div><div>     estimated scf accuracy    <       0.00000017 Ry</div><div>     estimated scf accuracy    <       0.00000003 Ry</div><div>     estimated scf accuracy    <       0.00000001 Ry</div><div>     estimated scf accuracy    <          5.6E-10 Ry</div><div>     estimated scf accuracy    <          4.6E-10 Ry</div><div>     estimated scf accuracy    <          3.3E-11 Ry</div><div>     estimated scf accuracy    <       0.00491362 Ry</div><div>.....</div><div> </div><div>When other system (magnetic) with second_order extrapolation for pot/wfc, it crashes after few steps (too many bands are not converged, at the beginning of SCF cycle), while without those options it works fine.</div><div> </div><div>I wonder if this is an expected behavior of MD simulations. I thought pot/wfc extrapolation was supposed to help.</div><div> </div><div>Thanks,</div><div> Sergey</div><div> </div><div>USF</div></div></div></div>