<div><div dir="auto">Dear Ernane,</div></div><div dir="auto">Have you thought of using a more sophisticated method (like GW) on [CO3]- to calculate its EA? This would give you the energy of [CO3]2- in vacuum. </div><div><div dir="auto">Best,</div><div dir="auto">Michal Krompiec </div><div dir="auto">University of Southampton & Merck KGaA</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 15 Mar 2019 at 18:22, Ernane de Freitas Martins <<a href="mailto:ernanefmg@gmail.com">ernanefmg@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div>Dear Giuseppe, </div><div dir="auto"><br></div><div dir="auto">I really appreciate your answer. Thank you very much for using your time to answer my question.</div><div dir="auto"><br></div><div dir="auto">I'll think on your suggestion about trying hybrid functionals. The point is that I need to estimate the solvation energy for carbonate ion using the environ module, then I'll need to run a vacuum calculation using the same functional I'm already using rVV-10).</div><div dir="auto"><br></div><div dir="auto">Thank you again for replying.<br><br><div data-smartmail="gmail_signature" dir="auto">Atenciosamente,</div></div></div><div dir="auto"><div dir="auto"><div data-smartmail="gmail_signature" dir="auto"><br><br>Dr. Ernane de Freitas Martins<br>Postdoctoral researcher<br>IF - USP<br>São Paulo, SP - Brazil</div><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">Em sex, 15 de mar de 2019 15:04, Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Ernane<br>
Your question contains part of the answer! Carbonate ion (CO3 2-) is  <br>
not stable outside water, and calculations of its properties in gas  <br>
phase are likely not so meaningful, but in the case of model  <br>
thermodynamics cycles (e.g. Born-Haber). The excess negative charge is  <br>
unbound when not stabilized by a strongly polar solvent, and this is  <br>
likely responsible for instabilities in the construction of the  <br>
Kohn-Sham potential along scf iterations. Moreover, this happens on  <br>
top of the strong delocalization error you experience when you use a  <br>
standard GGA exchange-correlation functional, when the  <br>
self-interaction of strongly localized electrons in the J[n] Coulomb  <br>
potential is not cancelled by a same term in the semi local exchange  <br>
potential. You may minimize this latter source of error by using a  <br>
hybrid GGA-EXX functional such as B3LYP, where the non local   <br>
Hartree-Fock part of the exchange functional can recover part of the  <br>
delocalization error, but you are not free yet from the instability of  <br>
carbonate in gas phase.<br>
HTH<br>
Giuseppe<br>
<br>
Ernane de Freitas Martins <<a href="mailto:ernanefmg@gmail.com" rel="noreferrer" target="_blank">ernanefmg@gmail.com</a>> ha scritto:<br>
<br>
> Hello,<br>
><br>
> I'm experiencing a problem to run a negatively charge molecule in quantum<br>
> espresso. The system is CO32-.<br>
><br>
> I try both vacuum and solvated (environ) calculations. The solvated one<br>
> works fine.<br>
><br>
> The problem is the calculation in vacuum. It never give the first ionic<br>
> step because the SCF accuracy never reaches the convence criterion.<br>
><br>
> I tried many different solutions (increase cutoffs and box size, use assume<br>
> isolated, decreasing and changing the mixing scheme and etc) and nothing<br>
> works.<br>
><br>
> The unique calculation that works fine for vacuum is the one with a box<br>
> size of 7.9 x 7.9 x 7.9 A. I really don't understand why it only works for<br>
> this specific box size.<br>
><br>
> I ran several other charged systems (+1, +2 and -1 total charge) and all of<br>
> them worked fine. The problem appears for -2 total charge in vacuum.<br>
><br>
> Would some of you kindly help me in this?<br>
><br>
> Cheers,<br>
><br>
> Dr. Ernane de Freitas Martins<br>
> Postdoctoral researcher<br>
> IF - USP<br>
> São Paulo, SP - Brazil<br>
<br>
<br>
<br>
GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) +39 373 7305625<br>
Tel + 39 06 90672342 - Fax +39 06 90672316<br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" rel="noreferrer" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div></div></div>
_______________________________________________<br>
users mailing list<br>
<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></blockquote></div></div>