mirror of https://github.com/abinit/abinit.git
1500 lines
72 KiB
Plaintext
1500 lines
72 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 19h08 )
|
|
|
|
- input file -> /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/TestBot_MPI1/v2_t36/t36.abi
|
|
- output file -> t36.abo
|
|
- root for input files -> t36i
|
|
- root for output files -> t36o
|
|
|
|
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 = 7 lmnmax = 3
|
|
lnmax = 3 mgfft = 24 mpssoang = 3 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 0 ntypat = 1
|
|
occopt = 4 xclevel = 1
|
|
- mband = 5 mffmem = 1 mkmem = 2
|
|
mpw = 182 nfft = 6144 nkpt = 2
|
|
================================================================================
|
|
P This job should need less than 2.978 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.030 Mbytes ; DEN or POT disk file : 0.049 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 = -2 lmnmax = 3
|
|
lnmax = 3 mgfft = 24 mpssoang = 3 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 0 ntypat = 1
|
|
occopt = 4 xclevel = 1
|
|
- mband = 5 mffmem = 1 mkmem = 2
|
|
mpw = 359 nfft = 6144 nkpt = 2
|
|
================================================================================
|
|
P This job should need less than 2.287 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.057 Mbytes ; DEN or POT disk file : 0.049 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 3 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 3 (RF).
|
|
intxc = 0 iscf = 7 lmnmax = 3 lnmax = 3
|
|
mgfft = 24 mpssoang = 3 mqgrid = 3001 natom = 2
|
|
nloc_mem = 1 nspden = 1 nspinor = 1 nsppol = 1
|
|
nsym = 1 n1xccc = 0 ntypat = 1 occopt = 4
|
|
xclevel = 1
|
|
- mband = 5 mffmem = 1 mkmem = 2
|
|
- mkqmem = 2 mk1mem = 2 mpw = 364
|
|
nfft = 6144 nkpt = 2
|
|
================================================================================
|
|
P This job should need less than 3.130 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.058 Mbytes ; DEN or POT disk file : 0.049 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 = 7 lmnmax = 3
|
|
lnmax = 3 mgfft = 48 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 0 ntypat = 1
|
|
occopt = 4 xclevel = 1
|
|
- mband = 10 mffmem = 1 mkmem = 1
|
|
mpw = 362 nfft = 12288 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 4.775 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.057 Mbytes ; DEN or POT disk file : 0.096 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 5 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 5.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 3
|
|
lnmax = 3 mgfft = 48 mpssoang = 3 mqgrid = 3001
|
|
natom = 4 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 0 ntypat = 1
|
|
occopt = 4 xclevel = 1
|
|
- mband = 10 mffmem = 1 mkmem = 1
|
|
mpw = 362 nfft = 12288 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 4.775 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.057 Mbytes ; DEN or POT disk file : 0.096 Mbytes.
|
|
================================================================================
|
|
|
|
--------------------------------------------------------------------------------
|
|
------------- Echo of variables that govern the present computation ------------
|
|
--------------------------------------------------------------------------------
|
|
-
|
|
- outvars: echo of selected default values
|
|
- iomode0 = 0 , fftalg0 =512 , wfoptalg0 = 0
|
|
-
|
|
- outvars: echo of global parameters not present in the input file
|
|
- max_nthreads = 0
|
|
-
|
|
-outvars: echo values of preprocessed input variables --------
|
|
acell1 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell2 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell3 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell4 1.1000000000E+01 1.1000000000E+01 3.2000000000E+01 Bohr
|
|
acell5 1.1000000000E+01 1.1000000000E+01 3.2000000000E+01 Bohr
|
|
amu 1.00000000E+00
|
|
asr 0
|
|
chneut 0
|
|
diemac 1.00000000E+00
|
|
diemix1 3.33333333E-01
|
|
diemix2 3.33333333E-01
|
|
diemix3 3.50000000E-01
|
|
diemix4 3.33333333E-01
|
|
diemix5 3.33333333E-01
|
|
ecut 2.50000000E+00 Hartree
|
|
- fftalg 512
|
|
getden1 0
|
|
getden2 1
|
|
getden3 0
|
|
getden4 0
|
|
getden5 0
|
|
getwfk1 0
|
|
getwfk2 1
|
|
getwfk3 1
|
|
getwfk4 0
|
|
getwfk5 4
|
|
getwfq1 0
|
|
getwfq2 0
|
|
getwfq3 2
|
|
getwfq4 0
|
|
getwfq5 0
|
|
iscf1 7
|
|
iscf2 -2
|
|
iscf3 7
|
|
iscf4 7
|
|
iscf5 7
|
|
istwfk1 2 4
|
|
istwfk2 4 0
|
|
istwfk3 1 1
|
|
istwfk4 2
|
|
istwfk5 2
|
|
jdtset 1 2 3 4 5
|
|
kpt1 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt3 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt4 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
kpt5 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
kptopt 0
|
|
P mkmem1 2
|
|
P mkmem2 2
|
|
P mkmem3 2
|
|
P mkmem4 1
|
|
P mkmem5 1
|
|
P mkqmem1 2
|
|
P mkqmem2 2
|
|
P mkqmem3 2
|
|
P mkqmem4 1
|
|
P mkqmem5 1
|
|
P mk1mem1 2
|
|
P mk1mem2 2
|
|
P mk1mem3 2
|
|
P mk1mem4 1
|
|
P mk1mem5 1
|
|
natom1 2
|
|
natom2 2
|
|
natom3 2
|
|
natom4 4
|
|
natom5 4
|
|
nband1 5
|
|
nband2 5
|
|
nband3 5
|
|
nband4 10
|
|
nband5 10
|
|
ndtset 5
|
|
ngfft1 16 16 24
|
|
ngfft2 16 16 24
|
|
ngfft3 16 16 24
|
|
ngfft4 16 16 48
|
|
ngfft5 16 16 48
|
|
nkpt1 2
|
|
nkpt2 2
|
|
nkpt3 2
|
|
nkpt4 1
|
|
nkpt5 1
|
|
nqpt1 0
|
|
nqpt2 1
|
|
nqpt3 1
|
|
nqpt4 0
|
|
nqpt5 0
|
|
nstep1 40
|
|
nstep2 40
|
|
nstep3 10
|
|
nstep4 40
|
|
nstep5 40
|
|
nsym 1
|
|
ntypat 1
|
|
occ1 2.000000 2.000000 2.000000 0.000000 0.000000
|
|
2.000000 2.000000 2.000000 0.000000 0.000000
|
|
occ3 2.000000 2.000000 2.000000 0.000000 0.000000
|
|
2.000000 2.000000 2.000000 0.000000 0.000000
|
|
occ4 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000
|
|
0.000000 0.000000 0.000000 0.000000
|
|
occ5 2.000000 2.000000 2.000000 2.000000 2.000000 2.000000
|
|
0.000000 0.000000 0.000000 0.000000
|
|
occopt 4
|
|
optdriver1 0
|
|
optdriver2 0
|
|
optdriver3 1
|
|
optdriver4 0
|
|
optdriver5 0
|
|
prtpot1 0
|
|
prtpot2 0
|
|
prtpot3 1
|
|
prtpot4 0
|
|
prtpot5 0
|
|
qpt1 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
qpt2 0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
qpt3 0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
qpt4 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
qpt5 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
rfatpol1 1 2
|
|
rfatpol2 1 2
|
|
rfatpol3 2 2
|
|
rfatpol4 1 4
|
|
rfatpol5 1 4
|
|
rfdir1 1 1 1
|
|
rfdir2 1 1 1
|
|
rfdir3 1 0 0
|
|
rfdir4 1 1 1
|
|
rfdir5 1 1 1
|
|
rfphon1 0
|
|
rfphon2 0
|
|
rfphon3 1
|
|
rfphon4 0
|
|
rfphon5 0
|
|
spgroup 1
|
|
toldfe1 0.00000000E+00 Hartree
|
|
toldfe2 0.00000000E+00 Hartree
|
|
toldfe3 1.00000000E-10 Hartree
|
|
toldfe4 0.00000000E+00 Hartree
|
|
toldfe5 0.00000000E+00 Hartree
|
|
tolvrs1 1.00000000E-20
|
|
tolvrs2 0.00000000E+00
|
|
tolvrs3 0.00000000E+00
|
|
tolvrs4 1.00000000E-20
|
|
tolvrs5 1.00000000E-20
|
|
tolwfr1 0.00000000E+00
|
|
tolwfr2 1.00000000E-26
|
|
tolwfr3 0.00000000E+00
|
|
tolwfr4 0.00000000E+00
|
|
tolwfr5 0.00000000E+00
|
|
tsmear 8.00000000E-02 Hartree
|
|
typat1 1 1
|
|
typat2 1 1
|
|
typat3 1 1
|
|
typat4 1 1 1 1
|
|
typat5 1 1 1 1
|
|
wtk1 0.50000 0.50000
|
|
wtk2 1.00000 1.00000
|
|
wtk3 0.50000 0.50000
|
|
wtk4 1.00000
|
|
wtk5 1.00000
|
|
xangst1 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst2 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst3 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst4 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 6.9428049767E+00
|
|
0.0000000000E+00 0.0000000000E+00 9.9908656982E+00
|
|
xangst5 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
5.8209492945E-04 0.0000000000E+00 1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 6.9428049767E+00
|
|
-5.8209492945E-04 0.0000000000E+00 9.9908656982E+00
|
|
xcart1 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart2 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart3 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart4 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.3120000000E+01
|
|
0.0000000000E+00 0.0000000000E+00 1.8880000000E+01
|
|
xcart5 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
1.1000000000E-03 0.0000000000E+00 2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.3120000000E+01
|
|
-1.1000000000E-03 0.0000000000E+00 1.8880000000E+01
|
|
xred1 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred2 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred3 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred4 0.0000000000E+00 0.0000000000E+00 -9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 4.1000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 5.9000000000E-01
|
|
xred5 0.0000000000E+00 0.0000000000E+00 -9.0000000000E-02
|
|
1.0000000000E-04 0.0000000000E+00 9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 4.1000000000E-01
|
|
-1.0000000000E-04 0.0000000000E+00 5.9000000000E-01
|
|
znucl 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.
|
|
|
|
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: 2, nkpt: 2, mband: 5, nsppol: 1, nspinor: 1, nspden: 1, mpw: 182, }
|
|
cutoff_energies: {ecut: 2.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 6.00000000E+00, charge: 0.00000000E+00, occopt: 4.00000000E+00, tsmear: 8.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
LDA: new Teter (4/93) with spin-polarized option - ixc=1
|
|
Citation for XC functional:
|
|
S. Goedecker, M. Teter, J. Huetter, PRB 54, 1703 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 11.0000000 0.0000000 0.0000000 G(1)= 0.0909091 0.0000000 0.0000000
|
|
R(2)= 0.0000000 11.0000000 0.0000000 G(2)= 0.0000000 0.0909091 0.0000000
|
|
R(3)= 0.0000000 0.0000000 16.0000000 G(3)= 0.0000000 0.0000000 0.0625000
|
|
Unit cell volume ucvol= 1.9360000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 24
|
|
ecut(hartree)= 2.500 => boxcut(ratio)= 2.04358
|
|
|
|
--- Pseudopotential description ------------------------------------------------
|
|
- pspini: atom type 1 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/13al.981214.fhi
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/13al.981214.fhi
|
|
- Aluminum, fhi98PP : Hamann-type, LDA CA PerdewWang, l=2 local
|
|
- 13.00000 3.00000 981214 znucl, zion, pspdat
|
|
6 7 2 2 493 0.00000 pspcod,pspxc,lmax,lloc,mmax,r2well
|
|
No XC core correction.
|
|
1.024700 amesh (Hamman grid)
|
|
pspatm : epsatm= 1.36305739
|
|
--- l ekb(1:nproj) -->
|
|
0 1.768744
|
|
1 0.900554
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
1.63566886E+01 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 361.500 361.491
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 1, }
|
|
solver: {iscf: 7, nstep: 40, nline: 4, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -3.9567382561827 -3.957E+00 4.908E-04 7.833E-01
|
|
ETOT 2 -3.9580271277641 -1.289E-03 3.373E-08 3.488E-01
|
|
ETOT 3 -3.9583095693872 -2.824E-04 1.182E-05 1.722E-02
|
|
ETOT 4 -3.9583469096491 -3.734E-05 2.881E-07 4.736E-03
|
|
ETOT 5 -3.9583612514512 -1.434E-05 2.851E-07 9.171E-04
|
|
ETOT 6 -3.9583625560626 -1.305E-06 9.178E-10 1.348E-04
|
|
ETOT 7 -3.9583629078058 -3.517E-07 4.149E-09 4.094E-06
|
|
ETOT 8 -3.9583629163964 -8.591E-09 6.097E-11 8.688E-07
|
|
ETOT 9 -3.9583629176411 -1.245E-09 5.212E-11 2.839E-07
|
|
ETOT 10 -3.9583629182393 -5.982E-10 1.330E-11 2.736E-08
|
|
ETOT 11 -3.9583629183015 -6.219E-11 1.236E-12 1.865E-09
|
|
ETOT 12 -3.9583629183058 -4.348E-12 5.692E-14 1.270E-11
|
|
ETOT 13 -3.9583629183058 -3.109E-15 1.247E-16 4.129E-12
|
|
ETOT 14 -3.9583629183058 -1.910E-14 3.834E-17 1.500E-13
|
|
ETOT 15 -3.9583629183058 7.994E-15 4.285E-18 2.472E-14
|
|
ETOT 16 -3.9583629183058 -7.994E-15 3.984E-19 6.887E-15
|
|
ETOT 17 -3.9583629183058 2.665E-15 1.279E-19 8.664E-16
|
|
ETOT 18 -3.9583629183058 -1.332E-15 3.927E-20 6.282E-17
|
|
ETOT 19 -3.9583629183058 -4.441E-16 3.686E-21 7.119E-18
|
|
ETOT 20 -3.9583629183059 -9.770E-15 8.956E-23 1.315E-18
|
|
ETOT 21 -3.9583629183058 1.688E-14 2.482E-23 1.198E-19
|
|
ETOT 22 -3.9583629183059 -1.554E-14 6.777E-24 1.615E-20
|
|
ETOT 23 -3.9583629183058 1.155E-14 7.478E-25 2.038E-21
|
|
|
|
At SCF step 23 vres2 = 2.04E-21 < tolvrs= 1.00E-20 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543351E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543351E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080290E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 1, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 11.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 11.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 16.0000000, ]
|
|
lattice_lengths: [ 11.00000, 11.00000, 16.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.9360000E+03
|
|
convergence: {deltae: 1.155E-14, res2: 2.038E-21, residm: 7.478E-25, diffor: null, }
|
|
etotal : -3.95836292E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -7.70666970E-02
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 5.63543351E-05, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 5.63543351E-05, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 7.57080290E-05, ]
|
|
pressure_GPa: -1.8478E+00
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, -1.8000E-01, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 1.8000E-01, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ 2.40184982E-14, 4.02301444E-15, 1.04287834E-02, ]
|
|
- [ -2.40184982E-14, -4.02301444E-15, -1.04287834E-02, ]
|
|
force_length_stats: {min: 1.04287834E-02, max: 1.04287834E-02, mean: 1.04287834E-02, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 0.74371586
|
|
2 2.00000 0.74371586
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 49.093E-26; max= 74.778E-26
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.000000000000 0.000000000000 -0.180000000000
|
|
0.000000000000 0.000000000000 0.180000000000
|
|
rms dE/dt= 9.6337E-02; max dE/dt= 1.6686E-01; dE/dt below (all hartree)
|
|
1 -0.000000000000 -0.000000000000 -0.166860533716
|
|
2 0.000000000000 0.000000000000 0.166860533716
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 -1.52403036073920
|
|
2 0.00000000000000 0.00000000000000 1.52403036073920
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 0.00000000000002 0.00000000000000 0.01042878335725
|
|
2 -0.00000000000002 -0.00000000000000 -0.01042878335725
|
|
frms,max,avg= 6.0210609E-03 1.0428783E-02 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 0.00000000000124 0.00000000000021 0.53626955681446
|
|
2 -0.00000000000124 -0.00000000000021 -0.53626955681446
|
|
frms,max,avg= 3.0961537E-01 5.3626956E-01 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 11.000000000000 11.000000000000 16.000000000000 bohr
|
|
= 5.820949294490 5.820949294490 8.466835337440 angstroms
|
|
prteigrs : about to open file t36o_DS1_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.07707 Average Vxc (hartree)= -0.13648
|
|
Eigenvalues (hartree) for nkpt= 2 k points:
|
|
kpt# 1, nband= 5, wtk= 0.50000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.25992 -0.19129 -0.09239 -0.05814 -0.05814
|
|
occupation numbers for kpt# 1
|
|
2.01380 2.12510 1.17566 0.47820 0.47820
|
|
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 : 1.12843373027589E+00
|
|
hartree : 9.62957767235848E-01
|
|
xc : -1.20599282023675E+00
|
|
Ewald energy : -2.26231457505637E+00
|
|
psp_core : 8.44870280674618E-03
|
|
local_psp : -3.62240544085170E+00
|
|
non_local_psp : 1.07072597218701E+00
|
|
internal : -3.92014666363932E+00
|
|
'-kT*entropy' : -3.82162546665197E-02
|
|
total_energy : -3.95836291830584E+00
|
|
total_energy_eV : -1.07712532827210E+02
|
|
band_energy : -1.06649401167774E+00
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543351E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543351E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080290E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -1.8478E+00 GPa]
|
|
- sigma(1 1)= 1.65800150E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 1.65800150E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 2.22740674E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
================================================================================
|
|
== DATASET 2 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 2, }
|
|
dimensions: {natom: 2, nkpt: 2, mband: 5, nsppol: 1, nspinor: 1, nspden: 1, mpw: 359, }
|
|
cutoff_energies: {ecut: 2.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 6.00000000E+00, charge: 0.00000000E+00, occopt: 4.00000000E+00, tsmear: 8.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: -2, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 1.
|
|
|
|
mkfilename : getden/=0, take file _DEN from output of DATASET 1.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
LDA: new Teter (4/93) with spin-polarized option - ixc=1
|
|
Citation for XC functional:
|
|
S. Goedecker, M. Teter, J. Huetter, PRB 54, 1703 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 11.0000000 0.0000000 0.0000000 G(1)= 0.0909091 0.0000000 0.0000000
|
|
R(2)= 0.0000000 11.0000000 0.0000000 G(2)= 0.0000000 0.0909091 0.0000000
|
|
R(3)= 0.0000000 0.0000000 16.0000000 G(3)= 0.0000000 0.0000000 0.0625000
|
|
Unit cell volume ucvol= 1.9360000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 24
|
|
ecut(hartree)= 2.500 => boxcut(ratio)= 2.04358
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t36o_DS1_WFK
|
|
|
|
================================================================================
|
|
prteigrs : about to open file t36o_DS2_EIG
|
|
Non-SCF case, kpt 1 ( 0.00000 0.00000 0.50000), residuals and eigenvalues=
|
|
2.52E-27 4.72E-28 2.24E-27 7.02E-27 7.64E-27
|
|
-2.5822E-01 -1.9808E-01 -7.0471E-02 -5.6327E-02 -5.6327E-02
|
|
Non-SCF case, kpt 2 ( 0.00000 0.00000 1.00000), residuals and eigenvalues=
|
|
1.45E-27 4.95E-28 1.12E-27 4.44E-27 4.59E-27
|
|
-2.5992E-01 -1.9129E-01 -9.2387E-02 -5.8135E-02 -5.8135E-02
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 2, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 11.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 11.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 16.0000000, ]
|
|
lattice_lengths: [ 11.00000, 11.00000, 16.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.9360000E+03
|
|
convergence: {deltae: 0.000E+00, res2: 0.000E+00, residm: 7.644E-27, diffor: 0.000E+00, }
|
|
etotal : -3.95836292E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -7.70666970E-02
|
|
cartesian_stress_tensor: null
|
|
pressure_GPa: null
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, -1.8000E-01, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 1.8000E-01, Al]
|
|
cartesian_forces: null
|
|
force_length_stats: {min: null, max: null, mean: null, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 0.74371586
|
|
2 2.00000 0.74371586
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 31.984E-28; max= 76.436E-28
|
|
reduced coordinates (array xred) for 2 atoms
|
|
0.000000000000 0.000000000000 -0.180000000000
|
|
0.000000000000 0.000000000000 0.180000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 -1.52403036073920
|
|
2 0.00000000000000 0.00000000000000 1.52403036073920
|
|
length scales= 11.000000000000 11.000000000000 16.000000000000 bohr
|
|
= 5.820949294490 5.820949294490 8.466835337440 angstroms
|
|
prteigrs : about to open file t36o_DS2_EIG
|
|
Eigenvalues (hartree) for nkpt= 2 k points:
|
|
kpt# 1, nband= 5, wtk= 1.00000, kpt= 0.0000 0.0000 0.5000 (reduced coord)
|
|
-0.25822 -0.19808 -0.07047 -0.05633 -0.05633
|
|
kpt# 2, nband= 5, wtk= 1.00000, kpt= 0.0000 0.0000 1.0000 (reduced coord)
|
|
-0.25992 -0.19129 -0.09239 -0.05814 -0.05814
|
|
|
|
================================================================================
|
|
== DATASET 3 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 3, }
|
|
dimensions: {natom: 2, nkpt: 2, mband: 5, nsppol: 1, nspinor: 1, nspden: 1, mpw: 364, }
|
|
cutoff_energies: {ecut: 2.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 6.00000000E+00, charge: 0.00000000E+00, occopt: 4.00000000E+00, tsmear: 8.00000000E-02, }
|
|
meta: {optdriver: 1, rfphon: 1, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 1.
|
|
|
|
mkfilename : getwfq/=0, take file _WFQ from output of DATASET 2.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
LDA: new Teter (4/93) with spin-polarized option - ixc=1
|
|
Citation for XC functional:
|
|
S. Goedecker, M. Teter, J. Huetter, PRB 54, 1703 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 11.0000000 0.0000000 0.0000000 G(1)= 0.0909091 0.0000000 0.0000000
|
|
R(2)= 0.0000000 11.0000000 0.0000000 G(2)= 0.0000000 0.0909091 0.0000000
|
|
R(3)= 0.0000000 0.0000000 16.0000000 G(3)= 0.0000000 0.0000000 0.0625000
|
|
Unit cell volume ucvol= 1.9360000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
setup1 : take into account q-point for computing boxcut.
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.5000 ngfft= 16 16 24
|
|
ecut(hartree)= 2.500 => boxcut(ratio)= 2.01963
|
|
--------------------------------------------------------------------------------
|
|
|
|
|
|
==> initialize data related to q vector <==
|
|
|
|
The list of irreducible perturbations for this q vector is:
|
|
1) idir= 1 ipert= 2
|
|
|
|
================================================================================
|
|
|
|
--------------------------------------------------------------------------------
|
|
Perturbation wavevector (in red.coord.) 0.000000 0.000000 0.500000
|
|
Perturbation : displacement of atom 2 along direction 1
|
|
The set of symmetries contains only one element for this perturbation.
|
|
symkpt : not enough symmetry to change the number of k points.
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
dfpt_looppert : total number of electrons, from k and k+q
|
|
fully or partially occupied states are 6.000000E+00 and 6.000000E+00.
|
|
Initialisation of the first-order wave-functions :
|
|
ireadwf= 0
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 3, }
|
|
solver: {iscf: 7, nstep: 10, nline: 4, wfoptalg: 0, }
|
|
tolerances: {toldfe: 1.00E-10, }
|
|
...
|
|
|
|
iter 2DEtotal(Ha) deltaE(Ha) residm vres2
|
|
-ETOT 1 17.999003819130 -1.637E+01 3.929E-02 1.626E+03
|
|
ETOT 2 2.0792146211899 -1.592E+01 2.808E-03 1.442E+02
|
|
ETOT 3 0.43600146452443 -1.643E+00 9.871E-04 1.122E+00
|
|
ETOT 4 0.42868326774327 -7.318E-03 1.164E-05 1.853E-01
|
|
ETOT 5 0.42820355312869 -4.797E-04 3.848E-07 1.172E-02
|
|
ETOT 6 0.42808306569045 -1.205E-04 1.412E-07 3.487E-04
|
|
ETOT 7 0.42808038850926 -2.677E-06 1.093E-09 9.644E-06
|
|
ETOT 8 0.42808037367573 -1.483E-08 7.159E-11 2.087E-06
|
|
ETOT 9 0.42808035172459 -2.195E-08 1.267E-11 6.708E-07
|
|
ETOT 10 0.42808035102285 -7.017E-10 8.497E-15 8.061E-09
|
|
|
|
scprqt: WARNING -
|
|
nstep= 10 was not enough SCF cycles to converge;
|
|
maximum energy difference= 7.017E-10 exceeds toldfe= 1.000E-10
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 24.260E-16; max= 84.973E-16
|
|
|
|
Thirteen components of 2nd-order total energy (hartree) are
|
|
1,2,3: 0th-order hamiltonian combined with 1st-order wavefunctions
|
|
kin0= 3.17425321E+01 eigvalue= 6.92550513E+00 local= -2.27242225E+01
|
|
4,5,6: 1st-order hamiltonian combined with 1st and 0th-order wfs
|
|
loc psp = -7.14668586E+01 Hartree= 1.59024238E+01 xc= -7.40150965E+00
|
|
note that "loc psp" includes a xc core correction that could be resolved
|
|
7,8,9: eventually, occupation + non-local contributions
|
|
edocc= 3.15887195E+00 enl0= 6.33723228E+00 enl1= 3.58522113E+00
|
|
1-9 gives the relaxation energy (to be shifted if some occ is /=2.0)
|
|
erelax= -3.39408044E+01
|
|
10,11,12 Non-relaxation contributions : frozen-wavefunctions and Ewald
|
|
fr.local= 3.93342827E+01 fr.nonlo= -1.04009160E+00 Ewald= -3.92530630E+00
|
|
13,14 Frozen wf xc core corrections (1) and (2)
|
|
frxc 1 = 0.00000000E+00 frxc 2 = 0.00000000E+00
|
|
Resulting in :
|
|
2DEtotal= 0.4280803510E+00 Ha. Also 2DEtotal= 0.116486587546E+02 eV
|
|
(2DErelax= -3.3940804401E+01 Ha. 2DEnonrelax= 3.4368884752E+01 Ha)
|
|
( non-var. 2DEtotal : 4.2806602060E-01 Ha)
|
|
|
|
================================================================================
|
|
|
|
---- first-order wavefunction calculations are completed ----
|
|
|
|
|
|
==> Compute Derivative Database <==
|
|
|
|
2nd-order matrix (non-cartesian coordinates, masses not included,
|
|
asr not included )
|
|
j1 j2 matrix element
|
|
dir pert dir pert real part imaginary part
|
|
|
|
1 1 1 2 -0.4126556577 -0.0000000000
|
|
|
|
2 1 1 2 -0.0000000727 -0.0000000000
|
|
|
|
3 1 1 2 -0.0000032454 -0.0000000000
|
|
|
|
1 2 1 1 -0.4126556577 0.0000000000
|
|
1 2 2 1 -0.0000000727 0.0000000000
|
|
1 2 3 1 -0.0000032454 0.0000000000
|
|
1 2 1 2 0.4280660206 0.0000000000
|
|
1 2 2 2 -0.0000001437 0.0000000000
|
|
1 2 3 2 0.0000032970 -0.0000000000
|
|
|
|
2 2 1 2 -0.0000001437 -0.0000000000
|
|
|
|
3 2 1 2 0.0000032970 0.0000000000
|
|
|
|
|
|
Dynamical matrix, in cartesian coordinates,
|
|
if specified in the inputs, asr has been imposed
|
|
j1 j2 matrix element
|
|
dir pert dir pert real part imaginary part
|
|
|
|
1 1 1 2 -0.0034103773 -0.0000000000
|
|
|
|
2 1 1 2 -0.0000000006 -0.0000000000
|
|
|
|
3 1 1 2 -0.0000000184 -0.0000000000
|
|
|
|
1 2 1 1 -0.0034103773 0.0000000000
|
|
1 2 2 1 -0.0000000006 0.0000000000
|
|
1 2 3 1 -0.0000000184 0.0000000000
|
|
1 2 1 2 0.0035377357 0.0000000000
|
|
1 2 2 2 -0.0000000012 0.0000000000
|
|
1 2 3 2 0.0000000187 -0.0000000000
|
|
|
|
2 2 1 2 -0.0000000012 -0.0000000000
|
|
|
|
3 2 1 2 0.0000000187 0.0000000000
|
|
|
|
Phonon wavevector (reduced coordinates) : 0.00000 0.00000 0.50000
|
|
Phonon energies in Hartree :
|
|
0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
|
|
1.393101E-03
|
|
Phonon frequencies in cm-1 :
|
|
- 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
|
|
- 3.057504E+02
|
|
chkph3 : WARNING -
|
|
Dynamical matrix incomplete, phonon frequencies may be wrong, see the log file for more explanations.
|
|
|
|
================================================================================
|
|
== DATASET 4 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 4, }
|
|
dimensions: {natom: 4, nkpt: 1, mband: 10, nsppol: 1, nspinor: 1, nspden: 1, mpw: 362, }
|
|
cutoff_energies: {ecut: 2.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 1.20000000E+01, charge: 0.00000000E+00, occopt: 4.00000000E+00, tsmear: 8.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
LDA: new Teter (4/93) with spin-polarized option - ixc=1
|
|
Citation for XC functional:
|
|
S. Goedecker, M. Teter, J. Huetter, PRB 54, 1703 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 11.0000000 0.0000000 0.0000000 G(1)= 0.0909091 0.0000000 0.0000000
|
|
R(2)= 0.0000000 11.0000000 0.0000000 G(2)= 0.0000000 0.0909091 0.0000000
|
|
R(3)= 0.0000000 0.0000000 32.0000000 G(3)= 0.0000000 0.0000000 0.0312500
|
|
Unit cell volume ucvol= 3.8720000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 48
|
|
ecut(hartree)= 2.500 => boxcut(ratio)= 2.04358
|
|
6.54267545E+01 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 723.000 723.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 4, }
|
|
solver: {iscf: 7, nstep: 40, nline: 4, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -7.9118316886265 -7.912E+00 6.282E-04 2.047E+00
|
|
ETOT 2 -7.9141376891438 -2.306E-03 6.990E-08 4.137E+00
|
|
ETOT 3 -7.9161622100846 -2.025E-03 5.160E-06 6.975E-01
|
|
ETOT 4 -7.9152821478752 8.801E-04 5.410E-06 2.696E+00
|
|
ETOT 5 -7.9163692506555 -1.087E-03 3.686E-06 5.715E-01
|
|
ETOT 6 -7.9167099484209 -3.407E-04 7.920E-07 1.018E-02
|
|
ETOT 7 -7.9167226778457 -1.273E-05 9.348E-08 2.890E-03
|
|
ETOT 8 -7.9167248279999 -2.150E-06 8.132E-09 1.324E-03
|
|
ETOT 9 -7.9167258006760 -9.727E-07 1.623E-09 3.728E-05
|
|
ETOT 10 -7.9167258327632 -3.209E-08 1.216E-10 3.261E-06
|
|
ETOT 11 -7.9167258364571 -3.694E-09 1.923E-11 8.065E-08
|
|
ETOT 12 -7.9167258365684 -1.113E-10 1.453E-12 6.347E-08
|
|
ETOT 13 -7.9167258366092 -4.080E-11 1.089E-13 1.326E-09
|
|
ETOT 14 -7.9167258366113 -2.048E-12 8.592E-15 5.624E-10
|
|
ETOT 15 -7.9167258366114 -1.537E-13 1.758E-15 2.399E-10
|
|
ETOT 16 -7.9167258366116 -1.341E-13 4.092E-16 2.295E-11
|
|
ETOT 17 -7.9167258366116 -4.441E-14 9.580E-17 3.996E-12
|
|
ETOT 18 -7.9167258366116 2.665E-14 1.434E-17 4.480E-14
|
|
ETOT 19 -7.9167258366116 1.599E-14 4.213E-19 3.506E-15
|
|
ETOT 20 -7.9167258366116 5.329E-15 1.852E-20 6.842E-16
|
|
ETOT 21 -7.9167258366116 -2.220E-14 4.522E-21 4.498E-16
|
|
ETOT 22 -7.9167258366116 -3.464E-14 6.076E-22 4.724E-17
|
|
ETOT 23 -7.9167258366115 1.128E-13 1.262E-22 1.901E-18
|
|
ETOT 24 -7.9167258366116 -1.288E-13 7.159E-24 1.222E-18
|
|
ETOT 25 -7.9167258366114 1.990E-13 7.219E-24 6.708E-19
|
|
ETOT 26 -7.9167258366115 -6.040E-14 1.523E-24 2.981E-20
|
|
ETOT 27 -7.9167258366116 -1.066E-13 1.181E-25 1.301E-21
|
|
|
|
At SCF step 27 vres2 = 1.30E-21 < tolvrs= 1.00E-20 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543351E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543351E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080290E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 4, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 11.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 11.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 32.0000000, ]
|
|
lattice_lengths: [ 11.00000, 11.00000, 32.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 3.8720000E+03
|
|
convergence: {deltae: -1.066E-13, res2: 1.301E-21, residm: 1.181E-25, diffor: null, }
|
|
etotal : -7.91672584E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -7.70666970E-02
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 5.63543351E-05, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 5.63543351E-05, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 7.57080290E-05, ]
|
|
pressure_GPa: -1.8478E+00
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, -9.0000E-02, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 9.0000E-02, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 4.1000E-01, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 5.9000E-01, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -7.43798460E-16, 5.60689876E-16, 1.04287834E-02, ]
|
|
- [ -3.32109299E-15, -1.80032212E-15, -1.04287834E-02, ]
|
|
- [ 4.33136402E-15, -1.16194816E-15, 1.04287834E-02, ]
|
|
- [ -2.66472572E-16, 2.40158041E-15, -1.04287834E-02, ]
|
|
force_length_stats: {min: 1.04287834E-02, max: 1.04287834E-02, mean: 1.04287834E-02, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 0.74371586
|
|
2 2.00000 0.74371586
|
|
3 2.00000 0.74371586
|
|
4 2.00000 0.74371586
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 61.713E-27; max= 11.807E-26
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 -0.090000000000
|
|
0.000000000000 0.000000000000 0.090000000000
|
|
0.000000000000 0.000000000000 0.410000000000
|
|
0.000000000000 0.000000000000 0.590000000000
|
|
rms dE/dt= 1.9267E-01; max dE/dt= 3.3372E-01; dE/dt below (all hartree)
|
|
1 0.000000000000 -0.000000000000 -0.333721067432
|
|
2 0.000000000000 0.000000000000 0.333721067432
|
|
3 -0.000000000000 0.000000000000 -0.333721067432
|
|
4 0.000000000000 -0.000000000000 0.333721067432
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 -1.52403036073920
|
|
2 0.00000000000000 0.00000000000000 1.52403036073920
|
|
3 0.00000000000000 0.00000000000000 6.94280497670080
|
|
4 0.00000000000000 0.00000000000000 9.99086569817920
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.00000000000000 0.00000000000000 0.01042878335725
|
|
2 -0.00000000000000 -0.00000000000000 -0.01042878335725
|
|
3 0.00000000000000 -0.00000000000000 0.01042878335725
|
|
4 -0.00000000000000 0.00000000000000 -0.01042878335725
|
|
frms,max,avg= 6.0210609E-03 1.0428783E-02 0.000E+00 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -0.00000000000004 0.00000000000003 0.53626955681443
|
|
2 -0.00000000000017 -0.00000000000009 -0.53626955681469
|
|
3 0.00000000000022 -0.00000000000006 0.53626955681458
|
|
4 -0.00000000000001 0.00000000000012 -0.53626955681432
|
|
frms,max,avg= 3.0961537E-01 5.3626956E-01 0.000E+00 0.000E+00 0.000E+00 e/A
|
|
length scales= 11.000000000000 11.000000000000 32.000000000000 bohr
|
|
= 5.820949294490 5.820949294490 16.933670674880 angstroms
|
|
prteigrs : about to open file t36o_DS4_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.07707 Average Vxc (hartree)= -0.13648
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 10, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.25992 -0.25822 -0.19808 -0.19129 -0.09239 -0.07047 -0.05814 -0.05814
|
|
-0.05633 -0.05633
|
|
occupation numbers for kpt# 1
|
|
2.01380 2.01494 2.11182 2.12510 1.17566 0.70527 0.47820 0.47820
|
|
0.44851 0.44851
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 4, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 2.25686746055175E+00
|
|
hartree : 1.92591553447194E+00
|
|
xc : -2.41198564047337E+00
|
|
Ewald energy : -4.52462915011266E+00
|
|
psp_core : 1.68974056134924E-02
|
|
local_psp : -7.24481088170569E+00
|
|
non_local_psp : 2.14145194437608E+00
|
|
internal : -7.84029332727847E+00
|
|
'-kT*entropy' : -7.64325093331471E-02
|
|
total_energy : -7.91672583661162E+00
|
|
total_energy_eV : -2.15425065654418E+02
|
|
band_energy : -2.13298802335214E+00
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543351E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543351E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080290E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -1.8478E+00 GPa]
|
|
- sigma(1 1)= 1.65800150E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 1.65800150E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 2.22740674E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
================================================================================
|
|
== DATASET 5 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 5, }
|
|
dimensions: {natom: 4, nkpt: 1, mband: 10, nsppol: 1, nspinor: 1, nspden: 1, mpw: 362, }
|
|
cutoff_energies: {ecut: 2.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 1.20000000E+01, charge: 0.00000000E+00, occopt: 4.00000000E+00, tsmear: 8.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 4.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
LDA: new Teter (4/93) with spin-polarized option - ixc=1
|
|
Citation for XC functional:
|
|
S. Goedecker, M. Teter, J. Huetter, PRB 54, 1703 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 11.0000000 0.0000000 0.0000000 G(1)= 0.0909091 0.0000000 0.0000000
|
|
R(2)= 0.0000000 11.0000000 0.0000000 G(2)= 0.0000000 0.0909091 0.0000000
|
|
R(3)= 0.0000000 0.0000000 32.0000000 G(3)= 0.0000000 0.0000000 0.0312500
|
|
Unit cell volume ucvol= 3.8720000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 16 16 48
|
|
ecut(hartree)= 2.500 => boxcut(ratio)= 2.04358
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t36o_DS4_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 723.000 723.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 5, }
|
|
solver: {iscf: 7, nstep: 40, nline: 4, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-20, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -7.9167256534328 -7.917E+00 9.494E-13 3.278E-05
|
|
ETOT 2 -7.9167258125588 -1.591E-07 5.155E-14 3.496E-06
|
|
ETOT 3 -7.9167258322346 -1.968E-08 1.937E-10 2.467E-08
|
|
ETOT 4 -7.9167258323255 -9.089E-11 6.792E-13 3.371E-09
|
|
ETOT 5 -7.9167258323307 -5.244E-12 2.420E-14 4.753E-11
|
|
ETOT 6 -7.9167258323308 -1.084E-13 1.261E-15 5.227E-12
|
|
ETOT 7 -7.9167258323308 3.020E-14 6.380E-17 6.219E-12
|
|
ETOT 8 -7.9167258323308 4.885E-14 1.167E-17 1.234E-12
|
|
ETOT 9 -7.9167258323309 -1.039E-13 2.129E-18 2.254E-14
|
|
ETOT 10 -7.9167258323308 3.819E-14 1.454E-18 8.877E-15
|
|
ETOT 11 -7.9167258323309 -3.819E-14 5.717E-20 1.417E-15
|
|
ETOT 12 -7.9167258323308 3.819E-14 8.923E-21 1.768E-16
|
|
ETOT 13 -7.9167258323308 7.994E-15 1.810E-21 1.761E-17
|
|
ETOT 14 -7.9167258323308 -2.132E-14 1.009E-22 4.724E-18
|
|
ETOT 15 -7.9167258323308 3.908E-14 1.255E-23 6.086E-19
|
|
ETOT 16 -7.9167258323308 3.375E-14 5.688E-24 1.467E-19
|
|
ETOT 17 -7.9167258323307 1.599E-14 3.182E-25 1.072E-19
|
|
ETOT 18 -7.9167258323308 -3.553E-14 1.659E-25 1.698E-20
|
|
ETOT 19 -7.9167258323308 -1.776E-14 8.178E-26 1.235E-21
|
|
|
|
At SCF step 19 vres2 = 1.24E-21 < tolvrs= 1.00E-20 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543362E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543353E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080295E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 5, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 11.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 11.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 32.0000000, ]
|
|
lattice_lengths: [ 11.00000, 11.00000, 32.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 3.8720000E+03
|
|
convergence: {deltae: -1.776E-14, res2: 1.235E-21, residm: 8.178E-26, diffor: null, }
|
|
etotal : -7.91672583E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -7.70666957E-02
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ 5.63543362E-05, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 5.63543353E-05, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 7.57080295E-05, ]
|
|
pressure_GPa: -1.8478E+00
|
|
xred :
|
|
- [ 0.0000E+00, 0.0000E+00, -9.0000E-02, Al]
|
|
- [ 1.0000E-04, 0.0000E+00, 9.0000E-02, Al]
|
|
- [ 0.0000E+00, 0.0000E+00, 4.1000E-01, Al]
|
|
- [ -1.0000E-04, 0.0000E+00, 5.9000E-01, Al]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ 3.75138892E-06, -1.95841050E-17, 1.04287836E-02, ]
|
|
- [ -3.89163923E-06, 3.46122071E-17, -1.04287836E-02, ]
|
|
- [ -3.75138892E-06, 7.38408944E-17, 1.04287836E-02, ]
|
|
- [ 3.89163923E-06, -8.88689965E-17, -1.04287836E-02, ]
|
|
force_length_stats: {min: 1.04287843E-02, max: 1.04287844E-02, mean: 1.04287843E-02, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 0.74371586
|
|
2 2.00000 0.74371575
|
|
3 2.00000 0.74371586
|
|
4 2.00000 0.74371575
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 37.440E-27; max= 81.778E-27
|
|
reduced coordinates (array xred) for 4 atoms
|
|
0.000000000000 0.000000000000 -0.090000000000
|
|
0.000100000000 0.000000000000 0.090000000000
|
|
0.000000000000 0.000000000000 0.410000000000
|
|
-0.000100000000 0.000000000000 0.590000000000
|
|
rms dE/dt= 1.9267E-01; max dE/dt= 3.3372E-01; dE/dt below (all hartree)
|
|
1 -0.000041265278 -0.000000000000 -0.333721076206
|
|
2 0.000042808031 -0.000000000000 0.333721076205
|
|
3 0.000041265278 -0.000000000000 -0.333721076206
|
|
4 -0.000042808031 -0.000000000000 0.333721076205
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 0.00000000000000 0.00000000000000 -1.52403036073920
|
|
2 0.00058209492945 0.00000000000000 1.52403036073920
|
|
3 0.00000000000000 0.00000000000000 6.94280497670080
|
|
4 -0.00058209492945 0.00000000000000 9.99086569817920
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 0.00000375138892 -0.00000000000000 0.01042878363142
|
|
2 -0.00000389163923 0.00000000000000 -0.01042878363143
|
|
3 -0.00000375138892 0.00000000000000 0.01042878363143
|
|
4 0.00000389163923 -0.00000000000000 -0.01042878363142
|
|
frms,max,avg= 6.0210614E-03 1.0428784E-02 0.000E+00 0.000E+00 1.126E-14 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 0.00019290415796 -0.00000000000000 0.53626957091268
|
|
2 -0.00020011611819 0.00000000000000 -0.53626957091324
|
|
3 -0.00019290415796 0.00000000000000 0.53626957091357
|
|
4 0.00020011611820 -0.00000000000000 -0.53626957091300
|
|
frms,max,avg= 3.0961540E-01 5.3626957E-01 0.000E+00 0.000E+00 5.790E-13 e/A
|
|
length scales= 11.000000000000 11.000000000000 32.000000000000 bohr
|
|
= 5.820949294490 5.820949294490 16.933670674880 angstroms
|
|
prteigrs : about to open file t36o_DS5_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.07707 Average Vxc (hartree)= -0.13648
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 10, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.25992 -0.25822 -0.19808 -0.19129 -0.09239 -0.07047 -0.05814 -0.05814
|
|
-0.05633 -0.05633
|
|
occupation numbers for kpt# 1
|
|
2.01380 2.01494 2.11182 2.12510 1.17566 0.70527 0.47820 0.47820
|
|
0.44851 0.44851
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 5, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 2.25686745103638E+00
|
|
hartree : 1.92591549884170E+00
|
|
xc : -2.41198563392494E+00
|
|
Ewald energy : -4.52462918936574E+00
|
|
psp_core : 1.68974056134924E-02
|
|
local_psp : -7.24481079992773E+00
|
|
non_local_psp : 2.14145194530581E+00
|
|
internal : -7.84029332242103E+00
|
|
'-kT*entropy' : -7.64325099097651E-02
|
|
total_energy : -7.91672583233080E+00
|
|
total_energy_eV : -2.15425065537931E+02
|
|
band_energy : -2.13298801285218E+00
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= 5.63543362E-05 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 5.63543353E-05 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 7.57080295E-05 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -1.8478E+00 GPa]
|
|
- sigma(1 1)= 1.65800153E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 1.65800151E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 2.22740675E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
== END DATASET(S) ==============================================================
|
|
================================================================================
|
|
|
|
-outvars: echo values of variables after computation --------
|
|
acell1 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell2 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell3 1.1000000000E+01 1.1000000000E+01 1.6000000000E+01 Bohr
|
|
acell4 1.1000000000E+01 1.1000000000E+01 3.2000000000E+01 Bohr
|
|
acell5 1.1000000000E+01 1.1000000000E+01 3.2000000000E+01 Bohr
|
|
amu 1.00000000E+00
|
|
asr 0
|
|
chneut 0
|
|
diemac 1.00000000E+00
|
|
diemix1 3.33333333E-01
|
|
diemix2 3.33333333E-01
|
|
diemix3 3.50000000E-01
|
|
diemix4 3.33333333E-01
|
|
diemix5 3.33333333E-01
|
|
ecut 2.50000000E+00 Hartree
|
|
etotal1 -3.9583629183E+00
|
|
etotal3 4.2808035102E-01
|
|
etotal4 -7.9167258366E+00
|
|
etotal5 -7.9167258323E+00
|
|
fcart1 2.4018498246E-14 4.0230144392E-15 1.0428783357E-02
|
|
-2.4018498246E-14 -4.0230144392E-15 -1.0428783357E-02
|
|
fcart3 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
fcart4 -7.4379846036E-16 5.6068987613E-16 1.0428783357E-02
|
|
-3.3210929867E-15 -1.8003221238E-15 -1.0428783357E-02
|
|
4.3313640192E-15 -1.1619481643E-15 1.0428783357E-02
|
|
-2.6647257211E-16 2.4015804120E-15 -1.0428783357E-02
|
|
fcart5 3.7513889172E-06 -1.9584105004E-17 1.0428783631E-02
|
|
-3.8916392258E-06 3.4612207114E-17 -1.0428783631E-02
|
|
-3.7513889173E-06 7.3840894397E-17 1.0428783631E-02
|
|
3.8916392259E-06 -8.8868996508E-17 -1.0428783631E-02
|
|
- fftalg 512
|
|
getden1 0
|
|
getden2 1
|
|
getden3 0
|
|
getden4 0
|
|
getden5 0
|
|
getwfk1 0
|
|
getwfk2 1
|
|
getwfk3 1
|
|
getwfk4 0
|
|
getwfk5 4
|
|
getwfq1 0
|
|
getwfq2 0
|
|
getwfq3 2
|
|
getwfq4 0
|
|
getwfq5 0
|
|
iscf1 7
|
|
iscf2 -2
|
|
iscf3 7
|
|
iscf4 7
|
|
iscf5 7
|
|
istwfk1 2 4
|
|
istwfk2 4 0
|
|
istwfk3 1 1
|
|
istwfk4 2
|
|
istwfk5 2
|
|
jdtset 1 2 3 4 5
|
|
kpt1 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt3 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
kpt4 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
kpt5 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
kptopt 0
|
|
P mkmem1 2
|
|
P mkmem2 2
|
|
P mkmem3 2
|
|
P mkmem4 1
|
|
P mkmem5 1
|
|
P mkqmem1 2
|
|
P mkqmem2 2
|
|
P mkqmem3 2
|
|
P mkqmem4 1
|
|
P mkqmem5 1
|
|
P mk1mem1 2
|
|
P mk1mem2 2
|
|
P mk1mem3 2
|
|
P mk1mem4 1
|
|
P mk1mem5 1
|
|
natom1 2
|
|
natom2 2
|
|
natom3 2
|
|
natom4 4
|
|
natom5 4
|
|
nband1 5
|
|
nband2 5
|
|
nband3 5
|
|
nband4 10
|
|
nband5 10
|
|
ndtset 5
|
|
ngfft1 16 16 24
|
|
ngfft2 16 16 24
|
|
ngfft3 16 16 24
|
|
ngfft4 16 16 48
|
|
ngfft5 16 16 48
|
|
nkpt1 2
|
|
nkpt2 2
|
|
nkpt3 2
|
|
nkpt4 1
|
|
nkpt5 1
|
|
nqpt1 0
|
|
nqpt2 1
|
|
nqpt3 1
|
|
nqpt4 0
|
|
nqpt5 0
|
|
nstep1 40
|
|
nstep2 40
|
|
nstep3 10
|
|
nstep4 40
|
|
nstep5 40
|
|
nsym 1
|
|
ntypat 1
|
|
occ1 2.013803 2.125104 1.175657 0.478195 0.478195
|
|
2.014940 2.111825 0.705269 0.448506 0.448506
|
|
occ3 2.013803 2.125104 1.175657 0.478195 0.478195
|
|
2.014940 2.111825 0.705269 0.448506 0.448506
|
|
occ4 2.013803 2.014940 2.111825 2.125104 1.175657 0.705269
|
|
0.478195 0.478195 0.448506 0.448506
|
|
occ5 2.013803 2.014940 2.111825 2.125104 1.175657 0.705269
|
|
0.478195 0.478195 0.448506 0.448506
|
|
occopt 4
|
|
optdriver1 0
|
|
optdriver2 0
|
|
optdriver3 1
|
|
optdriver4 0
|
|
optdriver5 0
|
|
prtpot1 0
|
|
prtpot2 0
|
|
prtpot3 1
|
|
prtpot4 0
|
|
prtpot5 0
|
|
qpt1 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
qpt2 0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
qpt3 0.00000000E+00 0.00000000E+00 5.00000000E-01
|
|
qpt4 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
qpt5 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
rfatpol1 1 2
|
|
rfatpol2 1 2
|
|
rfatpol3 2 2
|
|
rfatpol4 1 4
|
|
rfatpol5 1 4
|
|
rfdir1 1 1 1
|
|
rfdir2 1 1 1
|
|
rfdir3 1 0 0
|
|
rfdir4 1 1 1
|
|
rfdir5 1 1 1
|
|
rfphon1 0
|
|
rfphon2 0
|
|
rfphon3 1
|
|
rfphon4 0
|
|
rfphon5 0
|
|
spgroup 1
|
|
strten1 5.6354335095E-05 5.6354335091E-05 7.5708028999E-05
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten3 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten4 5.6354335095E-05 5.6354335091E-05 7.5708028999E-05
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten5 5.6354336207E-05 5.6354335267E-05 7.5708029480E-05
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
toldfe1 0.00000000E+00 Hartree
|
|
toldfe2 0.00000000E+00 Hartree
|
|
toldfe3 1.00000000E-10 Hartree
|
|
toldfe4 0.00000000E+00 Hartree
|
|
toldfe5 0.00000000E+00 Hartree
|
|
tolvrs1 1.00000000E-20
|
|
tolvrs2 0.00000000E+00
|
|
tolvrs3 0.00000000E+00
|
|
tolvrs4 1.00000000E-20
|
|
tolvrs5 1.00000000E-20
|
|
tolwfr1 0.00000000E+00
|
|
tolwfr2 1.00000000E-26
|
|
tolwfr3 0.00000000E+00
|
|
tolwfr4 0.00000000E+00
|
|
tolwfr5 0.00000000E+00
|
|
tsmear 8.00000000E-02 Hartree
|
|
typat1 1 1
|
|
typat2 1 1
|
|
typat3 1 1
|
|
typat4 1 1 1 1
|
|
typat5 1 1 1 1
|
|
wtk1 0.50000 0.50000
|
|
wtk2 1.00000 1.00000
|
|
wtk3 0.50000 0.50000
|
|
wtk4 1.00000
|
|
wtk5 1.00000
|
|
xangst1 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst2 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst3 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
xangst4 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 6.9428049767E+00
|
|
0.0000000000E+00 0.0000000000E+00 9.9908656982E+00
|
|
xangst5 0.0000000000E+00 0.0000000000E+00 -1.5240303607E+00
|
|
5.8209492945E-04 0.0000000000E+00 1.5240303607E+00
|
|
0.0000000000E+00 0.0000000000E+00 6.9428049767E+00
|
|
-5.8209492945E-04 0.0000000000E+00 9.9908656982E+00
|
|
xcart1 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart2 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart3 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
xcart4 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.3120000000E+01
|
|
0.0000000000E+00 0.0000000000E+00 1.8880000000E+01
|
|
xcart5 0.0000000000E+00 0.0000000000E+00 -2.8800000000E+00
|
|
1.1000000000E-03 0.0000000000E+00 2.8800000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 1.3120000000E+01
|
|
-1.1000000000E-03 0.0000000000E+00 1.8880000000E+01
|
|
xred1 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred2 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred3 0.0000000000E+00 0.0000000000E+00 -1.8000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 1.8000000000E-01
|
|
xred4 0.0000000000E+00 0.0000000000E+00 -9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 4.1000000000E-01
|
|
0.0000000000E+00 0.0000000000E+00 5.9000000000E-01
|
|
xred5 0.0000000000E+00 0.0000000000E+00 -9.0000000000E-02
|
|
1.0000000000E-04 0.0000000000E+00 9.0000000000E-02
|
|
0.0000000000E+00 0.0000000000E+00 4.1000000000E-01
|
|
-1.0000000000E-04 0.0000000000E+00 5.9000000000E-01
|
|
znucl 13.00000
|
|
|
|
================================================================================
|
|
|
|
The spacegroup number, the magnetic point group, and/or the number of symmetries
|
|
have changed between the initial recognition based on the input file
|
|
and a postprocessing based on the final acell, rprim, and xred.
|
|
More details in the log file.
|
|
|
|
|
|
- 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] 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
|
|
-
|
|
- [2] First-principles responses of solids to atomic displacements and homogeneous electric fields:,
|
|
- implementation of a conjugate-gradient algorithm. X. Gonze, Phys. Rev. B55, 10337 (1997).
|
|
- Comment: Non-vanishing rfphon and/or rfelfd, in the norm-conserving case.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze1997
|
|
-
|
|
- [3] Dynamical matrices, Born effective charges, dielectric permittivity tensors, and ,
|
|
- interatomic force constants from density-functional perturbation theory,
|
|
- X. Gonze and C. Lee, Phys. Rev. B55, 10355 (1997).
|
|
- Comment: Non-vanishing rfphon and/or rfelfd, in the norm-conserving case.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze1997a
|
|
-
|
|
- [4] Ab initio pseudopotentials for electronic structure calculations of poly-atomic systems,
|
|
- using density-functional theory.
|
|
- M. Fuchs and, M. Scheffler, Comput. Phys. Commun. 119, 67 (1999).
|
|
- Comment: Some pseudopotential generated using the FHI code were used.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#fuchs1999
|
|
-
|
|
- [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= 2.3 wall= 4.5
|
|
|
|
================================================================================
|
|
|
|
Calculation completed.
|
|
.Delivered 98 WARNINGs and 23 COMMENTs to log file.
|
|
+Overall time at end (sec) : cpu= 2.3 wall= 4.5
|