<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Thank you very much Natalie, Paolo and Giuseppe<div class="">You’ve been very helpful, cause I had the suspect SCAN wasn’t applicable to this subject but I was stuck trying to find tricks to make it work.</div><div class="">I’d have just a last question on the topic. I was able to calculate the single point energy of molecular oxygen (O2) using SCAN, hence calculate the absorption energy for a bridged oxygen on coronene at least with respect to O2. Shall I consider this result satisfactory, having used the same grid that I used for other functionals (PBE, PBE0, …)? The absorption energy is consistent with the others.</div><div class=""><br class=""></div><div class="">Thank you again for your help, have a great day!</div><div class=""><br class=""></div><div class="">Filippo</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Filippo Savazzi, PhD Student</div><div class="">Politecnico di Torino, Torino, Italy</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 21 Nov 2018, at 21:26, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" class="">p.giannozzi@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class="">I tend to agree with Natalie's conclusions. Unfortunately I am not aware of any simple solution. All meta-GGA functionals do strange things in systems with regions of vanishing charge, and SCAN is the one that gives less problems (!). SCF divergence is maybe extreme, but there are also problems like this one: <a href="https://gitlab.com/QEF/q-e/issues/32" class="">https://gitlab.com/QEF/q-e/issues/32</a>, indicating that one may need a surprisingly dense FFT grid to get reliable results<br class=""></div><div dir="ltr" class=""><br class=""></div><div class="">Paolo<br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="">On Wed, Nov 21, 2018 at 8:22 PM Holzwarth, Natalie <<a href="mailto:natalie@wfu.edu" class="">natalie@wfu.edu</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="">I have looked into the SCAN functional for an isolated atom, also finding that Vxc increases to unphysically large values as the density decreases exponentially. We were able to demonstrate for ourselves that the SCAN functional itself is responsible for this behavior due to the appearance of the density in denominators of various terms. It is my understanding that the SCAN developers are aware of this issue and are working on a remedy. In the meantime, perhaps it would be good to use the SCAN functional only for systems that do not have exponentially decreasing densities. At least that was my conclusion. Sincerely, Natalie Holzwarth<div class=""><br clear="all" class=""><div class=""><div dir="ltr" class="m_3234151201806222252gmail_signature" data-smartmail="gmail_signature">N. A. W. Holzwarth email: <a href="mailto:natalie@wfu.edu" target="_blank" class="">natalie@wfu.edu</a><div class="">Department of Physics web: <a href="http://www.wfu.edu/~natalie" target="_blank" class="">http://www.wfu.edu/~natalie</a></div><div class="">Wake Forest University phone: 1-336-758-5510 </div><div class="">Winston-Salem, NC 27109 USA office: Rm. 300 Olin Physical Lab</div></div></div><br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="">On Wed, Nov 21, 2018 at 1:26 PM Filippo Savazzi <<a href="mailto:filo.savazzi@gmail.com" target="_blank" class="">filo.savazzi@gmail.com</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thanks Giuseppe,<br class="">
I’ll give it a try, I didn’t think about setting the occupation from input.<br class="">
Thank you for your help.<br class="">
<br class="">
Best, <br class="">
<br class="">
Filippo<br class="">
<br class="">
<br class="">
<br class="">
Filippo Savazzi, PhD Student<br class="">
Politecnico di Torino, Torino, Italy<br class="">
<br class="">
> On 21 Nov 2018, at 19:00, Giuseppe Mattioli <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank" class="">giuseppe.mattioli@ism.cnr.it</a>> wrote:<br class="">
> <br class="">
> <br class="">
> Dear Filippo<br class="">
> Isolated open-shell atoms are tricky. Please look here for hints<br class="">
> <br class="">
> yourpathtoQE6.3/PW/examples/example05<br class="">
> <br class="">
> HTH<br class="">
> Giuseppe<br class="">
> <br class="">
> Quoting Filippo Savazzi <<a href="mailto:filo.savazzi@gmail.com" target="_blank" class="">filo.savazzi@gmail.com</a>>:<br class="">
> <br class="">
>> Dear QE users,<br class="">
>> <br class="">
>> I’m Filippo Savazzi, a PhD student from Politecnico di Torino, Italy.<br class="">
>> I’m using QE since a couple of years, and just browsing the archives of this mailing list or enquiring the mighty google I’ve been lucky enough that this is actually the first time I have to post a question.<br class="">
>> Thanks everybody for the indirect support you gave me in these years.<br class="">
>> <br class="">
>> I’m calculating absorption energies of oxygen on a coronene molecule, as part of a little evaluation of different available XC functionals that I’m doing. I’m calculating the absorption energy both with respect to the molecular oxygen and with atomic oxygen.<br class="">
>> My problems arise when I try to calculate the energy of a single atom of oxygen using SCAN. In this case the energy diverges until the SCF fails. I have no troubles running the calculation on the actual system (coronene + bridged O), neither on a molecular oxygen (triplet O2); the issue is just localised on the single oxygen atom. My first guess was that the implementation of SCAN has either troubles in dealing with isolated systems or with spin polarization, but I ruled out these issues when I was able to calculate the single point energy of an oxygen molecule. Are you aware of any other problems using SCAN in these conditions?<br class="">
>> In the following I attach the input for pw.x I’m using, as well as a brief example of the SCF part in the output. I already played with cell dimension (thinking too much vacuum could have been the problem), as well as with ecutwfc and ecutrho (so indirectly with the integration grid), and then, just to give it a chance, with the diagonalization method and beta-mixing.<br class="">
>> <br class="">
>> Thank you in advance for your attention.<br class="">
>> <br class="">
>> Best,<br class="">
>> Filippo<br class="">
>> <br class="">
>> <br class="">
>> Filippo Savazzi, PhD Student<br class="">
>> Politecnico di Torino, Torino, Italy<br class="">
>> <br class="">
>> <br class="">
>> <br class="">
>> <br class="">
>> &CONTROL<br class="">
>> calculation = ’scf',<br class="">
>> restart_mode='from_scratch',<br class="">
>> pseudo_dir = '/home/filippo/pseudo',<br class="">
>> outdir='./tmp_o',<br class="">
>> prefix='oxy',<br class="">
>> /<br class="">
>> &SYSTEM<br class="">
>> ibrav=6,<br class="">
>> celldm(1)=20,<br class="">
>> celldm(3)=1,<br class="">
>> nat=1,<br class="">
>> ntyp=1,<br class="">
>> input_dft='scan',<br class="">
>> ecutwfc=100,<br class="">
>> tot_magnetization=2,<br class="">
>> nspin=2,<br class="">
>> /<br class="">
>> &ELECTRONS<br class="">
>> diagonalization = 'cg',<br class="">
>> electron_maxstep = 300,<br class="">
>> mixing_mode = 'local-TF',<br class="">
>> mixing_beta = 0.7,<br class="">
>> conv_thr = 1.0d-7,<br class="">
>> /<br class="">
>> &IONS<br class="">
>> ion_dynamics='damp',<br class="">
>> upscale=1000,<br class="">
>> /<br class="">
>> ATOMIC_SPECIES<br class="">
>> O 15.99 O_ONCV_PBE-1.0.upf<br class="">
>> ATOMIC_POSITIONS {angstrom}<br class="">
>> O 0.0000 0.0000 0.0000<br class="">
>> K_POINTS{gamma}<br class="">
>> <br class="">
>> <br class="">
>> <br class="">
>> <br class="">
>> EXTRACT OF OUTPUT:<br class="">
>> <br class="">
>> Dense grid: 540376 G-vectors FFT dimensions: ( 128, 128, 128)<br class="">
>> <br class="">
>> Estimated max dynamical RAM per process > 116.06 MB<br class="">
>> <br class="">
>> Estimated total dynamical RAM > 1.36 GB<br class="">
>> Generating pointlists ...<br class="">
>> new r_m : 0.4125 (alat units) 8.2500 (a.u.) for type 1<br class="">
>> <br class="">
>> Initial potential from superposition of free atoms<br class="">
>> <br class="">
>> starting charge 5.99905, renormalised to 6.00000<br class="">
>> <br class="">
>> negative rho (up,down): 4.917E-05 4.917E-05<br class="">
>> Starting wfcs are random<br class="">
>> <br class="">
>> total cpu time spent up to now is 3.1 secs<br class="">
>> <br class="">
>> Self-consistent Calculation<br class="">
>> <br class="">
>> iteration # 1 ecut= 100.00 Ry beta= 0.70<br class="">
>> CG style diagonalization<br class="">
>> ethr = 1.00E-02, avg # of iterations = 9.6<br class="">
>> <br class="">
>> negative rho (up,down): 1.044E-03 1.363E-03<br class="">
>> <br class="">
>> total cpu time spent up to now is 240.5 secs<br class="">
>> <br class="">
>> total energy = 2429.63502350 Ry<br class="">
>> Harris-Foulkes estimate = -31.41616940 Ry<br class="">
>> estimated scf accuracy < 3.28037744 Ry<br class="">
>> <br class="">
>> total magnetization = 2.00 Bohr mag/cell<br class="">
>> absolute magnetization = 2.00 Bohr mag/cell<br class="">
>> <br class="">
>> iteration # 2 ecut= 100.00 Ry beta= 0.70<br class="">
>> CG style diagonalization<br class="">
>> c_bands: 4 eigenvalues not converged<br class="">
>> ethr = 1.00E-02, avg # of iterations = 13.6<br class="">
>> <br class="">
>> negative rho (up,down): 2.824E-03 3.815E-03<br class="">
>> <br class="">
>> total cpu time spent up to now is 331.7 secs<br class="">
>> <br class="">
>> total energy = 33638.43350785 Ry<br class="">
>> Harris-Foulkes estimate = 318.68713998 Ry<br class="">
>> estimated scf accuracy < 0.37936007 Ry<br class="">
>> <br class="">
>> total magnetization = 2.00 Bohr mag/cell<br class="">
>> absolute magnetization = 2.00 Bohr mag/cell<br class="">
> <br class="">
> <br class="">
> <br class="">
> GIUSEPPE MATTIOLI<br class="">
> CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br class="">
> Via Salaria Km 29,300 - C.P. 10<br class="">
> I-00015 - Monterotondo Scalo (RM)<br class="">
> Mob (*preferred*) +39 373 7305625<br class="">
> Tel + 39 06 90672342 - Fax +39 06 90672316<br class="">
> E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank" class="">giuseppe.mattioli@ism.cnr.it</a>><br class="">
> <br class="">
> _______________________________________________<br class="">
> users mailing list<br class="">
> <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class="">
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br class="">
<br class="">
_______________________________________________<br class="">
users mailing list<br class="">
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br class="">
users mailing list<br class="">
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all" class=""><br class="">-- <br class=""><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br class="">Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br class="">Phone +39-0432-558216, fax +39-0432-558222<br class=""><br class=""></div></div></div></div></div>
_______________________________________________<br class="">users mailing list<br class=""><a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class="">https://lists.quantum-espresso.org/mailman/listinfo/users</div></blockquote></div><br class=""></div></body></html>