<div dir="ltr">Dear Nicola,<div><br></div><div>Thanks a lot for the input - I will give it a read!</div><div><br></div><div>Chris </div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Tue, 25 Feb 2025 at 16:31, Nicola Marzari <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
<br>
Dear Chris,<br>
<br>
<br>
we discuss this problem quite extensively here: <br>
<a href="https://pubs.acs.org/doi/10.1021/acs.jctc.9b00552" rel="noreferrer" target="_blank">https://pubs.acs.org/doi/10.1021/acs.jctc.9b00552</a><br>
<br>
In addition, approximate DFT functionals do a poor job in describing <br>
these extra electrons, that are very delocalized (interestingly, CCSD(T) <br>
also struggles because the basis sets don't really extend far enough).<br>
<br>
So, not sure there are ideal solutions, and if what you get is actually <br>
accurate; Perdew-Zunger SIC applied only to the extra electrons could <br>
also help.<br>
<br>
Some of this is also discussed here, in terms of weakly dipole-bound <br>
electrons: <a href="https://pubs.acs.org/doi/10.1021/acs.jctc.6b00145" rel="noreferrer" target="_blank">https://pubs.acs.org/doi/10.1021/acs.jctc.6b00145</a> .<br>
<br>
But in general it's a tough problem, so don't expect a simple solution.<br>
<br>
                        nicola<br>
<br>
<br>
<br>
<br>
On 25/02/2025 02:44, Christoph Wolf wrote:<br>
> Dear all,<br>
> <br>
> I am facing a tricky system: I am trying to calculate the anion and di- <br>
> anion of a small polycyclic aromatic hydrocarbon system (~20-25 atoms).<br>
> <br>
> I am using ONCV PPs and B3LYP just for completeness.<br>
> <br>
> The neutral molecule is closed shell (66 electrons for example) and the <br>
> anion (67 e-) is a radical, both converge fast and to meaningful <br>
> results. However it seems almost impossible to me to stabilize the di- <br>
> anion (68 e-). A closed shell (total magnetic moment = 0) does not <br>
> converge at all, a bi-radical (mtot=2) converges numerically but the <br>
> resulting cell seems to have a lot of charge "outside the molecule" when <br>
> I calculate the charge difference between anion and di-anion (pp <br>
> plot_num =0 and then take the difference between two calculations). Also <br>
> the projwfc.x output for the di-anion has a lot of "nondescript" states:<br>
> <br>
> <br>
>   ==== e(  33) =    -4.37227 eV ====<br>
>      |psi|^2 = 0.985<br>
> ==== e(  34) =    -1.39249 eV ====<br>
>      |psi|^2 = 0.001<br>
> ==== e(  35) =    -1.07395 eV ====<br>
> ...<br>
> ==== e(  50) =    -0.36725 eV ====<br>
>      |psi|^2 = 0.933<br>
> <br>
> <br>
> Any trick to help the system localize the additional charge better?<br>
> <br>
> Any help or comment is as always much appreciated!<br>
> <br>
> Ciao!<br>
> Chris<br>
> -- <br>
> Christoph Wolf<br>
> IBS Center for Quantum Nanoscience<br>
> Seoul, South Korea<br>
> <br>
> <a href="https://stm25.org/" rel="noreferrer" target="_blank">https://stm25.org/</a> <<a href="https://stm25.org/" rel="noreferrer" target="_blank">https://stm25.org/</a>><br>
> <br>
> _______________________________________________________________________________<br>
> The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.<br>
> --------------------------------------------------------------------------------<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, SNSF<br>
Laboratory Head, PSI Center for Scientific Computing, Theory, and Data<br>
Contact info and websites at <a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a><br>
</blockquote></div><div><br clear="all"></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr">Group Leader "Theory of Quantum Systems at Surfaces"<br>IBS Center for Quantum Nanoscience<br>Seoul, South Korea</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>