mirror of https://github.com/abinit/abinit.git
3239 lines
185 KiB
Plaintext
3239 lines
185 KiB
Plaintext
|
|
.Version 10.1.4.5 of ABINIT, released Sep 2024.
|
|
.(MPI version, prepared for a x86_64_linux_gnu13.2 computer)
|
|
|
|
.Copyright (C) 1998-2025 ABINIT group .
|
|
ABINIT comes with ABSOLUTELY NO WARRANTY.
|
|
It is free software, and you are welcome to redistribute it
|
|
under certain conditions (GNU General Public License,
|
|
see ~abinit/COPYING or http://www.gnu.org/copyleft/gpl.txt).
|
|
|
|
ABINIT is a project of the Universite Catholique de Louvain,
|
|
Corning Inc. and other collaborators, see ~abinit/doc/developers/contributors.txt .
|
|
Please read https://docs.abinit.org/theory/acknowledgments for suggested
|
|
acknowledgments of the ABINIT effort.
|
|
For more information, see https://www.abinit.org .
|
|
|
|
.Starting date : Fri 13 Sep 2024.
|
|
- ( at 19h07 )
|
|
|
|
- input file -> /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/TestBot_MPI1/tutorial_tlruj_1-tlruj_2-tlruj_3/tlruj_2.abi
|
|
- output file -> /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/TestBot_MPI1/tutorial_tlruj_1-tlruj_2-tlruj_3/tlruj_2.abo
|
|
- root for input files -> tlruj_2i
|
|
- root for output files -> tlruj_2.o
|
|
|
|
- inpspheads : Reading pseudopotential header in XML form from
|
|
- /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- inpspheads : Reading pseudopotential header in XML form from
|
|
- /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- inpspheads : Reading pseudopotential header in XML form from
|
|
- /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/O.xml
|
|
DATASET 1 : space group R-3 m (#166); Bravais hR (rhombohedral)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 1.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 18
|
|
lnmax = 6 mgfft = 30 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 2 nspden = 2 nspinor = 1
|
|
nsppol = 2 nsym = 12 n1xccc = 1 ntypat = 3
|
|
occopt = 0 xclevel = 2
|
|
- mband = 26 mffmem = 1 mkmem = 10
|
|
mpw = 376 nfft = 27000 nkpt = 10
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 40 nfftf = 64000
|
|
================================================================================
|
|
P This job should need less than 41.508 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 2.985 Mbytes ; DEN or POT disk file : 0.979 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 2 : space group R-3 m (#166); Bravais hR (rhombohedral)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 2.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 18
|
|
lnmax = 6 mgfft = 30 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 2 nspden = 2 nspinor = 1
|
|
nsppol = 2 nsym = 12 n1xccc = 1 ntypat = 3
|
|
occopt = 0 xclevel = 2
|
|
- mband = 26 mffmem = 1 mkmem = 10
|
|
mpw = 376 nfft = 27000 nkpt = 10
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 40 nfftf = 64000
|
|
================================================================================
|
|
P This job should need less than 41.508 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 2.985 Mbytes ; DEN or POT disk file : 0.979 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 3 : space group R-3 m (#166); Bravais hR (rhombohedral)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 3.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 18
|
|
lnmax = 6 mgfft = 30 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 2 nspden = 2 nspinor = 1
|
|
nsppol = 2 nsym = 12 n1xccc = 1 ntypat = 3
|
|
occopt = 0 xclevel = 2
|
|
- mband = 26 mffmem = 1 mkmem = 10
|
|
mpw = 376 nfft = 27000 nkpt = 10
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 40 nfftf = 64000
|
|
================================================================================
|
|
P This job should need less than 41.508 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 2.985 Mbytes ; DEN or POT disk file : 0.979 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 4 : space group R-3 m (#166); Bravais hR (rhombohedral)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 4.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 18
|
|
lnmax = 6 mgfft = 30 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 2 nspden = 2 nspinor = 1
|
|
nsppol = 2 nsym = 12 n1xccc = 1 ntypat = 3
|
|
occopt = 0 xclevel = 2
|
|
- mband = 26 mffmem = 1 mkmem = 10
|
|
mpw = 376 nfft = 27000 nkpt = 10
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 40 nfftf = 64000
|
|
================================================================================
|
|
P This job should need less than 41.508 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 2.985 Mbytes ; DEN or POT disk file : 0.979 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 5 : space group R-3 m (#166); Bravais hR (rhombohedral)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 5.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 18
|
|
lnmax = 6 mgfft = 30 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 2 nspden = 2 nspinor = 1
|
|
nsppol = 2 nsym = 12 n1xccc = 1 ntypat = 3
|
|
occopt = 0 xclevel = 2
|
|
- mband = 26 mffmem = 1 mkmem = 10
|
|
mpw = 376 nfft = 27000 nkpt = 10
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 40 nfftf = 64000
|
|
================================================================================
|
|
P This job should need less than 41.508 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 2.985 Mbytes ; DEN or POT disk file : 0.979 Mbytes.
|
|
================================================================================
|
|
|
|
--------------------------------------------------------------------------------
|
|
------------- Echo of variables that govern the present computation ------------
|
|
--------------------------------------------------------------------------------
|
|
-
|
|
- outvars: echo of selected default values
|
|
- iomode0 = 0 , fftalg0 =512 , wfoptalg0 = 10
|
|
-
|
|
- outvars: echo of global parameters not present in the input file
|
|
- max_nthreads = 0
|
|
-
|
|
-outvars: echo values of preprocessed input variables --------
|
|
acell 7.8800000000E+00 7.8800000000E+00 7.8800000000E+00 Bohr
|
|
amu 5.86900000E+01 5.86900000E+01 1.59994000E+01
|
|
atvshift1 -0.00551 -0.00551 -0.00551 -0.00551 -0.00551
|
|
-0.00551 -0.00551 -0.00551 -0.00551 -0.00551
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift2 -0.00367 -0.00367 -0.00367 -0.00367 -0.00367
|
|
-0.00367 -0.00367 -0.00367 -0.00367 -0.00367
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift3 -0.00184 -0.00184 -0.00184 -0.00184 -0.00184
|
|
-0.00184 -0.00184 -0.00184 -0.00184 -0.00184
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift4 0.00184 0.00184 0.00184 0.00184 0.00184
|
|
0.00184 0.00184 0.00184 0.00184 0.00184
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift5 0.00367 0.00367 0.00367 0.00367 0.00367
|
|
0.00367 0.00367 0.00367 0.00367 0.00367
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
chkprim 0
|
|
chksymbreak 0
|
|
diemix 4.50000000E-01
|
|
dmatpuopt 3
|
|
ecut 1.00000000E+01 Hartree
|
|
- fftalg 512
|
|
ixc -101130
|
|
jdtset 1 2 3 4 5
|
|
kpt 1.25000000E-01 1.25000000E-01 1.25000000E-01
|
|
3.75000000E-01 1.25000000E-01 1.25000000E-01
|
|
-3.75000000E-01 1.25000000E-01 1.25000000E-01
|
|
-1.25000000E-01 1.25000000E-01 1.25000000E-01
|
|
3.75000000E-01 3.75000000E-01 1.25000000E-01
|
|
-3.75000000E-01 3.75000000E-01 1.25000000E-01
|
|
-1.25000000E-01 3.75000000E-01 1.25000000E-01
|
|
-3.75000000E-01 -3.75000000E-01 1.25000000E-01
|
|
3.75000000E-01 3.75000000E-01 3.75000000E-01
|
|
-3.75000000E-01 3.75000000E-01 3.75000000E-01
|
|
kptrlatt 4 0 0 0 4 0 0 0 4
|
|
kptrlen 2.22880057E+01
|
|
lpawu 2 2 1
|
|
macro_uj 1
|
|
P mkmem 10
|
|
natom 4
|
|
nband 26
|
|
natvshift 5
|
|
ndtset 5
|
|
ngfft 30 30 30
|
|
ngfftdg 40 40 40
|
|
nkpt 10
|
|
nspden 2
|
|
nsppol 2
|
|
nsym 12
|
|
ntypat 3
|
|
occ 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
0.000000 0.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
0.000000 0.000000
|
|
occopt 0
|
|
pawecutdg 2.00000000E+01 Hartree
|
|
pawujv1 -5.5123988097E-03 Hartree
|
|
pawujv2 -3.6749325398E-03 Hartree
|
|
pawujv3 -1.8374662699E-03 Hartree
|
|
pawujv4 1.8374662699E-03 Hartree
|
|
pawujv5 3.6749325398E-03 Hartree
|
|
prtden 0
|
|
prtebands 0
|
|
prteig 0
|
|
prtgsr 0
|
|
prtwf 0
|
|
rprim 5.0000000000E-01 5.0000000000E-01 1.0000000000E+00
|
|
1.0000000000E+00 5.0000000000E-01 5.0000000000E-01
|
|
5.0000000000E-01 1.0000000000E+00 5.0000000000E-01
|
|
shiftk 5.00000000E-01 5.00000000E-01 5.00000000E-01
|
|
spgroup 166
|
|
spinat 0.0000000000E+00 0.0000000000E+00 3.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 -3.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
symrel 1 0 0 0 1 0 0 0 1 -1 0 0 0 -1 0 0 0 -1
|
|
0 1 0 1 0 0 0 0 1 0 -1 0 -1 0 0 0 0 -1
|
|
0 0 1 1 0 0 0 1 0 0 0 -1 -1 0 0 0 -1 0
|
|
1 0 0 0 0 1 0 1 0 -1 0 0 0 0 -1 0 -1 0
|
|
0 1 0 0 0 1 1 0 0 0 -1 0 0 0 -1 -1 0 0
|
|
0 0 1 0 1 0 1 0 0 0 0 -1 0 -1 0 -1 0 0
|
|
tolvrs 1.00000000E-09
|
|
typat 1 2 3 3
|
|
usepawu 1
|
|
useylm 1
|
|
wtk 0.03125 0.09375 0.09375 0.09375 0.09375 0.18750
|
|
0.18750 0.09375 0.03125 0.09375
|
|
xangst 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
4.1699164037E+00 4.1699164037E+00 4.1699164037E+00
|
|
2.0849582018E+00 2.0849582018E+00 2.0849582018E+00
|
|
6.2548746055E+00 6.2548746055E+00 6.2548746055E+00
|
|
xcart 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
7.8800000000E+00 7.8800000000E+00 7.8800000000E+00
|
|
3.9400000000E+00 3.9400000000E+00 3.9400000000E+00
|
|
1.1820000000E+01 1.1820000000E+01 1.1820000000E+01
|
|
xred 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
5.0000000000E-01 5.0000000000E-01 5.0000000000E-01
|
|
2.5000000000E-01 2.5000000000E-01 2.5000000000E-01
|
|
7.5000000000E-01 7.5000000000E-01 7.5000000000E-01
|
|
znucl 28.00000 28.00000 8.00000
|
|
|
|
================================================================================
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 1.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 2.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 3.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 4.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 5.
|
|
|
|
================================================================================
|
|
== DATASET 1 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 1, }
|
|
dimensions: {natom: 4, nkpt: 10, mband: 26, nsppol: 2, nspinor: 1, nspden: 2, mpw: 376, }
|
|
cutoff_energies: {ecut: 10.0, pawecutdg: 20.0, }
|
|
electrons: {nelect: 4.80000000E+01, charge: 0.00000000E+00, occopt: 0.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename: getwfk from: ./tlruj_1.o_WFK
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 3.9400000 3.9400000 7.8800000 G(1)= -0.0634518 -0.0634518 0.1903553
|
|
R(2)= 7.8800000 3.9400000 3.9400000 G(2)= 0.1903553 -0.0634518 -0.0634518
|
|
R(3)= 3.9400000 7.8800000 3.9400000 G(3)= -0.0634518 0.1903553 -0.0634518
|
|
Unit cell volume ucvol= 2.4465194E+02 bohr^3
|
|
Angles (23,13,12)= 3.35573098E+01 3.35573098E+01 3.35573098E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 30 30 30
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.18911
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 40 40 40
|
|
ecut(hartree)= 20.000 => boxcut(ratio)= 2.06006
|
|
|
|
--- Pseudopotential description ------------------------------------------------
|
|
- pspini: atom type 1 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- pspatm : Reading pseudopotential header in XML form from /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 6 (lmn_size= 18), orbitals= 0 0 1 1 2 2
|
|
Spheres core radius: rc_sph= 1.81432190
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=1500 , AA= 0.29784E-03 BB= 0.83395E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60094920
|
|
mmax= 1500
|
|
Radial grid used for partial waves is grid 1
|
|
Radial grid used for projectors is grid 1
|
|
Radial grid used for (t)core density is grid 1
|
|
Radial grid used for Vloc is grid 1
|
|
Radial grid used for pseudo valence density is grid 1
|
|
Mesh size for Vloc has been set to 1333 to avoid numerical noise.
|
|
Compensation charge density is not taken into account in XC energy/potential
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
- pspini: atom type 2 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
- pspatm : Reading pseudopotential header in XML form from /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ni.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 6 (lmn_size= 18), orbitals= 0 0 1 1 2 2
|
|
Spheres core radius: rc_sph= 1.81432190
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=1500 , AA= 0.29784E-03 BB= 0.83395E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60094920
|
|
mmax= 1500
|
|
Radial grid used for partial waves is grid 1
|
|
Radial grid used for projectors is grid 1
|
|
Radial grid used for (t)core density is grid 1
|
|
Radial grid used for Vloc is grid 1
|
|
Radial grid used for pseudo valence density is grid 1
|
|
Mesh size for Vloc has been set to 1333 to avoid numerical noise.
|
|
Compensation charge density is not taken into account in XC energy/potential
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
- pspini: atom type 3 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/O.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/O.xml
|
|
- pspatm : Reading pseudopotential header in XML form from /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/O.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 4 (lmn_size= 8), orbitals= 0 0 1 1
|
|
Spheres core radius: rc_sph= 1.41465230
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=2001 , AA= 0.72565E-03 BB= 0.58052E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.20231231
|
|
mmax= 2001
|
|
Radial grid used for partial waves is grid 1
|
|
Radial grid used for projectors is grid 1
|
|
Radial grid used for (t)core density is grid 1
|
|
Radial grid used for Vloc is grid 1
|
|
Radial grid used for pseudo valence density is grid 1
|
|
Mesh size for Vloc has been set to 1762 to avoid numerical noise.
|
|
Compensation charge density is not taken into account in XC energy/potential
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
2.55203824E+03 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file ./tlruj_1.o_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 368.875 368.830
|
|
|
|
|
|
******************************************
|
|
DFT+U Method used: FLL
|
|
******************************************
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 1, }
|
|
solver: {iscf: 17, nstep: 30, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-09, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -363.75876619049 -3.638E+02 1.102E-07 4.793E-01
|
|
ETOT 2 -363.76534979982 -6.584E-03 1.367E-06 1.057E+01
|
|
ETOT 3 -363.75845928167 6.891E-03 4.526E-04 2.860E-02
|
|
ETOT 4 -363.75846309243 -3.811E-06 1.401E-06 6.914E-04
|
|
ETOT 5 -363.75846322059 -1.282E-07 6.874E-07 3.264E-04
|
|
ETOT 6 -363.75846373295 -5.124E-07 1.876E-07 5.345E-04
|
|
ETOT 7 -363.75846369933 3.361E-08 5.982E-08 6.175E-05
|
|
ETOT 8 -363.75846382840 -1.291E-07 3.339E-08 2.851E-05
|
|
ETOT 9 -363.75846387180 -4.339E-08 1.489E-08 2.468E-05
|
|
ETOT 10 -363.75846394355 -7.176E-08 2.713E-08 1.216E-05
|
|
ETOT 11 -363.75846396989 -2.633E-08 6.301E-09 1.491E-06
|
|
ETOT 12 -363.75846397999 -1.010E-08 4.374E-09 4.494E-07
|
|
ETOT 13 -363.75846399065 -1.067E-08 1.511E-09 1.081E-07
|
|
ETOT 14 -363.75846399159 -9.430E-10 9.341E-10 1.567E-08
|
|
ETOT 15 -363.75846399192 -3.292E-10 2.787E-10 2.745E-09
|
|
ETOT 16 -363.75846399199 -6.372E-11 1.391E-10 1.574E-09
|
|
ETOT 17 -363.75846399200 -8.356E-12 4.817E-11 2.589E-09
|
|
ETOT 18 -363.75846399200 -2.444E-12 1.901E-11 3.916E-11
|
|
|
|
At SCF step 18 nres2 = 3.92E-11 < tolvrs= 1.00E-09 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42483938E-01 sigma(3 2)= -4.85585375E-04
|
|
sigma(2 2)= 1.42483938E-01 sigma(3 1)= -4.85585375E-04
|
|
sigma(3 3)= 1.42483938E-01 sigma(2 1)= -4.85585375E-04
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 1, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 3.9400000, 3.9400000, 7.8800000, ]
|
|
- [ 7.8800000, 3.9400000, 3.9400000, ]
|
|
- [ 3.9400000, 7.8800000, 3.9400000, ]
|
|
lattice_lengths: [ 9.65099, 9.65099, 9.65099, ]
|
|
lattice_angles: [ 33.557, 33.557, 33.557, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.4465194E+02
|
|
convergence: {deltae: -2.444E-12, res2: 3.916E-11, residm: 1.901E-11, diffor: null, }
|
|
etotal : -3.63758464E+02
|
|
entropy : 0.00000000E+00
|
|
fermie : 4.06886672E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 1.42483938E-01, -4.85585375E-04, -4.85585375E-04, ]
|
|
- [ -4.85585375E-04, 1.42483938E-01, -4.85585375E-04, ]
|
|
- [ -4.85585375E-04, -4.85585375E-04, 1.42483938E-01, ]
|
|
pressure_GPa: -4.1920E+03
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Ni]
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ni]
|
|
- [ 2.5000E-01, 2.5000E-01, 2.5000E-01, O]
|
|
- [ 7.5000E-01, 7.5000E-01, 7.5000E-01, O]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -0.00000000E+00, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ -0.00000000E+00, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ -3.84191871E-04, -3.84191871E-04, -3.84191871E-04, ]
|
|
- [ 3.84191871E-04, 3.84191871E-04, 3.84191871E-04, ]
|
|
force_length_stats: {min: 0.00000000E+00, max: 6.65439841E-04, mean: 3.32719921E-04, }
|
|
...
|
|
|
|
Integrated electronic and magnetization densities in atomic spheres:
|
|
---------------------------------------------------------------------
|
|
Radius=ratsph(iatom), smearing ratsm= 0.0000. Diff(up-dn)=approximate z local magnetic moment.
|
|
Atom Radius up_density dn_density Total(up+dn) Diff(up-dn)
|
|
1 1.81432 8.027804 7.872077 15.899880 0.155727
|
|
2 1.81432 7.866218 8.022224 15.888442 -0.156007
|
|
3 1.41465 2.292181 2.291988 4.584169 0.000192
|
|
4 1.41465 2.292181 2.291988 4.584169 0.000192
|
|
---------------------------------------------------------------------
|
|
Sum: 20.478383 20.478278 40.956660 0.000105
|
|
Total magnetization (from the atomic spheres): 0.000105
|
|
Total magnetization (exact up - dn): -0.000000
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 12.721410689071114
|
|
Compensation charge over fine fft grid = 12.721698036683415
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1 - Spin component 1
|
|
-2.32864 0.37495 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00170 -0.00170 0.00000 -0.00170 ...
|
|
0.37495 0.04455 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00034 0.00034 0.00000 0.00034 ...
|
|
0.00000 0.00000 -2.24018 -0.00125 -0.00125 0.82007 0.00076 0.00076 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00125 -2.24018 -0.00125 0.00076 0.82007 0.00076 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00125 -0.00125 -2.24018 0.00076 0.00076 0.82007 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82007 0.00076 0.00076 0.37153 -0.00046 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00076 0.82007 0.00076 -0.00046 0.37153 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00076 0.00076 0.82007 -0.00046 -0.00046 0.37153 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91366 -0.00101 0.00093 -0.00101 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.91366 -0.00047 -0.00101 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00093 -0.00047 -0.91493 -0.00047 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.00101 -0.00047 -0.91366 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
-2.32487 0.37426 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00223 -0.00223 0.00000 -0.00223 ...
|
|
0.37426 0.04466 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00044 0.00044 0.00000 0.00044 ...
|
|
0.00000 0.00000 -2.23702 -0.00165 -0.00165 0.81830 0.00099 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -2.23702 -0.00165 0.00099 0.81830 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -0.00165 -2.23702 0.00099 0.00099 0.81830 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.81830 0.00099 0.00099 0.37258 -0.00060 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.81830 0.00099 -0.00060 0.37258 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.00099 0.81830 -0.00060 -0.00060 0.37258 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91093 -0.00114 0.00137 -0.00114 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.91093 -0.00068 -0.00114 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00137 -0.00068 -0.91201 -0.00068 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.00114 -0.00068 -0.91093 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
0.61343 -1.75573 0.00001 0.00001 0.00001 -0.00004 -0.00004 -0.00004
|
|
-1.75573 5.08565 -0.00002 -0.00002 -0.00002 0.00010 0.00010 0.00010
|
|
0.00001 -0.00002 -0.29408 -0.00002 -0.00002 0.88876 0.00011 0.00011
|
|
0.00001 -0.00002 -0.00002 -0.29408 -0.00002 0.00011 0.88876 0.00011
|
|
0.00001 -0.00002 -0.00002 -0.00002 -0.29408 0.00011 0.00011 0.88876
|
|
-0.00004 0.00010 0.88876 0.00011 0.00011 0.11666 -0.00056 -0.00056
|
|
-0.00004 0.00010 0.00011 0.88876 0.00011 -0.00056 0.11666 -0.00056
|
|
-0.00004 0.00010 0.00011 0.00011 0.88876 -0.00056 -0.00056 0.11666
|
|
Atom # 4 - Spin component 2
|
|
0.61343 -1.75572 -0.00002 -0.00002 -0.00002 0.00011 0.00011 0.00011
|
|
-1.75572 5.08562 0.00004 0.00004 0.00004 -0.00029 -0.00029 -0.00029
|
|
-0.00002 0.00004 -0.29407 -0.00002 -0.00002 0.88874 0.00011 0.00011
|
|
-0.00002 0.00004 -0.00002 -0.29407 -0.00002 0.00011 0.88874 0.00011
|
|
-0.00002 0.00004 -0.00002 -0.00002 -0.29407 0.00011 0.00011 0.88874
|
|
0.00011 -0.00029 0.88874 0.00011 0.00011 0.11677 -0.00057 -0.00057
|
|
0.00011 -0.00029 0.00011 0.88874 0.00011 -0.00057 0.11677 -0.00057
|
|
0.00011 -0.00029 0.00011 0.00011 0.88874 -0.00057 -0.00057 0.11677
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1 - Spin component 1
|
|
0.97699 -0.16656 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00012 -0.00012 0.00000 -0.00012 ...
|
|
-0.16656 0.45536 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00320 -0.00320 0.00000 -0.00320 ...
|
|
0.00000 0.00000 1.06702 0.00002 0.00002 0.13588 0.00008 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 1.06702 0.00002 0.00008 0.13588 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 0.00002 1.06702 0.00008 0.00008 0.13588 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13588 0.00008 0.00008 0.02463 0.00002 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.13588 0.00008 0.00002 0.02463 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.00008 0.13588 0.00002 0.00002 0.02463 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00012 -0.00320 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.15589 -0.02172 0.02663 -0.02172 ...
|
|
-0.00012 -0.00320 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02172 1.15589 -0.01331 -0.02172 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.02663 -0.01331 1.06102 -0.01331 ...
|
|
-0.00012 -0.00320 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02172 -0.02172 -0.01331 1.15589 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
0.97691 -0.16711 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00014 -0.00014 0.00000 -0.00014 ...
|
|
-0.16711 0.45290 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00182 -0.00182 0.00000 -0.00182 ...
|
|
0.00000 0.00000 1.06700 0.00003 0.00003 0.13578 0.00012 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 1.06700 0.00003 0.00012 0.13578 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 0.00003 1.06700 0.00012 0.00012 0.13578 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13578 0.00012 0.00012 0.02454 0.00003 0.00003 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.13578 0.00012 0.00003 0.02454 0.00003 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.00012 0.13578 0.00003 0.00003 0.02454 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00014 -0.00182 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.11188 -0.01376 0.00074 -0.01376 ...
|
|
-0.00014 -0.00182 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01376 1.11188 -0.00037 -0.01376 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00074 -0.00037 1.01378 -0.00037 ...
|
|
-0.00014 -0.00182 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01376 -0.01376 -0.00037 1.11188 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
1.01156 0.01395 -0.00082 -0.00082 -0.00082 0.00003 0.00003 0.00003
|
|
0.01395 0.00119 0.00015 0.00015 0.00015 0.00000 0.00000 0.00000
|
|
-0.00082 0.00015 1.06817 -0.00323 -0.00323 0.05050 -0.00033 -0.00033
|
|
-0.00082 0.00015 -0.00323 1.06817 -0.00323 -0.00033 0.05050 -0.00033
|
|
-0.00082 0.00015 -0.00323 -0.00323 1.06817 -0.00033 -0.00033 0.05050
|
|
0.00003 0.00000 0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
0.00003 0.00000 -0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
0.00003 0.00000 -0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - Spin component 2
|
|
1.01152 0.01396 0.00053 0.00053 0.00053 -0.00003 -0.00003 -0.00003
|
|
0.01396 0.00119 -0.00007 -0.00007 -0.00007 0.00000 0.00000 0.00000
|
|
0.00053 -0.00007 1.06810 -0.00321 -0.00321 0.05050 -0.00033 -0.00033
|
|
0.00053 -0.00007 -0.00321 1.06810 -0.00321 -0.00033 0.05050 -0.00033
|
|
0.00053 -0.00007 -0.00321 -0.00321 1.06810 -0.00033 -0.00033 0.05050
|
|
-0.00003 0.00000 0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00003 0.00000 -0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00003 0.00000 -0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
|
|
"PAW+U" part of augmentation waves occupancies Rhoij:
|
|
Atom # 1 - L=2 ONLY - Spin component 1
|
|
1.15589 -0.02172 0.02663 -0.02172 0.00000 0.42978 -0.00881 0.01034 -0.00881 0.00000
|
|
-0.02172 1.15589 -0.01331 -0.02172 0.02306 -0.00881 0.42978 -0.00517 -0.00881 0.00895
|
|
0.02663 -0.01331 1.06102 -0.01331 0.00000 0.01204 -0.00602 0.39247 -0.00602 0.00000
|
|
-0.02172 -0.02172 -0.01331 1.15589 -0.02306 -0.00881 -0.00881 -0.00517 0.42978 -0.00895
|
|
0.00000 0.02306 0.00000 -0.02306 1.06102 0.00000 0.01043 0.00000 -0.01043 0.39247
|
|
0.42978 -0.00881 0.01204 -0.00881 0.00000 0.16436 -0.00361 0.00473 -0.00361 0.00000
|
|
-0.00881 0.42978 -0.00602 -0.00881 0.01043 -0.00361 0.16436 -0.00236 -0.00361 0.00410
|
|
0.01034 -0.00517 0.39247 -0.00517 0.00000 0.00473 -0.00236 0.15215 -0.00236 0.00000
|
|
-0.00881 -0.00881 -0.00602 0.42978 -0.01043 -0.00361 -0.00361 -0.00236 0.16436 -0.00410
|
|
0.00000 0.00895 0.00000 -0.00895 0.39247 0.00000 0.00410 0.00000 -0.00410 0.15215
|
|
Atom # 1 - L=2 ONLY - Spin component 2
|
|
1.11188 -0.01376 0.00074 -0.01376 0.00000 0.40969 -0.00539 -0.00014 -0.00539 0.00000
|
|
-0.01376 1.11188 -0.00037 -0.01376 0.00064 -0.00539 0.40969 0.00007 -0.00539 -0.00012
|
|
0.00074 -0.00037 1.01378 -0.00037 0.00000 -0.00033 0.00017 0.37453 0.00017 0.00000
|
|
-0.01376 -0.01376 -0.00037 1.11188 -0.00064 -0.00539 -0.00539 0.00007 0.40969 0.00012
|
|
0.00000 0.00064 0.00000 -0.00064 1.01378 0.00000 -0.00029 0.00000 0.00029 0.37453
|
|
0.40969 -0.00539 -0.00033 -0.00539 0.00000 0.15518 -0.00214 -0.00026 -0.00214 0.00000
|
|
-0.00539 0.40969 0.00017 -0.00539 -0.00029 -0.00214 0.15518 0.00013 -0.00214 -0.00023
|
|
-0.00014 0.00007 0.37453 0.00007 0.00000 -0.00026 0.00013 0.14529 0.00013 0.00000
|
|
-0.00539 -0.00539 0.00017 0.40969 0.00029 -0.00214 -0.00214 0.00013 0.15518 0.00023
|
|
0.00000 -0.00012 0.00000 0.00012 0.37453 0.00000 -0.00023 0.00000 0.00023 0.14529
|
|
Atom # 2 - L=2 ONLY - Spin component 1
|
|
1.11007 -0.01371 0.00055 -0.01371 0.00000 0.40876 -0.00537 -0.00026 -0.00537 0.00000
|
|
-0.01371 1.11007 -0.00028 -0.01371 0.00048 -0.00537 0.40876 0.00013 -0.00537 -0.00022
|
|
0.00055 -0.00028 1.01105 -0.00028 0.00000 -0.00050 0.00025 0.37338 0.00025 0.00000
|
|
-0.01371 -0.01371 -0.00028 1.11007 -0.00048 -0.00537 -0.00537 0.00013 0.40876 0.00022
|
|
0.00000 0.00048 0.00000 -0.00048 1.01105 0.00000 -0.00043 0.00000 0.00043 0.37338
|
|
0.40876 -0.00537 -0.00050 -0.00537 0.00000 0.15474 -0.00213 -0.00034 -0.00213 0.00000
|
|
-0.00537 0.40876 0.00025 -0.00537 -0.00043 -0.00213 0.15474 0.00017 -0.00213 -0.00030
|
|
-0.00026 0.00013 0.37338 0.00013 0.00000 -0.00034 0.00017 0.14484 0.00017 0.00000
|
|
-0.00537 -0.00537 0.00025 0.40876 0.00043 -0.00213 -0.00213 0.00017 0.15474 0.00030
|
|
0.00000 -0.00022 0.00000 0.00022 0.37338 0.00000 -0.00030 0.00000 0.00030 0.14484
|
|
Atom # 2 - L=2 ONLY - Spin component 2
|
|
1.15432 -0.02174 0.02644 -0.02174 0.00000 0.42901 -0.00882 0.01022 -0.00882 0.00000
|
|
-0.02174 1.15432 -0.01322 -0.02174 0.02290 -0.00882 0.42901 -0.00511 -0.00882 0.00885
|
|
0.02644 -0.01322 1.05818 -0.01322 0.00000 0.01190 -0.00595 0.39128 -0.00595 0.00000
|
|
-0.02174 -0.02174 -0.01322 1.15432 -0.02290 -0.00882 -0.00882 -0.00511 0.42901 -0.00885
|
|
0.00000 0.02290 0.00000 -0.02290 1.05818 0.00000 0.01031 0.00000 -0.01031 0.39128
|
|
0.42901 -0.00882 0.01190 -0.00882 0.00000 0.16402 -0.00361 0.00466 -0.00361 0.00000
|
|
-0.00882 0.42901 -0.00595 -0.00882 0.01031 -0.00361 0.16402 -0.00233 -0.00361 0.00404
|
|
0.01022 -0.00511 0.39128 -0.00511 0.00000 0.00466 -0.00233 0.15168 -0.00233 0.00000
|
|
-0.00882 -0.00882 -0.00595 0.42901 -0.01031 -0.00361 -0.00361 -0.00233 0.16402 -0.00404
|
|
0.00000 0.00885 0.00000 -0.00885 0.39128 0.00000 0.00404 0.00000 -0.00404 0.15168
|
|
Atom # 3 - L=1 ONLY - Spin component 1
|
|
1.06817 -0.00323 -0.00323 0.05050 -0.00033 -0.00033
|
|
-0.00323 1.06817 -0.00323 -0.00033 0.05050 -0.00033
|
|
-0.00323 -0.00323 1.06817 -0.00033 -0.00033 0.05050
|
|
0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
Atom # 3 - L=1 ONLY - Spin component 2
|
|
1.06810 -0.00321 -0.00321 0.05050 -0.00033 -0.00033
|
|
-0.00321 1.06810 -0.00321 -0.00033 0.05050 -0.00033
|
|
-0.00321 -0.00321 1.06810 -0.00033 -0.00033 0.05050
|
|
0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 1
|
|
1.06817 -0.00323 -0.00323 0.05050 -0.00033 -0.00033
|
|
-0.00323 1.06817 -0.00323 -0.00033 0.05050 -0.00033
|
|
-0.00323 -0.00323 1.06817 -0.00033 -0.00033 0.05050
|
|
0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 2
|
|
1.06810 -0.00321 -0.00321 0.05050 -0.00033 -0.00033
|
|
-0.00321 1.06810 -0.00321 -0.00033 0.05050 -0.00033
|
|
-0.00321 -0.00321 1.06810 -0.00033 -0.00033 0.05050
|
|
0.05050 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05050 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05050 -0.00001 -0.00001 0.00287
|
|
|
|
---------- DFT+U DATA ---------------------------------------------------
|
|
|
|
====== For Atom 1, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 1. Occ. for lpawu and for spin 1 = 4.41365
|
|
Atom 1. Occ. for lpawu and for spin 2 = 4.23843
|
|
=> On atom 1, local Mag. for lpawu is -0.175220
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.91241 -0.01692 0.02063 -0.01692 0.00000
|
|
-0.01692 0.91241 -0.01032 -0.01692 0.01787
|
|
0.02063 -0.01032 0.83821 -0.01032 0.00000
|
|
-0.01692 -0.01692 -0.01032 0.91241 -0.01787
|
|
0.00000 0.01787 0.00000 -0.01787 0.83821
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.87878 -0.01078 0.00073 -0.01078 0.00000
|
|
-0.01078 0.87878 -0.00037 -0.01078 0.00064
|
|
0.00073 -0.00037 0.80104 -0.00037 0.00000
|
|
-0.01078 -0.01078 -0.00037 0.87878 -0.00064
|
|
0.00000 0.00064 0.00000 -0.00064 0.80104
|
|
|
|
|
|
====== For Atom 2, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 2. Occ. for lpawu and for spin 1 = 4.23015
|
|
Atom 2. Occ. for lpawu and for spin 2 = 4.40570
|
|
=> On atom 2, local Mag. for lpawu is 0.175555
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.87743 -0.01074 0.00061 -0.01074 0.00000
|
|
-0.01074 0.87743 -0.00031 -0.01074 0.00053
|
|
0.00061 -0.00031 0.79893 -0.00031 0.00000
|
|
-0.01074 -0.01074 -0.00031 0.87743 -0.00053
|
|
0.00000 0.00053 0.00000 -0.00053 0.79893
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.91123 -0.01694 0.02050 -0.01694 0.00000
|
|
-0.01694 0.91123 -0.01025 -0.01694 0.01775
|
|
0.02050 -0.01025 0.83601 -0.01025 0.00000
|
|
-0.01694 -0.01694 -0.01025 0.91123 -0.01775
|
|
0.00000 0.01775 0.00000 -0.01775 0.83601
|
|
|
|
|
|
====== For Atom 3, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 3. Occ. for lpawu and for spin 1 = 1.55460
|
|
Atom 3. Occ. for lpawu and for spin 2 = 1.55446
|
|
=> On atom 3, local Mag. for lpawu is -0.000140
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51820 -0.00084 -0.00084
|
|
-0.00084 0.51820 -0.00084
|
|
-0.00084 -0.00084 0.51820
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51815 -0.00082 -0.00082
|
|
-0.00082 0.51815 -0.00082
|
|
-0.00082 -0.00082 0.51815
|
|
|
|
|
|
====== For Atom 4, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 4. Occ. for lpawu and for spin 1 = 1.55460
|
|
Atom 4. Occ. for lpawu and for spin 2 = 1.55446
|
|
=> On atom 4, local Mag. for lpawu is -0.000140
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51820 -0.00084 -0.00084
|
|
-0.00084 0.51820 -0.00084
|
|
-0.00084 -0.00084 0.51820
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51815 -0.00082 -0.00082
|
|
-0.00082 0.51815 -0.00082
|
|
-0.00082 -0.00082 0.51815
|
|
|
|
|
|
|
|
|
|
---------- calculate U, (J) start ----------
|
|
|
|
All pawujat ok and equal to 1
|
|
macro_uj=1 and nspden=2: standard determination of Hubbard U-parameter
|
|
|
|
|
|
*********************************************************************
|
|
************************ Linear Response U ************************
|
|
|
|
Info printed for perturbed atom: 1
|
|
|
|
Perturbations Occupations
|
|
--------------- -----------------------------
|
|
alpha [eV] Unscreened Screened
|
|
--------------- -----------------------------
|
|
0.0000000000 8.6380182458 8.6380182458
|
|
-0.1500000000 8.6964981917 8.6520722001
|
|
|
|
Scalar response functions:
|
|
Chi0 [eV^-1]: -0.86637
|
|
Chi [eV^-1]: -0.09369
|
|
|
|
The scalar U from the two-point regression scheme is 9.51891 eV.
|
|
*********************************************************************
|
|
*********************************************************************
|
|
Note: For more reliable linear regressions of the response
|
|
matrices, it is advised that you have more than two points.
|
|
See the LRUJ protocol for more information.
|
|
|
|
|
|
URES ii nat r_max U(J)[eV] U_ASA[eV] U_inf[eV]
|
|
URES 1 1 0.00000 2.37973 1.91505 1.85249
|
|
URES 2 8 11.14400 8.34374 6.71448 6.49514
|
|
URES 3 27 12.45940 9.16681 7.37684 7.13586
|
|
URES 4 64 22.28800 9.37021 7.54052 7.29420
|
|
URES 5 125 24.28780 9.44276 7.59891 7.35068
|
|
URES 6 216 33.43200 9.47484 7.62472 7.37565
|
|
|
|
---------- calculate U, (J) end --------------
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 10.668E-14; max= 19.014E-12
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
0.250000000000 0.250000000000 0.250000000000
|
|
0.750000000000 0.750000000000 0.750000000000
|
|
rms dE/dt= 4.2814E-03; max dE/dt= 6.0549E-03; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
2 0.000000000000 0.000000000000 0.000000000000
|
|
3 0.006054863894 0.006054863894 0.006054863894
|
|
4 -0.006054863894 -0.006054863894 -0.006054863894
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 0.00000000000000
|
|
2 4.16991640368920 4.16991640368920 4.16991640368920
|
|
3 2.08495820184460 2.08495820184460 2.08495820184460
|
|
4 6.25487460553380 6.25487460553380 6.25487460553380
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 -0.00038419187146 -0.00038419187146 -0.00038419187146
|
|
4 0.00038419187146 0.00038419187146 0.00038419187146
|
|
frms,max,avg= 2.7166468E-04 3.8419187E-04 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 -0.01975593869218 -0.01975593869218 -0.01975593869218
|
|
4 0.01975593869218 0.01975593869218 0.01975593869218
|
|
frms,max,avg= 1.3969558E-02 1.9755939E-02 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 7.880000000000 7.880000000000 7.880000000000 bohr
|
|
= 4.169916403689 4.169916403689 4.169916403689 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.40689 Average Vxc (hartree)= -0.46897
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN UP:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.13802 -4.13610 -2.36951 -2.33717 -2.33037 -2.33037 -2.30834 -2.30834
|
|
-0.40354 -0.35639 0.07010 0.21050 0.21050 0.22570 0.22570 0.25257
|
|
0.26128 0.26128 0.31826 0.31826 0.37388 0.38649 0.38649 0.39600
|
|
0.40070 0.43571
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN DOWN:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14324 -4.13088 -2.37027 -2.33643 -2.33217 -2.33217 -2.30654 -2.30654
|
|
-0.40354 -0.35638 0.07009 0.21051 0.21051 0.22570 0.22570 0.25258
|
|
0.26129 0.26129 0.31825 0.31825 0.37396 0.38650 0.38650 0.39597
|
|
0.40063 0.43568
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 1, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.94961346946556E+01
|
|
hartree : 8.65906815458428E+01
|
|
xc : -2.18505986225962E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
local_psp : -2.68648474205162E+02
|
|
spherical_terms : -2.03982622980908E+01
|
|
total_energy : -3.63710771380203E+02
|
|
total_energy_eV : -9.89707341404342E+03
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 1, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -4.00263798498696E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
xc_dc : -8.00302998810592E+01
|
|
spherical_terms : -1.48015317662165E+01
|
|
total_energy_dc : -3.63758463991998E+02
|
|
total_energy_dc_eV : -9.89837119601025E+03
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42483938E-01 sigma(3 2)= -4.85585375E-04
|
|
sigma(2 2)= 1.42483938E-01 sigma(3 1)= -4.85585375E-04
|
|
sigma(3 3)= 1.42483938E-01 sigma(2 1)= -4.85585375E-04
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -4.1920E+03 GPa]
|
|
- sigma(1 1)= 4.19202148E+03 sigma(3 2)= -1.42864126E+01
|
|
- sigma(2 2)= 4.19202148E+03 sigma(3 1)= -1.42864126E+01
|
|
- sigma(3 3)= 4.19202148E+03 sigma(2 1)= -1.42864126E+01
|
|
|
|
================================================================================
|
|
== DATASET 2 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 2, }
|
|
dimensions: {natom: 4, nkpt: 10, mband: 26, nsppol: 2, nspinor: 1, nspden: 2, mpw: 376, }
|
|
cutoff_energies: {ecut: 10.0, pawecutdg: 20.0, }
|
|
electrons: {nelect: 4.80000000E+01, charge: 0.00000000E+00, occopt: 0.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename: getwfk from: ./tlruj_1.o_WFK
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 3.9400000 3.9400000 7.8800000 G(1)= -0.0634518 -0.0634518 0.1903553
|
|
R(2)= 7.8800000 3.9400000 3.9400000 G(2)= 0.1903553 -0.0634518 -0.0634518
|
|
R(3)= 3.9400000 7.8800000 3.9400000 G(3)= -0.0634518 0.1903553 -0.0634518
|
|
Unit cell volume ucvol= 2.4465194E+02 bohr^3
|
|
Angles (23,13,12)= 3.35573098E+01 3.35573098E+01 3.35573098E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 30 30 30
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.18911
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 40 40 40
|
|
ecut(hartree)= 20.000 => boxcut(ratio)= 2.06006
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file ./tlruj_1.o_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 368.875 368.830
|
|
|
|
|
|
******************************************
|
|
DFT+U Method used: FLL
|
|
******************************************
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: 17, nstep: 30, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-09, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -363.74270119920 -3.637E+02 5.081E-08 1.911E-01
|
|
ETOT 2 -363.74528981774 -2.589E-03 5.358E-07 5.178E+00
|
|
ETOT 3 -363.74256595595 2.724E-03 1.722E-04 6.078E-03
|
|
ETOT 4 -363.74257000372 -4.048E-06 4.625E-07 1.327E-04
|
|
ETOT 5 -363.74257042999 -4.263E-07 2.689E-07 5.670E-04
|
|
ETOT 6 -363.74257028143 1.486E-07 1.143E-08 1.175E-04
|
|
ETOT 7 -363.74257030307 -2.164E-08 4.787E-09 4.088E-05
|
|
ETOT 8 -363.74257031227 -9.204E-09 1.394E-10 3.213E-05
|
|
ETOT 9 -363.74257033177 -1.950E-08 1.592E-09 1.636E-05
|
|
ETOT 10 -363.74257034028 -8.510E-09 1.394E-10 1.578E-05
|
|
ETOT 11 -363.74257035047 -1.019E-08 5.512E-09 1.591E-06
|
|
ETOT 12 -363.74257035665 -6.176E-09 3.857E-11 5.781E-07
|
|
ETOT 13 -363.74257036049 -3.839E-09 1.402E-10 1.317E-07
|
|
ETOT 14 -363.74257036259 -2.102E-09 6.508E-11 1.433E-08
|
|
ETOT 15 -363.74257036288 -2.885E-10 2.732E-11 1.280E-08
|
|
ETOT 16 -363.74257036294 -6.094E-11 1.125E-11 4.256E-09
|
|
ETOT 17 -363.74257036294 -1.819E-12 1.076E-12 4.364E-10
|
|
|
|
At SCF step 17 nres2 = 4.36E-10 < tolvrs= 1.00E-09 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42479676E-01 sigma(3 2)= -4.86299688E-04
|
|
sigma(2 2)= 1.42479676E-01 sigma(3 1)= -4.86299688E-04
|
|
sigma(3 3)= 1.42479676E-01 sigma(2 1)= -4.86299688E-04
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 2, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 3.9400000, 3.9400000, 7.8800000, ]
|
|
- [ 7.8800000, 3.9400000, 3.9400000, ]
|
|
- [ 3.9400000, 7.8800000, 3.9400000, ]
|
|
lattice_lengths: [ 9.65099, 9.65099, 9.65099, ]
|
|
lattice_angles: [ 33.557, 33.557, 33.557, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.4465194E+02
|
|
convergence: {deltae: -1.819E-12, res2: 4.364E-10, residm: 1.076E-12, diffor: null, }
|
|
etotal : -3.63742570E+02
|
|
entropy : 0.00000000E+00
|
|
fermie : 4.07284737E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 1.42479676E-01, -4.86299688E-04, -4.86299688E-04, ]
|
|
- [ -4.86299688E-04, 1.42479676E-01, -4.86299688E-04, ]
|
|
- [ -4.86299688E-04, -4.86299688E-04, 1.42479676E-01, ]
|
|
pressure_GPa: -4.1919E+03
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Ni]
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ni]
|
|
- [ 2.5000E-01, 2.5000E-01, 2.5000E-01, O]
|
|
- [ 7.5000E-01, 7.5000E-01, 7.5000E-01, O]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -0.00000000E+00, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ -0.00000000E+00, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ -2.57269193E-04, -2.57269193E-04, -2.57269193E-04, ]
|
|
- [ 2.57269193E-04, 2.57269193E-04, 2.57269193E-04, ]
|
|
force_length_stats: {min: 0.00000000E+00, max: 4.45603314E-04, mean: 2.22801657E-04, }
|
|
...
|
|
|
|
Integrated electronic and magnetization densities in atomic spheres:
|
|
---------------------------------------------------------------------
|
|
Radius=ratsph(iatom), smearing ratsm= 0.0000. Diff(up-dn)=approximate z local magnetic moment.
|
|
Atom Radius up_density dn_density Total(up+dn) Diff(up-dn)
|
|
1 1.81432 8.026377 7.870255 15.896632 0.156122
|
|
2 1.81432 7.866323 8.022678 15.889001 -0.156355
|
|
3 1.41465 2.292227 2.292090 4.584317 0.000137
|
|
4 1.41465 2.292227 2.292090 4.584317 0.000137
|
|
---------------------------------------------------------------------
|
|
Sum: 20.477154 20.477113 40.954267 0.000041
|
|
Total magnetization (from the atomic spheres): 0.000041
|
|
Total magnetization (exact up - dn): -0.000000
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 12.720567589566176
|
|
Compensation charge over fine fft grid = 12.720854498626926
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1 - Spin component 1
|
|
-2.32959 0.37520 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00170 -0.00170 0.00000 -0.00170 ...
|
|
0.37520 0.04451 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00034 0.00034 0.00000 0.00034 ...
|
|
0.00000 0.00000 -2.24104 -0.00125 -0.00125 0.82074 0.00076 0.00076 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00125 -2.24104 -0.00125 0.00076 0.82074 0.00076 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00125 -0.00125 -2.24104 0.00076 0.00076 0.82074 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82074 0.00076 0.00076 0.37110 -0.00046 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00076 0.82074 0.00076 -0.00046 0.37110 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00076 0.00076 0.82074 -0.00046 -0.00046 0.37110 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91299 -0.00101 0.00094 -0.00101 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.91299 -0.00047 -0.00101 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00094 -0.00047 -0.91427 -0.00047 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.00101 -0.00047 -0.91299 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
-2.32581 0.37450 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00223 -0.00223 0.00000 -0.00223 ...
|
|
0.37450 0.04461 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00044 0.00044 0.00000 0.00044 ...
|
|
0.00000 0.00000 -2.23787 -0.00165 -0.00165 0.81896 0.00099 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -2.23787 -0.00165 0.00099 0.81896 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -0.00165 -2.23787 0.00099 0.00099 0.81896 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.81896 0.00099 0.00099 0.37216 -0.00060 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.81896 0.00099 -0.00060 0.37216 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.00099 0.81896 -0.00060 -0.00060 0.37216 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91025 -0.00114 0.00137 -0.00114 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.91025 -0.00068 -0.00114 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00137 -0.00068 -0.91133 -0.00068 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.00114 -0.00068 -0.91025 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
0.61342 -1.75567 0.00001 0.00001 0.00001 -0.00005 -0.00005 -0.00005
|
|
-1.75567 5.08546 -0.00002 -0.00002 -0.00002 0.00014 0.00014 0.00014
|
|
0.00001 -0.00002 -0.29405 -0.00002 -0.00002 0.88859 0.00011 0.00011
|
|
0.00001 -0.00002 -0.00002 -0.29405 -0.00002 0.00011 0.88859 0.00011
|
|
0.00001 -0.00002 -0.00002 -0.00002 -0.29405 0.00011 0.00011 0.88859
|
|
-0.00005 0.00014 0.88859 0.00011 0.00011 0.11755 -0.00056 -0.00056
|
|
-0.00005 0.00014 0.00011 0.88859 0.00011 -0.00056 0.11755 -0.00056
|
|
-0.00005 0.00014 0.00011 0.00011 0.88859 -0.00056 -0.00056 0.11755
|
|
Atom # 4 - Spin component 2
|
|
0.61341 -1.75566 -0.00002 -0.00002 -0.00002 0.00010 0.00010 0.00010
|
|
-1.75566 5.08543 0.00004 0.00004 0.00004 -0.00026 -0.00026 -0.00026
|
|
-0.00002 0.00004 -0.29405 -0.00002 -0.00002 0.88858 0.00011 0.00011
|
|
-0.00002 0.00004 -0.00002 -0.29405 -0.00002 0.00011 0.88858 0.00011
|
|
-0.00002 0.00004 -0.00002 -0.00002 -0.29405 0.00011 0.00011 0.88858
|
|
0.00010 -0.00026 0.88858 0.00011 0.00011 0.11763 -0.00057 -0.00057
|
|
0.00010 -0.00026 0.00011 0.88858 0.00011 -0.00057 0.11763 -0.00057
|
|
0.00010 -0.00026 0.00011 0.00011 0.88858 -0.00057 -0.00057 0.11763
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1 - Spin component 1
|
|
0.97696 -0.16669 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00013 -0.00013 0.00000 -0.00013 ...
|
|
-0.16669 0.45680 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00338 -0.00338 0.00000 -0.00338 ...
|
|
0.00000 0.00000 1.06703 0.00002 0.00002 0.13583 0.00008 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 1.06703 0.00002 0.00008 0.13583 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 0.00002 1.06703 0.00008 0.00008 0.13583 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13583 0.00008 0.00008 0.02463 0.00002 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.13583 0.00008 0.00002 0.02463 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.00008 0.13583 0.00002 0.00002 0.02463 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00013 -0.00338 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.15545 -0.02173 0.02663 -0.02173 ...
|
|
-0.00013 -0.00338 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02173 1.15545 -0.01332 -0.02173 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.02663 -0.01332 1.06033 -0.01332 ...
|
|
-0.00013 -0.00338 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02173 -0.02173 -0.01332 1.15545 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
0.97689 -0.16722 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00013 -0.00013 0.00000 -0.00013 ...
|
|
-0.16722 0.45473 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00157 -0.00157 0.00000 -0.00157 ...
|
|
0.00000 0.00000 1.06700 0.00003 0.00003 0.13572 0.00012 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 1.06700 0.00003 0.00012 0.13572 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 0.00003 1.06700 0.00012 0.00012 0.13572 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13572 0.00012 0.00012 0.02455 0.00003 0.00003 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.13572 0.00012 0.00003 0.02455 0.00003 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.00012 0.13572 0.00003 0.00003 0.02455 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00013 -0.00157 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.11131 -0.01381 0.00070 -0.01381 ...
|
|
-0.00013 -0.00157 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01381 1.11131 -0.00035 -0.01381 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00070 -0.00035 1.01297 -0.00035 ...
|
|
-0.00013 -0.00157 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01381 -0.01381 -0.00035 1.11131 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
1.01156 0.01395 -0.00079 -0.00079 -0.00079 0.00003 0.00003 0.00003
|
|
0.01395 0.00119 0.00014 0.00014 0.00014 0.00000 0.00000 0.00000
|
|
-0.00079 0.00014 1.06822 -0.00322 -0.00322 0.05051 -0.00033 -0.00033
|
|
-0.00079 0.00014 -0.00322 1.06822 -0.00322 -0.00033 0.05051 -0.00033
|
|
-0.00079 0.00014 -0.00322 -0.00322 1.06822 -0.00033 -0.00033 0.05051
|
|
0.00003 0.00000 0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
0.00003 0.00000 -0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
0.00003 0.00000 -0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - Spin component 2
|
|
1.01153 0.01396 0.00060 0.00060 0.00060 -0.00003 -0.00003 -0.00003
|
|
0.01396 0.00119 -0.00009 -0.00009 -0.00009 0.00000 0.00000 0.00000
|
|
0.00060 -0.00009 1.06817 -0.00321 -0.00321 0.05051 -0.00033 -0.00033
|
|
0.00060 -0.00009 -0.00321 1.06817 -0.00321 -0.00033 0.05051 -0.00033
|
|
0.00060 -0.00009 -0.00321 -0.00321 1.06817 -0.00033 -0.00033 0.05051
|
|
-0.00003 0.00000 0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00003 0.00000 -0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00003 0.00000 -0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
|
|
"PAW+U" part of augmentation waves occupancies Rhoij:
|
|
Atom # 1 - L=2 ONLY - Spin component 1
|
|
1.15545 -0.02173 0.02663 -0.02173 0.00000 0.42958 -0.00882 0.01033 -0.00882 0.00000
|
|
-0.02173 1.15545 -0.01332 -0.02173 0.02306 -0.00882 0.42958 -0.00517 -0.00882 0.00895
|
|
0.02663 -0.01332 1.06033 -0.01332 0.00000 0.01204 -0.00602 0.39217 -0.00602 0.00000
|
|
-0.02173 -0.02173 -0.01332 1.15545 -0.02306 -0.00882 -0.00882 -0.00517 0.42958 -0.00895
|
|
0.00000 0.02306 0.00000 -0.02306 1.06033 0.00000 0.01042 0.00000 -0.01042 0.39217
|
|
0.42958 -0.00882 0.01204 -0.00882 0.00000 0.16428 -0.00361 0.00473 -0.00361 0.00000
|
|
-0.00882 0.42958 -0.00602 -0.00882 0.01042 -0.00361 0.16428 -0.00236 -0.00361 0.00409
|
|
0.01033 -0.00517 0.39217 -0.00517 0.00000 0.00473 -0.00236 0.15203 -0.00236 0.00000
|
|
-0.00882 -0.00882 -0.00602 0.42958 -0.01042 -0.00361 -0.00361 -0.00236 0.16428 -0.00409
|
|
0.00000 0.00895 0.00000 -0.00895 0.39217 0.00000 0.00409 0.00000 -0.00409 0.15203
|
|
Atom # 1 - L=2 ONLY - Spin component 2
|
|
1.11131 -0.01381 0.00070 -0.01381 0.00000 0.40942 -0.00542 -0.00016 -0.00542 0.00000
|
|
-0.01381 1.11131 -0.00035 -0.01381 0.00061 -0.00542 0.40942 0.00008 -0.00542 -0.00014
|
|
0.00070 -0.00035 1.01297 -0.00035 0.00000 -0.00036 0.00018 0.37419 0.00018 0.00000
|
|
-0.01381 -0.01381 -0.00035 1.11131 -0.00061 -0.00542 -0.00542 0.00008 0.40942 0.00014
|
|
0.00000 0.00061 0.00000 -0.00061 1.01297 0.00000 -0.00031 0.00000 0.00031 0.37419
|
|
0.40942 -0.00542 -0.00036 -0.00542 0.00000 0.15506 -0.00215 -0.00028 -0.00215 0.00000
|
|
-0.00542 0.40942 0.00018 -0.00542 -0.00031 -0.00215 0.15506 0.00014 -0.00215 -0.00024
|
|
-0.00016 0.00008 0.37419 0.00008 0.00000 -0.00028 0.00014 0.14515 0.00014 0.00000
|
|
-0.00542 -0.00542 0.00018 0.40942 0.00031 -0.00215 -0.00215 0.00014 0.15506 0.00024
|
|
0.00000 -0.00014 0.00000 0.00014 0.37419 0.00000 -0.00024 0.00000 0.00024 0.14515
|
|
Atom # 2 - L=2 ONLY - Spin component 1
|
|
1.11008 -0.01381 0.00058 -0.01381 0.00000 0.40879 -0.00541 -0.00024 -0.00541 0.00000
|
|
-0.01381 1.11008 -0.00029 -0.01381 0.00050 -0.00541 0.40879 0.00012 -0.00541 -0.00021
|
|
0.00058 -0.00029 1.01115 -0.00029 0.00000 -0.00047 0.00024 0.37342 0.00024 0.00000
|
|
-0.01381 -0.01381 -0.00029 1.11008 -0.00050 -0.00541 -0.00541 0.00012 0.40879 0.00021
|
|
0.00000 0.00050 0.00000 -0.00050 1.01115 0.00000 -0.00041 0.00000 0.00041 0.37342
|
|
0.40879 -0.00541 -0.00047 -0.00541 0.00000 0.15476 -0.00215 -0.00033 -0.00215 0.00000
|
|
-0.00541 0.40879 0.00024 -0.00541 -0.00041 -0.00215 0.15476 0.00017 -0.00215 -0.00029
|
|
-0.00024 0.00012 0.37342 0.00012 0.00000 -0.00033 0.00017 0.14486 0.00017 0.00000
|
|
-0.00541 -0.00541 0.00024 0.40879 0.00041 -0.00215 -0.00215 0.00017 0.15476 0.00029
|
|
0.00000 -0.00021 0.00000 0.00021 0.37342 0.00000 -0.00029 0.00000 0.00029 0.14486
|
|
Atom # 2 - L=2 ONLY - Spin component 2
|
|
1.15441 -0.02174 0.02651 -0.02174 0.00000 0.42907 -0.00883 0.01026 -0.00883 0.00000
|
|
-0.02174 1.15441 -0.01325 -0.02174 0.02296 -0.00883 0.42907 -0.00513 -0.00883 0.00888
|
|
0.02651 -0.01325 1.05844 -0.01325 0.00000 0.01195 -0.00597 0.39138 -0.00597 0.00000
|
|
-0.02174 -0.02174 -0.01325 1.15441 -0.02296 -0.00883 -0.00883 -0.00513 0.42907 -0.00888
|
|
0.00000 0.02296 0.00000 -0.02296 1.05844 0.00000 0.01035 0.00000 -0.01035 0.39138
|
|
0.42907 -0.00883 0.01195 -0.00883 0.00000 0.16405 -0.00361 0.00468 -0.00361 0.00000
|
|
-0.00883 0.42907 -0.00597 -0.00883 0.01035 -0.00361 0.16405 -0.00234 -0.00361 0.00405
|
|
0.01026 -0.00513 0.39138 -0.00513 0.00000 0.00468 -0.00234 0.15172 -0.00234 0.00000
|
|
-0.00883 -0.00883 -0.00597 0.42907 -0.01035 -0.00361 -0.00361 -0.00234 0.16405 -0.00405
|
|
0.00000 0.00888 0.00000 -0.00888 0.39138 0.00000 0.00405 0.00000 -0.00405 0.15172
|
|
Atom # 3 - L=1 ONLY - Spin component 1
|
|
1.06822 -0.00322 -0.00322 0.05051 -0.00033 -0.00033
|
|
-0.00322 1.06822 -0.00322 -0.00033 0.05051 -0.00033
|
|
-0.00322 -0.00322 1.06822 -0.00033 -0.00033 0.05051
|
|
0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
Atom # 3 - L=1 ONLY - Spin component 2
|
|
1.06817 -0.00321 -0.00321 0.05051 -0.00033 -0.00033
|
|
-0.00321 1.06817 -0.00321 -0.00033 0.05051 -0.00033
|
|
-0.00321 -0.00321 1.06817 -0.00033 -0.00033 0.05051
|
|
0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 1
|
|
1.06822 -0.00322 -0.00322 0.05051 -0.00033 -0.00033
|
|
-0.00322 1.06822 -0.00322 -0.00033 0.05051 -0.00033
|
|
-0.00322 -0.00322 1.06822 -0.00033 -0.00033 0.05051
|
|
0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 2
|
|
1.06817 -0.00321 -0.00321 0.05051 -0.00033 -0.00033
|
|
-0.00321 1.06817 -0.00321 -0.00033 0.05051 -0.00033
|
|
-0.00321 -0.00321 1.06817 -0.00033 -0.00033 0.05051
|
|
0.05051 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05051 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05051 -0.00001 -0.00001 0.00287
|
|
|
|
---------- DFT+U DATA ---------------------------------------------------
|
|
|
|
====== For Atom 1, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 1. Occ. for lpawu and for spin 1 = 4.41157
|
|
Atom 1. Occ. for lpawu and for spin 2 = 4.23586
|
|
=> On atom 1, local Mag. for lpawu is -0.175705
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.91207 -0.01693 0.02064 -0.01693 0.00000
|
|
-0.01693 0.91207 -0.01032 -0.01693 0.01787
|
|
0.02064 -0.01032 0.83768 -0.01032 0.00000
|
|
-0.01693 -0.01693 -0.01032 0.91207 -0.01787
|
|
0.00000 0.01787 0.00000 -0.01787 0.83768
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.87835 -0.01082 0.00071 -0.01082 0.00000
|
|
-0.01082 0.87835 -0.00036 -0.01082 0.00062
|
|
0.00071 -0.00036 0.80041 -0.00036 0.00000
|
|
-0.01082 -0.01082 -0.00036 0.87835 -0.00062
|
|
0.00000 0.00062 0.00000 -0.00062 0.80041
|
|
|
|
|
|
====== For Atom 2, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 2. Occ. for lpawu and for spin 1 = 4.23030
|
|
Atom 2. Occ. for lpawu and for spin 2 = 4.40630
|
|
=> On atom 2, local Mag. for lpawu is 0.175995
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.87743 -0.01081 0.00063 -0.01081 0.00000
|
|
-0.01081 0.87743 -0.00031 -0.01081 0.00055
|
|
0.00063 -0.00031 0.79901 -0.00031 0.00000
|
|
-0.01081 -0.01081 -0.00031 0.87743 -0.00055
|
|
0.00000 0.00055 0.00000 -0.00055 0.79901
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.91129 -0.01694 0.02055 -0.01694 0.00000
|
|
-0.01694 0.91129 -0.01028 -0.01694 0.01780
|
|
0.02055 -0.01028 0.83622 -0.01028 0.00000
|
|
-0.01694 -0.01694 -0.01028 0.91129 -0.01780
|
|
0.00000 0.01780 0.00000 -0.01780 0.83622
|
|
|
|
|
|
====== For Atom 3, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 3. Occ. for lpawu and for spin 1 = 1.55463
|
|
Atom 3. Occ. for lpawu and for spin 2 = 1.55453
|
|
=> On atom 3, local Mag. for lpawu is -0.000094
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51821 -0.00083 -0.00083
|
|
-0.00083 0.51821 -0.00083
|
|
-0.00083 -0.00083 0.51821
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51818 -0.00082 -0.00082
|
|
-0.00082 0.51818 -0.00082
|
|
-0.00082 -0.00082 0.51818
|
|
|
|
|
|
====== For Atom 4, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 4. Occ. for lpawu and for spin 1 = 1.55463
|
|
Atom 4. Occ. for lpawu and for spin 2 = 1.55453
|
|
=> On atom 4, local Mag. for lpawu is -0.000094
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51821 -0.00083 -0.00083
|
|
-0.00083 0.51821 -0.00083
|
|
-0.00083 -0.00083 0.51821
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51818 -0.00082 -0.00082
|
|
-0.00082 0.51818 -0.00082
|
|
-0.00082 -0.00082 0.51818
|
|
|
|
|
|
|
|
|
|
---------- calculate U, (J) start ----------
|
|
|
|
All pawujat ok and equal to 1
|
|
macro_uj=1 and nspden=2: standard determination of Hubbard U-parameter
|
|
|
|
|
|
*********************************************************************
|
|
************************ Linear Response U ************************
|
|
|
|
Info printed for perturbed atom: 1
|
|
|
|
Perturbations Occupations
|
|
--------------- -----------------------------
|
|
alpha [eV] Unscreened Screened
|
|
--------------- -----------------------------
|
|
0.0000000000 8.6380182458 8.6380182458
|
|
-0.1000000000 8.6747201572 8.6474287224
|
|
|
|
Scalar response functions:
|
|
Chi0 [eV^-1]: -0.81560
|
|
Chi [eV^-1]: -0.09410
|
|
|
|
The scalar U from the two-point regression scheme is 9.40036 eV.
|
|
*********************************************************************
|
|
*********************************************************************
|
|
Note: For more reliable linear regressions of the response
|
|
matrices, it is advised that you have more than two points.
|
|
See the LRUJ protocol for more information.
|
|
|
|
|
|
URES ii nat r_max U(J)[eV] U_ASA[eV] U_inf[eV]
|
|
URES 1 1 0.00000 2.35009 1.89120 1.82942
|
|
URES 2 8 11.14400 8.23982 6.63086 6.41425
|
|
URES 3 27 12.45940 9.05264 7.28496 7.04699
|
|
URES 4 64 22.28800 9.25351 7.44661 7.20336
|
|
URES 5 125 24.28780 9.32516 7.50427 7.25913
|
|
URES 6 216 33.43200 9.35684 7.52976 7.28379
|
|
|
|
---------- calculate U, (J) end --------------
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 53.992E-15; max= 10.759E-13
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
0.250000000000 0.250000000000 0.250000000000
|
|
0.750000000000 0.750000000000 0.750000000000
|
|
rms dE/dt= 2.8670E-03; max dE/dt= 4.0546E-03; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
2 0.000000000000 0.000000000000 0.000000000000
|
|
3 0.004054562485 0.004054562485 0.004054562485
|
|
4 -0.004054562485 -0.004054562485 -0.004054562485
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 0.00000000000000
|
|
2 4.16991640368920 4.16991640368920 4.16991640368920
|
|
3 2.08495820184460 2.08495820184460 2.08495820184460
|
|
4 6.25487460553380 6.25487460553380 6.25487460553380
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 -0.00025726919321 -0.00025726919321 -0.00025726919321
|
|
4 0.00025726919321 0.00025726919321 0.00025726919321
|
|
frms,max,avg= 1.8191679E-04 2.5726919E-04 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 -0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 -0.01322931271084 -0.01322931271084 -0.01322931271084
|
|
4 0.01322931271084 0.01322931271084 0.01322931271084
|
|
frms,max,avg= 9.3545367E-03 1.3229313E-02 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 7.880000000000 7.880000000000 7.880000000000 bohr
|
|
= 4.169916403689 4.169916403689 4.169916403689 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.40728 Average Vxc (hartree)= -0.46901
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN UP:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.13800 -4.13750 -2.37010 -2.33787 -2.33104 -2.33104 -2.30892 -2.30892
|
|
-0.40344 -0.35627 0.07016 0.21073 0.21073 0.22613 0.22613 0.25268
|
|
0.26145 0.26145 0.31868 0.31868 0.37432 0.38694 0.38694 0.39642
|
|
0.40071 0.43613
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN DOWN:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14275 -4.13276 -2.37061 -2.33737 -2.33227 -2.33227 -2.30770 -2.30770
|
|
-0.40343 -0.35626 0.07015 0.21074 0.21074 0.22613 0.22613 0.25269
|
|
0.26146 0.26146 0.31868 0.31868 0.37437 0.38695 0.38695 0.39640
|
|
0.40066 0.43610
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 2, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.94890712607126E+01
|
|
hartree : 8.65750590804409E+01
|
|
xc : -2.18498794956413E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
local_psp : -2.68626248098313E+02
|
|
spherical_terms : -2.03985418272519E+01
|
|
total_energy : -3.63710791574906E+02
|
|
total_energy_eV : -9.89707396356921E+03
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 2, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -4.00276152129402E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
xc_dc : -8.00148744559162E+01
|
|
spherical_terms : -1.47998281992300E+01
|
|
total_energy_dc : -3.63742570362939E+02
|
|
total_energy_dc_eV : -9.89793870836900E+03
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42479676E-01 sigma(3 2)= -4.86299688E-04
|
|
sigma(2 2)= 1.42479676E-01 sigma(3 1)= -4.86299688E-04
|
|
sigma(3 3)= 1.42479676E-01 sigma(2 1)= -4.86299688E-04
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -4.1919E+03 GPa]
|
|
- sigma(1 1)= 4.19189608E+03 sigma(3 2)= -1.43074284E+01
|
|
- sigma(2 2)= 4.19189608E+03 sigma(3 1)= -1.43074284E+01
|
|
- sigma(3 3)= 4.19189608E+03 sigma(2 1)= -1.43074284E+01
|
|
|
|
================================================================================
|
|
== DATASET 3 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 3, }
|
|
dimensions: {natom: 4, nkpt: 10, mband: 26, nsppol: 2, nspinor: 1, nspden: 2, mpw: 376, }
|
|
cutoff_energies: {ecut: 10.0, pawecutdg: 20.0, }
|
|
electrons: {nelect: 4.80000000E+01, charge: 0.00000000E+00, occopt: 0.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename: getwfk from: ./tlruj_1.o_WFK
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 3.9400000 3.9400000 7.8800000 G(1)= -0.0634518 -0.0634518 0.1903553
|
|
R(2)= 7.8800000 3.9400000 3.9400000 G(2)= 0.1903553 -0.0634518 -0.0634518
|
|
R(3)= 3.9400000 7.8800000 3.9400000 G(3)= -0.0634518 0.1903553 -0.0634518
|
|
Unit cell volume ucvol= 2.4465194E+02 bohr^3
|
|
Angles (23,13,12)= 3.35573098E+01 3.35573098E+01 3.35573098E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 30 30 30
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.18911
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 40 40 40
|
|
ecut(hartree)= 20.000 => boxcut(ratio)= 2.06006
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file ./tlruj_1.o_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 368.875 368.830
|
|
|
|
|
|
******************************************
|
|
DFT+U Method used: FLL
|
|
******************************************
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 3, }
|
|
solver: {iscf: 17, nstep: 30, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-09, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -363.72671779804 -3.637E+02 1.318E-08 4.640E-02
|
|
ETOT 2 -363.72727648043 -5.587E-04 1.276E-07 1.331E+00
|
|
ETOT 3 -363.72668375830 5.927E-04 4.331E-05 4.581E-04
|
|
ETOT 4 -363.72668525412 -1.496E-06 8.870E-08 3.588E-05
|
|
ETOT 5 -363.72668532504 -7.092E-08 1.599E-08 1.232E-04
|
|
ETOT 6 -363.72668529166 3.337E-08 1.532E-09 3.746E-05
|
|
ETOT 7 -363.72668528522 6.450E-09 1.557E-09 9.362E-06
|
|
ETOT 8 -363.72668528953 -4.312E-09 1.007E-10 6.943E-06
|
|
ETOT 9 -363.72668529286 -3.335E-09 4.035E-10 3.772E-06
|
|
ETOT 10 -363.72668529549 -2.626E-09 2.279E-10 1.805E-06
|
|
ETOT 11 -363.72668529606 -5.758E-10 8.386E-10 1.256E-07
|
|
ETOT 12 -363.72668529729 -1.225E-09 8.607E-11 2.855E-08
|
|
ETOT 13 -363.72668529776 -4.758E-10 1.252E-10 2.086E-08
|
|
ETOT 14 -363.72668529788 -1.151E-10 1.226E-11 1.655E-09
|
|
ETOT 15 -363.72668529791 -2.569E-11 1.272E-11 3.676E-10
|
|
|
|
At SCF step 15 nres2 = 3.68E-10 < tolvrs= 1.00E-09 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42475182E-01 sigma(3 2)= -4.87136612E-04
|
|
sigma(2 2)= 1.42475182E-01 sigma(3 1)= -4.87136612E-04
|
|
sigma(3 3)= 1.42475182E-01 sigma(2 1)= -4.87136612E-04
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 3, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 3.9400000, 3.9400000, 7.8800000, ]
|
|
- [ 7.8800000, 3.9400000, 3.9400000, ]
|
|
- [ 3.9400000, 7.8800000, 3.9400000, ]
|
|
lattice_lengths: [ 9.65099, 9.65099, 9.65099, ]
|
|
lattice_angles: [ 33.557, 33.557, 33.557, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.4465194E+02
|
|
convergence: {deltae: -2.569E-11, res2: 3.676E-10, residm: 1.272E-11, diffor: null, }
|
|
etotal : -3.63726685E+02
|
|
entropy : 0.00000000E+00
|
|
fermie : 4.07678431E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 1.42475182E-01, -4.87136612E-04, -4.87136612E-04, ]
|
|
- [ -4.87136612E-04, 1.42475182E-01, -4.87136612E-04, ]
|
|
- [ -4.87136612E-04, -4.87136612E-04, 1.42475182E-01, ]
|
|
pressure_GPa: -4.1918E+03
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Ni]
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ni]
|
|
- [ 2.5000E-01, 2.5000E-01, 2.5000E-01, O]
|
|
- [ 7.5000E-01, 7.5000E-01, 7.5000E-01, O]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ 2.03287907E-20, -3.38813179E-20, 6.77626358E-21, ]
|
|
- [ 2.03287907E-20, -3.38813179E-20, 6.77626358E-21, ]
|
|
- [ -1.30232881E-04, -1.30232881E-04, -1.30232881E-04, ]
|
|
- [ 1.30232881E-04, 1.30232881E-04, 1.30232881E-04, ]
|
|
force_length_stats: {min: 4.00889160E-20, max: 2.25569967E-04, mean: 1.12784984E-04, }
|
|
...
|
|
|
|
Integrated electronic and magnetization densities in atomic spheres:
|
|
---------------------------------------------------------------------
|
|
Radius=ratsph(iatom), smearing ratsm= 0.0000. Diff(up-dn)=approximate z local magnetic moment.
|
|
Atom Radius up_density dn_density Total(up+dn) Diff(up-dn)
|
|
1 1.81432 8.024974 7.868389 15.893363 0.156585
|
|
2 1.81432 7.866410 8.023142 15.889551 -0.156732
|
|
3 1.41465 2.292268 2.292194 4.584462 0.000074
|
|
4 1.41465 2.292268 2.292194 4.584462 0.000074
|
|
---------------------------------------------------------------------
|
|
Sum: 20.475920 20.475919 40.951839 0.000001
|
|
Total magnetization (from the atomic spheres): 0.000001
|
|
Total magnetization (exact up - dn): -0.000000
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 12.719705992644119
|
|
Compensation charge over fine fft grid = 12.719993068917242
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1 - Spin component 1
|
|
-2.33054 0.37546 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00170 -0.00170 0.00000 -0.00170 ...
|
|
0.37546 0.04446 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00034 0.00034 0.00000 0.00034 ...
|
|
0.00000 0.00000 -2.24190 -0.00126 -0.00126 0.82140 0.00077 0.00077 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00126 -2.24190 -0.00126 0.00077 0.82140 0.00077 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00126 -0.00126 -2.24190 0.00077 0.00077 0.82140 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82140 0.00077 0.00077 0.37067 -0.00046 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00077 0.82140 0.00077 -0.00046 0.37067 -0.00046 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00077 0.00077 0.82140 -0.00046 -0.00046 0.37067 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91232 -0.00101 0.00094 -0.00101 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.91232 -0.00047 -0.00101 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00094 -0.00047 -0.91360 -0.00047 ...
|
|
-0.00170 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00101 -0.00101 -0.00047 -0.91232 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
-2.32675 0.37475 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00223 -0.00223 0.00000 -0.00223 ...
|
|
0.37475 0.04457 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00044 0.00044 0.00000 0.00044 ...
|
|
0.00000 0.00000 -2.23873 -0.00165 -0.00165 0.81962 0.00099 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -2.23873 -0.00165 0.00099 0.81962 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -0.00165 -2.23873 0.00099 0.00099 0.81962 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.81962 0.00099 0.00099 0.37173 -0.00060 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.81962 0.00099 -0.00060 0.37173 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.00099 0.81962 -0.00060 -0.00060 0.37173 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.90957 -0.00114 0.00137 -0.00114 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.90957 -0.00068 -0.00114 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00137 -0.00068 -0.91066 -0.00068 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.00114 -0.00068 -0.90957 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
0.61340 -1.75562 0.00001 0.00001 0.00001 -0.00007 -0.00007 -0.00007
|
|
-1.75562 5.08527 -0.00003 -0.00003 -0.00003 0.00017 0.00017 0.00017
|
|
0.00001 -0.00003 -0.29402 -0.00002 -0.00002 0.88843 0.00011 0.00011
|
|
0.00001 -0.00003 -0.00002 -0.29402 -0.00002 0.00011 0.88843 0.00011
|
|
0.00001 -0.00003 -0.00002 -0.00002 -0.29402 0.00011 0.00011 0.88843
|
|
-0.00007 0.00017 0.88843 0.00011 0.00011 0.11845 -0.00056 -0.00056
|
|
-0.00007 0.00017 0.00011 0.88843 0.00011 -0.00056 0.11845 -0.00056
|
|
-0.00007 0.00017 0.00011 0.00011 0.88843 -0.00056 -0.00056 0.11845
|
|
Atom # 4 - Spin component 2
|
|
0.61340 -1.75561 -0.00001 -0.00001 -0.00001 0.00009 0.00009 0.00009
|
|
-1.75561 5.08525 0.00004 0.00004 0.00004 -0.00024 -0.00024 -0.00024
|
|
-0.00001 0.00004 -0.29402 -0.00002 -0.00002 0.88842 0.00011 0.00011
|
|
-0.00001 0.00004 -0.00002 -0.29402 -0.00002 0.00011 0.88842 0.00011
|
|
-0.00001 0.00004 -0.00002 -0.00002 -0.29402 0.00011 0.00011 0.88842
|
|
0.00009 -0.00024 0.88842 0.00011 0.00011 0.11849 -0.00056 -0.00056
|
|
0.00009 -0.00024 0.00011 0.88842 0.00011 -0.00056 0.11849 -0.00056
|
|
0.00009 -0.00024 0.00011 0.00011 0.88842 -0.00056 -0.00056 0.11849
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1 - Spin component 1
|
|
0.97694 -0.16681 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00014 -0.00014 0.00000 -0.00014 ...
|
|
-0.16681 0.45824 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00358 -0.00358 0.00000 -0.00358 ...
|
|
0.00000 0.00000 1.06704 0.00002 0.00002 0.13577 0.00008 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 1.06704 0.00002 0.00008 0.13577 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 0.00002 1.06704 0.00008 0.00008 0.13577 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13577 0.00008 0.00008 0.02463 0.00002 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.13577 0.00008 0.00002 0.02463 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.00008 0.13577 0.00002 0.00002 0.02463 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00014 -0.00358 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.15501 -0.02174 0.02664 -0.02174 ...
|
|
-0.00014 -0.00358 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02174 1.15501 -0.01332 -0.02174 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.02664 -0.01332 1.05965 -0.01332 ...
|
|
-0.00014 -0.00358 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02174 -0.02174 -0.01332 1.15501 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
0.97687 -0.16733 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00013 -0.00013 0.00000 -0.00013 ...
|
|
-0.16733 0.45666 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00128 -0.00128 0.00000 -0.00128 ...
|
|
0.00000 0.00000 1.06701 0.00003 0.00003 0.13566 0.00012 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 1.06701 0.00003 0.00012 0.13566 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 0.00003 1.06701 0.00012 0.00012 0.13566 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13566 0.00012 0.00012 0.02455 0.00004 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.13566 0.00012 0.00004 0.02455 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.00012 0.13566 0.00004 0.00004 0.02455 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00013 -0.00128 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.11070 -0.01390 0.00067 -0.01390 ...
|
|
-0.00013 -0.00128 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01390 1.11070 -0.00033 -0.01390 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00067 -0.00033 1.01216 -0.00033 ...
|
|
-0.00013 -0.00128 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01390 -0.01390 -0.00033 1.11070 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
1.01157 0.01395 -0.00077 -0.00077 -0.00077 0.00003 0.00003 0.00003
|
|
0.01395 0.00119 0.00013 0.00013 0.00013 -0.00000 -0.00000 -0.00000
|
|
-0.00077 0.00013 1.06827 -0.00321 -0.00321 0.05052 -0.00033 -0.00033
|
|
-0.00077 0.00013 -0.00321 1.06827 -0.00321 -0.00033 0.05052 -0.00033
|
|
-0.00077 0.00013 -0.00321 -0.00321 1.06827 -0.00033 -0.00033 0.05052
|
|
0.00003 -0.00000 0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
0.00003 -0.00000 -0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
0.00003 -0.00000 -0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - Spin component 2
|
|
1.01155 0.01396 0.00067 0.00067 0.00067 -0.00003 -0.00003 -0.00003
|
|
0.01396 0.00119 -0.00011 -0.00011 -0.00011 0.00000 0.00000 0.00000
|
|
0.00067 -0.00011 1.06824 -0.00320 -0.00320 0.05052 -0.00033 -0.00033
|
|
0.00067 -0.00011 -0.00320 1.06824 -0.00320 -0.00033 0.05052 -0.00033
|
|
0.00067 -0.00011 -0.00320 -0.00320 1.06824 -0.00033 -0.00033 0.05052
|
|
-0.00003 0.00000 0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00003 0.00000 -0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00003 0.00000 -0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
|
|
"PAW+U" part of augmentation waves occupancies Rhoij:
|
|
Atom # 1 - L=2 ONLY - Spin component 1
|
|
1.15501 -0.02174 0.02664 -0.02174 0.00000 0.42939 -0.00883 0.01033 -0.00883 0.00000
|
|
-0.02174 1.15501 -0.01332 -0.02174 0.02307 -0.00883 0.42939 -0.00516 -0.00883 0.00894
|
|
0.02664 -0.01332 1.05965 -0.01332 0.00000 0.01204 -0.00602 0.39188 -0.00602 0.00000
|
|
-0.02174 -0.02174 -0.01332 1.15501 -0.02307 -0.00883 -0.00883 -0.00516 0.42939 -0.00894
|
|
0.00000 0.02307 0.00000 -0.02307 1.05965 0.00000 0.01042 0.00000 -0.01042 0.39188
|
|
0.42939 -0.00883 0.01204 -0.00883 0.00000 0.16420 -0.00362 0.00472 -0.00362 0.00000
|
|
-0.00883 0.42939 -0.00602 -0.00883 0.01042 -0.00362 0.16420 -0.00236 -0.00362 0.00409
|
|
0.01033 -0.00516 0.39188 -0.00516 0.00000 0.00472 -0.00236 0.15191 -0.00236 0.00000
|
|
-0.00883 -0.00883 -0.00602 0.42939 -0.01042 -0.00362 -0.00362 -0.00236 0.16420 -0.00409
|
|
0.00000 0.00894 0.00000 -0.00894 0.39188 0.00000 0.00409 0.00000 -0.00409 0.15191
|
|
Atom # 1 - L=2 ONLY - Spin component 2
|
|
1.11070 -0.01390 0.00067 -0.01390 0.00000 0.40914 -0.00545 -0.00019 -0.00545 0.00000
|
|
-0.01390 1.11070 -0.00033 -0.01390 0.00058 -0.00545 0.40914 0.00010 -0.00545 -0.00017
|
|
0.00067 -0.00033 1.01216 -0.00033 0.00000 -0.00039 0.00020 0.37384 0.00020 0.00000
|
|
-0.01390 -0.01390 -0.00033 1.11070 -0.00058 -0.00545 -0.00545 0.00010 0.40914 0.00017
|
|
0.00000 0.00058 0.00000 -0.00058 1.01216 0.00000 -0.00034 0.00000 0.00034 0.37384
|
|
0.40914 -0.00545 -0.00039 -0.00545 0.00000 0.15494 -0.00217 -0.00030 -0.00217 0.00000
|
|
-0.00545 0.40914 0.00020 -0.00545 -0.00034 -0.00217 0.15494 0.00015 -0.00217 -0.00026
|
|
-0.00019 0.00010 0.37384 0.00010 0.00000 -0.00030 0.00015 0.14502 0.00015 0.00000
|
|
-0.00545 -0.00545 0.00020 0.40914 0.00034 -0.00217 -0.00217 0.00015 0.15494 0.00026
|
|
0.00000 -0.00017 0.00000 0.00017 0.37384 0.00000 -0.00026 0.00000 0.00026 0.14502
|
|
Atom # 2 - L=2 ONLY - Spin component 1
|
|
1.11008 -0.01391 0.00060 -0.01391 0.00000 0.40882 -0.00546 -0.00023 -0.00546 0.00000
|
|
-0.01391 1.11008 -0.00030 -0.01391 0.00052 -0.00546 0.40882 0.00012 -0.00546 -0.00020
|
|
0.00060 -0.00030 1.01125 -0.00030 0.00000 -0.00045 0.00022 0.37346 0.00022 0.00000
|
|
-0.01391 -0.01391 -0.00030 1.11008 -0.00052 -0.00546 -0.00546 0.00012 0.40882 0.00020
|
|
0.00000 0.00052 0.00000 -0.00052 1.01125 0.00000 -0.00039 0.00000 0.00039 0.37346
|
|
0.40882 -0.00546 -0.00045 -0.00546 0.00000 0.15479 -0.00217 -0.00032 -0.00217 0.00000
|
|
-0.00546 0.40882 0.00022 -0.00546 -0.00039 -0.00217 0.15479 0.00016 -0.00217 -0.00028
|
|
-0.00023 0.00012 0.37346 0.00012 0.00000 -0.00032 0.00016 0.14487 0.00016 0.00000
|
|
-0.00546 -0.00546 0.00022 0.40882 0.00039 -0.00217 -0.00217 0.00016 0.15479 0.00028
|
|
0.00000 -0.00020 0.00000 0.00020 0.37346 0.00000 -0.00028 0.00000 0.00028 0.14487
|
|
Atom # 2 - L=2 ONLY - Spin component 2
|
|
1.15450 -0.02174 0.02658 -0.02174 0.00000 0.42913 -0.00883 0.01029 -0.00883 0.00000
|
|
-0.02174 1.15450 -0.01329 -0.02174 0.02302 -0.00883 0.42913 -0.00514 -0.00883 0.00891
|
|
0.02658 -0.01329 1.05871 -0.01329 0.00000 0.01199 -0.00600 0.39148 -0.00600 0.00000
|
|
-0.02174 -0.02174 -0.01329 1.15450 -0.02302 -0.00883 -0.00883 -0.00514 0.42913 -0.00891
|
|
0.00000 0.02302 0.00000 -0.02302 1.05871 0.00000 0.01038 0.00000 -0.01038 0.39148
|
|
0.42913 -0.00883 0.01199 -0.00883 0.00000 0.16409 -0.00362 0.00470 -0.00362 0.00000
|
|
-0.00883 0.42913 -0.00600 -0.00883 0.01038 -0.00362 0.16409 -0.00235 -0.00362 0.00407
|
|
0.01029 -0.00514 0.39148 -0.00514 0.00000 0.00470 -0.00235 0.15176 -0.00235 0.00000
|
|
-0.00883 -0.00883 -0.00600 0.42913 -0.01038 -0.00362 -0.00362 -0.00235 0.16409 -0.00407
|
|
0.00000 0.00891 0.00000 -0.00891 0.39148 0.00000 0.00407 0.00000 -0.00407 0.15176
|
|
Atom # 3 - L=1 ONLY - Spin component 1
|
|
1.06827 -0.00321 -0.00321 0.05052 -0.00033 -0.00033
|
|
-0.00321 1.06827 -0.00321 -0.00033 0.05052 -0.00033
|
|
-0.00321 -0.00321 1.06827 -0.00033 -0.00033 0.05052
|
|
0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
Atom # 3 - L=1 ONLY - Spin component 2
|
|
1.06824 -0.00320 -0.00320 0.05052 -0.00033 -0.00033
|
|
-0.00320 1.06824 -0.00320 -0.00033 0.05052 -0.00033
|
|
-0.00320 -0.00320 1.06824 -0.00033 -0.00033 0.05052
|
|
0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 1
|
|
1.06827 -0.00321 -0.00321 0.05052 -0.00033 -0.00033
|
|
-0.00321 1.06827 -0.00321 -0.00033 0.05052 -0.00033
|
|
-0.00321 -0.00321 1.06827 -0.00033 -0.00033 0.05052
|
|
0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 2
|
|
1.06824 -0.00320 -0.00320 0.05052 -0.00033 -0.00033
|
|
-0.00320 1.06824 -0.00320 -0.00033 0.05052 -0.00033
|
|
-0.00320 -0.00320 1.06824 -0.00033 -0.00033 0.05052
|
|
0.05052 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05052 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05052 -0.00001 -0.00001 0.00287
|
|
|
|
---------- DFT+U DATA ---------------------------------------------------
|
|
|
|
====== For Atom 1, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 1. Occ. for lpawu and for spin 1 = 4.40952
|
|
Atom 1. Occ. for lpawu and for spin 2 = 4.23323
|
|
=> On atom 1, local Mag. for lpawu is -0.176283
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.91173 -0.01694 0.02064 -0.01694 0.00000
|
|
-0.01694 0.91173 -0.01032 -0.01694 0.01788
|
|
0.02064 -0.01032 0.83716 -0.01032 0.00000
|
|
-0.01694 -0.01694 -0.01032 0.91173 -0.01788
|
|
0.00000 0.01788 0.00000 -0.01788 0.83716
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.87789 -0.01089 0.00069 -0.01089 0.00000
|
|
-0.01089 0.87789 -0.00034 -0.01089 0.00059
|
|
0.00069 -0.00034 0.79979 -0.00034 0.00000
|
|
-0.01089 -0.01089 -0.00034 0.87789 -0.00059
|
|
0.00000 0.00059 0.00000 -0.00059 0.79979
|
|
|
|
|
|
====== For Atom 2, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 2. Occ. for lpawu and for spin 1 = 4.23043
|
|
Atom 2. Occ. for lpawu and for spin 2 = 4.40690
|
|
=> On atom 2, local Mag. for lpawu is 0.176472
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.87742 -0.01089 0.00065 -0.01089 0.00000
|
|
-0.01089 0.87742 -0.00032 -0.01089 0.00056
|
|
0.00065 -0.00032 0.79909 -0.00032 0.00000
|
|
-0.01089 -0.01089 -0.00032 0.87742 -0.00056
|
|
0.00000 0.00056 0.00000 -0.00056 0.79909
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.91135 -0.01694 0.02060 -0.01694 0.00000
|
|
-0.01694 0.91135 -0.01030 -0.01694 0.01784
|
|
0.02060 -0.01030 0.83643 -0.01030 0.00000
|
|
-0.01694 -0.01694 -0.01030 0.91135 -0.01784
|
|
0.00000 0.01784 0.00000 -0.01784 0.83643
|
|
|
|
|
|
====== For Atom 3, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 3. Occ. for lpawu and for spin 1 = 1.55465
|
|
Atom 3. Occ. for lpawu and for spin 2 = 1.55461
|
|
=> On atom 3, local Mag. for lpawu is -0.000048
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51822 -0.00082 -0.00082
|
|
-0.00082 0.51822 -0.00082
|
|
-0.00082 -0.00082 0.51822
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51820 -0.00081 -0.00081
|
|
-0.00081 0.51820 -0.00081
|
|
-0.00081 -0.00081 0.51820
|
|
|
|
|
|
====== For Atom 4, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 4. Occ. for lpawu and for spin 1 = 1.55465
|
|
Atom 4. Occ. for lpawu and for spin 2 = 1.55461
|
|
=> On atom 4, local Mag. for lpawu is -0.000048
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51822 -0.00082 -0.00082
|
|
-0.00082 0.51822 -0.00082
|
|
-0.00082 -0.00082 0.51822
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51820 -0.00081 -0.00081
|
|
-0.00081 0.51820 -0.00081
|
|
-0.00081 -0.00081 0.51820
|
|
|
|
|
|
|
|
|
|
---------- calculate U, (J) start ----------
|
|
|
|
All pawujat ok and equal to 1
|
|
macro_uj=1 and nspden=2: standard determination of Hubbard U-parameter
|
|
|
|
|
|
*********************************************************************
|
|
************************ Linear Response U ************************
|
|
|
|
Info printed for perturbed atom: 1
|
|
|
|
Perturbations Occupations
|
|
--------------- -----------------------------
|
|
alpha [eV] Unscreened Screened
|
|
--------------- -----------------------------
|
|
0.0000000000 8.6380182458 8.6380182458
|
|
-0.0500000000 8.6560314461 8.6427490201
|
|
|
|
Scalar response functions:
|
|
Chi0 [eV^-1]: -0.80059
|
|
Chi [eV^-1]: -0.09462
|
|
|
|
The scalar U from the two-point regression scheme is 9.32001 eV.
|
|
*********************************************************************
|
|
*********************************************************************
|
|
Note: For more reliable linear regressions of the response
|
|
matrices, it is advised that you have more than two points.
|
|
See the LRUJ protocol for more information.
|
|
|
|
|
|
URES ii nat r_max U(J)[eV] U_ASA[eV] U_inf[eV]
|
|
URES 1 1 0.00000 2.33000 1.87503 1.81378
|
|
URES 2 8 11.14400 8.16939 6.57418 6.35943
|
|
URES 3 27 12.45940 8.97527 7.22270 6.98675
|
|
URES 4 64 22.28800 9.17442 7.38296 7.14179
|
|
URES 5 125 24.28780 9.24546 7.44013 7.19708
|
|
URES 6 216 33.43200 9.27686 7.46540 7.22153
|
|
|
|
---------- calculate U, (J) end --------------
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 10.175E-14; max= 12.724E-12
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
0.250000000000 0.250000000000 0.250000000000
|
|
0.750000000000 0.750000000000 0.750000000000
|
|
rms dE/dt= 1.4513E-03; max dE/dt= 2.0525E-03; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
2 0.000000000000 0.000000000000 0.000000000000
|
|
3 0.002052470207 0.002052470207 0.002052470207
|
|
4 -0.002052470207 -0.002052470207 -0.002052470207
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 0.00000000000000
|
|
2 4.16991640368920 4.16991640368920 4.16991640368920
|
|
3 2.08495820184460 2.08495820184460 2.08495820184460
|
|
4 6.25487460553380 6.25487460553380 6.25487460553380
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
3 -0.00013023288118 -0.00013023288118 -0.00013023288118
|
|
4 0.00013023288118 0.00013023288118 0.00013023288118
|
|
frms,max,avg= 9.2088553E-05 1.3023288E-04 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
3 -0.00669684344569 -0.00669684344569 -0.00669684344569
|
|
4 0.00669684344569 0.00669684344569 0.00669684344569
|
|
frms,max,avg= 4.7353834E-03 6.6968434E-03 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 7.880000000000 7.880000000000 7.880000000000 bohr
|
|
= 4.169916403689 4.169916403689 4.169916403689 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.40768 Average Vxc (hartree)= -0.46905
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN UP:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.13990 -4.13702 -2.37077 -2.33851 -2.33184 -2.33184 -2.30940 -2.30940
|
|
-0.40334 -0.35615 0.07022 0.21096 0.21096 0.22655 0.22655 0.25279
|
|
0.26163 0.26163 0.31910 0.31910 0.37475 0.38738 0.38738 0.39683
|
|
0.40071 0.43654
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN DOWN:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14228 -4.13465 -2.37103 -2.33825 -2.33246 -2.33246 -2.30878 -2.30878
|
|
-0.40333 -0.35614 0.07022 0.21097 0.21097 0.22655 0.22655 0.25279
|
|
0.26163 0.26163 0.31910 0.31910 0.37478 0.38739 0.38739 0.39683
|
|
0.40069 0.43653
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 3, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.94818960907360E+01
|
|
hartree : 8.65592597754744E+01
|
|
xc : -2.18491558914781E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
local_psp : -2.68603774711125E+02
|
|
spherical_terms : -2.03987767341144E+01
|
|
total_energy : -3.63710803965360E+02
|
|
total_energy_eV : -9.89707430073063E+03
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 3, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -4.00290379971530E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
xc_dc : -7.99992731699368E+01
|
|
spherical_terms : -1.47981216359626E+01
|
|
total_energy_dc : -3.63726685297905E+02
|
|
total_energy_dc_eV : -9.89750645376673E+03
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42475182E-01 sigma(3 2)= -4.87136612E-04
|
|
sigma(2 2)= 1.42475182E-01 sigma(3 1)= -4.87136612E-04
|
|
sigma(3 3)= 1.42475182E-01 sigma(2 1)= -4.87136612E-04
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -4.1918E+03 GPa]
|
|
- sigma(1 1)= 4.19176387E+03 sigma(3 2)= -1.43320515E+01
|
|
- sigma(2 2)= 4.19176387E+03 sigma(3 1)= -1.43320515E+01
|
|
- sigma(3 3)= 4.19176387E+03 sigma(2 1)= -1.43320515E+01
|
|
|
|
================================================================================
|
|
== DATASET 4 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 4, }
|
|
dimensions: {natom: 4, nkpt: 10, mband: 26, nsppol: 2, nspinor: 1, nspden: 2, mpw: 376, }
|
|
cutoff_energies: {ecut: 10.0, pawecutdg: 20.0, }
|
|
electrons: {nelect: 4.80000000E+01, charge: 0.00000000E+00, occopt: 0.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename: getwfk from: ./tlruj_1.o_WFK
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 3.9400000 3.9400000 7.8800000 G(1)= -0.0634518 -0.0634518 0.1903553
|
|
R(2)= 7.8800000 3.9400000 3.9400000 G(2)= 0.1903553 -0.0634518 -0.0634518
|
|
R(3)= 3.9400000 7.8800000 3.9400000 G(3)= -0.0634518 0.1903553 -0.0634518
|
|
Unit cell volume ucvol= 2.4465194E+02 bohr^3
|
|
Angles (23,13,12)= 3.35573098E+01 3.35573098E+01 3.35573098E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 30 30 30
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.18911
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 40 40 40
|
|
ecut(hartree)= 20.000 => boxcut(ratio)= 2.06006
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file ./tlruj_1.o_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 368.875 368.830
|
|
|
|
|
|
******************************************
|
|
DFT+U Method used: FLL
|
|
******************************************
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 4, }
|
|
solver: {iscf: 17, nstep: 30, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-09, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -363.69497467037 -3.637E+02 1.422E-08 5.078E-02
|
|
ETOT 2 -363.69557245537 -5.978E-04 1.316E-07 1.168E+00
|
|
ETOT 3 -363.69494015885 6.323E-04 4.555E-05 1.072E-03
|
|
ETOT 4 -363.69494111039 -9.515E-07 9.840E-08 3.917E-05
|
|
ETOT 5 -363.69494117247 -6.208E-08 3.072E-08 6.116E-05
|
|
ETOT 6 -363.69494118720 -1.473E-08 1.710E-09 3.891E-05
|
|
ETOT 7 -363.69494117653 1.067E-08 2.341E-09 7.586E-06
|
|
ETOT 8 -363.69494117635 1.818E-10 3.748E-11 7.713E-06
|
|
ETOT 9 -363.69494118033 -3.976E-09 9.942E-11 8.342E-06
|
|
ETOT 10 -363.69494119081 -1.048E-08 9.396E-10 3.229E-06
|
|
ETOT 11 -363.69494119855 -7.740E-09 1.473E-09 7.514E-07
|
|
ETOT 12 -363.69494120048 -1.934E-09 1.152E-10 1.973E-07
|
|
ETOT 13 -363.69494120328 -2.797E-09 4.707E-10 1.220E-07
|
|
ETOT 14 -363.69494120364 -3.582E-10 3.796E-11 9.369E-08
|
|
ETOT 15 -363.69494120392 -2.811E-10 2.794E-11 1.453E-08
|
|
ETOT 16 -363.69494120404 -1.172E-10 1.072E-11 9.965E-09
|
|
ETOT 17 -363.69494120409 -5.002E-11 5.391E-12 3.755E-10
|
|
|
|
At SCF step 17 nres2 = 3.76E-10 < tolvrs= 1.00E-09 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42464714E-01 sigma(3 2)= -4.89478733E-04
|
|
sigma(2 2)= 1.42464714E-01 sigma(3 1)= -4.89478733E-04
|
|
sigma(3 3)= 1.42464714E-01 sigma(2 1)= -4.89478733E-04
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 4, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 3.9400000, 3.9400000, 7.8800000, ]
|
|
- [ 7.8800000, 3.9400000, 3.9400000, ]
|
|
- [ 3.9400000, 7.8800000, 3.9400000, ]
|
|
lattice_lengths: [ 9.65099, 9.65099, 9.65099, ]
|
|
lattice_angles: [ 33.557, 33.557, 33.557, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.4465194E+02
|
|
convergence: {deltae: -5.002E-11, res2: 3.755E-10, residm: 5.391E-12, diffor: null, }
|
|
etotal : -3.63694941E+02
|
|
entropy : 0.00000000E+00
|
|
fermie : 4.08504528E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 1.42464714E-01, -4.89478733E-04, -4.89478733E-04, ]
|
|
- [ -4.89478733E-04, 1.42464714E-01, -4.89478733E-04, ]
|
|
- [ -4.89478733E-04, -4.89478733E-04, 1.42464714E-01, ]
|
|
pressure_GPa: -4.1915E+03
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Ni]
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ni]
|
|
- [ 2.5000E-01, 2.5000E-01, 2.5000E-01, O]
|
|
- [ 7.5000E-01, 7.5000E-01, 7.5000E-01, O]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ 6.77626358E-21, -2.03287907E-20, 6.77626358E-21, ]
|
|
- [ 6.77626358E-21, -2.03287907E-20, 6.77626358E-21, ]
|
|
- [ 1.32502130E-04, 1.32502130E-04, 1.32502130E-04, ]
|
|
- [ -1.32502130E-04, -1.32502130E-04, -1.32502130E-04, ]
|
|
force_length_stats: {min: 2.24743238E-20, max: 2.29500421E-04, mean: 1.14750211E-04, }
|
|
...
|
|
|
|
Integrated electronic and magnetization densities in atomic spheres:
|
|
---------------------------------------------------------------------
|
|
Radius=ratsph(iatom), smearing ratsm= 0.0000. Diff(up-dn)=approximate z local magnetic moment.
|
|
Atom Radius up_density dn_density Total(up+dn) Diff(up-dn)
|
|
1 1.81432 8.022336 7.864412 15.886747 0.157924
|
|
2 1.81432 7.866440 8.024121 15.890561 -0.157681
|
|
3 1.41465 2.292330 2.292422 4.584752 -0.000092
|
|
4 1.41465 2.292330 2.292422 4.584752 -0.000092
|
|
---------------------------------------------------------------------
|
|
Sum: 20.473436 20.473377 40.946813 0.000059
|
|
Total magnetization (from the atomic spheres): 0.000059
|
|
Total magnetization (exact up - dn): 0.000000
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 12.717893955413377
|
|
Compensation charge over fine fft grid = 12.718180256555661
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1 - Spin component 1
|
|
-2.33247 0.37596 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00171 -0.00171 0.00000 -0.00171 ...
|
|
0.37596 0.04437 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00034 0.00034 0.00000 0.00034 ...
|
|
0.00000 0.00000 -2.24365 -0.00127 -0.00127 0.82275 0.00077 0.00077 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00127 -2.24365 -0.00127 0.00077 0.82275 0.00077 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00127 -0.00127 -2.24365 0.00077 0.00077 0.82275 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82275 0.00077 0.00077 0.36979 -0.00047 -0.00047 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00077 0.82275 0.00077 -0.00047 0.36979 -0.00047 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00077 0.00077 0.82275 -0.00047 -0.00047 0.36979 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00171 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91101 -0.00102 0.00094 -0.00102 ...
|
|
-0.00171 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00102 -0.91101 -0.00047 -0.00102 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00094 -0.00047 -0.91229 -0.00047 ...
|
|
-0.00171 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00102 -0.00102 -0.00047 -0.91101 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
-2.32864 0.37525 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00223 -0.00223 0.00000 -0.00223 ...
|
|
0.37525 0.04448 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00044 0.00044 0.00000 0.00044 ...
|
|
0.00000 0.00000 -2.24045 -0.00165 -0.00165 0.82095 0.00099 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -2.24045 -0.00165 0.00099 0.82095 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -0.00165 -2.24045 0.00099 0.00099 0.82095 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82095 0.00099 0.00099 0.37086 -0.00060 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.82095 0.00099 -0.00060 0.37086 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.00099 0.82095 -0.00060 -0.00060 0.37086 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.90823 -0.00114 0.00137 -0.00114 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.90823 -0.00069 -0.00114 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00137 -0.00069 -0.90933 -0.00069 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.00114 -0.00069 -0.90823 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
0.61336 -1.75550 0.00002 0.00002 0.00002 -0.00010 -0.00010 -0.00010
|
|
-1.75550 5.08489 -0.00004 -0.00004 -0.00004 0.00025 0.00025 0.00025
|
|
0.00002 -0.00004 -0.29395 -0.00002 -0.00002 0.88808 0.00011 0.00011
|
|
0.00002 -0.00004 -0.00002 -0.29395 -0.00002 0.00011 0.88808 0.00011
|
|
0.00002 -0.00004 -0.00002 -0.00002 -0.29395 0.00011 0.00011 0.88808
|
|
-0.00010 0.00025 0.88808 0.00011 0.00011 0.12028 -0.00056 -0.00056
|
|
-0.00010 0.00025 0.00011 0.88808 0.00011 -0.00056 0.12028 -0.00056
|
|
-0.00010 0.00025 0.00011 0.00011 0.88808 -0.00056 -0.00056 0.12028
|
|
Atom # 4 - Spin component 2
|
|
0.61337 -1.75551 -0.00001 -0.00001 -0.00001 0.00007 0.00007 0.00007
|
|
-1.75551 5.08492 0.00003 0.00003 0.00003 -0.00019 -0.00019 -0.00019
|
|
-0.00001 0.00003 -0.29395 -0.00002 -0.00002 0.88809 0.00011 0.00011
|
|
-0.00001 0.00003 -0.00002 -0.29395 -0.00002 0.00011 0.88809 0.00011
|
|
-0.00001 0.00003 -0.00002 -0.00002 -0.29395 0.00011 0.00011 0.88809
|
|
0.00007 -0.00019 0.88809 0.00011 0.00011 0.12023 -0.00055 -0.00055
|
|
0.00007 -0.00019 0.00011 0.88809 0.00011 -0.00055 0.12023 -0.00055
|
|
0.00007 -0.00019 0.00011 0.00011 0.88809 -0.00055 -0.00055 0.12023
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1 - Spin component 1
|
|
0.97690 -0.16706 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00016 -0.00016 0.00000 -0.00016 ...
|
|
-0.16706 0.46117 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00409 -0.00409 0.00000 -0.00409 ...
|
|
0.00000 0.00000 1.06706 0.00002 0.00002 0.13566 0.00008 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 1.06706 0.00002 0.00008 0.13566 0.00008 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 0.00002 1.06706 0.00008 0.00008 0.13566 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13566 0.00008 0.00008 0.02463 0.00002 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.13566 0.00008 0.00002 0.02463 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00008 0.00008 0.13566 0.00002 0.00002 0.02463 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00016 -0.00409 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.15418 -0.02175 0.02666 -0.02175 ...
|
|
-0.00016 -0.00409 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02175 1.15418 -0.01333 -0.02175 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.02666 -0.01333 1.05835 -0.01333 ...
|
|
-0.00016 -0.00409 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02175 -0.02175 -0.01333 1.15418 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
0.97682 -0.16754 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00011 -0.00011 0.00000 -0.00011 ...
|
|
-0.16754 0.46102 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00055 -0.00055 0.00000 -0.00055 ...
|
|
0.00000 0.00000 1.06703 0.00003 0.00003 0.13555 0.00012 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 1.06703 0.00003 0.00012 0.13555 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 0.00003 1.06703 0.00012 0.00012 0.13555 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13555 0.00012 0.00012 0.02455 0.00004 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.13555 0.00012 0.00004 0.02455 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.00012 0.13555 0.00004 0.00004 0.02455 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00011 -0.00055 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.10934 -0.01422 0.00057 -0.01422 ...
|
|
-0.00011 -0.00055 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01422 1.10934 -0.00029 -0.01422 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00057 -0.00029 1.01054 -0.00029 ...
|
|
-0.00011 -0.00055 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01422 -0.01422 -0.00029 1.10934 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
1.01160 0.01395 -0.00073 -0.00073 -0.00073 0.00003 0.00003 0.00003
|
|
0.01395 0.00120 0.00012 0.00012 0.00012 -0.00000 -0.00000 -0.00000
|
|
-0.00073 0.00012 1.06834 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
-0.00073 0.00012 -0.00318 1.06834 -0.00318 -0.00033 0.05053 -0.00033
|
|
-0.00073 0.00012 -0.00318 -0.00318 1.06834 -0.00033 -0.00033 0.05053
|
|
0.00003 -0.00000 0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
0.00003 -0.00000 -0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
0.00003 -0.00000 -0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - Spin component 2
|
|
1.01162 0.01395 0.00083 0.00083 0.00083 -0.00003 -0.00003 -0.00003
|
|
0.01395 0.00120 -0.00015 -0.00015 -0.00015 -0.00000 -0.00000 -0.00000
|
|
0.00083 -0.00015 1.06837 -0.00319 -0.00319 0.05053 -0.00033 -0.00033
|
|
0.00083 -0.00015 -0.00319 1.06837 -0.00319 -0.00033 0.05053 -0.00033
|
|
0.00083 -0.00015 -0.00319 -0.00319 1.06837 -0.00033 -0.00033 0.05053
|
|
-0.00003 -0.00000 0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00003 -0.00000 -0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00003 -0.00000 -0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
|
|
"PAW+U" part of augmentation waves occupancies Rhoij:
|
|
Atom # 1 - L=2 ONLY - Spin component 1
|
|
1.15418 -0.02175 0.02666 -0.02175 0.00000 0.42901 -0.00884 0.01032 -0.00884 0.00000
|
|
-0.02175 1.15418 -0.01333 -0.02175 0.02309 -0.00884 0.42901 -0.00516 -0.00884 0.00894
|
|
0.02666 -0.01333 1.05835 -0.01333 0.00000 0.01204 -0.00602 0.39131 -0.00602 0.00000
|
|
-0.02175 -0.02175 -0.01333 1.15418 -0.02309 -0.00884 -0.00884 -0.00516 0.42901 -0.00894
|
|
0.00000 0.02309 0.00000 -0.02309 1.05835 0.00000 0.01042 0.00000 -0.01042 0.39131
|
|
0.42901 -0.00884 0.01204 -0.00884 0.00000 0.16405 -0.00362 0.00472 -0.00362 0.00000
|
|
-0.00884 0.42901 -0.00602 -0.00884 0.01042 -0.00362 0.16405 -0.00236 -0.00362 0.00408
|
|
0.01032 -0.00516 0.39131 -0.00516 0.00000 0.00472 -0.00236 0.15168 -0.00236 0.00000
|
|
-0.00884 -0.00884 -0.00602 0.42901 -0.01042 -0.00362 -0.00362 -0.00236 0.16405 -0.00408
|
|
0.00000 0.00894 0.00000 -0.00894 0.39131 0.00000 0.00408 0.00000 -0.00408 0.15168
|
|
Atom # 1 - L=2 ONLY - Spin component 2
|
|
1.10934 -0.01422 0.00057 -0.01422 0.00000 0.40850 -0.00560 -0.00025 -0.00560 0.00000
|
|
-0.01422 1.10934 -0.00029 -0.01422 0.00050 -0.00560 0.40850 0.00013 -0.00560 -0.00022
|
|
0.00057 -0.00029 1.01054 -0.00029 0.00000 -0.00046 0.00023 0.37315 0.00023 0.00000
|
|
-0.01422 -0.01422 -0.00029 1.10934 -0.00050 -0.00560 -0.00560 0.00013 0.40850 0.00022
|
|
0.00000 0.00050 0.00000 -0.00050 1.01054 0.00000 -0.00040 0.00000 0.00040 0.37315
|
|
0.40850 -0.00560 -0.00046 -0.00560 0.00000 0.15465 -0.00223 -0.00033 -0.00223 0.00000
|
|
-0.00560 0.40850 0.00023 -0.00560 -0.00040 -0.00223 0.15465 0.00017 -0.00223 -0.00029
|
|
-0.00025 0.00013 0.37315 0.00013 0.00000 -0.00033 0.00017 0.14474 0.00017 0.00000
|
|
-0.00560 -0.00560 0.00023 0.40850 0.00040 -0.00223 -0.00223 0.00017 0.15465 0.00029
|
|
0.00000 -0.00022 0.00000 0.00022 0.37315 0.00000 -0.00029 0.00000 0.00029 0.14474
|
|
Atom # 2 - L=2 ONLY - Spin component 1
|
|
1.11000 -0.01416 0.00063 -0.01416 0.00000 0.40884 -0.00557 -0.00021 -0.00557 0.00000
|
|
-0.01416 1.11000 -0.00032 -0.01416 0.00055 -0.00557 0.40884 0.00011 -0.00557 -0.00019
|
|
0.00063 -0.00032 1.01144 -0.00032 0.00000 -0.00041 0.00021 0.37353 0.00021 0.00000
|
|
-0.01416 -0.01416 -0.00032 1.11000 -0.00055 -0.00557 -0.00557 0.00011 0.40884 0.00019
|
|
0.00000 0.00055 0.00000 -0.00055 1.01144 0.00000 -0.00036 0.00000 0.00036 0.37353
|
|
0.40884 -0.00557 -0.00041 -0.00557 0.00000 0.15481 -0.00222 -0.00031 -0.00222 0.00000
|
|
-0.00557 0.40884 0.00021 -0.00557 -0.00036 -0.00222 0.15481 0.00015 -0.00222 -0.00027
|
|
-0.00021 0.00011 0.37353 0.00011 0.00000 -0.00031 0.00015 0.14489 0.00015 0.00000
|
|
-0.00557 -0.00557 0.00021 0.40884 0.00036 -0.00222 -0.00222 0.00015 0.15481 0.00027
|
|
0.00000 -0.00019 0.00000 0.00019 0.37353 0.00000 -0.00027 0.00000 0.00027 0.14489
|
|
Atom # 2 - L=2 ONLY - Spin component 2
|
|
1.15467 -0.02177 0.02672 -0.02177 0.00000 0.42925 -0.00885 0.01036 -0.00885 0.00000
|
|
-0.02177 1.15467 -0.01336 -0.02177 0.02314 -0.00885 0.42925 -0.00518 -0.00885 0.00897
|
|
0.02672 -0.01336 1.05928 -0.01336 0.00000 0.01208 -0.00604 0.39170 -0.00604 0.00000
|
|
-0.02177 -0.02177 -0.01336 1.15467 -0.02314 -0.00885 -0.00885 -0.00518 0.42925 -0.00897
|
|
0.00000 0.02314 0.00000 -0.02314 1.05928 0.00000 0.01046 0.00000 -0.01046 0.39170
|
|
0.42925 -0.00885 0.01208 -0.00885 0.00000 0.16415 -0.00363 0.00474 -0.00363 0.00000
|
|
-0.00885 0.42925 -0.00604 -0.00885 0.01046 -0.00363 0.16415 -0.00237 -0.00363 0.00410
|
|
0.01036 -0.00518 0.39170 -0.00518 0.00000 0.00474 -0.00237 0.15184 -0.00237 0.00000
|
|
-0.00885 -0.00885 -0.00604 0.42925 -0.01046 -0.00363 -0.00363 -0.00237 0.16415 -0.00410
|
|
0.00000 0.00897 0.00000 -0.00897 0.39170 0.00000 0.00410 0.00000 -0.00410 0.15184
|
|
Atom # 3 - L=1 ONLY - Spin component 1
|
|
1.06834 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
-0.00318 1.06834 -0.00318 -0.00033 0.05053 -0.00033
|
|
-0.00318 -0.00318 1.06834 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 3 - L=1 ONLY - Spin component 2
|
|
1.06837 -0.00319 -0.00319 0.05053 -0.00033 -0.00033
|
|
-0.00319 1.06837 -0.00319 -0.00033 0.05053 -0.00033
|
|
-0.00319 -0.00319 1.06837 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 1
|
|
1.06834 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
-0.00318 1.06834 -0.00318 -0.00033 0.05053 -0.00033
|
|
-0.00318 -0.00318 1.06834 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 2
|
|
1.06837 -0.00319 -0.00319 0.05053 -0.00033 -0.00033
|
|
-0.00319 1.06837 -0.00319 -0.00033 0.05053 -0.00033
|
|
-0.00319 -0.00319 1.06837 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
|
|
---------- DFT+U DATA ---------------------------------------------------
|
|
|
|
====== For Atom 1, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 1. Occ. for lpawu and for spin 1 = 4.40560
|
|
Atom 1. Occ. for lpawu and for spin 2 = 4.22762
|
|
=> On atom 1, local Mag. for lpawu is -0.177977
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.91110 -0.01694 0.02066 -0.01694 0.00000
|
|
-0.01694 0.91110 -0.01033 -0.01694 0.01789
|
|
0.02066 -0.01033 0.83615 -0.01033 0.00000
|
|
-0.01694 -0.01694 -0.01033 0.91110 -0.01789
|
|
0.00000 0.01789 0.00000 -0.01789 0.83615
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.87685 -0.01113 0.00062 -0.01113 0.00000
|
|
-0.01113 0.87685 -0.00031 -0.01113 0.00054
|
|
0.00062 -0.00031 0.79854 -0.00031 0.00000
|
|
-0.01113 -0.01113 -0.00031 0.87685 -0.00054
|
|
0.00000 0.00054 0.00000 -0.00054 0.79854
|
|
|
|
|
|
====== For Atom 2, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 2. Occ. for lpawu and for spin 1 = 4.23050
|
|
Atom 2. Occ. for lpawu and for spin 2 = 4.40815
|
|
=> On atom 2, local Mag. for lpawu is 0.177653
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.87734 -0.01108 0.00066 -0.01108 0.00000
|
|
-0.01108 0.87734 -0.00033 -0.01108 0.00057
|
|
0.00066 -0.00033 0.79923 -0.00033 0.00000
|
|
-0.01108 -0.01108 -0.00033 0.87734 -0.00057
|
|
0.00000 0.00057 0.00000 -0.00057 0.79923
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.91147 -0.01696 0.02070 -0.01696 0.00000
|
|
-0.01696 0.91147 -0.01035 -0.01696 0.01793
|
|
0.02070 -0.01035 0.83688 -0.01035 0.00000
|
|
-0.01696 -0.01696 -0.01035 0.91147 -0.01793
|
|
0.00000 0.01793 0.00000 -0.01793 0.83688
|
|
|
|
|
|
====== For Atom 3, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 3. Occ. for lpawu and for spin 1 = 1.55467
|
|
Atom 3. Occ. for lpawu and for spin 2 = 1.55472
|
|
=> On atom 3, local Mag. for lpawu is 0.000050
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51822 -0.00081 -0.00081
|
|
-0.00081 0.51822 -0.00081
|
|
-0.00081 -0.00081 0.51822
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51824 -0.00081 -0.00081
|
|
-0.00081 0.51824 -0.00081
|
|
-0.00081 -0.00081 0.51824
|
|
|
|
|
|
====== For Atom 4, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 4. Occ. for lpawu and for spin 1 = 1.55467
|
|
Atom 4. Occ. for lpawu and for spin 2 = 1.55472
|
|
=> On atom 4, local Mag. for lpawu is 0.000050
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51822 -0.00081 -0.00081
|
|
-0.00081 0.51822 -0.00081
|
|
-0.00081 -0.00081 0.51822
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51824 -0.00081 -0.00081
|
|
-0.00081 0.51824 -0.00081
|
|
-0.00081 -0.00081 0.51824
|
|
|
|
|
|
|
|
|
|
---------- calculate U, (J) start ----------
|
|
|
|
All pawujat ok and equal to 1
|
|
macro_uj=1 and nspden=2: standard determination of Hubbard U-parameter
|
|
|
|
|
|
*********************************************************************
|
|
************************ Linear Response U ************************
|
|
|
|
Info printed for perturbed atom: 1
|
|
|
|
Perturbations Occupations
|
|
--------------- -----------------------------
|
|
alpha [eV] Unscreened Screened
|
|
--------------- -----------------------------
|
|
0.0000000000 8.6380182458 8.6380182458
|
|
0.0500000000 8.6192514477 8.6332174981
|
|
|
|
Scalar response functions:
|
|
Chi0 [eV^-1]: -0.83408
|
|
Chi [eV^-1]: -0.09601
|
|
|
|
The scalar U from the two-point regression scheme is 9.21612 eV.
|
|
*********************************************************************
|
|
*********************************************************************
|
|
Note: For more reliable linear regressions of the response
|
|
matrices, it is advised that you have more than two points.
|
|
See the LRUJ protocol for more information.
|
|
|
|
|
|
URES ii nat r_max U(J)[eV] U_ASA[eV] U_inf[eV]
|
|
URES 1 1 0.00000 2.30403 1.85413 1.79356
|
|
URES 2 8 11.14400 8.07833 6.50090 6.28854
|
|
URES 3 27 12.45940 8.87522 7.14218 6.90887
|
|
URES 4 64 22.28800 9.07215 7.30066 7.06217
|
|
URES 5 125 24.28780 9.14239 7.35719 7.11686
|
|
URES 6 216 33.43200 9.17345 7.38218 7.14103
|
|
|
|
---------- calculate U, (J) end --------------
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 32.694E-14; max= 53.913E-13
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
0.250000000000 0.250000000000 0.250000000000
|
|
0.750000000000 0.750000000000 0.750000000000
|
|
rms dE/dt= 1.4766E-03; max dE/dt= 2.0882E-03; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
2 0.000000000000 0.000000000000 0.000000000000
|
|
3 -0.002088233570 -0.002088233570 -0.002088233570
|
|
4 0.002088233570 0.002088233570 0.002088233570
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 0.00000000000000
|
|
2 4.16991640368920 4.16991640368920 4.16991640368920
|
|
3 2.08495820184460 2.08495820184460 2.08495820184460
|
|
4 6.25487460553380 6.25487460553380 6.25487460553380
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
3 0.00013250213009 0.00013250213009 0.00013250213009
|
|
4 -0.00013250213009 -0.00013250213009 -0.00013250213009
|
|
frms,max,avg= 9.3693155E-05 1.3250213E-04 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 0.00000000000000
|
|
3 0.00681353290635 0.00681353290635 0.00681353290635
|
|
4 -0.00681353290635 -0.00681353290635 -0.00681353290635
|
|
frms,max,avg= 4.8178953E-03 6.8135329E-03 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 7.880000000000 7.880000000000 7.880000000000 bohr
|
|
= 4.169916403689 4.169916403689 4.169916403689 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.40850 Average Vxc (hartree)= -0.46912
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN UP:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14374 -4.13607 -2.37239 -2.33959 -2.33377 -2.33377 -2.31007 -2.31007
|
|
-0.40313 -0.35590 0.07034 0.21141 0.21141 0.22739 0.22739 0.25301
|
|
0.26198 0.26198 0.31991 0.31991 0.37558 0.38825 0.38825 0.39763
|
|
0.40073 0.43735
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN DOWN:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14136 -4.13845 -2.37212 -2.33984 -2.33314 -2.33314 -2.31069 -2.31069
|
|
-0.40313 -0.35590 0.07034 0.21141 0.21141 0.22739 0.22739 0.25301
|
|
0.26197 0.26197 0.31991 0.31991 0.37556 0.38825 0.38825 0.39763
|
|
0.40077 0.43736
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 4, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.94669717791619E+01
|
|
hartree : 8.65267496867786E+01
|
|
xc : -2.18476794820045E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
local_psp : -2.68557569777781E+02
|
|
spherical_terms : -2.03990180069380E+01
|
|
total_energy : -3.63710798295636E+02
|
|
total_energy_eV : -9.89707414644957E+03
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 4, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -4.00328213969163E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
xc_dc : -7.99671660289640E+01
|
|
spherical_terms : -1.47947012833539E+01
|
|
total_energy_dc : -3.63694941204087E+02
|
|
total_energy_dc_eV : -9.89664265304455E+03
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42464714E-01 sigma(3 2)= -4.89478733E-04
|
|
sigma(2 2)= 1.42464714E-01 sigma(3 1)= -4.89478733E-04
|
|
sigma(3 3)= 1.42464714E-01 sigma(2 1)= -4.89478733E-04
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -4.1915E+03 GPa]
|
|
- sigma(1 1)= 4.19145589E+03 sigma(3 2)= -1.44009591E+01
|
|
- sigma(2 2)= 4.19145589E+03 sigma(3 1)= -1.44009591E+01
|
|
- sigma(3 3)= 4.19145589E+03 sigma(2 1)= -1.44009591E+01
|
|
|
|
================================================================================
|
|
== DATASET 5 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 5, }
|
|
dimensions: {natom: 4, nkpt: 10, mband: 26, nsppol: 2, nspinor: 1, nspden: 2, mpw: 376, }
|
|
cutoff_energies: {ecut: 10.0, pawecutdg: 20.0, }
|
|
electrons: {nelect: 4.80000000E+01, charge: 0.00000000E+00, occopt: 0.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename: getwfk from: ./tlruj_1.o_WFK
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 3.9400000 3.9400000 7.8800000 G(1)= -0.0634518 -0.0634518 0.1903553
|
|
R(2)= 7.8800000 3.9400000 3.9400000 G(2)= 0.1903553 -0.0634518 -0.0634518
|
|
R(3)= 3.9400000 7.8800000 3.9400000 G(3)= -0.0634518 0.1903553 -0.0634518
|
|
Unit cell volume ucvol= 2.4465194E+02 bohr^3
|
|
Angles (23,13,12)= 3.35573098E+01 3.35573098E+01 3.35573098E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 30 30 30
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.18911
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 40 40 40
|
|
ecut(hartree)= 20.000 => boxcut(ratio)= 2.06006
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file ./tlruj_1.o_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 368.875 368.830
|
|
|
|
|
|
******************************************
|
|
DFT+U Method used: FLL
|
|
******************************************
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 5, }
|
|
solver: {iscf: 17, nstep: 30, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-09, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -363.67922154433 -3.637E+02 5.914E-08 2.311E-01
|
|
ETOT 2 -363.68213345805 -2.912E-03 5.813E-07 4.661E+00
|
|
ETOT 3 -363.67908298786 3.050E-03 1.930E-04 1.244E-02
|
|
ETOT 4 -363.67908195238 1.035E-06 5.630E-07 6.073E-04
|
|
ETOT 5 -363.67908203699 -8.461E-08 3.281E-07 1.308E-04
|
|
ETOT 6 -363.67908228985 -2.529E-07 2.570E-08 1.486E-04
|
|
ETOT 7 -363.67908227738 1.248E-08 1.661E-08 3.523E-05
|
|
ETOT 8 -363.67908229015 -1.277E-08 2.671E-09 3.097E-05
|
|
ETOT 9 -363.67908231464 -2.449E-08 2.987E-09 2.349E-05
|
|
ETOT 10 -363.67908238539 -7.075E-08 5.235E-09 2.201E-05
|
|
ETOT 11 -363.67908241952 -3.413E-08 1.977E-08 1.513E-06
|
|
ETOT 12 -363.67908243720 -1.768E-08 1.259E-09 2.766E-07
|
|
ETOT 13 -363.67908244056 -3.362E-09 1.427E-09 1.039E-07
|
|
ETOT 14 -363.67908244288 -2.320E-09 2.806E-10 1.292E-08
|
|
ETOT 15 -363.67908244305 -1.687E-10 2.330E-10 5.456E-09
|
|
ETOT 16 -363.67908244312 -7.236E-11 5.206E-11 1.050E-09
|
|
ETOT 17 -363.67908244313 -8.356E-12 3.529E-11 6.783E-10
|
|
|
|
At SCF step 17 nres2 = 6.78E-10 < tolvrs= 1.00E-09 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42458350E-01 sigma(3 2)= -4.91229279E-04
|
|
sigma(2 2)= 1.42458350E-01 sigma(3 1)= -4.91229279E-04
|
|
sigma(3 3)= 1.42458350E-01 sigma(2 1)= -4.91229279E-04
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 5, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 3.9400000, 3.9400000, 7.8800000, ]
|
|
- [ 7.8800000, 3.9400000, 3.9400000, ]
|
|
- [ 3.9400000, 7.8800000, 3.9400000, ]
|
|
lattice_lengths: [ 9.65099, 9.65099, 9.65099, ]
|
|
lattice_angles: [ 33.557, 33.557, 33.557, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.4465194E+02
|
|
convergence: {deltae: -8.356E-12, res2: 6.783E-10, residm: 3.529E-11, diffor: null, }
|
|
etotal : -3.63679082E+02
|
|
entropy : 0.00000000E+00
|
|
fermie : 4.08931266E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 1.42458350E-01, -4.91229279E-04, -4.91229279E-04, ]
|
|
- [ -4.91229279E-04, 1.42458350E-01, -4.91229279E-04, ]
|
|
- [ -4.91229279E-04, -4.91229279E-04, 1.42458350E-01, ]
|
|
pressure_GPa: -4.1913E+03
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Ni]
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ni]
|
|
- [ 2.5000E-01, 2.5000E-01, 2.5000E-01, O]
|
|
- [ 7.5000E-01, 7.5000E-01, 7.5000E-01, O]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ 2.71050543E-20, -2.71050543E-20, -0.00000000E+00, ]
|
|
- [ 2.71050543E-20, -2.71050543E-20, -0.00000000E+00, ]
|
|
- [ 2.65775977E-04, 2.65775977E-04, 2.65775977E-04, ]
|
|
- [ -2.65775977E-04, -2.65775977E-04, -2.65775977E-04, ]
|
|
force_length_stats: {min: 3.83323354E-20, max: 4.60337496E-04, mean: 2.30168748E-04, }
|
|
...
|
|
|
|
Integrated electronic and magnetization densities in atomic spheres:
|
|
---------------------------------------------------------------------
|
|
Radius=ratsph(iatom), smearing ratsm= 0.0000. Diff(up-dn)=approximate z local magnetic moment.
|
|
Atom Radius up_density dn_density Total(up+dn) Diff(up-dn)
|
|
1 1.81432 8.021144 7.862221 15.883365 0.158923
|
|
2 1.81432 7.866357 8.024650 15.891007 -0.158293
|
|
3 1.41465 2.292341 2.292552 4.584894 -0.000211
|
|
4 1.41465 2.292341 2.292552 4.584894 -0.000211
|
|
---------------------------------------------------------------------
|
|
Sum: 20.472184 20.471975 40.944159 0.000208
|
|
Total magnetization (from the atomic spheres): 0.000208
|
|
Total magnetization (exact up - dn): 0.000000
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 12.716914199061700
|
|
Compensation charge over fine fft grid = 12.717201808624242
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1 - Spin component 1
|
|
-2.33345 0.37622 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00172 -0.00172 0.00000 -0.00172 ...
|
|
0.37622 0.04433 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00034 0.00034 0.00000 0.00034 ...
|
|
0.00000 0.00000 -2.24455 -0.00127 -0.00127 0.82343 0.00078 0.00078 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00127 -2.24455 -0.00127 0.00078 0.82343 0.00078 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00127 -0.00127 -2.24455 0.00078 0.00078 0.82343 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82343 0.00078 0.00078 0.36934 -0.00047 -0.00047 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00078 0.82343 0.00078 -0.00047 0.36934 -0.00047 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00078 0.00078 0.82343 -0.00047 -0.00047 0.36934 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00172 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.91037 -0.00103 0.00094 -0.00103 ...
|
|
-0.00172 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00103 -0.91037 -0.00047 -0.00103 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00094 -0.00047 -0.91166 -0.00047 ...
|
|
-0.00172 0.00034 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00103 -0.00103 -0.00047 -0.91037 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
-2.32960 0.37551 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00223 -0.00223 0.00000 -0.00223 ...
|
|
0.37551 0.04443 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00044 0.00044 0.00000 0.00044 ...
|
|
0.00000 0.00000 -2.24132 -0.00165 -0.00165 0.82163 0.00099 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -2.24132 -0.00165 0.00099 0.82163 0.00099 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 -0.00165 -0.00165 -2.24132 0.00099 0.00099 0.82163 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.82163 0.00099 0.00099 0.37042 -0.00060 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.82163 0.00099 -0.00060 0.37042 -0.00060 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00099 0.00099 0.82163 -0.00060 -0.00060 0.37042 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.90758 -0.00114 0.00137 -0.00114 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.90758 -0.00069 -0.00114 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00137 -0.00069 -0.90867 -0.00069 ...
|
|
-0.00223 0.00044 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00114 -0.00114 -0.00069 -0.90758 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
0.61335 -1.75544 0.00002 0.00002 0.00002 -0.00011 -0.00011 -0.00011
|
|
-1.75544 5.08470 -0.00004 -0.00004 -0.00004 0.00029 0.00029 0.00029
|
|
0.00002 -0.00004 -0.29392 -0.00002 -0.00002 0.88790 0.00011 0.00011
|
|
0.00002 -0.00004 -0.00002 -0.29392 -0.00002 0.00011 0.88790 0.00011
|
|
0.00002 -0.00004 -0.00002 -0.00002 -0.29392 0.00011 0.00011 0.88790
|
|
-0.00011 0.00029 0.88790 0.00011 0.00011 0.12123 -0.00056 -0.00056
|
|
-0.00011 0.00029 0.00011 0.88790 0.00011 -0.00056 0.12123 -0.00056
|
|
-0.00011 0.00029 0.00011 0.00011 0.88790 -0.00056 -0.00056 0.12123
|
|
Atom # 4 - Spin component 2
|
|
0.61335 -1.75547 -0.00001 -0.00001 -0.00001 0.00006 0.00006 0.00006
|
|
-1.75547 5.08478 0.00003 0.00003 0.00003 -0.00017 -0.00017 -0.00017
|
|
-0.00001 0.00003 -0.29392 -0.00002 -0.00002 0.88792 0.00011 0.00011
|
|
-0.00001 0.00003 -0.00002 -0.29392 -0.00002 0.00011 0.88792 0.00011
|
|
-0.00001 0.00003 -0.00002 -0.00002 -0.29392 0.00011 0.00011 0.88792
|
|
0.00006 -0.00017 0.88792 0.00011 0.00011 0.12112 -0.00055 -0.00055
|
|
0.00006 -0.00017 0.00011 0.88792 0.00011 -0.00055 0.12112 -0.00055
|
|
0.00006 -0.00017 0.00011 0.00011 0.88792 -0.00055 -0.00055 0.12112
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1 - Spin component 1
|
|
0.97687 -0.16719 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00017 -0.00017 0.00000 -0.00017 ...
|
|
-0.16719 0.46268 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00442 -0.00442 0.00000 -0.00442 ...
|
|
0.00000 0.00000 1.06707 0.00002 0.00002 0.13560 0.00007 0.00007 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 1.06707 0.00002 0.00007 0.13560 0.00007 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00002 0.00002 1.06707 0.00007 0.00007 0.13560 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13560 0.00007 0.00007 0.02463 0.00002 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00007 0.13560 0.00007 0.00002 0.02463 0.00002 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00007 0.00007 0.13560 0.00002 0.00002 0.02463 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00017 -0.00442 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.15379 -0.02175 0.02667 -0.02175 ...
|
|
-0.00017 -0.00442 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02175 1.15379 -0.01334 -0.02175 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.02667 -0.01334 1.05774 -0.01334 ...
|
|
-0.00017 -0.00442 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.02175 -0.02175 -0.01334 1.15379 ...
|
|
... only 12 components have been written...
|
|
Atom # 1 - Spin component 2
|
|
0.97680 -0.16763 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00009 -0.00009 0.00000 -0.00009 ...
|
|
-0.16763 0.46353 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00011 -0.00011 0.00000 -0.00011 ...
|
|
0.00000 0.00000 1.06704 0.00003 0.00003 0.13549 0.00012 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 1.06704 0.00003 0.00012 0.13549 0.00012 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00003 0.00003 1.06704 0.00012 0.00012 0.13549 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.13549 0.00012 0.00012 0.02455 0.00004 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.13549 0.00012 0.00004 0.02455 0.00004 0.00000 0.00000 0.00000 0.00000 ...
|
|
0.00000 0.00000 0.00012 0.00012 0.13549 0.00004 0.00004 0.02455 0.00000 0.00000 0.00000 0.00000 ...
|
|
-0.00009 -0.00011 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 1.10853 -0.01449 0.00051 -0.01449 ...
|
|
-0.00009 -0.00011 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01449 1.10853 -0.00026 -0.01449 ...
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00051 -0.00026 1.00973 -0.00026 ...
|
|
-0.00009 -0.00011 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.01449 -0.01449 -0.00026 1.10853 ...
|
|
... only 12 components have been written...
|
|
Atom # 4 - Spin component 1
|
|
1.01162 0.01395 -0.00072 -0.00072 -0.00072 0.00003 0.00003 0.00003
|
|
0.01395 0.00120 0.00011 0.00011 0.00011 -0.00000 -0.00000 -0.00000
|
|
-0.00072 0.00011 1.06836 -0.00317 -0.00317 0.05053 -0.00033 -0.00033
|
|
-0.00072 0.00011 -0.00317 1.06836 -0.00317 -0.00033 0.05053 -0.00033
|
|
-0.00072 0.00011 -0.00317 -0.00317 1.06836 -0.00033 -0.00033 0.05053
|
|
0.00003 -0.00000 0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
0.00003 -0.00000 -0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
0.00003 -0.00000 -0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - Spin component 2
|
|
1.01169 0.01394 0.00092 0.00092 0.00092 -0.00003 -0.00003 -0.00003
|
|
0.01394 0.00120 -0.00017 -0.00017 -0.00017 -0.00000 -0.00000 -0.00000
|
|
0.00092 -0.00017 1.06841 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
0.00092 -0.00017 -0.00318 1.06841 -0.00318 -0.00033 0.05053 -0.00033
|
|
0.00092 -0.00017 -0.00318 -0.00318 1.06841 -0.00033 -0.00033 0.05053
|
|
-0.00003 -0.00000 0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00003 -0.00000 -0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00003 -0.00000 -0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
|
|
"PAW+U" part of augmentation waves occupancies Rhoij:
|
|
Atom # 1 - L=2 ONLY - Spin component 1
|
|
1.15379 -0.02175 0.02667 -0.02175 0.00000 0.42883 -0.00885 0.01032 -0.00885 0.00000
|
|
-0.02175 1.15379 -0.01334 -0.02175 0.02310 -0.00885 0.42883 -0.00516 -0.00885 0.00894
|
|
0.02667 -0.01334 1.05774 -0.01334 0.00000 0.01204 -0.00602 0.39104 -0.00602 0.00000
|
|
-0.02175 -0.02175 -0.01334 1.15379 -0.02310 -0.00885 -0.00885 -0.00516 0.42883 -0.00894
|
|
0.00000 0.02310 0.00000 -0.02310 1.05774 0.00000 0.01043 0.00000 -0.01043 0.39104
|
|
0.42883 -0.00885 0.01204 -0.00885 0.00000 0.16397 -0.00363 0.00472 -0.00363 0.00000
|
|
-0.00885 0.42883 -0.00602 -0.00885 0.01043 -0.00363 0.16397 -0.00236 -0.00363 0.00408
|
|
0.01032 -0.00516 0.39104 -0.00516 0.00000 0.00472 -0.00236 0.15157 -0.00236 0.00000
|
|
-0.00885 -0.00885 -0.00602 0.42883 -0.01043 -0.00363 -0.00363 -0.00236 0.16397 -0.00408
|
|
0.00000 0.00894 0.00000 -0.00894 0.39104 0.00000 0.00408 0.00000 -0.00408 0.15157
|
|
Atom # 1 - L=2 ONLY - Spin component 2
|
|
1.10853 -0.01449 0.00051 -0.01449 0.00000 0.40813 -0.00572 -0.00029 -0.00572 0.00000
|
|
-0.01449 1.10853 -0.00026 -0.01449 0.00044 -0.00572 0.40813 0.00014 -0.00572 -0.00025
|
|
0.00051 -0.00026 1.00973 -0.00026 0.00000 -0.00051 0.00025 0.37281 0.00025 0.00000
|
|
-0.01449 -0.01449 -0.00026 1.10853 -0.00044 -0.00572 -0.00572 0.00014 0.40813 0.00025
|
|
0.00000 0.00044 0.00000 -0.00044 1.00973 0.00000 -0.00044 0.00000 0.00044 0.37281
|
|
0.40813 -0.00572 -0.00051 -0.00572 0.00000 0.15448 -0.00229 -0.00035 -0.00229 0.00000
|
|
-0.00572 0.40813 0.00025 -0.00572 -0.00044 -0.00229 0.15448 0.00018 -0.00229 -0.00031
|
|
-0.00029 0.00014 0.37281 0.00014 0.00000 -0.00035 0.00018 0.14461 0.00018 0.00000
|
|
-0.00572 -0.00572 0.00025 0.40813 0.00044 -0.00229 -0.00229 0.00018 0.15448 0.00031
|
|
0.00000 -0.00025 0.00000 0.00025 0.37281 0.00000 -0.00031 0.00000 0.00031 0.14461
|
|
Atom # 2 - L=2 ONLY - Spin component 1
|
|
1.10991 -0.01431 0.00063 -0.01431 0.00000 0.40882 -0.00564 -0.00021 -0.00564 0.00000
|
|
-0.01431 1.10991 -0.00031 -0.01431 0.00054 -0.00564 0.40882 0.00011 -0.00564 -0.00018
|
|
0.00063 -0.00031 1.01152 -0.00031 0.00000 -0.00040 0.00020 0.37356 0.00020 0.00000
|
|
-0.01431 -0.01431 -0.00031 1.10991 -0.00054 -0.00564 -0.00564 0.00011 0.40882 0.00018
|
|
0.00000 0.00054 0.00000 -0.00054 1.01152 0.00000 -0.00035 0.00000 0.00035 0.37356
|
|
0.40882 -0.00564 -0.00040 -0.00564 0.00000 0.15481 -0.00225 -0.00030 -0.00225 0.00000
|
|
-0.00564 0.40882 0.00020 -0.00564 -0.00035 -0.00225 0.15481 0.00015 -0.00225 -0.00026
|
|
-0.00021 0.00011 0.37356 0.00011 0.00000 -0.00030 0.00015 0.14490 0.00015 0.00000
|
|
-0.00564 -0.00564 0.00020 0.40882 0.00035 -0.00225 -0.00225 0.00015 0.15481 0.00026
|
|
0.00000 -0.00018 0.00000 0.00018 0.37356 0.00000 -0.00026 0.00000 0.00026 0.14490
|
|
Atom # 2 - L=2 ONLY - Spin component 2
|
|
1.15474 -0.02181 0.02679 -0.02181 0.00000 0.42930 -0.00887 0.01039 -0.00887 0.00000
|
|
-0.02181 1.15474 -0.01340 -0.02181 0.02320 -0.00887 0.42930 -0.00520 -0.00887 0.00900
|
|
0.02679 -0.01340 1.05960 -0.01340 0.00000 0.01213 -0.00606 0.39182 -0.00606 0.00000
|
|
-0.02181 -0.02181 -0.01340 1.15474 -0.02320 -0.00887 -0.00887 -0.00520 0.42930 -0.00900
|
|
0.00000 0.02320 0.00000 -0.02320 1.05960 0.00000 0.01050 0.00000 -0.01050 0.39182
|
|
0.42930 -0.00887 0.01213 -0.00887 0.00000 0.16418 -0.00364 0.00476 -0.00364 0.00000
|
|
-0.00887 0.42930 -0.00606 -0.00887 0.01050 -0.00364 0.16418 -0.00238 -0.00364 0.00412
|
|
0.01039 -0.00520 0.39182 -0.00520 0.00000 0.00476 -0.00238 0.15188 -0.00238 0.00000
|
|
-0.00887 -0.00887 -0.00606 0.42930 -0.01050 -0.00364 -0.00364 -0.00238 0.16418 -0.00412
|
|
0.00000 0.00900 0.00000 -0.00900 0.39182 0.00000 0.00412 0.00000 -0.00412 0.15188
|
|
Atom # 3 - L=1 ONLY - Spin component 1
|
|
1.06836 -0.00317 -0.00317 0.05053 -0.00033 -0.00033
|
|
-0.00317 1.06836 -0.00317 -0.00033 0.05053 -0.00033
|
|
-0.00317 -0.00317 1.06836 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 3 - L=1 ONLY - Spin component 2
|
|
1.06841 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
-0.00318 1.06841 -0.00318 -0.00033 0.05053 -0.00033
|
|
-0.00318 -0.00318 1.06841 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 1
|
|
1.06836 -0.00317 -0.00317 0.05053 -0.00033 -0.00033
|
|
-0.00317 1.06836 -0.00317 -0.00033 0.05053 -0.00033
|
|
-0.00317 -0.00317 1.06836 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
Atom # 4 - L=1 ONLY - Spin component 2
|
|
1.06841 -0.00318 -0.00318 0.05053 -0.00033 -0.00033
|
|
-0.00318 1.06841 -0.00318 -0.00033 0.05053 -0.00033
|
|
-0.00318 -0.00318 1.06841 -0.00033 -0.00033 0.05053
|
|
0.05053 -0.00033 -0.00033 0.00287 -0.00001 -0.00001
|
|
-0.00033 0.05053 -0.00033 -0.00001 0.00287 -0.00001
|
|
-0.00033 -0.00033 0.05053 -0.00001 -0.00001 0.00287
|
|
|
|
---------- DFT+U DATA ---------------------------------------------------
|
|
|
|
====== For Atom 1, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 1. Occ. for lpawu and for spin 1 = 4.40378
|
|
Atom 1. Occ. for lpawu and for spin 2 = 4.22453
|
|
=> On atom 1, local Mag. for lpawu is -0.179246
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.91080 -0.01694 0.02067 -0.01694 0.00000
|
|
-0.01694 0.91080 -0.01034 -0.01694 0.01790
|
|
0.02067 -0.01034 0.83569 -0.01034 0.00000
|
|
-0.01694 -0.01694 -0.01034 0.91080 -0.01790
|
|
0.00000 0.01790 0.00000 -0.01790 0.83569
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.87624 -0.01134 0.00058 -0.01134 0.00000
|
|
-0.01134 0.87624 -0.00029 -0.01134 0.00050
|
|
0.00058 -0.00029 0.79791 -0.00029 0.00000
|
|
-0.01134 -0.01134 -0.00029 0.87624 -0.00050
|
|
0.00000 0.00050 0.00000 -0.00050 0.79791
|
|
|
|
|
|
====== For Atom 2, occupations for correlated orbitals. lpawu = 2
|
|
|
|
Atom 2. Occ. for lpawu and for spin 1 = 4.23040
|
|
Atom 2. Occ. for lpawu and for spin 2 = 4.40879
|
|
=> On atom 2, local Mag. for lpawu is 0.178393
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.87727 -0.01120 0.00066 -0.01120 0.00000
|
|
-0.01120 0.87727 -0.00033 -0.01120 0.00057
|
|
0.00066 -0.00033 0.79929 -0.00033 0.00000
|
|
-0.01120 -0.01120 -0.00033 0.87727 -0.00057
|
|
0.00000 0.00057 0.00000 -0.00057 0.79929
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.91151 -0.01699 0.02076 -0.01699 0.00000
|
|
-0.01699 0.91151 -0.01038 -0.01699 0.01798
|
|
0.02076 -0.01038 0.83713 -0.01038 0.00000
|
|
-0.01699 -0.01699 -0.01038 0.91151 -0.01798
|
|
0.00000 0.01798 0.00000 -0.01798 0.83713
|
|
|
|
|
|
====== For Atom 3, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 3. Occ. for lpawu and for spin 1 = 1.55464
|
|
Atom 3. Occ. for lpawu and for spin 2 = 1.55475
|
|
=> On atom 3, local Mag. for lpawu is 0.000105
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51821 -0.00080 -0.00080
|
|
-0.00080 0.51821 -0.00080
|
|
-0.00080 -0.00080 0.51821
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51825 -0.00081 -0.00081
|
|
-0.00081 0.51825 -0.00081
|
|
-0.00081 -0.00081 0.51825
|
|
|
|
|
|
====== For Atom 4, occupations for correlated orbitals. lpawu = 1
|
|
|
|
Atom 4. Occ. for lpawu and for spin 1 = 1.55464
|
|
Atom 4. Occ. for lpawu and for spin 2 = 1.55475
|
|
=> On atom 4, local Mag. for lpawu is 0.000105
|
|
|
|
== Occupation matrix for correlated orbitals:
|
|
|
|
Occupation matrix for spin 1
|
|
0.51821 -0.00080 -0.00080
|
|
-0.00080 0.51821 -0.00080
|
|
-0.00080 -0.00080 0.51821
|
|
|
|
|
|
Occupation matrix for spin 2
|
|
0.51825 -0.00081 -0.00081
|
|
-0.00081 0.51825 -0.00081
|
|
-0.00081 -0.00081 0.51825
|
|
|
|
|
|
|
|
|
|
---------- calculate U, (J) start ----------
|
|
|
|
All pawujat ok and equal to 1
|
|
macro_uj=1 and nspden=2: standard determination of Hubbard U-parameter
|
|
|
|
|
|
*********************************************************************
|
|
************************ Linear Response U ************************
|
|
|
|
Info printed for perturbed atom: 1
|
|
|
|
Perturbations Occupations
|
|
--------------- -----------------------------
|
|
alpha [eV] Unscreened Screened
|
|
--------------- -----------------------------
|
|
0.0000000000 8.6380182458 8.6380182458
|
|
0.1000000000 8.5980252002 8.6283148749
|
|
|
|
Scalar response functions:
|
|
Chi0 [eV^-1]: -0.88873
|
|
Chi [eV^-1]: -0.09703
|
|
|
|
The scalar U from the two-point regression scheme is 9.18050 eV.
|
|
*********************************************************************
|
|
*********************************************************************
|
|
Note: For more reliable linear regressions of the response
|
|
matrices, it is advised that you have more than two points.
|
|
See the LRUJ protocol for more information.
|
|
|
|
|
|
URES ii nat r_max U(J)[eV] U_ASA[eV] U_inf[eV]
|
|
URES 1 1 0.00000 2.29513 1.84696 1.78663
|
|
URES 2 8 11.14400 8.04711 6.47577 6.26423
|
|
URES 3 27 12.45940 8.84092 7.11458 6.88217
|
|
URES 4 64 22.28800 9.03709 7.27245 7.03488
|
|
URES 5 125 24.28780 9.10706 7.32876 7.08935
|
|
URES 6 216 33.43200 9.13800 7.35365 7.11343
|
|
|
|
---------- calculate U, (J) end --------------
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 85.855E-15; max= 35.292E-12
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
0.250000000000 0.250000000000 0.250000000000
|
|
0.750000000000 0.750000000000 0.750000000000
|
|
rms dE/dt= 2.9618E-03; max dE/dt= 4.1886E-03; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
2 0.000000000000 0.000000000000 0.000000000000
|
|
3 -0.004188629398 -0.004188629398 -0.004188629398
|
|
4 0.004188629398 0.004188629398 0.004188629398
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 0.00000000000000
|
|
2 4.16991640368920 4.16991640368920 4.16991640368920
|
|
3 2.08495820184460 2.08495820184460 2.08495820184460
|
|
4 6.25487460553380 6.25487460553380 6.25487460553380
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 0.00026577597706 0.00026577597706 0.00026577597706
|
|
4 -0.00026577597706 -0.00026577597706 -0.00026577597706
|
|
frms,max,avg= 1.8793200E-04 2.6577598E-04 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
2 0.00000000000000 -0.00000000000000 -0.00000000000000
|
|
3 0.01366674908723 0.01366674908723 0.01366674908723
|
|
4 -0.01366674908723 -0.01366674908723 -0.01366674908723
|
|
frms,max,avg= 9.6638510E-03 1.3666749E-02 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 7.880000000000 7.880000000000 7.880000000000 bohr
|
|
= 4.169916403689 4.169916403689 4.169916403689 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.40893 Average Vxc (hartree)= -0.46916
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN UP:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14570 -4.13562 -2.37334 -2.34006 -2.33491 -2.33491 -2.31028 -2.31028
|
|
-0.40302 -0.35577 0.07041 0.21163 0.21163 0.22778 0.22778 0.25312
|
|
0.26215 0.26215 0.32030 0.32030 0.37597 0.38867 0.38867 0.39799
|
|
0.40075 0.43774
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
Eigenvalues (hartree) for nkpt= 10 k points, SPIN DOWN:
|
|
kpt# 1, nband= 26, wtk= 0.03125, kpt= 0.1250 0.1250 0.1250 (reduced coord)
|
|
-4.14094 -4.14036 -2.37280 -2.34056 -2.33366 -2.33366 -2.31151 -2.31151
|
|
-0.40303 -0.35578 0.07040 0.21162 0.21162 0.22778 0.22778 0.25311
|
|
0.26214 0.26214 0.32031 0.32031 0.37594 0.38867 0.38867 0.39799
|
|
0.40082 0.43777
|
|
prteigrs : prtvol=0 or 1, do not print more k-points.
|
|
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 5, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.94590398971003E+01
|
|
hartree : 8.65097584018907E+01
|
|
xc : -2.18469212836194E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
local_psp : -2.68533458846765E+02
|
|
spherical_terms : -2.03989579211726E+01
|
|
total_energy : -3.63710792247419E+02
|
|
total_energy_eV : -9.89707398186923E+03
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 5, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -4.00354704420935E+01
|
|
Ewald energy : -2.39331554544165E+02
|
|
psp_core : 1.04313020493122E+01
|
|
xc_dc : -7.99503807961414E+01
|
|
spherical_terms : -1.47929787100448E+01
|
|
total_energy_dc : -3.63679082443133E+02
|
|
total_energy_dc_eV : -9.89621111421267E+03
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 1.42458350E-01 sigma(3 2)= -4.91229279E-04
|
|
sigma(2 2)= 1.42458350E-01 sigma(3 1)= -4.91229279E-04
|
|
sigma(3 3)= 1.42458350E-01 sigma(2 1)= -4.91229279E-04
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -4.1913E+03 GPa]
|
|
- sigma(1 1)= 4.19126865E+03 sigma(3 2)= -1.44524619E+01
|
|
- sigma(2 2)= 4.19126865E+03 sigma(3 1)= -1.44524619E+01
|
|
- sigma(3 3)= 4.19126865E+03 sigma(2 1)= -1.44524619E+01
|
|
|
|
== END DATASET(S) ==============================================================
|
|
================================================================================
|
|
|
|
-outvars: echo values of variables after computation --------
|
|
acell 7.8800000000E+00 7.8800000000E+00 7.8800000000E+00 Bohr
|
|
amu 5.86900000E+01 5.86900000E+01 1.59994000E+01
|
|
atvshift1 -0.00551 -0.00551 -0.00551 -0.00551 -0.00551
|
|
-0.00551 -0.00551 -0.00551 -0.00551 -0.00551
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift2 -0.00367 -0.00367 -0.00367 -0.00367 -0.00367
|
|
-0.00367 -0.00367 -0.00367 -0.00367 -0.00367
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift3 -0.00184 -0.00184 -0.00184 -0.00184 -0.00184
|
|
-0.00184 -0.00184 -0.00184 -0.00184 -0.00184
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift4 0.00184 0.00184 0.00184 0.00184 0.00184
|
|
0.00184 0.00184 0.00184 0.00184 0.00184
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
atvshift5 0.00367 0.00367 0.00367 0.00367 0.00367
|
|
0.00367 0.00367 0.00367 0.00367 0.00367
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000
|
|
chkprim 0
|
|
chksymbreak 0
|
|
diemix 4.50000000E-01
|
|
dmatpuopt 3
|
|
ecut 1.00000000E+01 Hartree
|
|
etotal1 -3.6375846399E+02
|
|
etotal2 -3.6374257036E+02
|
|
etotal3 -3.6372668530E+02
|
|
etotal4 -3.6369494120E+02
|
|
etotal5 -3.6367908244E+02
|
|
fcart1 -0.0000000000E+00 -0.0000000000E+00 -0.0000000000E+00
|
|
-0.0000000000E+00 -0.0000000000E+00 -0.0000000000E+00
|
|
-3.8419187146E-04 -3.8419187146E-04 -3.8419187146E-04
|
|
3.8419187146E-04 3.8419187146E-04 3.8419187146E-04
|
|
fcart2 -0.0000000000E+00 -0.0000000000E+00 -0.0000000000E+00
|
|
-0.0000000000E+00 -0.0000000000E+00 -0.0000000000E+00
|
|
-2.5726919321E-04 -2.5726919321E-04 -2.5726919321E-04
|
|
2.5726919321E-04 2.5726919321E-04 2.5726919321E-04
|
|
fcart3 2.0328790734E-20 -3.3881317890E-20 6.7762635780E-21
|
|
2.0328790734E-20 -3.3881317890E-20 6.7762635780E-21
|
|
-1.3023288118E-04 -1.3023288118E-04 -1.3023288118E-04
|
|
1.3023288118E-04 1.3023288118E-04 1.3023288118E-04
|
|
fcart4 6.7762635780E-21 -2.0328790734E-20 6.7762635780E-21
|
|
6.7762635780E-21 -2.0328790734E-20 6.7762635780E-21
|
|
1.3250213009E-04 1.3250213009E-04 1.3250213009E-04
|
|
-1.3250213009E-04 -1.3250213009E-04 -1.3250213009E-04
|
|
fcart5 2.7105054312E-20 -2.7105054312E-20 -0.0000000000E+00
|
|
2.7105054312E-20 -2.7105054312E-20 -0.0000000000E+00
|
|
2.6577597706E-04 2.6577597706E-04 2.6577597706E-04
|
|
-2.6577597706E-04 -2.6577597706E-04 -2.6577597706E-04
|
|
- fftalg 512
|
|
ixc -101130
|
|
jdtset 1 2 3 4 5
|
|
kpt 1.25000000E-01 1.25000000E-01 1.25000000E-01
|
|
3.75000000E-01 1.25000000E-01 1.25000000E-01
|
|
-3.75000000E-01 1.25000000E-01 1.25000000E-01
|
|
-1.25000000E-01 1.25000000E-01 1.25000000E-01
|
|
3.75000000E-01 3.75000000E-01 1.25000000E-01
|
|
-3.75000000E-01 3.75000000E-01 1.25000000E-01
|
|
-1.25000000E-01 3.75000000E-01 1.25000000E-01
|
|
-3.75000000E-01 -3.75000000E-01 1.25000000E-01
|
|
3.75000000E-01 3.75000000E-01 3.75000000E-01
|
|
-3.75000000E-01 3.75000000E-01 3.75000000E-01
|
|
kptrlatt 4 0 0 0 4 0 0 0 4
|
|
kptrlen 2.22880057E+01
|
|
lpawu 2 2 1
|
|
macro_uj 1
|
|
P mkmem 10
|
|
natom 4
|
|
nband 26
|
|
natvshift 5
|
|
ndtset 5
|
|
ngfft 30 30 30
|
|
ngfftdg 40 40 40
|
|
nkpt 10
|
|
nspden 2
|
|
nsppol 2
|
|
nsym 12
|
|
ntypat 3
|
|
occ 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
0.000000 0.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
|
|
0.000000 0.000000
|
|
occopt 0
|
|
pawecutdg 2.00000000E+01 Hartree
|
|
pawujv1 -5.5123988097E-03 Hartree
|
|
pawujv2 -3.6749325398E-03 Hartree
|
|
pawujv3 -1.8374662699E-03 Hartree
|
|
pawujv4 1.8374662699E-03 Hartree
|
|
pawujv5 3.6749325398E-03 Hartree
|
|
prtden 0
|
|
prtebands 0
|
|
prteig 0
|
|
prtgsr 0
|
|
prtwf 0
|
|
rprim 5.0000000000E-01 5.0000000000E-01 1.0000000000E+00
|
|
1.0000000000E+00 5.0000000000E-01 5.0000000000E-01
|
|
5.0000000000E-01 1.0000000000E+00 5.0000000000E-01
|
|
shiftk 5.00000000E-01 5.00000000E-01 5.00000000E-01
|
|
spgroup 166
|
|
spinat 0.0000000000E+00 0.0000000000E+00 3.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 -3.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten1 1.4248393786E-01 1.4248393786E-01 1.4248393786E-01
|
|
-4.8558537481E-04 -4.8558537481E-04 -4.8558537481E-04
|
|
strten2 1.4247967568E-01 1.4247967568E-01 1.4247967568E-01
|
|
-4.8629968752E-04 -4.8629968752E-04 -4.8629968752E-04
|
|
strten3 1.4247518211E-01 1.4247518211E-01 1.4247518211E-01
|
|
-4.8713661179E-04 -4.8713661179E-04 -4.8713661179E-04
|
|
strten4 1.4246471388E-01 1.4246471388E-01 1.4246471388E-01
|
|
-4.8947873269E-04 -4.8947873269E-04 -4.8947873269E-04
|
|
strten5 1.4245835000E-01 1.4245835000E-01 1.4245835000E-01
|
|
-4.9122927857E-04 -4.9122927857E-04 -4.9122927857E-04
|
|
symrel 1 0 0 0 1 0 0 0 1 -1 0 0 0 -1 0 0 0 -1
|
|
0 1 0 1 0 0 0 0 1 0 -1 0 -1 0 0 0 0 -1
|
|
0 0 1 1 0 0 0 1 0 0 0 -1 -1 0 0 0 -1 0
|
|
1 0 0 0 0 1 0 1 0 -1 0 0 0 0 -1 0 -1 0
|
|
0 1 0 0 0 1 1 0 0 0 -1 0 0 0 -1 -1 0 0
|
|
0 0 1 0 1 0 1 0 0 0 0 -1 0 -1 0 -1 0 0
|
|
tolvrs 1.00000000E-09
|
|
typat 1 2 3 3
|
|
usepawu 1
|
|
useylm 1
|
|
wtk 0.03125 0.09375 0.09375 0.09375 0.09375 0.18750
|
|
0.18750 0.09375 0.03125 0.09375
|
|
xangst 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
4.1699164037E+00 4.1699164037E+00 4.1699164037E+00
|
|
2.0849582018E+00 2.0849582018E+00 2.0849582018E+00
|
|
6.2548746055E+00 6.2548746055E+00 6.2548746055E+00
|
|
xcart 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
7.8800000000E+00 7.8800000000E+00 7.8800000000E+00
|
|
3.9400000000E+00 3.9400000000E+00 3.9400000000E+00
|
|
1.1820000000E+01 1.1820000000E+01 1.1820000000E+01
|
|
xred 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
5.0000000000E-01 5.0000000000E-01 5.0000000000E-01
|
|
2.5000000000E-01 2.5000000000E-01 2.5000000000E-01
|
|
7.5000000000E-01 7.5000000000E-01 7.5000000000E-01
|
|
znucl 28.00000 28.00000 8.00000
|
|
|
|
================================================================================
|
|
|
|
|
|
- Timing analysis has been suppressed with timopt=0
|
|
|
|
|
|
|
|
================================================================================
|
|
|
|
Suggested references for the acknowledgment of ABINIT usage.
|
|
|
|
The users of ABINIT have little formal obligations with respect to the ABINIT group
|
|
(those specified in the GNU General Public License, http://www.gnu.org/copyleft/gpl.txt).
|
|
However, it is common practice in the scientific literature,
|
|
to acknowledge the efforts of people that have made the research possible.
|
|
In this spirit, please find below suggested citations of work written by ABINIT developers,
|
|
corresponding to implementations inside of ABINIT that you have used in the present run.
|
|
Note also that it will be of great value to readers of publications presenting these results,
|
|
to read papers enabling them to understand the theoretical formalism and details
|
|
of the ABINIT implementation.
|
|
For information on why they are suggested, see also https://docs.abinit.org/theory/acknowledgments.
|
|
-
|
|
- [1] Facilities and practices for linear response Hubbard parameters U and J in Abinit.
|
|
- L. MacEnulty, M. Giantomassi, B. Amadon, G.-M. Rignanese and D.D. O'Regan Electron. Struct. 6 037003 (2024).
|
|
- Comment: to be cited in case the Hubbard U or Hund's J are calculated using the lrUJ utility
|
|
- or the renovated functionalities of UJdet, i.e., macro_uj>0.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#MacEnulty2024
|
|
-
|
|
- [2] Gamma and beta cerium: DFT+U calculations of ground-state parameters.
|
|
- B. Amadon, F. Jollet and M. Torrent, Phys. Rev. B 77, 155104 (2008).
|
|
- Comment: DFT+U calculations, usepawu/=0. Strong suggestion to cite this paper.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#amadon2008a
|
|
-
|
|
- [3] Implementation of the Projector Augmented-Wave Method in the ABINIT code.
|
|
- M. Torrent, F. Jollet, F. Bottin, G. Zerah, and X. Gonze Comput. Mat. Science 42, 337, (2008).
|
|
- Comment: PAW calculations. Strong suggestion to cite this paper.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#torrent2008
|
|
-
|
|
- [4] Libxc: A library of exchange and correlation functionals for density functional theory.
|
|
- M.A.L. Marques, M.J.T. Oliveira, T. Burnus, Computer Physics Communications 183, 2227 (2012).
|
|
- Comment: to be cited when LibXC is used (negative value of ixc)
|
|
- Strong suggestion to cite this paper.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#marques2012
|
|
-
|
|
- [5] The Abinit project: Impact, environment and recent developments.
|
|
- Computer Phys. Comm. 248, 107042 (2020).
|
|
- X.Gonze, B. Amadon, G. Antonius, F.Arnardi, L.Baguet, J.-M.Beuken,
|
|
- J.Bieder, F.Bottin, J.Bouchet, E.Bousquet, N.Brouwer, F.Bruneval,
|
|
- G.Brunin, T.Cavignac, J.-B. Charraud, Wei Chen, M.Cote, S.Cottenier,
|
|
- J.Denier, G.Geneste, Ph.Ghosez, M.Giantomassi, Y.Gillet, O.Gingras,
|
|
- D.R.Hamann, G.Hautier, Xu He, N.Helbig, N.Holzwarth, Y.Jia, F.Jollet,
|
|
- W.Lafargue-Dit-Hauret, K.Lejaeghere, M.A.L.Marques, A.Martin, C.Martins,
|
|
- H.P.C. Miranda, F.Naccarato, K. Persson, G.Petretto, V.Planes, Y.Pouillon,
|
|
- S.Prokhorenko, F.Ricci, G.-M.Rignanese, A.H.Romero, M.M.Schmitt, M.Torrent,
|
|
- M.J.van Setten, B.Van Troeye, M.J.Verstraete, G.Zerah and J.W.Zwanzig
|
|
- Comment: the fifth generic paper describing the ABINIT project.
|
|
- Note that a version of this paper, that is not formatted for Computer Phys. Comm.
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT20.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze2020
|
|
-
|
|
- [6] ABINIT: Overview, and focus on selected capabilities
|
|
- J. Chem. Phys. 152, 124102 (2020).
|
|
- A. Romero, D.C. Allan, B. Amadon, G. Antonius, T. Applencourt, L.Baguet,
|
|
- J.Bieder, F.Bottin, J.Bouchet, E.Bousquet, F.Bruneval,
|
|
- G.Brunin, D.Caliste, M.Cote,
|
|
- J.Denier, C. Dreyer, Ph.Ghosez, M.Giantomassi, Y.Gillet, O.Gingras,
|
|
- D.R.Hamann, G.Hautier, F.Jollet, G. Jomard,
|
|
- A.Martin,
|
|
- H.P.C. Miranda, F.Naccarato, G.Petretto, N.A. Pike, V.Planes,
|
|
- S.Prokhorenko, T. Rangel, F.Ricci, G.-M.Rignanese, M.Royo, M.Stengel, M.Torrent,
|
|
- M.J.van Setten, B.Van Troeye, M.J.Verstraete, J.Wiktor, J.W.Zwanziger, and X.Gonze.
|
|
- Comment: a global overview of ABINIT, with focus on selected capabilities .
|
|
- Note that a version of this paper, that is not formatted for J. Chem. Phys
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT20_JPC.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#romero2020
|
|
-
|
|
- [7] Recent developments in the ABINIT software package.
|
|
- Computer Phys. Comm. 205, 106 (2016).
|
|
- X.Gonze, F.Jollet, F.Abreu Araujo, D.Adams, B.Amadon, T.Applencourt,
|
|
- C.Audouze, J.-M.Beuken, J.Bieder, A.Bokhanchuk, E.Bousquet, F.Bruneval
|
|
- D.Caliste, M.Cote, F.Dahm, F.Da Pieve, M.Delaveau, M.Di Gennaro,
|
|
- B.Dorado, C.Espejo, G.Geneste, L.Genovese, A.Gerossier, M.Giantomassi,
|
|
- Y.Gillet, D.R.Hamann, L.He, G.Jomard, J.Laflamme Janssen, S.Le Roux,
|
|
- A.Levitt, A.Lherbier, F.Liu, I.Lukacevic, A.Martin, C.Martins,
|
|
- M.J.T.Oliveira, S.Ponce, Y.Pouillon, T.Rangel, G.-M.Rignanese,
|
|
- A.H.Romero, B.Rousseau, O.Rubel, A.A.Shukri, M.Stankovski, M.Torrent,
|
|
- M.J.Van Setten, B.Van Troeye, M.J.Verstraete, D.Waroquier, J.Wiktor,
|
|
- B.Xu, A.Zhou, J.W.Zwanziger.
|
|
- Comment: the fourth generic paper describing the ABINIT project.
|
|
- Note that a version of this paper, that is not formatted for Computer Phys. Comm.
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT16.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze2016
|
|
-
|
|
- Proc. 0 individual time (sec): cpu= 103.9 wall= 105.6
|
|
|
|
================================================================================
|
|
|
|
Calculation completed.
|
|
.Delivered 20 WARNINGs and 6 COMMENTs to log file.
|
|
+Overall time at end (sec) : cpu= 103.9 wall= 105.6
|