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