<div dir="ltr"><div><span class="gmail-im"><br></span></div><div>I tried to compute energy of antiferromagnetic (AFM) V3Al using the below script:</div><div><br></div><div><br></div><div><span class="gmail-im">&CONTROL<br>  calculation = 'scf'<br>  prefix = 'v3al'<br>  outdir = './'<br>  pseudo_dir = './'<br>  tprnfor = .true.<br>  tstress = .true.<br>  etot_conv_thr =   4.0000000000d-05<br>  forc_conv_thr =   1.0000000000d-04 <br>  verbosity = 'high'<br>/<br><br>&SYSTEM<br></span>  ibrav = 0<br>! celldm(1) = 11.512210722<br>! space_group = 225<br>  nat = 16<br>  ntyp = 3<span class="gmail-im"><br>  ecutwfc = 60<br>  ecutrho = 400<br>  occupations = 'smearing'<br>  smearing = 'cold'<br>  degauss = 0.005<br>  nspin = 2 <br>  starting_magnetization(1)=1<br>  starting_magnetization(2)=-1<br></span>  starting_magnetization(3)=1<br>  starting_magnetization(4)=-1<br>  starting_magnetization(5)=1<br>  starting_magnetization(6)=-1<br>  starting_magnetization(7)=1<br>  starting_magnetization(8)=-1<span class="gmail-im"><br>/<br><br>&ELECTRONS<br>  conv_thr =   8.0000000000d-10<br>  mixing_beta =   4.0000000000d-01<br>/<br><br>ATOMIC_SPECIES<br>Al 26.9815 Al.pbe-nl-rrkjus_psl.1.0.0.UPF<br>V1 50.9415 V.pbe-spnl-rrkjus_psl.1.0.0.UPF<br></span>V2 50.9415 V.pbe-spnl-rrkjus_psl.1.0.0.UPF<br> <br>ATOMIC_POSITIONS {crystal}<br>V2            0.2500000000       0.2500000000       0.7500000000<br>V2            0.2500000000       0.2500000000       0.2500000000 <br>V2            0.7500000000       0.7500000000       0.7500000000 <br>V2            0.7500000000       0.7500000000       0.2500000000 <br>V2            0.7500000000       0.2500000000       0.7500000000 <br>V2            0.2500000000       0.7500000000       0.2500000000 <br>V2            0.2500000000       0.7500000000       0.7500000000 <br>V2            0.7500000000       0.2500000000       0.2500000000 <br>Al            0.0000000000       0.0000000000       0.0000000000<br>Al            0.0000000000       0.5000000000       0.5000000000<br>Al            0.5000000000       0.0000000000       0.5000000000<br>Al            0.5000000000       0.5000000000       0.0000000000<br>V1            0.5000000000       0.5000000000       0.5000000000<br>V1            0.5000000000       0.0000000000       0.0000000000<br>V1            0.0000000000       0.5000000000       0.0000000000<br>V1            0.0000000000       0.0000000000       0.5000000000<br><br>K_POINTS automatic<br>  10 10 10 0 0 0<br><br>CELL_PARAMETERS angstrom<br>6.092 0.000 0.000<br>0.000 6.092 0.000<br>0.000 0.000 6.092</div><div><br></div><div>-------------------------------------------------------------------------</div><div><br></div><div>But as the iterations proceed, I see that the values of total and absolute magnetisation both approaches zero value.</div><div>     total magnetization       =     0.03 Bohr mag/cell<br> 
    absolute magnetization    =     0.14 Bohr mag/cell (This is at 27th 
iteration and calculation is still going as energy is not yet converged.
 So, at final iteration, probably both will become zero)</div><div><br></div><div>Shouldn't
 atleast the absolute magnetisation be nonzero value? I guess for AFM, 
only total magnetisation at the end should be zero.</div><div>I'm saying
 this because when I had performed AFM calculation for chromium, total 
magnetisation was 0.0 but absolute was 2.28 bohr mag/cell.</div><div>are things going correct? Attached below is the cif file. <br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 18, 2024 at 4:14 PM Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it">paolo.giannozzi@uniud.it</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">On 6/18/24 08:59, Pooja Vyas wrote:<br>
<br>
>    space_group = 225<br>
<br>
> ATOMIC_POSITIONS {crystal_sg}<br>
> Al 4a<br>
> V1 4b<br>
> V2 8c<br>
> V3 0.25 0.25 0.75<br>
<br>
>       Error in routine check_atoms (1):<br>
>       atoms #   3 and #   6 overlap!<br>
<br>
they do: atom V3 overlaps with one of the two positions for atoms V2.<br>
See <a href="https://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-wp-list" rel="noreferrer" target="_blank">https://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-wp-list</a> for <br>
group 225.<br>
<br>
4a = 4 atoms in the conventional cell (cubic), 1 atom in the primitive <br>
cell (fcc). Same for 4b and 8c. So you have nat=3 (number of Wyckoff <br>
positions, not of atoms) and 4 atoms, 1 Al and 3 V, if you set<br>
ATOMIC_POSITIONS {crystal_sg}<br>
Al 4a<br>
V1 4b<br>
V2 8c<br>
Note that ntyp=3: there are two different types of V atoms, not three.<br>
Disclaimer: just guessing<br>
<br>
Paolo<br>
<br>
>       stopping ...<br>
> <br>
> <br>
> <br>
> ----------------------------------------------------------------------------<br>
> <br>
> What is causing the issue and what should be changed in the input file?<br>
> Any kind of help is appreciated.<br>
> <br>
> <br>
> *Pooja,<br>
> *<br>
> *Research Scholar,*<br>
> *Department of Physics,*<br>
> *MKBU,*<br>
> *Bhavnagar,*<br>
> *Gujarat.*<br>
> *<br>
> *<br>
> <br>
> _______________________________________________<br>
> The Quantum ESPRESSO community stands by the Ukrainian<br>
> people and expresses its concerns about the devastating<br>
> effects that the Russian military offensive has on their<br>
> country and on the free and peaceful scientific, cultural,<br>
> and economic cooperation amongst peoples<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>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 206, 33100 Udine Italy, +39-0432-558216<br>
X SCUOLA ESTIVA DI ENERGIE RINNOVABILI <a href="https://scuola-rinnovabili.uniud.it/" rel="noreferrer" target="_blank">https://scuola-rinnovabili.uniud.it/</a><br>
</blockquote></div>