mirror of https://github.com/abinit/abinit.git
1365 lines
67 KiB
Plaintext
1365 lines
67 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 19h17 )
|
|
|
|
- input file -> /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/TestBot_MPI1/v9_t206/t206.abi
|
|
- output file -> t206.abo
|
|
- root for input files -> t206i
|
|
- root for output files -> t206o
|
|
|
|
- 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/P.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/Al.xml
|
|
DATASET 1 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 1.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 8
|
|
lnmax = 4 mgfft = 16 mpssoang = 2 mqgrid = 3001
|
|
natom = 2 nloc_mem = 2 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 1 ntypat = 2
|
|
occopt = 1 xclevel = 2
|
|
- mband = 8 mffmem = 1 mkmem = 8
|
|
mpw = 138 nfft = 4096 nkpt = 8
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 24 nfftf = 13824
|
|
================================================================================
|
|
P This job should need less than 6.159 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.137 Mbytes ; DEN or POT disk file : 0.107 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 2 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 2.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 8
|
|
lnmax = 4 mgfft = 16 mpssoang = 2 mqgrid = 3001
|
|
natom = 2 nloc_mem = 2 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 1 ntypat = 2
|
|
occopt = 1 xclevel = 2
|
|
- mband = 8 mffmem = 1 mkmem = 8
|
|
mpw = 138 nfft = 4096 nkpt = 8
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 24 nfftf = 13824
|
|
================================================================================
|
|
P This job should need less than 6.159 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.137 Mbytes ; DEN or POT disk file : 0.107 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 3 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 3.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 8
|
|
lnmax = 4 mgfft = 16 mpssoang = 2 mqgrid = 3001
|
|
natom = 2 nloc_mem = 2 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 1 ntypat = 2
|
|
occopt = 1 xclevel = 2
|
|
- mband = 8 mffmem = 1 mkmem = 8
|
|
mpw = 138 nfft = 4096 nkpt = 8
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 24 nfftf = 13824
|
|
================================================================================
|
|
P This job should need less than 6.159 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.137 Mbytes ; DEN or POT disk file : 0.107 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 4 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 4.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 8
|
|
lnmax = 4 mgfft = 16 mpssoang = 2 mqgrid = 3001
|
|
natom = 2 nloc_mem = 2 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 1 ntypat = 2
|
|
occopt = 1 xclevel = 2
|
|
- mband = 8 mffmem = 1 mkmem = 8
|
|
mpw = 138 nfft = 4096 nkpt = 8
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 24 nfftf = 13824
|
|
================================================================================
|
|
P This job should need less than 6.159 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.137 Mbytes ; DEN or POT disk file : 0.107 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 1.0000000000E+01 1.0000000000E+01 1.0000000000E+01 Bohr
|
|
amu 3.09737620E+01 2.69815390E+01
|
|
bandpp1 1
|
|
bandpp2 8
|
|
bandpp3 8
|
|
bandpp4 8
|
|
ecut 5.00000000E+00 Hartree
|
|
- fftalg 512
|
|
istwfk 1 1 1 1 1 1 1 1
|
|
ixc -101130
|
|
jdtset 1 2 3 4
|
|
kpt 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
5.00000000E-01 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 5.00000000E-01 0.00000000E+00
|
|
5.00000000E-01 5.00000000E-01 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
5.00000000E-01 0.00000000E+00 5.00000000E-01
|
|
0.00000000E+00 5.00000000E-01 5.00000000E-01
|
|
5.00000000E-01 5.00000000E-01 5.00000000E-01
|
|
kptopt 3
|
|
kptrlatt 2 0 0 0 2 0 0 0 2
|
|
kptrlen 1.41421356E+01
|
|
P mkmem 8
|
|
natom 2
|
|
nband 8
|
|
nbdbuf 2
|
|
ndtset 4
|
|
ngfft 16 16 16
|
|
ngfftdg 24 24 24
|
|
nkpt 8
|
|
nline 8
|
|
nstep 15
|
|
nsym 1
|
|
ntypat 2
|
|
occ 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000
|
|
0.000000 0.000000
|
|
pawecutdg 1.00000000E+01 Hartree
|
|
prtden 0
|
|
prteig 0
|
|
prtwf 0
|
|
rprim 0.0000000000E+00 5.0000000000E-01 5.0000000000E-01
|
|
5.0000000000E-01 0.0000000000E+00 5.0000000000E-01
|
|
5.0000000000E-01 5.0000000000E-01 0.0000000000E+00
|
|
spgroup 1
|
|
tolvrs 1.00000000E-20
|
|
typat 1 2
|
|
useylm 1
|
|
wfoptalg1 0
|
|
wfoptalg2 14
|
|
wfoptalg3 114
|
|
wfoptalg4 111
|
|
wtk 0.12500 0.12500 0.12500 0.12500 0.12500 0.12500
|
|
0.12500 0.12500
|
|
xangst 1.3335265656E+00 1.2938382750E+00 1.3097135913E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
xcart 2.5200000000E+00 2.4450000000E+00 2.4750000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
xred 2.4000000000E-01 2.5500000000E-01 2.4900000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
znucl 15.00000 13.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.
|
|
|
|
================================================================================
|
|
== DATASET 1 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 1, }
|
|
dimensions: {natom: 2, nkpt: 8, mband: 8, nsppol: 1, nspinor: 1, nspden: 1, mpw: 138, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 10.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 0.0000000 5.0000000 5.0000000 G(1)= -0.1000000 0.1000000 0.1000000
|
|
R(2)= 5.0000000 0.0000000 5.0000000 G(2)= 0.1000000 -0.1000000 0.1000000
|
|
R(3)= 5.0000000 5.0000000 0.0000000 G(3)= 0.1000000 0.1000000 -0.1000000
|
|
Unit cell volume ucvol= 2.5000000E+02 bohr^3
|
|
Angles (23,13,12)= 6.00000000E+01 6.00000000E+01 6.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 16
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.24794
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 6.316547 Hartrees makes boxcut=2
|
|
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 24 24 24
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.38430
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 14.212230 Hartrees makes boxcut=2
|
|
|
|
|
|
--- 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/P.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/P.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/P.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 4 (lmn_size= 8), orbitals= 0 0 1 1
|
|
Spheres core radius: rc_sph= 1.90690075
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=2001 , AA= 0.40634E-03 BB= 0.60952E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60765221
|
|
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 1773 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/Al.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Al.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/Al.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 4 (lmn_size= 8), orbitals= 0 0 1 1
|
|
Spheres core radius: rc_sph= 1.90363307
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=2001 , AA= 0.46377E-03 BB= 0.60291E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60786206
|
|
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 1771 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.65993774E+01 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 133.375 133.250
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 1, }
|
|
solver: {iscf: 17, nstep: 15, nline: 8, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -8.6353324191942 -8.635E+00 1.456E-01 1.010E+00
|
|
ETOT 2 -8.7305549888055 -9.522E-02 1.400E-05 3.969E-01
|
|
ETOT 3 -8.7215850875914 8.970E-03 3.611E-04 2.222E-02
|
|
ETOT 4 -8.7217777052295 -1.926E-04 1.247E-06 8.501E-03
|
|
ETOT 5 -8.7218787832869 -1.011E-04 7.978E-06 3.493E-04
|
|
ETOT 6 -8.7218981452984 -1.936E-05 1.010E-07 6.408E-05
|
|
ETOT 7 -8.7219049054492 -6.760E-06 4.218E-08 1.008E-06
|
|
ETOT 8 -8.7219049894897 -8.404E-08 4.501E-10 2.189E-08
|
|
ETOT 9 -8.7219049941177 -4.628E-09 2.466E-11 2.313E-09
|
|
ETOT 10 -8.7219049944094 -2.917E-10 5.835E-12 2.469E-10
|
|
ETOT 11 -8.7219049944313 -2.188E-11 1.926E-13 2.394E-11
|
|
ETOT 12 -8.7219049944333 -1.997E-12 3.064E-14 1.215E-12
|
|
ETOT 13 -8.7219049944336 -3.766E-13 2.058E-15 9.558E-14
|
|
ETOT 14 -8.7219049944337 -9.059E-14 5.095E-17 1.440E-14
|
|
ETOT 15 -8.7219049944337 5.507E-14 2.785E-17 9.768E-16
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910451E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639483E-04 sigma(3 1)= -5.15338770E-05
|
|
sigma(3 3)= -6.02738603E-04 sigma(2 1)= -2.25449262E-05
|
|
|
|
|
|
scprqt: WARNING -
|
|
nstep= 15 was not enough SCF cycles to converge;
|
|
density residual= 9.768E-16 exceeds tolvrs= 1.000E-20
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 1, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 0.0000000, 5.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 0.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 5.0000000, 0.0000000, ]
|
|
lattice_lengths: [ 7.07107, 7.07107, 7.07107, ]
|
|
lattice_angles: [ 60.000, 60.000, 60.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.5000000E+02
|
|
convergence: {deltae: 5.507E-14, res2: 9.768E-16, residm: 2.785E-17, diffor: null, }
|
|
etotal : -8.72190499E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : 2.53959489E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -6.02910451E-04, -2.25449262E-05, -5.15338770E-05, ]
|
|
- [ -2.25449262E-05, -6.00639483E-04, 1.75404114E-05, ]
|
|
- [ -5.15338770E-05, 1.75404114E-05, -6.02738603E-04, ]
|
|
pressure_GPa: 1.7714E+01
|
|
xred :
|
|
- [ 2.4000E-01, 2.5500E-01, 2.4900E-01, P]
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.56720153E-03, 7.78450884E-03, 3.34988651E-03, ]
|
|
- [ 2.56720153E-03, -7.78450884E-03, -3.34988651E-03, ]
|
|
force_length_stats: {min: 8.85498962E-03, max: 8.85498962E-03, mean: 8.85498962E-03, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 1.90690 2.68246148
|
|
2 1.90363 0.77420350
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = -0.524738063033980
|
|
Compensation charge over fine fft grid = -0.524720183571195
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1
|
|
1.36052 0.01008 -0.00007 -0.00003 0.00002 -0.00009 -0.00004 0.00003
|
|
0.01008 28.05725 -0.00011 -0.00004 0.00003 -0.00248 -0.00108 0.00083
|
|
-0.00007 -0.00011 0.26741 -0.00000 0.00000 0.01302 -0.00002 0.00003
|
|
-0.00003 -0.00004 -0.00000 0.26741 0.00001 -0.00002 0.01302 0.00007
|
|
0.00002 0.00003 0.00000 0.00001 0.26741 0.00003 0.00007 0.01302
|
|
-0.00009 -0.00248 0.01302 -0.00002 0.00003 17.66220 -0.00030 0.00042
|
|
-0.00004 -0.00108 -0.00002 0.01302 0.00007 -0.00030 17.66219 0.00107
|
|
0.00003 0.00083 0.00003 0.00007 0.01302 0.00042 0.00107 17.66218
|
|
Atom # 2
|
|
0.32629 -0.05739 0.00003 0.00001 -0.00001 -0.00024 -0.00010 0.00008
|
|
-0.05739 39.68611 -0.00018 -0.00008 0.00006 0.00698 0.00304 -0.00235
|
|
0.00003 -0.00018 0.07346 -0.00000 0.00000 -0.15662 0.00002 -0.00002
|
|
0.00001 -0.00008 -0.00000 0.07346 0.00001 0.00002 -0.15662 -0.00006
|
|
-0.00001 0.00006 0.00000 0.00001 0.07346 -0.00002 -0.00006 -0.15662
|
|
-0.00024 0.00698 -0.15662 0.00002 -0.00002 24.40602 -0.00183 0.00230
|
|
-0.00010 0.00304 0.00002 -0.15662 -0.00006 -0.00183 24.40618 0.00512
|
|
0.00008 -0.00235 -0.00002 -0.00006 -0.15662 0.00230 0.00512 24.40619
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1
|
|
1.50794 0.00728 -0.00593 -0.00262 0.00206 -0.00006 -0.00003 0.00002
|
|
0.00728 0.00004 -0.00001 -0.00000 -0.00000 -0.00000 -0.00000 0.00000
|
|
-0.00593 -0.00001 1.31374 0.00020 -0.00021 0.00718 -0.00002 0.00002
|
|
-0.00262 -0.00000 0.00020 1.31382 -0.00031 -0.00002 0.00718 0.00005
|
|
0.00206 -0.00000 -0.00021 -0.00031 1.31383 0.00002 0.00005 0.00718
|
|
-0.00006 -0.00000 0.00718 -0.00002 0.00002 0.00004 -0.00000 0.00000
|
|
-0.00003 -0.00000 -0.00002 0.00718 0.00005 -0.00000 0.00004 0.00000
|
|
0.00002 0.00000 0.00002 0.00005 0.00718 0.00000 0.00000 0.00004
|
|
Atom # 2
|
|
1.03892 -0.00154 0.01727 0.00749 -0.00579 -0.00001 -0.00000 0.00000
|
|
-0.00154 0.00000 0.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
0.01727 0.00000 0.99942 -0.01315 0.01656 -0.00073 0.00000 -0.00000
|
|
0.00749 -0.00000 -0.01315 1.00049 0.03683 0.00000 -0.00073 -0.00000
|
|
-0.00579 0.00000 0.01656 0.03683 1.00059 -0.00000 -0.00000 -0.00073
|
|
-0.00001 0.00000 -0.00073 0.00000 -0.00000 0.00000 -0.00000 0.00000
|
|
-0.00000 0.00000 0.00000 -0.00073 -0.00000 -0.00000 0.00000 0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00073 0.00000 0.00000 0.00000
|
|
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 79.121E-19; max= 27.847E-18
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.240000000000 0.255000000000 0.249000000000
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 3.5569E-02; max dE/dt= 5.5266E-02; dE/dt below (all hartree)
|
|
1 -0.056077753739 -0.003943882416 -0.026236800521
|
|
2 0.055266199733 0.003882967376 0.025936272591
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 1.33352656564680 1.29383827500255 1.30971359126025
|
|
2 0.00000000000000 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00256720152836 0.00778450883965 0.00334988650753
|
|
2 0.00256720152836 -0.00778450883965 -0.00334988650753
|
|
frms,max,avg= 5.1124306E-03 7.7845088E-03 -2.251E-05 5.256E-05 2.860E-05 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.13201079921863 0.40029550547284 0.17225807567724
|
|
2 0.13201079921863 -0.40029550547284 -0.17225807567724
|
|
frms,max,avg= 2.6289173E-01 4.0029551E-01 -1.157E-03 2.703E-03 1.471E-03 e/A
|
|
length scales= 10.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 5.291772085900 5.291772085900 5.291772085900 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.25396 Average Vxc (hartree)= -0.34016
|
|
Eigenvalues (hartree) for nkpt= 8 k points:
|
|
kpt# 1, nband= 8, wtk= 0.12500, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.19950 0.24708 0.24941 0.25396 0.39574 0.41680 0.41829 0.41922
|
|
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 : 3.88297384689434E+00
|
|
hartree : 9.47041929716052E-01
|
|
xc : -3.18572949447927E+00
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
local_psp : -2.96817232966421E+00
|
|
spherical_terms : 1.70505178467474E+00
|
|
total_energy : -8.72190498423356E+00
|
|
total_energy_eV : -2.37335104516427E+02
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 1, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : 6.63514815980366E-01
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
xc_dc : -3.61188812617761E-01
|
|
spherical_terms : 7.88397235789140E-02
|
|
total_energy_dc : -8.72190499443368E+00
|
|
total_energy_dc_eV : -2.37335104793986E+02
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910451E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639483E-04 sigma(3 1)= -5.15338770E-05
|
|
sigma(3 3)= -6.02738603E-04 sigma(2 1)= -2.25449262E-05
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= 1.7714E+01 GPa]
|
|
- sigma(1 1)= -1.77382349E+01 sigma(3 2)= 5.16056634E-01
|
|
- sigma(2 2)= -1.76714207E+01 sigma(3 1)= -1.51617875E+00
|
|
- sigma(3 3)= -1.77331790E+01 sigma(2 1)= -6.63294518E-01
|
|
|
|
================================================================================
|
|
== DATASET 2 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 2, }
|
|
dimensions: {natom: 2, nkpt: 8, mband: 8, nsppol: 1, nspinor: 1, nspden: 1, mpw: 138, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 10.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 0.0000000 5.0000000 5.0000000 G(1)= -0.1000000 0.1000000 0.1000000
|
|
R(2)= 5.0000000 0.0000000 5.0000000 G(2)= 0.1000000 -0.1000000 0.1000000
|
|
R(3)= 5.0000000 5.0000000 0.0000000 G(3)= 0.1000000 0.1000000 -0.1000000
|
|
Unit cell volume ucvol= 2.5000000E+02 bohr^3
|
|
Angles (23,13,12)= 6.00000000E+01 6.00000000E+01 6.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 16
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.24794
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 6.316547 Hartrees makes boxcut=2
|
|
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 24 24 24
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.38430
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 14.212230 Hartrees makes boxcut=2
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 133.375 133.250
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: 17, nstep: 15, nline: 8, wfoptalg: 14, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -8.7465072911079 -8.747E+00 6.058E-10 1.449E+00
|
|
ETOT 2 -8.7227908923171 2.372E-02 2.371E-17 1.863E-01
|
|
ETOT 3 -8.7214810043079 1.310E-03 9.419E-10 1.265E-02
|
|
ETOT 4 -8.7217969146142 -3.159E-04 1.083E-10 9.954E-04
|
|
ETOT 5 -8.7219026309223 -1.057E-04 1.606E-11 1.762E-05
|
|
ETOT 6 -8.7219049753334 -2.344E-06 1.930E-17 2.389E-07
|
|
ETOT 7 -8.7219049907672 -1.543E-08 5.715E-20 3.042E-08
|
|
ETOT 8 -8.7219049940139 -3.247E-09 5.564E-21 2.663E-09
|
|
ETOT 9 -8.7219049944344 -4.205E-10 8.818E-22 1.046E-10
|
|
ETOT 10 -8.7219049944338 6.537E-13 2.788E-23 6.032E-12
|
|
ETOT 11 -8.7219049944338 -1.776E-15 2.402E-25 2.224E-14
|
|
ETOT 12 -8.7219049944336 1.279E-13 3.516E-27 2.564E-15
|
|
ETOT 13 -8.7219049944337 -7.283E-14 1.099E-27 5.605E-17
|
|
ETOT 14 -8.7219049944337 -3.908E-14 2.421E-29 6.344E-18
|
|
ETOT 15 -8.7219049944337 7.105E-15 1.321E-29 2.139E-19
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449274E-05
|
|
|
|
|
|
scprqt: WARNING -
|
|
nstep= 15 was not enough SCF cycles to converge;
|
|
density residual= 2.139E-19 exceeds tolvrs= 1.000E-20
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 2, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 0.0000000, 5.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 0.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 5.0000000, 0.0000000, ]
|
|
lattice_lengths: [ 7.07107, 7.07107, 7.07107, ]
|
|
lattice_angles: [ 60.000, 60.000, 60.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.5000000E+02
|
|
convergence: {deltae: 7.105E-15, res2: 2.139E-19, residm: 1.321E-29, diffor: null, }
|
|
etotal : -8.72190499E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : 2.53959488E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -6.02910429E-04, -2.25449274E-05, -5.15338774E-05, ]
|
|
- [ -2.25449274E-05, -6.00639463E-04, 1.75404114E-05, ]
|
|
- [ -5.15338774E-05, 1.75404114E-05, -6.02738583E-04, ]
|
|
pressure_GPa: 1.7714E+01
|
|
xred :
|
|
- [ 2.4000E-01, 2.5500E-01, 2.4900E-01, P]
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.56720139E-03, 7.78450885E-03, 3.34988677E-03, ]
|
|
- [ 2.56720139E-03, -7.78450885E-03, -3.34988677E-03, ]
|
|
force_length_stats: {min: 8.85498969E-03, max: 8.85498969E-03, mean: 8.85498969E-03, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 1.90690 2.68246148
|
|
2 1.90363 0.77420350
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = -0.524738057081647
|
|
Compensation charge over fine fft grid = -0.524720184220758
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1
|
|
1.36052 0.01008 -0.00007 -0.00003 0.00002 -0.00009 -0.00004 0.00003
|
|
0.01008 28.05725 -0.00011 -0.00004 0.00003 -0.00248 -0.00108 0.00083
|
|
-0.00007 -0.00011 0.26741 -0.00000 0.00000 0.01302 -0.00002 0.00003
|
|
-0.00003 -0.00004 -0.00000 0.26741 0.00001 -0.00002 0.01302 0.00007
|
|
0.00002 0.00003 0.00000 0.00001 0.26741 0.00003 0.00007 0.01302
|
|
-0.00009 -0.00248 0.01302 -0.00002 0.00003 17.66220 -0.00030 0.00042
|
|
-0.00004 -0.00108 -0.00002 0.01302 0.00007 -0.00030 17.66219 0.00107
|
|
0.00003 0.00083 0.00003 0.00007 0.01302 0.00042 0.00107 17.66218
|
|
Atom # 2
|
|
0.32629 -0.05739 0.00003 0.00001 -0.00001 -0.00024 -0.00010 0.00008
|
|
-0.05739 39.68611 -0.00018 -0.00008 0.00006 0.00698 0.00304 -0.00235
|
|
0.00003 -0.00018 0.07346 -0.00000 0.00000 -0.15662 0.00002 -0.00002
|
|
0.00001 -0.00008 -0.00000 0.07346 0.00001 0.00002 -0.15662 -0.00006
|
|
-0.00001 0.00006 0.00000 0.00001 0.07346 -0.00002 -0.00006 -0.15662
|
|
-0.00024 0.00698 -0.15662 0.00002 -0.00002 24.40602 -0.00183 0.00230
|
|
-0.00010 0.00304 0.00002 -0.15662 -0.00006 -0.00183 24.40618 0.00512
|
|
0.00008 -0.00235 -0.00002 -0.00006 -0.15662 0.00230 0.00512 24.40619
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1
|
|
1.50794 0.00728 -0.00593 -0.00262 0.00206 -0.00006 -0.00003 0.00002
|
|
0.00728 0.00004 -0.00001 -0.00000 -0.00000 -0.00000 -0.00000 0.00000
|
|
-0.00593 -0.00001 1.31374 0.00020 -0.00021 0.00718 -0.00002 0.00002
|
|
-0.00262 -0.00000 0.00020 1.31382 -0.00031 -0.00002 0.00718 0.00005
|
|
0.00206 -0.00000 -0.00021 -0.00031 1.31383 0.00002 0.00005 0.00718
|
|
-0.00006 -0.00000 0.00718 -0.00002 0.00002 0.00004 -0.00000 0.00000
|
|
-0.00003 -0.00000 -0.00002 0.00718 0.00005 -0.00000 0.00004 0.00000
|
|
0.00002 0.00000 0.00002 0.00005 0.00718 0.00000 0.00000 0.00004
|
|
Atom # 2
|
|
1.03892 -0.00154 0.01727 0.00749 -0.00579 -0.00001 -0.00000 0.00000
|
|
-0.00154 0.00000 0.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
0.01727 0.00000 0.99942 -0.01315 0.01656 -0.00073 0.00000 -0.00000
|
|
0.00749 -0.00000 -0.01315 1.00049 0.03683 0.00000 -0.00073 -0.00000
|
|
-0.00579 0.00000 0.01656 0.03683 1.00059 -0.00000 -0.00000 -0.00073
|
|
-0.00001 0.00000 -0.00073 0.00000 -0.00000 0.00000 -0.00000 0.00000
|
|
-0.00000 0.00000 0.00000 -0.00073 -0.00000 -0.00000 0.00000 0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00073 0.00000 0.00000 0.00000
|
|
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 12.156E-31; max= 13.212E-30
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.240000000000 0.255000000000 0.249000000000
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 3.5569E-02; max dE/dt= 5.5266E-02; dE/dt below (all hartree)
|
|
1 -0.056077754318 -0.003943883782 -0.026236800633
|
|
2 0.055266201901 0.003882970028 0.025936273901
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 1.33352656564680 1.29383827500255 1.30971359126025
|
|
2 0.00000000000000 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00256720139377 0.00778450884713 0.00334988677474
|
|
2 0.00256720139377 -0.00778450884713 -0.00334988677474
|
|
frms,max,avg= 5.1124307E-03 7.7845088E-03 -2.251E-05 5.256E-05 2.860E-05 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.13201079229773 0.40029550585735 0.17225808941733
|
|
2 0.13201079229773 -0.40029550585735 -0.17225808941733
|
|
frms,max,avg= 2.6289173E-01 4.0029551E-01 -1.157E-03 2.703E-03 1.471E-03 e/A
|
|
length scales= 10.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 5.291772085900 5.291772085900 5.291772085900 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.25396 Average Vxc (hartree)= -0.34016
|
|
Eigenvalues (hartree) for nkpt= 8 k points:
|
|
kpt# 1, nband= 8, wtk= 0.12500, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.19950 0.24708 0.24941 0.25396 0.39574 0.41680 0.41829 0.41922
|
|
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 : 3.88297384962404E+00
|
|
hartree : 9.47041924568268E-01
|
|
xc : -3.18572948868150E+00
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
local_psp : -2.96817232810210E+00
|
|
spherical_terms : 1.70505176945914E+00
|
|
total_energy : -8.72190499450735E+00
|
|
total_energy_eV : -2.37335104795991E+02
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 2, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : 6.63514810037876E-01
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
xc_dc : -3.61188808905360E-01
|
|
spherical_terms : 7.88397258089488E-02
|
|
total_energy_dc : -8.72190499443373E+00
|
|
total_energy_dc_eV : -2.37335104793988E+02
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449274E-05
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= 1.7714E+01 GPa]
|
|
- sigma(1 1)= -1.77382343E+01 sigma(3 2)= 5.16056633E-01
|
|
- sigma(2 2)= -1.76714201E+01 sigma(3 1)= -1.51617876E+00
|
|
- sigma(3 3)= -1.77331784E+01 sigma(2 1)= -6.63294554E-01
|
|
|
|
================================================================================
|
|
== DATASET 3 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 3, }
|
|
dimensions: {natom: 2, nkpt: 8, mband: 8, nsppol: 1, nspinor: 1, nspden: 1, mpw: 138, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 10.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 0.0000000 5.0000000 5.0000000 G(1)= -0.1000000 0.1000000 0.1000000
|
|
R(2)= 5.0000000 0.0000000 5.0000000 G(2)= 0.1000000 -0.1000000 0.1000000
|
|
R(3)= 5.0000000 5.0000000 0.0000000 G(3)= 0.1000000 0.1000000 -0.1000000
|
|
Unit cell volume ucvol= 2.5000000E+02 bohr^3
|
|
Angles (23,13,12)= 6.00000000E+01 6.00000000E+01 6.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 16
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.24794
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 6.316547 Hartrees makes boxcut=2
|
|
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 24 24 24
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.38430
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 14.212230 Hartrees makes boxcut=2
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 133.375 133.250
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 3, }
|
|
solver: {iscf: 17, nstep: 15, nline: 8, wfoptalg: 114, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -8.7465072911124 -8.747E+00 3.385E-17 1.449E+00
|
|
ETOT 2 -8.7227908927729 2.372E-02 7.748E-21 1.863E-01
|
|
ETOT 3 -8.7214810091268 1.310E-03 7.363E-15 1.264E-02
|
|
ETOT 4 -8.7217968845093 -3.159E-04 2.556E-16 9.961E-04
|
|
ETOT 5 -8.7219026250630 -1.057E-04 7.931E-17 1.770E-05
|
|
ETOT 6 -8.7219049761928 -2.351E-06 2.275E-17 2.395E-07
|
|
ETOT 7 -8.7219049907517 -1.456E-08 7.272E-20 3.053E-08
|
|
ETOT 8 -8.7219049940110 -3.259E-09 4.995E-21 2.726E-09
|
|
ETOT 9 -8.7219049944345 -4.234E-10 9.444E-21 1.074E-10
|
|
ETOT 10 -8.7219049944337 7.194E-13 6.308E-21 6.207E-12
|
|
ETOT 11 -8.7219049944337 3.553E-15 7.634E-21 2.194E-14
|
|
ETOT 12 -8.7219049944337 7.638E-14 4.209E-21 2.554E-15
|
|
ETOT 13 -8.7219049944337 -6.750E-14 7.965E-21 5.454E-17
|
|
ETOT 14 -8.7219049944337 4.441E-14 7.650E-21 5.924E-18
|
|
ETOT 15 -8.7219049944337 -3.020E-14 3.397E-21 1.661E-19
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449274E-05
|
|
|
|
|
|
scprqt: WARNING -
|
|
nstep= 15 was not enough SCF cycles to converge;
|
|
density residual= 1.661E-19 exceeds tolvrs= 1.000E-20
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 3, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 0.0000000, 5.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 0.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 5.0000000, 0.0000000, ]
|
|
lattice_lengths: [ 7.07107, 7.07107, 7.07107, ]
|
|
lattice_angles: [ 60.000, 60.000, 60.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.5000000E+02
|
|
convergence: {deltae: -3.020E-14, res2: 1.661E-19, residm: 3.397E-21, diffor: null, }
|
|
etotal : -8.72190499E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : 2.53959488E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -6.02910429E-04, -2.25449274E-05, -5.15338774E-05, ]
|
|
- [ -2.25449274E-05, -6.00639463E-04, 1.75404114E-05, ]
|
|
- [ -5.15338774E-05, 1.75404114E-05, -6.02738583E-04, ]
|
|
pressure_GPa: 1.7714E+01
|
|
xred :
|
|
- [ 2.4000E-01, 2.5500E-01, 2.4900E-01, P]
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.56720140E-03, 7.78450885E-03, 3.34988678E-03, ]
|
|
- [ 2.56720140E-03, -7.78450885E-03, -3.34988678E-03, ]
|
|
force_length_stats: {min: 8.85498970E-03, max: 8.85498970E-03, mean: 8.85498970E-03, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 1.90690 2.68246148
|
|
2 1.90363 0.77420350
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = -0.524738056987817
|
|
Compensation charge over fine fft grid = -0.524720184218282
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1
|
|
1.36052 0.01008 -0.00007 -0.00003 0.00002 -0.00009 -0.00004 0.00003
|
|
0.01008 28.05725 -0.00011 -0.00004 0.00003 -0.00248 -0.00108 0.00083
|
|
-0.00007 -0.00011 0.26741 -0.00000 0.00000 0.01302 -0.00002 0.00003
|
|
-0.00003 -0.00004 -0.00000 0.26741 0.00001 -0.00002 0.01302 0.00007
|
|
0.00002 0.00003 0.00000 0.00001 0.26741 0.00003 0.00007 0.01302
|
|
-0.00009 -0.00248 0.01302 -0.00002 0.00003 17.66220 -0.00030 0.00042
|
|
-0.00004 -0.00108 -0.00002 0.01302 0.00007 -0.00030 17.66219 0.00107
|
|
0.00003 0.00083 0.00003 0.00007 0.01302 0.00042 0.00107 17.66218
|
|
Atom # 2
|
|
0.32629 -0.05739 0.00003 0.00001 -0.00001 -0.00024 -0.00010 0.00008
|
|
-0.05739 39.68611 -0.00018 -0.00008 0.00006 0.00698 0.00304 -0.00235
|
|
0.00003 -0.00018 0.07346 -0.00000 0.00000 -0.15662 0.00002 -0.00002
|
|
0.00001 -0.00008 -0.00000 0.07346 0.00001 0.00002 -0.15662 -0.00006
|
|
-0.00001 0.00006 0.00000 0.00001 0.07346 -0.00002 -0.00006 -0.15662
|
|
-0.00024 0.00698 -0.15662 0.00002 -0.00002 24.40602 -0.00183 0.00230
|
|
-0.00010 0.00304 0.00002 -0.15662 -0.00006 -0.00183 24.40618 0.00512
|
|
0.00008 -0.00235 -0.00002 -0.00006 -0.15662 0.00230 0.00512 24.40619
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1
|
|
1.50794 0.00728 -0.00593 -0.00262 0.00206 -0.00006 -0.00003 0.00002
|
|
0.00728 0.00004 -0.00001 -0.00000 -0.00000 -0.00000 -0.00000 0.00000
|
|
-0.00593 -0.00001 1.31374 0.00020 -0.00021 0.00718 -0.00002 0.00002
|
|
-0.00262 -0.00000 0.00020 1.31382 -0.00031 -0.00002 0.00718 0.00005
|
|
0.00206 -0.00000 -0.00021 -0.00031 1.31383 0.00002 0.00005 0.00718
|
|
-0.00006 -0.00000 0.00718 -0.00002 0.00002 0.00004 -0.00000 0.00000
|
|
-0.00003 -0.00000 -0.00002 0.00718 0.00005 -0.00000 0.00004 0.00000
|
|
0.00002 0.00000 0.00002 0.00005 0.00718 0.00000 0.00000 0.00004
|
|
Atom # 2
|
|
1.03892 -0.00154 0.01727 0.00749 -0.00579 -0.00001 -0.00000 0.00000
|
|
-0.00154 0.00000 0.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
0.01727 0.00000 0.99942 -0.01315 0.01656 -0.00073 0.00000 -0.00000
|
|
0.00749 -0.00000 -0.01315 1.00049 0.03683 0.00000 -0.00073 -0.00000
|
|
-0.00579 0.00000 0.01656 0.03683 1.00059 -0.00000 -0.00000 -0.00073
|
|
-0.00001 0.00000 -0.00073 0.00000 -0.00000 0.00000 -0.00000 0.00000
|
|
-0.00000 0.00000 0.00000 -0.00073 -0.00000 -0.00000 0.00000 0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00073 0.00000 0.00000 0.00000
|
|
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 68.655E-23; max= 33.970E-22
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.240000000000 0.255000000000 0.249000000000
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 3.5569E-02; max dE/dt= 5.5266E-02; dE/dt below (all hartree)
|
|
1 -0.056077754453 -0.003943883792 -0.026236800694
|
|
2 0.055266201870 0.003882970026 0.025936273887
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 1.33352656564680 1.29383827500255 1.30971359126025
|
|
2 0.00000000000000 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00256720139619 0.00778450885431 0.00334988677799
|
|
2 0.00256720139619 -0.00778450885431 -0.00334988677799
|
|
frms,max,avg= 5.1124307E-03 7.7845089E-03 -2.251E-05 5.256E-05 2.860E-05 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.13201079242206 0.40029550622667 0.17225808958469
|
|
2 0.13201079242206 -0.40029550622667 -0.17225808958469
|
|
frms,max,avg= 2.6289173E-01 4.0029551E-01 -1.157E-03 2.703E-03 1.471E-03 e/A
|
|
length scales= 10.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 5.291772085900 5.291772085900 5.291772085900 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.25396 Average Vxc (hartree)= -0.34016
|
|
Eigenvalues (hartree) for nkpt= 8 k points:
|
|
kpt# 1, nband= 8, wtk= 0.12500, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.19950 0.24708 0.24941 0.25396 0.39574 0.41680 0.41829 0.41922
|
|
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 : 3.88297384962551E+00
|
|
hartree : 9.47041924568493E-01
|
|
xc : -3.18572948863175E+00
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
local_psp : -2.96817232809363E+00
|
|
spherical_terms : 1.70505176924363E+00
|
|
total_energy : -8.72190499466295E+00
|
|
total_energy_eV : -2.37335104800225E+02
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 3, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : 6.63514810037323E-01
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
xc_dc : -3.61188808916867E-01
|
|
spherical_terms : 7.88397258210294E-02
|
|
total_energy_dc : -8.72190499443371E+00
|
|
total_energy_dc_eV : -2.37335104793987E+02
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449274E-05
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= 1.7714E+01 GPa]
|
|
- sigma(1 1)= -1.77382343E+01 sigma(3 2)= 5.16056633E-01
|
|
- sigma(2 2)= -1.76714201E+01 sigma(3 1)= -1.51617876E+00
|
|
- sigma(3 3)= -1.77331784E+01 sigma(2 1)= -6.63294554E-01
|
|
|
|
================================================================================
|
|
== DATASET 4 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 4, }
|
|
dimensions: {natom: 2, nkpt: 8, mband: 8, nsppol: 1, nspinor: 1, nspden: 1, mpw: 138, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 10.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 0.0000000 5.0000000 5.0000000 G(1)= -0.1000000 0.1000000 0.1000000
|
|
R(2)= 5.0000000 0.0000000 5.0000000 G(2)= 0.1000000 -0.1000000 0.1000000
|
|
R(3)= 5.0000000 5.0000000 0.0000000 G(3)= 0.1000000 0.1000000 -0.1000000
|
|
Unit cell volume ucvol= 2.5000000E+02 bohr^3
|
|
Angles (23,13,12)= 6.00000000E+01 6.00000000E+01 6.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 16
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.24794
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 6.316547 Hartrees makes boxcut=2
|
|
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 24 24 24
|
|
ecut(hartree)= 10.000 => boxcut(ratio)= 2.38430
|
|
|
|
getcut : COMMENT -
|
|
Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2
|
|
is sufficient for exact treatment of convolution.
|
|
Such a large boxcut is a waste : you could raise ecut
|
|
e.g. ecut= 14.212230 Hartrees makes boxcut=2
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 133.375 133.250
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 4, }
|
|
solver: {iscf: 17, nstep: 15, nline: 8, wfoptalg: 111, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -8.7465038108563 -8.747E+00 5.436E-04 1.448E+00
|
|
ETOT 2 -8.7227954860066 2.371E-02 3.377E-05 1.864E-01
|
|
ETOT 3 -8.7214807082120 1.315E-03 2.774E-05 1.260E-02
|
|
ETOT 4 -8.7217971631571 -3.165E-04 1.629E-05 9.951E-04
|
|
ETOT 5 -8.7219026986909 -1.055E-04 1.301E-05 1.735E-05
|
|
ETOT 6 -8.7219049783788 -2.280E-06 8.753E-06 2.850E-07
|
|
ETOT 7 -8.7219049901113 -1.173E-08 7.108E-06 3.323E-08
|
|
ETOT 8 -8.7219049940364 -3.925E-09 5.230E-06 2.445E-09
|
|
ETOT 9 -8.7219049944302 -3.938E-10 4.320E-06 1.006E-10
|
|
ETOT 10 -8.7219049944340 -3.794E-12 3.368E-06 6.014E-12
|
|
ETOT 11 -8.7219049944337 3.020E-13 2.807E-06 2.120E-14
|
|
ETOT 12 -8.7219049944337 1.776E-14 2.262E-06 2.728E-15
|
|
ETOT 13 -8.7219049944338 -4.796E-14 1.898E-06 3.076E-17
|
|
ETOT 14 -8.7219049944337 4.796E-14 1.558E-06 2.543E-18
|
|
ETOT 15 -8.7219049944337 -1.776E-14 1.315E-06 1.038E-19
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449275E-05
|
|
|
|
|
|
scprqt: WARNING -
|
|
nstep= 15 was not enough SCF cycles to converge;
|
|
density residual= 1.038E-19 exceeds tolvrs= 1.000E-20
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 4, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 0.0000000, 5.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 0.0000000, 5.0000000, ]
|
|
- [ 5.0000000, 5.0000000, 0.0000000, ]
|
|
lattice_lengths: [ 7.07107, 7.07107, 7.07107, ]
|
|
lattice_angles: [ 60.000, 60.000, 60.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 2.5000000E+02
|
|
convergence: {deltae: -1.776E-14, res2: 1.038E-19, residm: 1.315E-06, diffor: null, }
|
|
etotal : -8.72190499E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : 2.53959488E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -6.02910429E-04, -2.25449275E-05, -5.15338774E-05, ]
|
|
- [ -2.25449275E-05, -6.00639463E-04, 1.75404114E-05, ]
|
|
- [ -5.15338774E-05, 1.75404114E-05, -6.02738583E-04, ]
|
|
pressure_GPa: 1.7714E+01
|
|
xred :
|
|
- [ 2.4000E-01, 2.5500E-01, 2.4900E-01, P]
|
|
- [ 0.0000E+00, 0.0000E+00, 0.0000E+00, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.56720139E-03, 7.78450885E-03, 3.34988678E-03, ]
|
|
- [ 2.56720139E-03, -7.78450885E-03, -3.34988678E-03, ]
|
|
force_length_stats: {min: 8.85498969E-03, max: 8.85498969E-03, mean: 8.85498969E-03, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 1.90690 2.68246148
|
|
2 1.90363 0.77420350
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = -0.524738056937835
|
|
Compensation charge over fine fft grid = -0.524720184229335
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
Atom # 1
|
|
1.36052 0.01008 -0.00007 -0.00003 0.00002 -0.00009 -0.00004 0.00003
|
|
0.01008 28.05725 -0.00011 -0.00004 0.00003 -0.00248 -0.00108 0.00083
|
|
-0.00007 -0.00011 0.26741 -0.00000 0.00000 0.01302 -0.00002 0.00003
|
|
-0.00003 -0.00004 -0.00000 0.26741 0.00001 -0.00002 0.01302 0.00007
|
|
0.00002 0.00003 0.00000 0.00001 0.26741 0.00003 0.00007 0.01302
|
|
-0.00009 -0.00248 0.01302 -0.00002 0.00003 17.66220 -0.00030 0.00042
|
|
-0.00004 -0.00108 -0.00002 0.01302 0.00007 -0.00030 17.66219 0.00107
|
|
0.00003 0.00083 0.00003 0.00007 0.01302 0.00042 0.00107 17.66218
|
|
Atom # 2
|
|
0.32629 -0.05739 0.00003 0.00001 -0.00001 -0.00024 -0.00010 0.00008
|
|
-0.05739 39.68611 -0.00018 -0.00008 0.00006 0.00698 0.00304 -0.00235
|
|
0.00003 -0.00018 0.07346 -0.00000 0.00000 -0.15662 0.00002 -0.00002
|
|
0.00001 -0.00008 -0.00000 0.07346 0.00001 0.00002 -0.15662 -0.00006
|
|
-0.00001 0.00006 0.00000 0.00001 0.07346 -0.00002 -0.00006 -0.15662
|
|
-0.00024 0.00698 -0.15662 0.00002 -0.00002 24.40602 -0.00183 0.00230
|
|
-0.00010 0.00304 0.00002 -0.15662 -0.00006 -0.00183 24.40618 0.00512
|
|
0.00008 -0.00235 -0.00002 -0.00006 -0.15662 0.00230 0.00512 24.40619
|
|
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
Atom # 1
|
|
1.50794 0.00728 -0.00593 -0.00262 0.00206 -0.00006 -0.00003 0.00002
|
|
0.00728 0.00004 -0.00001 -0.00000 -0.00000 -0.00000 -0.00000 0.00000
|
|
-0.00593 -0.00001 1.31374 0.00020 -0.00021 0.00718 -0.00002 0.00002
|
|
-0.00262 -0.00000 0.00020 1.31382 -0.00031 -0.00002 0.00718 0.00005
|
|
0.00206 -0.00000 -0.00021 -0.00031 1.31383 0.00002 0.00005 0.00718
|
|
-0.00006 -0.00000 0.00718 -0.00002 0.00002 0.00004 -0.00000 0.00000
|
|
-0.00003 -0.00000 -0.00002 0.00718 0.00005 -0.00000 0.00004 0.00000
|
|
0.00002 0.00000 0.00002 0.00005 0.00718 0.00000 0.00000 0.00004
|
|
Atom # 2
|
|
1.03892 -0.00154 0.01727 0.00749 -0.00579 -0.00001 -0.00000 0.00000
|
|
-0.00154 0.00000 0.00000 -0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
0.01727 0.00000 0.99942 -0.01315 0.01656 -0.00073 0.00000 -0.00000
|
|
0.00749 -0.00000 -0.01315 1.00049 0.03683 0.00000 -0.00073 -0.00000
|
|
-0.00579 0.00000 0.01656 0.03683 1.00059 -0.00000 -0.00000 -0.00073
|
|
-0.00001 0.00000 -0.00073 0.00000 -0.00000 0.00000 -0.00000 0.00000
|
|
-0.00000 0.00000 0.00000 -0.00073 -0.00000 -0.00000 0.00000 0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00073 0.00000 0.00000 0.00000
|
|
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 27.401E-09; max= 13.152E-07
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.240000000000 0.255000000000 0.249000000000
|
|
0.000000000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 3.5569E-02; max dE/dt= 5.5266E-02; dE/dt below (all hartree)
|
|
1 -0.056077754325 -0.003943883781 -0.026236800633
|
|
2 0.055266201912 0.003882970027 0.025936273902
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 1.33352656564680 1.29383827500255 1.30971359126025
|
|
2 0.00000000000000 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00256720139472 0.00778450884819 0.00334988677551
|
|
2 0.00256720139472 -0.00778450884819 -0.00334988677551
|
|
frms,max,avg= 5.1124307E-03 7.7845088E-03 -2.251E-05 5.256E-05 2.860E-05 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.13201079234628 0.40029550591208 0.17225808945698
|
|
2 0.13201079234628 -0.40029550591208 -0.17225808945698
|
|
frms,max,avg= 2.6289173E-01 4.0029551E-01 -1.157E-03 2.703E-03 1.471E-03 e/A
|
|
length scales= 10.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 5.291772085900 5.291772085900 5.291772085900 angstroms
|
|
Fermi (or HOMO) energy (hartree) = 0.25396 Average Vxc (hartree)= -0.34016
|
|
Eigenvalues (hartree) for nkpt= 8 k points:
|
|
kpt# 1, nband= 8, wtk= 0.12500, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.19950 0.24708 0.24941 0.25396 0.39574 0.41681 0.41838 0.41992
|
|
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 : 3.88297384965792E+00
|
|
hartree : 9.47041924555123E-01
|
|
xc : -3.18572948860555E+00
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
local_psp : -2.96817232805007E+00
|
|
spherical_terms : 1.70505176916819E+00
|
|
total_energy : -8.72190499464960E+00
|
|
total_energy_eV : -2.37335104799861E+02
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 4, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : 6.63514810025726E-01
|
|
Ewald energy : -8.99667321169082E+00
|
|
psp_core : -1.06397509684384E-01
|
|
xc_dc : -3.61188808912380E-01
|
|
spherical_terms : 7.88397258281326E-02
|
|
total_energy_dc : -8.72190499443372E+00
|
|
total_energy_dc_eV : -2.37335104793987E+02
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -6.02910429E-04 sigma(3 2)= 1.75404114E-05
|
|
sigma(2 2)= -6.00639463E-04 sigma(3 1)= -5.15338774E-05
|
|
sigma(3 3)= -6.02738583E-04 sigma(2 1)= -2.25449275E-05
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= 1.7714E+01 GPa]
|
|
- sigma(1 1)= -1.77382343E+01 sigma(3 2)= 5.16056633E-01
|
|
- sigma(2 2)= -1.76714201E+01 sigma(3 1)= -1.51617876E+00
|
|
- sigma(3 3)= -1.77331784E+01 sigma(2 1)= -6.63294554E-01
|
|
|
|
== END DATASET(S) ==============================================================
|
|
================================================================================
|
|
|
|
-outvars: echo values of variables after computation --------
|
|
acell 1.0000000000E+01 1.0000000000E+01 1.0000000000E+01 Bohr
|
|
amu 3.09737620E+01 2.69815390E+01
|
|
bandpp1 1
|
|
bandpp2 8
|
|
bandpp3 8
|
|
bandpp4 8
|
|
ecut 5.00000000E+00 Hartree
|
|
etotal1 -8.7219049944E+00
|
|
etotal2 -8.7219049944E+00
|
|
etotal3 -8.7219049944E+00
|
|
etotal4 -8.7219049944E+00
|
|
fcart1 -2.5672015284E-03 7.7845088396E-03 3.3498865075E-03
|
|
2.5672015284E-03 -7.7845088396E-03 -3.3498865075E-03
|
|
fcart2 -2.5672013938E-03 7.7845088471E-03 3.3498867747E-03
|
|
2.5672013938E-03 -7.7845088471E-03 -3.3498867747E-03
|
|
fcart3 -2.5672013962E-03 7.7845088543E-03 3.3498867780E-03
|
|
2.5672013962E-03 -7.7845088543E-03 -3.3498867780E-03
|
|
fcart4 -2.5672013947E-03 7.7845088482E-03 3.3498867755E-03
|
|
2.5672013947E-03 -7.7845088482E-03 -3.3498867755E-03
|
|
- fftalg 512
|
|
istwfk 1 1 1 1 1 1 1 1
|
|
ixc -101130
|
|
jdtset 1 2 3 4
|
|
kpt 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
5.00000000E-01 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 5.00000000E-01 0.00000000E+00
|
|
5.00000000E-01 5.00000000E-01 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
5.00000000E-01 0.00000000E+00 5.00000000E-01
|
|
0.00000000E+00 5.00000000E-01 5.00000000E-01
|
|
5.00000000E-01 5.00000000E-01 5.00000000E-01
|
|
kptopt 3
|
|
kptrlatt 2 0 0 0 2 0 0 0 2
|
|
kptrlen 1.41421356E+01
|
|
P mkmem 8
|
|
natom 2
|
|
nband 8
|
|
nbdbuf 2
|
|
ndtset 4
|
|
ngfft 16 16 16
|
|
ngfftdg 24 24 24
|
|
nkpt 8
|
|
nline 8
|
|
nstep 15
|
|
nsym 1
|
|
ntypat 2
|
|
occ 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000
|
|
0.000000 0.000000
|
|
pawecutdg 1.00000000E+01 Hartree
|
|
prtden 0
|
|
prteig 0
|
|
prtwf 0
|
|
rprim 0.0000000000E+00 5.0000000000E-01 5.0000000000E-01
|
|
5.0000000000E-01 0.0000000000E+00 5.0000000000E-01
|
|
5.0000000000E-01 5.0000000000E-01 0.0000000000E+00
|
|
spgroup 1
|
|
strten1 -6.0291045063E-04 -6.0063948326E-04 -6.0273860346E-04
|
|
1.7540411438E-05 -5.1533877001E-05 -2.2544926224E-05
|
|
strten2 -6.0291042928E-04 -6.0063946271E-04 -6.0273858340E-04
|
|
1.7540411389E-05 -5.1533877419E-05 -2.2544927447E-05
|
|
strten3 -6.0291042926E-04 -6.0063946269E-04 -6.0273858337E-04
|
|
1.7540411384E-05 -5.1533877409E-05 -2.2544927442E-05
|
|
strten4 -6.0291042922E-04 -6.0063946265E-04 -6.0273858333E-04
|
|
1.7540411392E-05 -5.1533877424E-05 -2.2544927452E-05
|
|
tolvrs 1.00000000E-20
|
|
typat 1 2
|
|
useylm 1
|
|
wfoptalg1 0
|
|
wfoptalg2 14
|
|
wfoptalg3 114
|
|
wfoptalg4 111
|
|
wtk 0.12500 0.12500 0.12500 0.12500 0.12500 0.12500
|
|
0.12500 0.12500
|
|
xangst 1.3335265656E+00 1.2938382750E+00 1.3097135913E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
xcart 2.5200000000E+00 2.4450000000E+00 2.4750000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
xred 2.4000000000E-01 2.5500000000E-01 2.4900000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
znucl 15.00000 13.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] Parallel eigensolvers in plane-wave Density Functional Theory
|
|
- A. Levitt and M. Torrent, Computer Phys. Comm. 187, 98-105 (2015).
|
|
- Comment: in case Chebyshev Filtering algorithm is used (wfoptalg=1/111).
|
|
- Strong suggestion to cite this paper in your publications.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#levitt2015
|
|
-
|
|
- [2] 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
|
|
-
|
|
- [3] 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
|
|
-
|
|
- [4] 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
|
|
-
|
|
- [5] 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
|
|
-
|
|
- [6] 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= 5.0 wall= 5.0
|
|
|
|
================================================================================
|
|
|
|
Calculation completed.
|
|
.Delivered 4 WARNINGs and 33 COMMENTs to log file.
|
|
+Overall time at end (sec) : cpu= 5.0 wall= 5.0
|