<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hi Mehrdad,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">If you have a large energy difference between the end of the vc-relax algorithm and the final scf step, you may be using a fairly low plane wave cutoff energy for your system, though that may be fine depending on what quantities you are looking at. An important question here is: how are you choosing your convergence criteria? In my experience stresses (vc-relax) require a higher cutoff compared to total energy or forces (just ionic relaxation).</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">This all depends of course on what you are trying to get out of your calculation. In my case I was comparing different magnetic ground states that could be very close in energy, where very small changes in geometry can make a difference. For my specific systems I would converge plane wave cutoff, k-points and convergence threshold with respect to unit cell stress using a reference calculation with very high cutoff, low threshold (10^-9 Ry), and dense k-point mesh.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Just some food for thought. I've definitely seen papers in the literature where they claim state A has lower energy than state B. I've reproduced such results using the somewhat lax cutoffs reported, and then found when you actually do CONVERGED calculations, state B is actually lower in energy (whoops!). Convergence is important.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Best,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Kevin May, PhD<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Postdoctoral Associate</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Department of Materials Science and Engineering</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Massachusetts Institute of Technology</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><span style="font-family:arial,helvetica,sans-serif"><font size="2"><font size="1"><span style="color:rgb(102,102,102)"></span></font></font></span><span style="font-family:arial,helvetica,sans-serif"><font size="2"><font size="1"><span style="color:rgb(102,102,102)"></span></font></font></span></div></div></div></div></div></div></div>