<div dir="ltr"><div dir="ltr">On Sat, Jan 11, 2020 at 11:00 AM Victor Bermudez <<a href="mailto:bermudez@alum.mit.edu">bermudez@alum.mit.edu</a>> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Recent versions of QE print a table of "atom charges" upon convergence of each SCF cycle. What do these numbers mean ? </blockquote><div><br></div><div>they are meant
to give a  vague idea of where the magnetization is located, I think. They are electronic charges (number of electrons, I guess) and magnetization - difference between up and down charges, integrated in a spherical region of radius r_m (printed in output).<br></div><div><br></div><div>Paolo</div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I'm aware that accurate estimates of atom charges require something like a Bader analysis, but I'd still like to know how to interpret these "atom charge" values. <br>
<br>
         For example, here's output (input given below) for a free NO2 molecule. This is just a "simple" example of a general result. <br>
<br>
Magnetic moment per site:<br>
     atom:    1    charge:    1.9288    magn:    0.1688    constr:    0.0000<br>
     atom:    2    charge:    2.9353    magn:    0.1673    constr:    0.0000<br>
     atom:    3    charge:    2.9339    magn:    0.1665    constr:    0.0000 <br>
<br>
Atom 1 is N, and atoms 2 and 3 are O. The core charges (Zval) are +5 for N and +6 for O, giving a total of +17. For a charge-neutral molecule I would therefore expect the "charge" number to sum to 17 if these numbers represent electron densities. Instead the sum is 7.798 |e|. If the charges are in Rydberg atomic units (1 au = |e|/rt(2)) then the sum is actually 5.515 |e|. <br>
<br>
Another question I have concerns the "magnetization" values. What quantities do they represent, and what are the units ? <br>
<br>
         Here's the input (for a cubic lattice of molecules with a 20-Angstrom lattice constant): <br>
*******************************************************************************<br>
&CONTROL<br>
calculation='relax',<br>
restart_mode='from_scratch',<br>
title='free NO2; start with #237406 opt.; 200/800 cut-offs',<br>
pseudo_dir="/p/home/bermudvm/Pseudo",<br>
outdir="/p/work1/workspace/bermudvm/2753418.THUNDER/tmp",<br>
! outdir='/lustre/cmf/scratch/b/bermudez/171133/',<br>
etot_conv_thr=1.0D-5,<br>
forc_conv_thr=1.0D-4,<br>
nstep=500,<br>
verbosity='default',<br>
max_seconds=21500.0<br>
/<br>
<br>
&SYSTEM<br>
ibrav=1,<br>
a=20.0,<br>
nat=3,<br>
ntyp=2,<br>
ecutwfc=200.0,<br>
ecutrho=800.0,<br>
occupations='smearing',<br>
smearing='gaussian',<br>
degauss=0.001,<br>
nspin=2,<br>
starting_magnetization(1)=0.20,<br>
starting_magnetization(2)=0.0,<br>
! nosym=.true.,<br>
! nbnd=14,<br>
vdw_corr='TS'<br>
/<br>
<br>
&ELECTRONS<br>
electron_maxstep=100,<br>
conv_thr=1.0D-9,<br>
mixing_beta=0.7,<br>
mixing_mode='plain'<br>
! mixing_ndim=8<br>
/<br>
<br>
&IONS<br>
ion_dynamics='bfgs',<br>
upscale=10000.D0<br>
/<br>
<br>
ATOMIC_SPECIES<br>
N 14.0067 N_ONCV_PBE-1.0.upf <br>
O 15.9994 O_ONCV_PBE-1.0.upf<br>
<br>
ATOMIC_POSITIONS angstrom<br>
N       -0.016847220  -0.000000025  -0.014422960<br>
O        1.092527105   0.000000576   0.456691560<br>
O       -1.125679885  -0.000000551   0.457731400<br>
<br>
K_POINTS gamma<br>
<br>
*************************************************************<br>
<br>
Any advice would be much appreciated.<br>
<br>
Vic Bermudez<br>
US Naval Research Lab. (retired) <br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</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>
</blockquote></div><div class="gmail_quote"><br></div><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>