<div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I am facing a tricky system: I am trying to calculate the anion and di-anion of a small polycyclic aromatic hydrocarbon system (~20-25 atoms).</div><div><br></div><div>I am using ONCV PPs and B3LYP just for completeness.</div><div><br></div><div>The neutral molecule is closed shell (66 electrons for example) and the anion (67 e-) is a radical, both converge fast and to meaningful results. However it seems almost impossible to me to stabilize the di-anion (68 e-). A closed shell (total magnetic moment = 0) does not converge at all, a bi-radical (mtot=2) converges numerically but the resulting cell seems to have a lot of charge "outside the molecule" when I calculate the charge difference between anion and di-anion (pp plot_num =0 and then take the difference between two calculations). Also the projwfc.x output for the di-anion has a lot of "nondescript" states:</div><div><br></div><div><br></div><div> ==== e(  33) =    -4.37227 eV ==== </div><div>    |psi|^2 = 0.985</div><div>==== e(  34) =    -1.39249 eV ==== </div><div> </div><div>    |psi|^2 = 0.001</div><div>==== e(  35) =    -1.07395 eV ==== </div><div>...</div><div><div>==== e(  50) =    -0.36725 eV ==== </div><div>    |psi|^2 = 0.933</div></div><div><br></div><div><br></div><div>Any trick to help the system localize the additional charge better?</div><div><br></div><div>Any help or comment is as always much appreciated!</div><div><br></div><div>Ciao!</div><div>Chris </div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr">Christoph Wolf<br>IBS Center for Quantum Nanoscience<br>Seoul, South Korea</div></div></div></div></div>

<br>
<a href="https://stm25.org/" target="_blank"><img src="https://qns.science/wp-content/uploads/2024/12/stm25-email-footer.jpg" alt="https://stm25.org/"></a><br>