mirror of https://github.com/abinit/abinit.git
1510 lines
72 KiB
Plaintext
1510 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/v3_t07/t07.abi
|
|
- output file -> t07.abo
|
|
- root for input files -> t07i
|
|
- root for output files -> t07o
|
|
|
|
DATASET 1 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 1.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 1
|
|
lnmax = 1 mgfft = 24 mpssoang = 1 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 4 n1xccc = 0 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
mpw = 280 nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 3.330 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.006 Mbytes ; DEN or POT disk file : 0.075 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 2 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 2.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 1
|
|
lnmax = 1 mgfft = 24 mpssoang = 1 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 4 n1xccc = 0 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
mpw = 280 nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 3.330 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.006 Mbytes ; DEN or POT disk file : 0.075 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 3 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 3.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 1
|
|
lnmax = 1 mgfft = 24 mpssoang = 1 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 4 n1xccc = 0 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
mpw = 280 nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 3.330 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.006 Mbytes ; DEN or POT disk file : 0.075 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 4 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 4.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 1
|
|
lnmax = 1 mgfft = 24 mpssoang = 1 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 4 n1xccc = 0 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
mpw = 280 nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 3.330 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.006 Mbytes ; DEN or POT disk file : 0.075 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 5 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 5.
|
|
intxc = 0 ionmov = 0 iscf = 7 lmnmax = 1
|
|
lnmax = 1 mgfft = 24 mpssoang = 1 mqgrid = 3001
|
|
natom = 2 nloc_mem = 1 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 4 n1xccc = 0 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
mpw = 280 nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 3.330 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.006 Mbytes ; DEN or POT disk file : 0.075 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 6 : space group Pm m 2 (# 25); Bravais oP (primitive ortho.)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 6 (RF).
|
|
intxc = 0 iscf = 7 lmnmax = 1 lnmax = 1
|
|
mgfft = 24 mpssoang = 1 mqgrid = 3001 natom = 2
|
|
nloc_mem = 1 nspden = 1 nspinor = 1 nsppol = 1
|
|
nsym = 4 n1xccc = 0 ntypat = 1 occopt = 1
|
|
xclevel = 2
|
|
- mband = 1 mffmem = 1 mkmem = 1
|
|
- mkqmem = 1 mk1mem = 1 mpw = 559
|
|
nfft = 9600 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 2.538 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.011 Mbytes ; DEN or POT disk file : 0.075 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 --------
|
|
acell 1.2000000000E+01 1.0000000000E+01 1.0000000000E+01 Bohr
|
|
amu 1.00800000E+00
|
|
asr 0
|
|
chneut 0
|
|
diemac 1.00000000E+00
|
|
diemix 5.00000000E-01
|
|
ecut 4.50000000E+00 Hartree
|
|
- fftalg 512
|
|
getwfk1 -1
|
|
getwfk2 -1
|
|
getwfk3 -1
|
|
getwfk4 -1
|
|
getwfk5 -1
|
|
getwfk6 3
|
|
istwfk1 2
|
|
istwfk2 2
|
|
istwfk3 2
|
|
istwfk4 2
|
|
istwfk5 2
|
|
istwfk6 1
|
|
ixc 12
|
|
jdtset 1 2 3 4 5 6
|
|
kptopt 0
|
|
P mkmem 1
|
|
P mkqmem 1
|
|
P mk1mem 1
|
|
natom 2
|
|
nband 1
|
|
ndtset 6
|
|
ngfft 24 20 20
|
|
nkpt 1
|
|
nline 3
|
|
nqpt1 0
|
|
nqpt2 0
|
|
nqpt3 0
|
|
nqpt4 0
|
|
nqpt5 0
|
|
nqpt6 1
|
|
nstep1 30
|
|
nstep2 30
|
|
nstep3 30
|
|
nstep4 30
|
|
nstep5 30
|
|
nstep6 20
|
|
nsym 4
|
|
ntypat 1
|
|
occ 2.000000
|
|
optdriver1 0
|
|
optdriver2 0
|
|
optdriver3 0
|
|
optdriver4 0
|
|
optdriver5 0
|
|
optdriver6 1
|
|
prtpot1 0
|
|
prtpot2 0
|
|
prtpot3 0
|
|
prtpot4 0
|
|
prtpot5 0
|
|
prtpot6 1
|
|
rfatpol1 1 2
|
|
rfatpol2 1 2
|
|
rfatpol3 1 2
|
|
rfatpol4 1 2
|
|
rfatpol5 1 2
|
|
rfatpol6 2 2
|
|
rfdir1 1 1 1
|
|
rfdir2 1 1 1
|
|
rfdir3 1 1 1
|
|
rfdir4 1 1 1
|
|
rfdir5 1 1 1
|
|
rfdir6 1 0 0
|
|
rfphon1 0
|
|
rfphon2 0
|
|
rfphon3 0
|
|
rfphon4 0
|
|
rfphon5 0
|
|
rfphon6 1
|
|
spgroup 25
|
|
symrel 1 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 0 -1
|
|
1 0 0 0 -1 0 0 0 1 1 0 0 0 -1 0 0 0 -1
|
|
tolvrs 1.00000000E-08
|
|
typat 1 1
|
|
xangst1 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9782093299E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst2 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9813843932E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst3 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst4 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9877345197E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst5 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9909095830E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst6 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart1 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6280000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart2 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6340000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart3 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart4 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6460000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart5 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6520000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart6 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xred1 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.6900000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred2 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.6950000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred3 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred4 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7050000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred5 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7100000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred6 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
znucl 1.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.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 6.
|
|
|
|
================================================================================
|
|
== DATASET 1 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 1, }
|
|
dimensions: {natom: 2, nkpt: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 280, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
|
|
--- Pseudopotential description ------------------------------------------------
|
|
- pspini: atom type 1 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/PseudosTM_pwteter/1h.pspnc
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/PseudosTM_pwteter/1h.pspnc
|
|
- Troullier-Martins psp for element H Thu Oct 27 17:28:54 EDT 1994
|
|
- 1.00000 1.00000 940714 znucl, zion, pspdat
|
|
1 1 0 0 2001 0.00000 pspcod,pspxc,lmax,lloc,mmax,r2well
|
|
0 7.740 11.990 0 1.5855604 l,e99.0,e99.9,nproj,rcpsp
|
|
0.00000000 0.00000000 0.00000000 0.00000000 rms, ekb1, ekb2, epsatm
|
|
0.00000000000000 0.00000000000000 0.00000000000000 rchrg,fchrg,qchrg
|
|
Note: local psp for atom with Z= 1.0
|
|
pspatm : epsatm= 0.04198703
|
|
--- l ekb(1:nproj) -->
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
1.67948119E-01 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 559.000 559.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 1, }
|
|
solver: {iscf: 7, nstep: 30, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -1.0373131705935 -1.037E+00 4.047E-04 5.791E+00
|
|
ETOT 2 -1.0395244613655 -2.211E-03 4.373E-10 2.019E+00
|
|
ETOT 3 -1.0398352725554 -3.108E-04 3.929E-07 5.057E-01
|
|
ETOT 4 -1.0398374436416 -2.171E-06 1.755E-09 1.116E-01
|
|
ETOT 5 -1.0398380228789 -5.792E-07 1.088E-09 4.217E-04
|
|
ETOT 6 -1.0398380251713 -2.292E-09 4.824E-12 3.535E-05
|
|
ETOT 7 -1.0398380252948 -1.235E-10 2.167E-13 1.696E-07
|
|
ETOT 8 -1.0398380253006 -5.883E-12 7.638E-15 1.814E-08
|
|
ETOT 9 -1.0398380253007 -4.885E-14 1.736E-16 7.281E-11
|
|
|
|
At SCF step 9 vres2 = 7.28E-11 < tolvrs= 1.00E-08 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.03108870E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15673558E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15673558E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 1, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 12.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 10.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 10.0000000, ]
|
|
lattice_lengths: [ 12.00000, 10.00000, 10.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.2000000E+03
|
|
convergence: {deltae: -4.885E-14, res2: 7.281E-11, residm: 1.736E-16, diffor: null, }
|
|
etotal : -1.03983803E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -3.75839724E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -2.03108870E-04, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 1.15673558E-04, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 1.15673558E-04, ]
|
|
pressure_GPa: -2.7693E-01
|
|
xred :
|
|
- [ -4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
- [ 4.6900E-02, 0.0000E+00, 0.0000E+00, H]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.58137395E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ 2.58137395E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
force_length_stats: {min: 2.58137395E-01, max: 2.58137395E-01, mean: 2.58137395E-01, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 1.58516000
|
|
2 2.00000 1.58314716
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 17.359E-17; max= 17.359E-17
|
|
reduced coordinates (array xred) for 2 atoms
|
|
-0.047000000000 0.000000000000 0.000000000000
|
|
0.046900000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 1.7884E+00; max dE/dt= 3.0976E+00; dE/dt below (all hartree)
|
|
1 3.097649212049 0.000000000000 0.000000000000
|
|
2 -3.097648264787 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 -0.29845594564476 0.00000000000000 0.00000000000000
|
|
2 0.29782093299445 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.25813739486820 -0.00000000000000 -0.00000000000000
|
|
2 0.25813739486820 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 1.4903569E-01 2.5813739E-01 -3.947E-08 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -13.27395743118894 -0.00000000000000 -0.00000000000000
|
|
2 13.27395743118894 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 7.6637229E+00 1.3273957E+01 -2.030E-06 0.000E+00 0.000E+00 e/A
|
|
length scales= 12.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 6.350126503080 5.291772085900 5.291772085900 angstroms
|
|
prteigrs : about to open file t07o_DS1_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.37584 Average Vxc (hartree)= -0.06014
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 1, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.37584
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 1, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 1.10059241157634E+00
|
|
hartree : 8.32640618571442E-01
|
|
xc : -6.52122311976919E-01
|
|
Ewald energy : 3.60683130499491E-01
|
|
psp_core : 1.39956766170961E-04
|
|
local_psp : -2.68177183073720E+00
|
|
non_local_psp : 0.00000000000000E+00
|
|
total_energy : -1.03983802530068E+00
|
|
total_energy_eV : -2.82954316586813E+01
|
|
band_energy : -7.51679448383716E-01
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.03108870E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15673558E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15673558E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -2.7693E-01 GPa]
|
|
- sigma(1 1)= -5.97566827E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 3.40323300E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 3.40323300E+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: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 280, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 1.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t07o_DS1_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 559.000 559.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: 7, nstep: 30, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -1.0399926560218 -1.040E+00 6.323E-13 5.171E-04
|
|
ETOT 2 -1.0399926691333 -1.311E-08 2.093E-14 6.776E-05
|
|
ETOT 3 -1.0399926705258 -1.392E-09 8.236E-12 1.092E-06
|
|
ETOT 4 -1.0399926705313 -5.513E-12 7.905E-15 3.284E-07
|
|
ETOT 5 -1.0399926705339 -2.674E-12 1.494E-14 8.280E-08
|
|
ETOT 6 -1.0399926705339 1.221E-14 2.574E-17 7.043E-09
|
|
|
|
At SCF step 6 vres2 = 7.04E-09 < tolvrs= 1.00E-08 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.02569732E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15593581E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15593581E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 2, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 12.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 10.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 10.0000000, ]
|
|
lattice_lengths: [ 12.00000, 10.00000, 10.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.2000000E+03
|
|
convergence: {deltae: 1.221E-14, res2: 7.043E-09, residm: 2.574E-17, diffor: null, }
|
|
etotal : -1.03999267E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -3.75765495E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -2.02569732E-04, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 1.15593581E-04, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 1.15593581E-04, ]
|
|
pressure_GPa: -2.8065E-01
|
|
xred :
|
|
- [ -4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
- [ 4.6950E-02, 0.0000E+00, 0.0000E+00, H]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.57347031E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ 2.57347031E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
force_length_stats: {min: 2.57347031E-01, max: 2.57347031E-01, mean: 2.57347031E-01, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 1.58489860
|
|
2 2.00000 1.58304309
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 25.737E-18; max= 25.737E-18
|
|
reduced coordinates (array xred) for 2 atoms
|
|
-0.047000000000 0.000000000000 0.000000000000
|
|
0.046950000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 1.7830E+00; max dE/dt= 3.0882E+00; dE/dt below (all hartree)
|
|
1 3.088164615073 0.000000000000 0.000000000000
|
|
2 -3.088164121443 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 -0.29845594564476 0.00000000000000 0.00000000000000
|
|
2 0.29813843931961 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.25734703068817 -0.00000000000000 -0.00000000000000
|
|
2 0.25734703068817 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 1.4857938E-01 2.5734703E-01 -2.057E-08 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -13.23331527437827 -0.00000000000000 -0.00000000000000
|
|
2 13.23331527437827 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 7.6402581E+00 1.3233315E+01 -1.058E-06 0.000E+00 0.000E+00 e/A
|
|
length scales= 12.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 6.350126503080 5.291772085900 5.291772085900 angstroms
|
|
prteigrs : about to open file t07o_DS2_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.37577 Average Vxc (hartree)= -0.06015
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 1, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.37577
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 2, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 1.10027684434829E+00
|
|
hartree : 8.32444685852306E-01
|
|
xc : -6.52028599392947E-01
|
|
Ewald energy : 3.60211990582588E-01
|
|
psp_core : 1.39956766170961E-04
|
|
local_psp : -2.68103754869034E+00
|
|
non_local_psp : 0.00000000000000E+00
|
|
total_energy : -1.03999267053394E+00
|
|
total_energy_eV : -2.82996397694855E+01
|
|
band_energy : -7.51530989588326E-01
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.02569732E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15593581E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15593581E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -2.8065E-01 GPa]
|
|
- sigma(1 1)= -5.95980627E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 3.40087999E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 3.40087999E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
================================================================================
|
|
== DATASET 3 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 3, }
|
|
dimensions: {natom: 2, nkpt: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 280, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 2.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t07o_DS2_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 559.000 559.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 3, }
|
|
solver: {iscf: 7, nstep: 30, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -1.0401468275192 -1.040E+00 6.341E-13 5.172E-04
|
|
ETOT 2 -1.0401468406329 -1.311E-08 2.100E-14 6.790E-05
|
|
ETOT 3 -1.0401468420246 -1.392E-09 8.252E-12 1.098E-06
|
|
ETOT 4 -1.0401468420302 -5.663E-12 8.242E-15 3.305E-07
|
|
ETOT 5 -1.0401468420329 -2.683E-12 1.501E-14 8.458E-08
|
|
ETOT 6 -1.0401468420329 -9.992E-15 2.441E-17 9.101E-09
|
|
|
|
At SCF step 6 vres2 = 9.10E-09 < tolvrs= 1.00E-08 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.02031183E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15513603E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15513603E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 3, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 12.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 10.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 10.0000000, ]
|
|
lattice_lengths: [ 12.00000, 10.00000, 10.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.2000000E+03
|
|
convergence: {deltae: -9.992E-15, res2: 9.101E-09, residm: 2.441E-17, diffor: null, }
|
|
etotal : -1.04014684E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -3.75691325E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -2.02031183E-04, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 1.15513603E-04, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 1.15513603E-04, ]
|
|
pressure_GPa: -2.8436E-01
|
|
xred :
|
|
- [ -4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
- [ 4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.56558260E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ 2.56558260E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
force_length_stats: {min: 2.56558260E-01, max: 2.56558260E-01, mean: 2.56558260E-01, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 1.58463713
|
|
2 2.00000 1.58463713
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 24.405E-18; max= 24.405E-18
|
|
reduced coordinates (array xred) for 2 atoms
|
|
-0.047000000000 0.000000000000 0.000000000000
|
|
0.047000000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 1.7775E+00; max dE/dt= 3.0787E+00; dE/dt below (all hartree)
|
|
1 3.078699120236 0.000000000000 0.000000000000
|
|
2 -3.078699114878 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 -0.29845594564476 0.00000000000000 0.00000000000000
|
|
2 0.29845594564476 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.25655825979639 -0.00000000000000 -0.00000000000000
|
|
2 0.25655825979639 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 1.4812398E-01 2.5655826E-01 -2.232E-10 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -13.19275504773714 -0.00000000000000 -0.00000000000000
|
|
2 13.19275504773714 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 7.6168407E+00 1.3192755E+01 -1.148E-08 0.000E+00 0.000E+00 e/A
|
|
length scales= 12.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 6.350126503080 5.291772085900 5.291772085900 angstroms
|
|
prteigrs : about to open file t07o_DS3_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.37569 Average Vxc (hartree)= -0.06016
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 1, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.37569
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 3, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 1.09996145374418E+00
|
|
hartree : 8.32248811449797E-01
|
|
xc : -6.51934916146786E-01
|
|
Ewald energy : 3.59741353824815E-01
|
|
psp_core : 1.39956766170961E-04
|
|
local_psp : -2.68030350167111E+00
|
|
non_local_psp : 0.00000000000000E+00
|
|
total_energy : -1.04014684203294E+00
|
|
total_energy_eV : -2.83038349893250E+01
|
|
band_energy : -7.51382649728407E-01
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.02031183E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15513603E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15513603E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -2.8436E-01 GPa]
|
|
- sigma(1 1)= -5.94396161E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 3.39852696E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 3.39852696E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
================================================================================
|
|
== DATASET 4 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 4, }
|
|
dimensions: {natom: 2, nkpt: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 280, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 7, paral_kgb: 0, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 3.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t07o_DS3_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 559.000 559.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 4, }
|
|
solver: {iscf: 7, nstep: 30, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -1.0403005262448 -1.040E+00 6.360E-13 5.109E-04
|
|
ETOT 2 -1.0403005393645 -1.312E-08 2.108E-14 6.718E-05
|
|
ETOT 3 -1.0403005407561 -1.392E-09 8.263E-12 1.090E-06
|
|
ETOT 4 -1.0403005407619 -5.823E-12 8.620E-15 3.264E-07
|
|
ETOT 5 -1.0403005407646 -2.632E-12 1.487E-14 8.228E-08
|
|
ETOT 6 -1.0403005407645 4.885E-15 2.503E-17 7.165E-09
|
|
|
|
At SCF step 6 vres2 = 7.16E-09 < tolvrs= 1.00E-08 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.01493217E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15433671E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15433671E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 4, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 12.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 10.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 10.0000000, ]
|
|
lattice_lengths: [ 12.00000, 10.00000, 10.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.2000000E+03
|
|
convergence: {deltae: 4.885E-15, res2: 7.165E-09, residm: 2.503E-17, diffor: null, }
|
|
etotal : -1.04030054E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -3.75617181E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -2.01493217E-04, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 1.15433671E-04, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 1.15433671E-04, ]
|
|
pressure_GPa: -2.8807E-01
|
|
xred :
|
|
- [ -4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
- [ 4.7050E-02, 0.0000E+00, 0.0000E+00, H]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.55771098E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ 2.55771098E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
force_length_stats: {min: 2.55771098E-01, max: 2.55771098E-01, mean: 2.55771098E-01, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 1.58437552
|
|
2 2.00000 1.58453540
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 25.025E-18; max= 25.025E-18
|
|
reduced coordinates (array xred) for 2 atoms
|
|
-0.047000000000 0.000000000000 0.000000000000
|
|
0.047050000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 1.7720E+00; max dE/dt= 3.0693E+00; dE/dt below (all hartree)
|
|
1 3.069252928379 0.000000000000 0.000000000000
|
|
2 -3.069253419286 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 -0.29845594564476 0.00000000000000 0.00000000000000
|
|
2 0.29877345196991 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.25577109781936 -0.00000000000000 -0.00000000000000
|
|
2 0.25577109781936 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 1.4766951E-01 2.5577110E-01 2.045E-08 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -13.15227755481171 -0.00000000000000 -0.00000000000000
|
|
2 13.15227755481171 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 7.5934710E+00 1.3152278E+01 1.052E-06 0.000E+00 0.000E+00 e/A
|
|
length scales= 12.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 6.350126503080 5.291772085900 5.291772085900 angstroms
|
|
prteigrs : about to open file t07o_DS4_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.37562 Average Vxc (hartree)= -0.06017
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 1, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.37562
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 4, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 1.09964611531279E+00
|
|
hartree : 8.32052920505909E-01
|
|
xc : -6.51841225832714E-01
|
|
Ewald energy : 3.59271219424895E-01
|
|
psp_core : 1.39956766170961E-04
|
|
local_psp : -2.67956952694159E+00
|
|
non_local_psp : 0.00000000000000E+00
|
|
total_energy : -1.04030054076455E+00
|
|
total_energy_eV : -2.83080173445097E+01
|
|
band_energy : -7.51234361161007E-01
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.01493217E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15433671E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15433671E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -2.8807E-01 GPa]
|
|
- sigma(1 1)= -5.92813411E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 3.39617528E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 3.39617528E+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: 2, nkpt: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 280, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.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:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
--------------------------------------------------------------------------------
|
|
|
|
-inwffil : will read wavefunctions from disk file t07o_DS4_WFK
|
|
_setup2: Arith. and geom. avg. npw (full set) are 559.000 559.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 5, }
|
|
solver: {iscf: 7, nstep: 30, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm vres2
|
|
ETOT 1 -1.0404537531678 -1.040E+00 6.379E-13 4.980E-04
|
|
ETOT 2 -1.0404537662937 -1.313E-08 2.115E-14 6.547E-05
|
|
ETOT 3 -1.0404537676852 -1.392E-09 8.266E-12 1.067E-06
|
|
ETOT 4 -1.0404537676910 -5.848E-12 8.896E-15 3.156E-07
|
|
ETOT 5 -1.0404537676936 -2.562E-12 1.450E-14 8.017E-08
|
|
ETOT 6 -1.0404537676936 -1.577E-14 2.461E-17 6.978E-09
|
|
|
|
At SCF step 6 vres2 = 6.98E-09 < tolvrs= 1.00E-08 =>converged.
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.00955833E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15353787E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15353787E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 5, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 12.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 10.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 10.0000000, ]
|
|
lattice_lengths: [ 12.00000, 10.00000, 10.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 1.2000000E+03
|
|
convergence: {deltae: -1.577E-14, res2: 6.978E-09, residm: 2.461E-17, diffor: null, }
|
|
etotal : -1.04045377E+00
|
|
entropy : 0.00000000E+00
|
|
fermie : -3.75543060E-01
|
|
cartesian_stress_tensor: # hartree/bohr^3
|
|
- [ -2.00955833E-04, 0.00000000E+00, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 1.15353787E-04, 0.00000000E+00, ]
|
|
- [ 0.00000000E+00, 0.00000000E+00, 1.15353787E-04, ]
|
|
pressure_GPa: -2.9178E-01
|
|
xred :
|
|
- [ -4.7000E-02, 0.0000E+00, 0.0000E+00, H]
|
|
- [ 4.7100E-02, 0.0000E+00, 0.0000E+00, H]
|
|
cartesian_forces: # hartree/bohr
|
|
- [ -2.54985543E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
- [ 2.54985543E-01, -0.00000000E+00, -0.00000000E+00, ]
|
|
force_length_stats: {min: 2.54985543E-01, max: 2.54985543E-01, mean: 2.54985543E-01, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 2.00000 1.58411377
|
|
2 2.00000 1.58443358
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 24.613E-18; max= 24.613E-18
|
|
reduced coordinates (array xred) for 2 atoms
|
|
-0.047000000000 0.000000000000 0.000000000000
|
|
0.047100000000 0.000000000000 0.000000000000
|
|
rms dE/dt= 1.7666E+00; max dE/dt= 3.0598E+00; dE/dt below (all hartree)
|
|
1 3.059826018937 0.000000000000 0.000000000000
|
|
2 -3.059827005363 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 -0.29845594564476 0.00000000000000 0.00000000000000
|
|
2 0.29909095829507 0.00000000000000 0.00000000000000
|
|
|
|
cartesian forces (hartree/bohr) at end:
|
|
1 -0.25498554267916 -0.00000000000000 -0.00000000000000
|
|
2 0.25498554267916 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 1.4721597E-01 2.5498554E-01 4.110E-08 0.000E+00 0.000E+00 h/b
|
|
|
|
cartesian forces (eV/Angstrom) at end:
|
|
1 -13.11188268875130 -0.00000000000000 -0.00000000000000
|
|
2 13.11188268875130 -0.00000000000000 -0.00000000000000
|
|
frms,max,avg= 7.5701490E+00 1.3111883E+01 2.114E-06 0.000E+00 0.000E+00 e/A
|
|
length scales= 12.000000000000 10.000000000000 10.000000000000 bohr
|
|
= 6.350126503080 5.291772085900 5.291772085900 angstroms
|
|
prteigrs : about to open file t07o_DS5_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.37554 Average Vxc (hartree)= -0.06018
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 1, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.37554
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 5, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 1.09933082212860E+00
|
|
hartree : 8.31857008408353E-01
|
|
xc : -6.51747526268826E-01
|
|
Ewald energy : 3.58801586583330E-01
|
|
psp_core : 1.39956766170961E-04
|
|
local_psp : -2.67883561531124E+00
|
|
non_local_psp : 0.00000000000000E+00
|
|
total_energy : -1.04045376769361E+00
|
|
total_energy_eV : -2.83121868612941E+01
|
|
band_energy : -7.51086119091980E-01
|
|
...
|
|
|
|
|
|
Cartesian components of stress tensor (hartree/bohr^3)
|
|
sigma(1 1)= -2.00955833E-04 sigma(3 2)= 0.00000000E+00
|
|
sigma(2 2)= 1.15353787E-04 sigma(3 1)= 0.00000000E+00
|
|
sigma(3 3)= 1.15353787E-04 sigma(2 1)= 0.00000000E+00
|
|
|
|
-Cartesian components of stress tensor (GPa) [Pressure= -2.9178E-01 GPa]
|
|
- sigma(1 1)= -5.91232373E+00 sigma(3 2)= 0.00000000E+00
|
|
- sigma(2 2)= 3.39382503E+00 sigma(3 1)= 0.00000000E+00
|
|
- sigma(3 3)= 3.39382503E+00 sigma(2 1)= 0.00000000E+00
|
|
|
|
================================================================================
|
|
== DATASET 6 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 6, }
|
|
dimensions: {natom: 2, nkpt: 1, mband: 1, nsppol: 1, nspinor: 1, nspden: 1, mpw: 559, }
|
|
cutoff_energies: {ecut: 4.5, pawecutdg: -1.0, }
|
|
electrons: {nelect: 2.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 1, rfphon: 1, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 3.
|
|
|
|
Exchange-correlation functional for the present dataset will be:
|
|
GGA: x-only Perdew-Burke-Ernzerhof functional - ixc=12
|
|
Citation for XC functional:
|
|
J.P.Perdew, K.Burke, M.Ernzerhof, PRL 77, 3865 (1996)
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 12.0000000 0.0000000 0.0000000 G(1)= 0.0833333 0.0000000 0.0000000
|
|
R(2)= 0.0000000 10.0000000 0.0000000 G(2)= 0.0000000 0.1000000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 10.0000000 G(3)= 0.0000000 0.0000000 0.1000000
|
|
Unit cell volume ucvol= 1.2000000E+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.0000 ngfft= 24 20 20
|
|
ecut(hartree)= 4.500 => boxcut(ratio)= 2.09440
|
|
--------------------------------------------------------------------------------
|
|
|
|
|
|
==> 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.000000
|
|
Perturbation : displacement of atom 2 along direction 1
|
|
Found 4 symmetries that leave the perturbation invariant.
|
|
symkpt : not enough symmetry to change the number of k points.
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
Initialisation of the first-order wave-functions :
|
|
ireadwf= 0
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 6, }
|
|
solver: {iscf: 7, nstep: 20, nline: 3, wfoptalg: 0, }
|
|
tolerances: {tolvrs: 1.00E-08, }
|
|
...
|
|
|
|
iter 2DEtotal(Ha) deltaE(Ha) residm vres2
|
|
-ETOT 1 200.51102648047 -9.477E+01 3.630E-01 1.497E+05
|
|
ETOT 2 190.26251197478 -1.025E+01 1.533E-02 4.728E+04
|
|
ETOT 3 189.23359892371 -1.029E+00 5.137E-03 1.295E+03
|
|
ETOT 4 189.10792902270 -1.257E-01 3.380E-04 1.301E+02
|
|
ETOT 5 189.10693078117 -9.982E-04 2.680E-06 1.818E+01
|
|
ETOT 6 189.10690809039 -2.269E-05 1.529E-08 4.599E+00
|
|
ETOT 7 189.10690404282 -4.048E-06 9.512E-09 3.687E-01
|
|
ETOT 8 189.10690392845 -1.144E-07 3.785E-10 9.271E-05
|
|
ETOT 9 189.10690392794 -5.026E-10 1.311E-12 2.284E-03
|
|
ETOT 10 189.10690392774 -2.003E-10 2.525E-13 2.012E-05
|
|
ETOT 11 189.10690392774 -2.970E-12 2.414E-15 1.498E-07
|
|
ETOT 12 189.10690392774 4.263E-13 1.186E-16 7.871E-09
|
|
|
|
At SCF step 12 vres2 = 7.87E-09 < tolvrs= 1.00E-08 =>converged.
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 11.860E-17; max= 11.860E-17
|
|
|
|
Thirteen components of 2nd-order total energy (hartree) are
|
|
1,2,3: 0th-order hamiltonian combined with 1st-order wavefunctions
|
|
kin0= 1.11225728E+02 eigvalue= 3.17522165E+01 local= -5.86118957E+01
|
|
4,5,6: 1st-order hamiltonian combined with 1st and 0th-order wfs
|
|
loc psp = -2.12352208E+02 Hartree= 4.01684966E+01 xc= -1.83584408E+01
|
|
note that "loc psp" includes a xc core correction that could be resolved
|
|
7,8,9: eventually, occupation + non-local contributions
|
|
edocc= 0.00000000E+00 enl0= 0.00000000E+00 enl1= 0.00000000E+00
|
|
1-9 gives the relaxation energy (to be shifted if some occ is /=2.0)
|
|
erelax= -1.06176104E+02
|
|
10,11,12 Non-relaxation contributions : frozen-wavefunctions and Ewald
|
|
fr.local= 9.43399143E+01 fr.nonlo= 0.00000000E+00 Ewald= 2.00943093E+02
|
|
13,14 Frozen wf xc core corrections (1) and (2)
|
|
frxc 1 = 0.00000000E+00 frxc 2 = 0.00000000E+00
|
|
Resulting in :
|
|
2DEtotal= 0.1891069039E+03 Ha. Also 2DEtotal= 0.514586055335E+04 eV
|
|
(2DErelax= -1.0617610369E+02 Ha. 2DEnonrelax= 2.9528300762E+02 Ha)
|
|
( non-var. 2DEtotal : 1.8910690359E+02 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 2 1 0.0000000000 0.0000000000
|
|
1 1 3 1 -0.0000000000 0.0000000000
|
|
1 1 1 2 -189.1167560782 0.0000000000
|
|
1 1 2 2 -0.0000000000 0.0000000000
|
|
1 1 3 2 0.0000000000 -0.0000000000
|
|
1 1 2 4 0.0000000000 0.0000000000
|
|
1 1 3 4 0.0000000000 0.0000000000
|
|
|
|
2 1 1 1 0.0000000000 0.0000000000
|
|
2 1 3 1 -0.0000000000 0.0000000000
|
|
2 1 1 2 -0.0000000000 0.0000000000
|
|
2 1 3 2 0.0000000000 0.0000000000
|
|
2 1 1 4 0.0000000000 0.0000000000
|
|
2 1 3 4 0.0000000000 0.0000000000
|
|
|
|
3 1 1 1 -0.0000000000 0.0000000000
|
|
3 1 2 1 -0.0000000000 0.0000000000
|
|
3 1 1 2 -0.0000000000 -0.0000000000
|
|
3 1 2 2 0.0000000000 0.0000000000
|
|
3 1 1 4 0.0000000000 0.0000000000
|
|
3 1 2 4 0.0000000000 0.0000000000
|
|
|
|
1 2 1 1 -189.1167560782 -0.0000000000
|
|
1 2 2 1 -0.0000000000 -0.0000000000
|
|
1 2 3 1 -0.0000000000 0.0000000000
|
|
1 2 1 2 189.1069035883 0.0000000000
|
|
1 2 2 2 -0.0000000000 0.0000000000
|
|
1 2 3 2 0.0000000000 0.0000000000
|
|
1 2 2 4 0.0000000000 0.0000000000
|
|
1 2 3 4 0.0000000000 0.0000000000
|
|
|
|
2 2 1 1 -0.0000000000 -0.0000000000
|
|
2 2 3 1 0.0000000000 -0.0000000000
|
|
2 2 1 2 -0.0000000000 0.0000000000
|
|
2 2 3 2 -0.0000000000 0.0000000000
|
|
2 2 1 4 0.0000000000 0.0000000000
|
|
2 2 3 4 0.0000000000 0.0000000000
|
|
|
|
3 2 1 1 0.0000000000 0.0000000000
|
|
3 2 2 1 0.0000000000 -0.0000000000
|
|
3 2 1 2 0.0000000000 0.0000000000
|
|
3 2 2 2 -0.0000000000 0.0000000000
|
|
3 2 1 4 0.0000000000 0.0000000000
|
|
3 2 2 4 0.0000000000 0.0000000000
|
|
|
|
1 4 2 1 0.0000000000 0.0000000000
|
|
1 4 3 1 0.0000000000 0.0000000000
|
|
1 4 2 2 0.0000000000 0.0000000000
|
|
1 4 3 2 0.0000000000 0.0000000000
|
|
1 4 2 4 0.0000000000 0.0000000000
|
|
1 4 3 4 0.0000000000 0.0000000000
|
|
|
|
2 4 1 1 0.0000000000 0.0000000000
|
|
2 4 3 1 0.0000000000 0.0000000000
|
|
2 4 1 2 0.0000000000 0.0000000000
|
|
2 4 3 2 0.0000000000 0.0000000000
|
|
2 4 1 4 0.0000000000 0.0000000000
|
|
2 4 3 4 0.0000000000 0.0000000000
|
|
|
|
3 4 1 1 0.0000000000 0.0000000000
|
|
3 4 2 1 0.0000000000 0.0000000000
|
|
3 4 1 2 0.0000000000 0.0000000000
|
|
3 4 2 2 0.0000000000 0.0000000000
|
|
3 4 1 4 0.0000000000 0.0000000000
|
|
3 4 2 4 0.0000000000 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 2 1 0.0000000000 0.0000000000
|
|
1 1 3 1 -0.0000000000 0.0000000000
|
|
1 1 1 2 -1.3133108061 0.0000000000
|
|
1 1 2 2 -0.0000000000 0.0000000000
|
|
1 1 3 2 0.0000000000 -0.0000000000
|
|
|
|
2 1 1 1 0.0000000000 0.0000000000
|
|
2 1 3 1 -0.0000000000 0.0000000000
|
|
2 1 1 2 -0.0000000000 0.0000000000
|
|
2 1 3 2 0.0000000000 0.0000000000
|
|
|
|
3 1 1 1 -0.0000000000 0.0000000000
|
|
3 1 2 1 -0.0000000000 0.0000000000
|
|
3 1 1 2 -0.0000000000 -0.0000000000
|
|
3 1 2 2 0.0000000000 0.0000000000
|
|
|
|
1 2 1 1 -1.3133108061 -0.0000000000
|
|
1 2 2 1 -0.0000000000 -0.0000000000
|
|
1 2 3 1 -0.0000000000 0.0000000000
|
|
1 2 1 2 1.3132423860 0.0000000000
|
|
1 2 2 2 -0.0000000000 0.0000000000
|
|
1 2 3 2 0.0000000000 0.0000000000
|
|
|
|
2 2 1 1 -0.0000000000 -0.0000000000
|
|
2 2 3 1 0.0000000000 -0.0000000000
|
|
2 2 1 2 -0.0000000000 0.0000000000
|
|
2 2 3 2 -0.0000000000 0.0000000000
|
|
|
|
3 2 1 1 0.0000000000 0.0000000000
|
|
3 2 2 1 0.0000000000 -0.0000000000
|
|
3 2 1 2 0.0000000000 0.0000000000
|
|
3 2 2 2 -0.0000000000 0.0000000000
|
|
|
|
Phonon wavevector (reduced coordinates) : 0.00000 0.00000 0.00000
|
|
Phonon energies in Hartree :
|
|
0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
|
|
2.673389E-02
|
|
Phonon frequencies in cm-1 :
|
|
- 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
|
|
- 5.867410E+03
|
|
chkph3 : WARNING -
|
|
Dynamical matrix incomplete, phonon frequencies may be wrong, see the log file for more explanations.
|
|
|
|
== END DATASET(S) ==============================================================
|
|
================================================================================
|
|
|
|
-outvars: echo values of variables after computation --------
|
|
acell 1.2000000000E+01 1.0000000000E+01 1.0000000000E+01 Bohr
|
|
amu 1.00800000E+00
|
|
asr 0
|
|
chneut 0
|
|
diemac 1.00000000E+00
|
|
diemix 5.00000000E-01
|
|
ecut 4.50000000E+00 Hartree
|
|
etotal1 -1.0398380253E+00
|
|
etotal2 -1.0399926705E+00
|
|
etotal3 -1.0401468420E+00
|
|
etotal4 -1.0403005408E+00
|
|
etotal5 -1.0404537677E+00
|
|
etotal6 1.8910690393E+02
|
|
fcart1 -2.5813739487E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
2.5813739487E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
fcart2 -2.5734703069E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
2.5734703069E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
fcart3 -2.5655825980E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
2.5655825980E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
fcart4 -2.5577109782E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
2.5577109782E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
fcart5 -2.5498554268E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
2.5498554268E-01 -0.0000000000E+00 -0.0000000000E+00
|
|
fcart6 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
- fftalg 512
|
|
getwfk1 -1
|
|
getwfk2 -1
|
|
getwfk3 -1
|
|
getwfk4 -1
|
|
getwfk5 -1
|
|
getwfk6 3
|
|
istwfk1 2
|
|
istwfk2 2
|
|
istwfk3 2
|
|
istwfk4 2
|
|
istwfk5 2
|
|
istwfk6 1
|
|
ixc 12
|
|
jdtset 1 2 3 4 5 6
|
|
kptopt 0
|
|
P mkmem 1
|
|
P mkqmem 1
|
|
P mk1mem 1
|
|
natom 2
|
|
nband 1
|
|
ndtset 6
|
|
ngfft 24 20 20
|
|
nkpt 1
|
|
nline 3
|
|
nqpt1 0
|
|
nqpt2 0
|
|
nqpt3 0
|
|
nqpt4 0
|
|
nqpt5 0
|
|
nqpt6 1
|
|
nstep1 30
|
|
nstep2 30
|
|
nstep3 30
|
|
nstep4 30
|
|
nstep5 30
|
|
nstep6 20
|
|
nsym 4
|
|
ntypat 1
|
|
occ 2.000000
|
|
optdriver1 0
|
|
optdriver2 0
|
|
optdriver3 0
|
|
optdriver4 0
|
|
optdriver5 0
|
|
optdriver6 1
|
|
prtpot1 0
|
|
prtpot2 0
|
|
prtpot3 0
|
|
prtpot4 0
|
|
prtpot5 0
|
|
prtpot6 1
|
|
rfatpol1 1 2
|
|
rfatpol2 1 2
|
|
rfatpol3 1 2
|
|
rfatpol4 1 2
|
|
rfatpol5 1 2
|
|
rfatpol6 2 2
|
|
rfdir1 1 1 1
|
|
rfdir2 1 1 1
|
|
rfdir3 1 1 1
|
|
rfdir4 1 1 1
|
|
rfdir5 1 1 1
|
|
rfdir6 1 0 0
|
|
rfphon1 0
|
|
rfphon2 0
|
|
rfphon3 0
|
|
rfphon4 0
|
|
rfphon5 0
|
|
rfphon6 1
|
|
spgroup 25
|
|
strten1 -2.0310887023E-04 1.1567355794E-04 1.1567355790E-04
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten2 -2.0256973183E-04 1.1559358085E-04 1.1559358084E-04
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten3 -2.0203118287E-04 1.1551360298E-04 1.1551360297E-04
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten4 -2.0149321696E-04 1.1543367100E-04 1.1543367100E-04
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten5 -2.0095583284E-04 1.1535378745E-04 1.1535378744E-04
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
strten6 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
symrel 1 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 0 -1
|
|
1 0 0 0 -1 0 0 0 1 1 0 0 0 -1 0 0 0 -1
|
|
tolvrs 1.00000000E-08
|
|
typat 1 1
|
|
xangst1 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9782093299E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst2 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9813843932E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst3 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst4 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9877345197E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst5 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9909095830E-01 0.0000000000E+00 0.0000000000E+00
|
|
xangst6 -2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
2.9845594564E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart1 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6280000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart2 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6340000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart3 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart4 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6460000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart5 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6520000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xcart6 -5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
5.6400000000E-01 0.0000000000E+00 0.0000000000E+00
|
|
xred1 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.6900000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred2 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.6950000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred3 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred4 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7050000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred5 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7100000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
xred6 -4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
4.7000000000E-02 0.0000000000E+00 0.0000000000E+00
|
|
znucl 1.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] 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
|
|
-
|
|
- [5] 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= 1.3 wall= 1.4
|
|
|
|
================================================================================
|
|
|
|
Calculation completed.
|
|
.Delivered 40 WARNINGs and 25 COMMENTs to log file.
|
|
+Overall time at end (sec) : cpu= 1.3 wall= 1.4
|