<div class="iw_mail" dir="auto" style="font-size: 13px;"><div>Dear QE Community,</div><div><br></div><div>Hello.</div><div>I am interested to know how can one find the best Hubbard occupations (or the best in a given model subspace ) in DFT+U on the fly?</div><div>The question arises from the fact that in an AIMD run which the system stays on BO-PES, within a DFT+U an update to occupation matrices is needed at each time-step. </div><div>Any comments is highly appreciated.</div><div><br></div><div>Best regards, and</div><div>BEST WISHES in this NEW YEAR<br></div><div><br><div>Mahmoud Payami</div><div>NSTRI, AEOI</div><div>Tehran, Iran</div><div>Email: mpayami@aeoi.org.ir</div><div>Phone: +98(0)2182066504</div><div>------------------------------------<br></div></div></div>