mirror of https://github.com/abinit/abinit.git
911 lines
44 KiB
Plaintext
911 lines
44 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 19h15 )
|
|
|
|
- input file -> /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/TestBot_MPI1/v9_t44/t44.abi
|
|
- output file -> t44.abo
|
|
- root for input files -> t44i
|
|
- root for output files -> t44o
|
|
|
|
- inpspheads : Reading pseudopotential header in XML form from
|
|
- /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
DATASET 1 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 1.
|
|
intxc = 0 ionmov = 0 iscf = 17 lmnmax = 8
|
|
lnmax = 4 mgfft = 45 mpssoang = 2 mqgrid = 3001
|
|
natom = 1 nloc_mem = 2 nspden = 1 nspinor = 1
|
|
nsppol = 1 nsym = 1 n1xccc = 1 ntypat = 1
|
|
occopt = 1 xclevel = 2
|
|
- mband = 4 mffmem = 1 mkmem = 1
|
|
mpw = 4337 nfft = 91125 nkpt = 1
|
|
PAW method is used; the additional fine FFT grid is defined by:
|
|
mgfftf= 48 nfftf = 110592
|
|
================================================================================
|
|
P This job should need less than 38.922 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.267 Mbytes ; DEN or POT disk file : 0.846 Mbytes.
|
|
================================================================================
|
|
|
|
DATASET 2 : space group P1 (# 1); Bravais aP (primitive triclinic)
|
|
================================================================================
|
|
Values of the parameters that define the memory need for DATASET 2 (RF).
|
|
intxc = 0 iscf = -3 lmnmax = 8 lnmax = 4
|
|
mgfft = 45 mpssoang = 2 mqgrid = 3001 natom = 1
|
|
nloc_mem = 2 nspden = 1 nspinor = 1 nsppol = 1
|
|
nsym = 1 n1xccc = 1 ntypat = 1 occopt = 1
|
|
xclevel = 2
|
|
- mband = 4 mffmem = 1 mkmem = 1
|
|
- mkqmem = 1 mk1mem = 1 mpw = 4337
|
|
nfft = 91125 nkpt = 1
|
|
================================================================================
|
|
P This job should need less than 22.004 Mbytes of memory.
|
|
Rough estimation (10% accuracy) of disk space for files :
|
|
_ WF disk file : 0.267 Mbytes ; DEN or POT disk file : 0.697 Mbytes.
|
|
================================================================================
|
|
|
|
--------------------------------------------------------------------------------
|
|
------------- Echo of variables that govern the present computation ------------
|
|
--------------------------------------------------------------------------------
|
|
-
|
|
- outvars: echo of selected default values
|
|
- iomode0 = 0 , fftalg0 =512 , wfoptalg0 = 10
|
|
-
|
|
- outvars: echo of global parameters not present in the input file
|
|
- max_nthreads = 0
|
|
-
|
|
-outvars: echo values of preprocessed input variables --------
|
|
acell 2.0000000000E+01 2.0000000000E+01 2.0000000000E+01 Bohr
|
|
amu 2.01797000E+01
|
|
ecut 5.00000000E+00 Hartree
|
|
enunit 2
|
|
- fftalg 512
|
|
getwfk1 0
|
|
getwfk2 -1
|
|
iscf1 17
|
|
iscf2 -3
|
|
istwfk 1
|
|
ixc -101130
|
|
jdtset 1 2
|
|
kptopt 0
|
|
lambsig1 0.00000000E+00
|
|
lambsig2 3.41860000E-04
|
|
P mkmem 1
|
|
P mkqmem 1
|
|
P mk1mem 1
|
|
natom 1
|
|
nband 4
|
|
ndtset 2
|
|
ngfft 45 45 45
|
|
ngfftdg 48 48 48
|
|
nkpt 1
|
|
nstep 500
|
|
nsym 1
|
|
ntypat 1
|
|
nucdipmom 1.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
occ 2.000000 2.000000 2.000000 2.000000
|
|
optdriver1 0
|
|
optdriver2 1
|
|
optforces 0
|
|
optstress 0
|
|
orbmag1 0
|
|
orbmag2 2
|
|
pawcpxocc 2
|
|
pawecutdg 6.00000000E+00 Hartree
|
|
prtpot1 0
|
|
prtpot2 1
|
|
rfddk1 0
|
|
rfddk2 1
|
|
spgroup 1
|
|
symmorphi 0
|
|
tolvrs1 1.00000000E-18
|
|
tolvrs2 0.00000000E+00
|
|
tolwfr1 0.00000000E+00
|
|
tolwfr2 1.00000000E-20
|
|
typat 1
|
|
usexcnhat 0
|
|
useylm 1
|
|
xangst 5.2917720859E+00 5.2917720859E+00 5.2917720859E+00
|
|
xcart 1.0000000000E+01 1.0000000000E+01 1.0000000000E+01
|
|
xred 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01
|
|
znucl 10.00000
|
|
|
|
================================================================================
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 1.
|
|
|
|
chkinp: Checking input parameters for consistency, jdtset= 2.
|
|
|
|
================================================================================
|
|
== DATASET 1 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 1, }
|
|
dimensions: {natom: 1, nkpt: 1, mband: 4, nsppol: 1, nspinor: 1, nspden: 1, mpw: 4337, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 6.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 0, ionmov: 0, optcell: 0, iscf: 17, paral_kgb: 0, }
|
|
...
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 20.0000000 0.0000000 0.0000000 G(1)= 0.0500000 0.0000000 0.0000000
|
|
R(2)= 0.0000000 20.0000000 0.0000000 G(2)= 0.0000000 0.0500000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 20.0000000 G(3)= 0.0000000 0.0000000 0.0500000
|
|
Unit cell volume ucvol= 8.0000000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 45 45 45
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.18561
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 48 48 48
|
|
ecut(hartree)= 6.000 => boxcut(ratio)= 2.17656
|
|
|
|
--- Pseudopotential description ------------------------------------------------
|
|
- pspini: atom type 1 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
- pspatm : Reading pseudopotential header in XML form from /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 4 (lmn_size= 8), orbitals= 0 0 1 1
|
|
Spheres core radius: rc_sph= 1.80173543
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=2001 , AA= 0.59080E-03 BB= 0.59080E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60087216
|
|
mmax= 2001
|
|
Radial grid used for partial waves is grid 1
|
|
Radial grid used for projectors is grid 1
|
|
Radial grid used for (t)core density is grid 1
|
|
Radial grid used for Vloc is grid 1
|
|
Radial grid used for pseudo valence density is grid 1
|
|
Mesh size for Vloc has been set to 1766 to avoid numerical noise.
|
|
Compensation charge density is not taken into account in XC energy/potential
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
5.19616736E+01 ecore*ucvol(ha*bohr**3)
|
|
--------------------------------------------------------------------------------
|
|
|
|
_setup2: Arith. and geom. avg. npw (full set) are 4337.000 4337.000
|
|
|
|
================================================================================
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 1, }
|
|
solver: {iscf: 17, nstep: 500, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolvrs: 1.00E-18, }
|
|
...
|
|
|
|
iter Etot(hartree) deltaE(h) residm nres2
|
|
ETOT 1 -33.255360784645 -3.326E+01 3.625E-02 1.466E+01
|
|
ETOT 2 -34.012550385011 -7.572E-01 1.076E-03 3.284E+00
|
|
ETOT 3 -33.708275266402 3.043E-01 1.749E-05 7.110E-01
|
|
ETOT 4 -33.627593781821 8.068E-02 1.004E-05 8.482E-02
|
|
ETOT 5 -33.616724848324 1.087E-02 2.073E-07 2.908E-02
|
|
ETOT 6 -33.610624008243 6.101E-03 2.946E-07 1.279E-03
|
|
ETOT 7 -33.610728999119 -1.050E-04 1.781E-07 1.105E-03
|
|
ETOT 8 -33.611787299872 -1.058E-03 8.726E-08 7.600E-04
|
|
ETOT 9 -33.611643053345 1.442E-04 1.535E-08 1.455E-04
|
|
ETOT 10 -33.611689906363 -4.685E-05 3.655E-09 1.793E-05
|
|
ETOT 11 -33.611882497270 -1.926E-04 7.970E-10 5.996E-06
|
|
ETOT 12 -33.611968263688 -8.577E-05 5.557E-11 3.121E-06
|
|
ETOT 13 -33.612033358844 -6.510E-05 1.141E-10 1.270E-06
|
|
ETOT 14 -33.612058073493 -2.471E-05 4.102E-11 6.423E-07
|
|
ETOT 15 -33.612069821769 -1.175E-05 2.993E-12 3.724E-07
|
|
ETOT 16 -33.612081642562 -1.182E-05 8.517E-12 1.477E-07
|
|
ETOT 17 -33.612086839467 -5.197E-06 3.402E-12 4.985E-08
|
|
ETOT 18 -33.612089088658 -2.249E-06 2.090E-12 1.167E-08
|
|
ETOT 19 -33.612089814006 -7.253E-07 3.069E-13 3.263E-09
|
|
ETOT 20 -33.612090100530 -2.865E-07 4.768E-14 4.633E-10
|
|
ETOT 21 -33.612090128703 -2.817E-08 2.582E-14 1.367E-10
|
|
ETOT 22 -33.612090132216 -3.513E-09 8.123E-16 6.254E-11
|
|
ETOT 23 -33.612090133771 -1.554E-09 6.431E-16 2.362E-11
|
|
ETOT 24 -33.612090134858 -1.088E-09 3.074E-15 2.921E-12
|
|
ETOT 25 -33.612090134971 -1.127E-10 9.335E-17 9.778E-13
|
|
ETOT 26 -33.612090135016 -4.530E-11 1.249E-17 1.207E-13
|
|
ETOT 27 -33.612090135022 -5.677E-12 2.353E-17 1.098E-14
|
|
ETOT 28 -33.612090135022 4.761E-13 1.921E-18 4.593E-15
|
|
ETOT 29 -33.612090135022 -6.892E-13 2.453E-19 1.742E-15
|
|
ETOT 30 -33.612090135023 -7.603E-13 2.099E-19 2.863E-16
|
|
ETOT 31 -33.612090135022 7.603E-13 2.766E-20 1.485E-16
|
|
ETOT 32 -33.612090135023 -5.613E-13 2.279E-20 9.778E-17
|
|
ETOT 33 -33.612090135023 3.411E-13 4.969E-20 5.204E-17
|
|
ETOT 34 -33.612090135022 1.421E-13 5.979E-20 3.753E-17
|
|
ETOT 35 -33.612090135022 -7.105E-14 2.858E-20 1.422E-17
|
|
ETOT 36 -33.612090135023 -3.908E-13 3.522E-21 9.072E-18
|
|
ETOT 37 -33.612090135022 7.674E-13 1.713E-20 1.611E-18
|
|
ETOT 38 -33.612090135023 -9.877E-13 4.591E-22 4.902E-19
|
|
|
|
At SCF step 38 nres2 = 4.90E-19 < tolvrs= 1.00E-18 =>converged.
|
|
|
|
--- !ResultsGS
|
|
iteration_state: {dtset: 1, }
|
|
comment : Summary of ground state results
|
|
lattice_vectors:
|
|
- [ 20.0000000, 0.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 20.0000000, 0.0000000, ]
|
|
- [ 0.0000000, 0.0000000, 20.0000000, ]
|
|
lattice_lengths: [ 20.00000, 20.00000, 20.00000, ]
|
|
lattice_angles: [ 90.000, 90.000, 90.000, ] # degrees, (23, 13, 12)
|
|
lattice_volume: 8.0000000E+03
|
|
convergence: {deltae: -9.877E-13, res2: 4.902E-19, residm: 4.591E-22, diffor: 0.000E+00, }
|
|
etotal : -3.36120901E+01
|
|
entropy : 0.00000000E+00
|
|
fermie : -5.70833084E-01
|
|
cartesian_stress_tensor: null
|
|
pressure_GPa: null
|
|
xred :
|
|
- [ 5.0000E-01, 5.0000E-01, 5.0000E-01, Ne]
|
|
cartesian_forces: null
|
|
force_length_stats: {min: null, max: null, mean: null, }
|
|
...
|
|
|
|
Integrated electronic density in atomic spheres:
|
|
------------------------------------------------
|
|
Atom Sphere_radius Integrated_density
|
|
1 1.80174 7.53346492
|
|
|
|
PAW TEST:
|
|
==== Compensation charge inside spheres ============
|
|
The following values must be close to each other ...
|
|
Compensation charge over spherical meshes = 2.562890881227881
|
|
Compensation charge over fine fft grid = 2.562232923990812
|
|
|
|
==== Results concerning PAW augmentation regions ====
|
|
|
|
Total pseudopotential strength Dij (hartree):
|
|
=== REAL PART:
|
|
-0.07569 0.01106 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.01106 0.14608 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 -0.84377 0.00000 0.00000 0.21729 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 -0.84377 0.00000 0.00000 0.21729 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 -0.84377 0.00000 0.00000 0.21729
|
|
0.00000 0.00000 0.21729 0.00000 0.00000 0.42377 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.21729 0.00000 0.00000 0.42377 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.21729 0.00000 0.00000 0.42377
|
|
=== IMAGINARY PART:
|
|
-0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 -0.00000 0.00057 -0.00000 -0.00000 -0.00052 -0.00000
|
|
0.00000 0.00000 -0.00057 -0.00000 -0.00000 0.00052 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 -0.00052 0.00000 -0.00000 0.00046 -0.00000
|
|
0.00000 0.00000 0.00052 0.00000 0.00000 -0.00046 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
|
|
Total pseudopotential strength Dij (eV):
|
|
=== REAL PART:
|
|
-2.05975 0.30108 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.30108 3.97507 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 -22.96021 0.00000 0.00000 5.91272 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 -22.96021 0.00000 0.00000 5.91272 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 -22.96021 0.00000 0.00000 5.91272
|
|
0.00000 0.00000 5.91272 0.00000 0.00000 11.53125 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 5.91272 0.00000 0.00000 11.53125 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 5.91272 0.00000 0.00000 11.53125
|
|
=== IMAGINARY PART:
|
|
-0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 -0.00000 0.01545 -0.00000 -0.00000 -0.01405 -0.00000
|
|
0.00000 0.00000 -0.01545 -0.00000 -0.00000 0.01405 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 -0.01405 0.00000 -0.00000 0.01259 -0.00000
|
|
0.00000 0.00000 0.01405 0.00000 0.00000 -0.01259 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
|
|
Augmentation waves occupancies Rhoij:
|
|
=== REAL PART:
|
|
2.34657 0.45073 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.45073 0.08657 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
|
|
0.00000 0.00000 2.19709 -0.00000 0.00000 0.45218 0.00000 0.00000
|
|
0.00000 0.00000 -0.00000 2.19709 0.00000 0.00000 0.45218 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 2.19709 0.00000 0.00000 0.45218
|
|
0.00000 0.00000 0.45218 0.00000 0.00000 0.09306 0.00000 0.00000
|
|
0.00000 0.00000 0.00000 0.45218 0.00000 0.00000 0.09306 0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.45218 0.00000 0.00000 0.09306
|
|
=== IMAGINARY PART:
|
|
0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00018 -0.00000 0.00000 -0.00026 -0.00000
|
|
0.00000 0.00000 -0.00018 0.00000 -0.00000 0.00026 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 -0.00000 -0.00000 -0.00000 -0.00000
|
|
0.00000 0.00000 -0.00000 -0.00026 0.00000 -0.00000 -0.00011 -0.00000
|
|
0.00000 0.00000 0.00026 0.00000 0.00000 0.00011 -0.00000 -0.00000
|
|
0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 -0.00000
|
|
|
|
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 18.191E-23; max= 45.909E-23
|
|
reduced coordinates (array xred) for 1 atoms
|
|
0.500000000000 0.500000000000 0.500000000000
|
|
rms dE/dt= 0.0000E+00; max dE/dt= 0.0000E+00; dE/dt below (all hartree)
|
|
1 0.000000000000 0.000000000000 0.000000000000
|
|
|
|
cartesian coordinates (angstrom) at end:
|
|
1 5.29177208590000 5.29177208590000 5.29177208590000
|
|
length scales= 20.000000000000 20.000000000000 20.000000000000 bohr
|
|
= 10.583544171800 10.583544171800 10.583544171800 angstroms
|
|
prteigrs : about to open file t44o_DS1_EIG
|
|
Fermi (or HOMO) energy (hartree) = -0.57083 Average Vxc (hartree)= -0.03454
|
|
Eigenvalues (hartree) for nkpt= 1 k points:
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-1.54900 -0.57173 -0.57128 -0.57083
|
|
Fermi (or HOMO) energy (eV) = -15.53316 Average Vxc (eV)= -0.93982
|
|
Eigenvalues ( eV ) for nkpt= 1 k points:
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-42.15052 -15.55750 -15.54532 -15.53316
|
|
|
|
--- !EnergyTerms
|
|
iteration_state : {dtset: 1, }
|
|
comment : Components of total free energy in Hartree
|
|
kinetic : 8.82280629898322E+00
|
|
hartree : 2.25816005187074E+01
|
|
xc : -2.88349598085601E+00
|
|
Ewald energy : -4.53967596716905E+00
|
|
psp_core : 6.49520920517537E-03
|
|
local_psp : -5.77829347833744E+01
|
|
spherical_terms : 1.83114575780156E-01
|
|
nucl. magn. dipoles : -5.84931565005995E-09
|
|
total_energy : -3.36120901345728E+01
|
|
total_energy_eV : -9.14631486988780E+02
|
|
...
|
|
|
|
|
|
--- !EnergyTermsDC
|
|
iteration_state : {dtset: 1, }
|
|
comment : '"Double-counting" decomposition of free energy'
|
|
band_energy : -6.52568789974358E+00
|
|
Ewald energy : -4.53967596716905E+00
|
|
psp_core : 6.49520920517537E-03
|
|
xc_dc : -2.17230570683999E+01
|
|
spherical_terms : -8.30164408915766E-01
|
|
total_energy_dc : -3.36120901350231E+01
|
|
total_energy_dc_eV : -9.14631487001032E+02
|
|
...
|
|
|
|
|
|
================================================================================
|
|
== DATASET 2 ==================================================================
|
|
- mpi_nproc: 1, omp_nthreads: -1 (-1 if OMP is not activated)
|
|
|
|
|
|
--- !DatasetInfo
|
|
iteration_state: {dtset: 2, }
|
|
dimensions: {natom: 1, nkpt: 1, mband: 4, nsppol: 1, nspinor: 1, nspden: 1, mpw: 4337, }
|
|
cutoff_energies: {ecut: 5.0, pawecutdg: 6.0, }
|
|
electrons: {nelect: 8.00000000E+00, charge: 0.00000000E+00, occopt: 1.00000000E+00, tsmear: 1.00000000E-02, }
|
|
meta: {optdriver: 1, rfddk: 1, }
|
|
...
|
|
|
|
mkfilename : getwfk/=0, take file _WFK from output of DATASET 1.
|
|
|
|
Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1):
|
|
R(1)= 20.0000000 0.0000000 0.0000000 G(1)= 0.0500000 0.0000000 0.0000000
|
|
R(2)= 0.0000000 20.0000000 0.0000000 G(2)= 0.0000000 0.0500000 0.0000000
|
|
R(3)= 0.0000000 0.0000000 20.0000000 G(3)= 0.0000000 0.0000000 0.0500000
|
|
Unit cell volume ucvol= 8.0000000E+03 bohr^3
|
|
Angles (23,13,12)= 9.00000000E+01 9.00000000E+01 9.00000000E+01 degrees
|
|
|
|
Coarse grid specifications (used for wave-functions):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 45 45 45
|
|
ecut(hartree)= 5.000 => boxcut(ratio)= 2.18561
|
|
|
|
Fine grid specifications (used for densities):
|
|
|
|
getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 48 48 48
|
|
ecut(hartree)= 6.000 => boxcut(ratio)= 2.17656
|
|
|
|
--- Pseudopotential description ------------------------------------------------
|
|
- pspini: atom type 1 psp file is /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
- pspatm: opening atomic psp file /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
- pspatm : Reading pseudopotential header in XML form from /home/buildbot/ABINIT3/eos_gnu_13.2_mpich/trunk_merge-10.0/tests/Pspdir/Psdj_paw_pbe_std/Ne.xml
|
|
Pseudopotential format is: paw10
|
|
basis_size (lnmax)= 4 (lmn_size= 8), orbitals= 0 0 1 1
|
|
Spheres core radius: rc_sph= 1.80173543
|
|
1 radial meshes are used:
|
|
- mesh 1: r(i)=AA*[exp(BB*(i-1))-1], size=2001 , AA= 0.59080E-03 BB= 0.59080E-02
|
|
Shapefunction is SIN type: shapef(r)=[sin(pi*r/rshp)/(pi*r/rshp)]**2
|
|
Radius for shape functions = 1.60087216
|
|
mmax= 2001
|
|
Radial grid used for partial waves is grid 1
|
|
Radial grid used for projectors is grid 1
|
|
Radial grid used for (t)core density is grid 1
|
|
Radial grid used for Vloc is grid 1
|
|
Radial grid used for pseudo valence density is grid 1
|
|
Mesh size for Vloc has been set to 1766 to avoid numerical noise.
|
|
Compensation charge density is not taken into account in XC energy/potential
|
|
pspatm: atomic psp has been read and splines computed
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
|
|
==> initialize data related to q vector <==
|
|
|
|
The list of irreducible perturbations for this q vector is:
|
|
1) idir= 1 ipert= 2
|
|
2) idir= 2 ipert= 2
|
|
3) idir= 3 ipert= 2
|
|
|
|
================================================================================
|
|
|
|
--------------------------------------------------------------------------------
|
|
Perturbation wavevector (in red.coord.) 0.000000 0.000000 0.000000
|
|
Perturbation : derivative vs k 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.
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
Initialisation of the first-order wave-functions :
|
|
ireadwf= 0
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: -3, nstep: 500, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolwfr: 1.00E-20, }
|
|
...
|
|
|
|
iter 2DEtotal(Ha) deltaE(Ha) residm vres2
|
|
-ETOT 1 -0.46224064118488 -4.622E-01 8.891E-05 0.000E+00
|
|
ETOT 2 -0.46226092441753 -2.028E-05 4.787E-08 0.000E+00
|
|
ETOT 3 -0.46226092849039 -4.073E-09 4.509E-12 0.000E+00
|
|
ETOT 4 -0.46226092849184 -1.455E-12 3.795E-15 0.000E+00
|
|
ETOT 5 -0.46226092849183 7.327E-15 7.274E-19 0.000E+00
|
|
ETOT 6 -0.46226092849184 -4.163E-15 8.861E-21 0.000E+00
|
|
|
|
At SCF step 6 max residual= 8.86E-21 < tolwfr= 1.00E-20 =>converged.
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 41.500E-22; max= 88.614E-22
|
|
dfpt_looppert : ek2= 1.5791367042E+00
|
|
f-sum rule ratio= 5.6448323357E-01
|
|
prteigrs : about to open file t44t_1WF1_EIG
|
|
Expectation of eigenvalue derivatives (hartree) for nkpt= 1 k points:
|
|
(in case of degenerate eigenvalues, averaged derivative)
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 -0.00000 0.00000 0.00000
|
|
Expectation of eigenvalue derivatives ( eV ) for nkpt= 1 k points:
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 -0.00000 0.00000 0.00000
|
|
|
|
Nine components of 2nd-order total energy (hartree) are
|
|
1,2,3: 0th-order hamiltonian combined with 1st-order wavefunctions
|
|
kin0= 4.07363314E-01 eigvalue= 1.52974216E-01 local= -1.98583979E-01 nclr dpl0= 1.18196114E-09
|
|
4,5,6: 1st-order hamiltonian combined with 1st and 0th-order wfs
|
|
kin1= -8.91396193E-01 Hartree= 0.00000000E+00 xc= 0.00000000E+00 nclr dpl1= 0.00000000E+00
|
|
7,8,9: eventually, occupation + non-local contributions
|
|
edocc= 0.00000000E+00 enl0= -1.09482618E-02 enl1= 7.83299735E-02
|
|
10: eventually, PAW "on-site" Hxc contribution: epaw1= 0.00000000E+00
|
|
1-10 gives the relaxation energy (to be shifted if some occ is /=2.0)
|
|
erelax= -4.62260928E-01
|
|
11 Contribution from 1st-order change of wavefunctions overlap
|
|
eovl1 = -5.57278188E-02
|
|
No Ewald or frozen-wf contrib.: the relaxation energy is the total one
|
|
2DEtotal= -0.4622609285E+00 Ha. Also 2DEtotal= -0.125787595687E+02 eV
|
|
( non-var. 2DEtotal : -4.6226092853E-01 Ha)
|
|
|
|
--------------------------------------------------------------------------------
|
|
Perturbation wavevector (in red.coord.) 0.000000 0.000000 0.000000
|
|
Perturbation : derivative vs k along direction 2
|
|
The set of symmetries contains only one element for this perturbation.
|
|
symkpt : not enough symmetry to change the number of k points.
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
Initialisation of the first-order wave-functions :
|
|
ireadwf= 0
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: -3, nstep: 500, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolwfr: 1.00E-20, }
|
|
...
|
|
|
|
iter 2DEtotal(Ha) deltaE(Ha) residm vres2
|
|
-ETOT 1 -0.46223821341521 -4.622E-01 7.283E-05 0.000E+00
|
|
ETOT 2 -0.46226093150842 -2.272E-05 2.705E-08 0.000E+00
|
|
ETOT 3 -0.46226093584214 -4.334E-09 2.958E-12 0.000E+00
|
|
ETOT 4 -0.46226093584369 -1.552E-12 1.963E-15 0.000E+00
|
|
ETOT 5 -0.46226093584369 -2.165E-15 4.582E-19 0.000E+00
|
|
ETOT 6 -0.46226093584369 7.216E-16 4.996E-21 0.000E+00
|
|
|
|
At SCF step 6 max residual= 5.00E-21 < tolwfr= 1.00E-20 =>converged.
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 38.467E-22; max= 49.960E-22
|
|
dfpt_looppert : ek2= 1.5791367042E+00
|
|
f-sum rule ratio= 5.6448323383E-01
|
|
prteigrs : about to open file t44t_1WF1_EIG
|
|
Expectation of eigenvalue derivatives (hartree) for nkpt= 1 k points:
|
|
(in case of degenerate eigenvalues, averaged derivative)
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 -0.00000 -0.00000 0.00000
|
|
Expectation of eigenvalue derivatives ( eV ) for nkpt= 1 k points:
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 -0.00000 -0.00000 0.00000
|
|
|
|
Nine components of 2nd-order total energy (hartree) are
|
|
1,2,3: 0th-order hamiltonian combined with 1st-order wavefunctions
|
|
kin0= 4.07363355E-01 eigvalue= 1.52974247E-01 local= -1.98584005E-01 nclr dpl0= 2.22455655E-09
|
|
4,5,6: 1st-order hamiltonian combined with 1st and 0th-order wfs
|
|
kin1= -8.91396193E-01 Hartree= 0.00000000E+00 xc= 0.00000000E+00 nclr dpl1= -1.18098012E-09
|
|
7,8,9: eventually, occupation + non-local contributions
|
|
edocc= 0.00000000E+00 enl0= -1.09482608E-02 enl1= 7.83299186E-02
|
|
10: eventually, PAW "on-site" Hxc contribution: epaw1= 0.00000000E+00
|
|
1-10 gives the relaxation energy (to be shifted if some occ is /=2.0)
|
|
erelax= -4.62260936E-01
|
|
11 Contribution from 1st-order change of wavefunctions overlap
|
|
eovl1 = -5.57277985E-02
|
|
No Ewald or frozen-wf contrib.: the relaxation energy is the total one
|
|
2DEtotal= -0.4622609358E+00 Ha. Also 2DEtotal= -0.125787597687E+02 eV
|
|
( non-var. 2DEtotal : -4.6226093589E-01 Ha)
|
|
|
|
--------------------------------------------------------------------------------
|
|
Perturbation wavevector (in red.coord.) 0.000000 0.000000 0.000000
|
|
Perturbation : derivative vs k along direction 3
|
|
The set of symmetries contains only one element for this perturbation.
|
|
symkpt : not enough symmetry to change the number of k points.
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
--------------------------------------------------------------------------------
|
|
|
|
Initialisation of the first-order wave-functions :
|
|
ireadwf= 0
|
|
|
|
--- !BeginCycle
|
|
iteration_state: {dtset: 2, }
|
|
solver: {iscf: -3, nstep: 500, nline: 4, wfoptalg: 10, }
|
|
tolerances: {tolwfr: 1.00E-20, }
|
|
...
|
|
|
|
iter 2DEtotal(Ha) deltaE(Ha) residm vres2
|
|
-ETOT 1 -0.46223821340888 -4.622E-01 7.283E-05 0.000E+00
|
|
ETOT 2 -0.46226093150209 -2.272E-05 2.705E-08 0.000E+00
|
|
ETOT 3 -0.46226093583580 -4.334E-09 2.958E-12 0.000E+00
|
|
ETOT 4 -0.46226093583735 -1.555E-12 1.963E-15 0.000E+00
|
|
ETOT 5 -0.46226093583735 -1.443E-15 4.582E-19 0.000E+00
|
|
ETOT 6 -0.46226093583735 2.220E-16 4.996E-21 0.000E+00
|
|
|
|
At SCF step 6 max residual= 5.00E-21 < tolwfr= 1.00E-20 =>converged.
|
|
================================================================================
|
|
|
|
----iterations are completed or convergence reached----
|
|
|
|
Mean square residual over all n,k,spin= 38.467E-22; max= 49.960E-22
|
|
dfpt_looppert : ek2= 1.5791367042E+00
|
|
f-sum rule ratio= 5.6448323382E-01
|
|
prteigrs : about to open file t44t_1WF1_EIG
|
|
Expectation of eigenvalue derivatives (hartree) for nkpt= 1 k points:
|
|
(in case of degenerate eigenvalues, averaged derivative)
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 0.00000 -0.00000 -0.00000
|
|
Expectation of eigenvalue derivatives ( eV ) for nkpt= 1 k points:
|
|
kpt# 1, nband= 4, wtk= 1.00000, kpt= 0.0000 0.0000 0.0000 (reduced coord)
|
|
-0.00000 0.00000 -0.00000 -0.00000
|
|
|
|
Nine components of 2nd-order total energy (hartree) are
|
|
1,2,3: 0th-order hamiltonian combined with 1st-order wavefunctions
|
|
kin0= 4.07363355E-01 eigvalue= 1.52974247E-01 local= -1.98584005E-01 nclr dpl0= 2.22455655E-09
|
|
4,5,6: 1st-order hamiltonian combined with 1st and 0th-order wfs
|
|
kin1= -8.91396193E-01 Hartree= 0.00000000E+00 xc= 0.00000000E+00 nclr dpl1= -1.18098012E-09
|
|
7,8,9: eventually, occupation + non-local contributions
|
|
edocc= 0.00000000E+00 enl0= -1.09482608E-02 enl1= 7.83299186E-02
|
|
10: eventually, PAW "on-site" Hxc contribution: epaw1= 0.00000000E+00
|
|
1-10 gives the relaxation energy (to be shifted if some occ is /=2.0)
|
|
erelax= -4.62260936E-01
|
|
11 Contribution from 1st-order change of wavefunctions overlap
|
|
eovl1 = -5.57277985E-02
|
|
No Ewald or frozen-wf contrib.: the relaxation energy is the total one
|
|
2DEtotal= -0.4622609358E+00 Ha. Also 2DEtotal= -0.125787597686E+02 eV
|
|
( non-var. 2DEtotal : -4.6226093588E-01 Ha)
|
|
|
|
====================================================
|
|
|
|
Orbital magnetic moment computed with DFPT derivative wavefunctions
|
|
|
|
Orbital magnetic moment, Cartesian directions :
|
|
-2.70058261E-04 -2.60006285E-17 2.93706745E-17
|
|
|
|
|
|
Chern vector, Cartesian directions :
|
|
-6.60484183E-08 1.08253354E-20 -2.23440100E-20
|
|
|
|
|
|
Orbital magnetic moment, term-by-term breakdown :
|
|
rho(1) CC : -6.07370136E-04 1.81789488E-17 -4.77561140E-17
|
|
rho(1) VV1 : -1.82528822E-04 1.82263893E-17 -4.11711228E-17
|
|
rho(1) VV2 : -3.95173368E-05 4.22889226E-18 -2.27989199E-17
|
|
rho(0) NL : 7.22869026E-04 -1.11501758E-17 2.72896672E-17
|
|
<L_R> terms : 2.69604004E-04 -5.79365334E-17 1.20267791E-16
|
|
<A0.An> terms : -9.12549967E-05 2.45185043E-18 -6.46062728E-18
|
|
Lamb terms : -3.41860000E-04 0.00000000E+00 0.00000000E+00
|
|
Chern vector, term-by-term breakdown :
|
|
Ch CC : -6.21436311E-08 2.07506015E-20 -4.42037389E-20
|
|
Ch vv1 : -3.90478725E-09 -9.83500781E-21 2.18798388E-20
|
|
Ch vv2 : 0.00000000E+00 -9.02582865E-23 -2.01098994E-23
|
|
|
|
|
|
Term-by-term breakdowns for each band :
|
|
|
|
|
|
band 1 of 4
|
|
Orbital magnetic moment : 6.50502060E-04 -2.44272359E-17 6.56601895E-17
|
|
rho(1) CC : 1.81220508E-05 -1.06825640E-18 2.93513043E-18
|
|
rho(1) VV1 : -1.75662774E-04 7.11701104E-18 -9.49849755E-18
|
|
rho(1) VV2 : 8.80459165E-05 -1.79826489E-17 4.17353689E-17
|
|
rho(0) NL : 7.28880988E-04 -1.24933416E-17 3.04881877E-17
|
|
<L_R> terms : 6.27931563E-26 -6.51568356E-27 1.18853588E-26
|
|
<A0.An> terms : -8.88412145E-06 2.15625455E-28 5.95940252E-30
|
|
|
|
|
|
Chern vector : 1.82314389E-08 4.43219430E-22 -1.63590160E-21
|
|
Ch CC : 7.03491306E-09 -9.45296412E-22 2.48363627E-21
|
|
Ch VV1 : 2.24478809E-08 -9.09480209E-22 1.21380949E-21
|
|
Ch VV2 : -1.12513551E-08 2.29799605E-21 -5.33334736E-21
|
|
|
|
|
|
band 2 of 4
|
|
Orbital magnetic moment : 8.06456375E-01 -8.21428425E-11 2.40013748E-10
|
|
rho(1) CC : 3.67746652E-01 -3.74551802E-11 1.09440315E-10
|
|
rho(1) VV1 : -6.86567224E-06 7.19161530E-16 -2.10050830E-15
|
|
rho(1) VV2 : -1.27563253E-04 1.29947793E-14 -3.79687789E-14
|
|
rho(0) NL : 3.25209590E-02 -3.31008336E-12 9.67267078E-12
|
|
<L_R> terms : 4.06356153E-01 -4.13921323E-11 1.20943284E-10
|
|
<A0.An> terms : -3.29599365E-05 8.39372449E-16 -2.45254176E-15
|
|
|
|
|
|
Chern vector : 6.73390190E-05 -6.85957324E-15 2.00429546E-14
|
|
Ch CC : 4.85509237E-05 -4.94288139E-15 1.44427791E-14
|
|
Ch VV1 : 9.59562116E-07 -1.00511667E-16 2.93571863E-16
|
|
Ch VV2 : 1.78285332E-05 -1.81618019E-15 5.30660371E-15
|
|
|
|
|
|
band 3 of 4
|
|
Orbital magnetic moment : -2.54628624E-05 1.64165279E-10 -4.79670739E-10
|
|
rho(1) CC : -8.98831084E-06 7.48969352E-11 -2.18839624E-10
|
|
rho(1) VV1 : -3.75873890E-10 -7.08052152E-16 2.06883568E-15
|
|
rho(1) VV2 : -1.99108918E-27 -1.29725677E-14 3.79042446E-14
|
|
rho(0) NL : 5.80263957E-24 6.61843138E-12 -1.93382411E-11
|
|
<L_R> terms : 1.76384438E-23 8.26635928E-11 -2.41532848E-10
|
|
<A0.An> terms : -1.64741757E-05 -9.42265728E-20 9.77180889E-19
|
|
|
|
|
|
Chern vector : -4.43461989E-10 1.37116138E-14 -4.00636452E-14
|
|
Ch CC : -5.48548344E-10 9.88680308E-15 -2.88880210E-14
|
|
Ch VV1 : 1.05086356E-10 1.97956342E-16 -5.78402510E-16
|
|
Ch VV2 : 5.56666242E-28 3.62685438E-15 -1.05972217E-14
|
|
|
|
|
|
band 4 of 4
|
|
Orbital magnetic moment : -8.07009613E-01 -8.20224378E-11 2.39656955E-10
|
|
rho(1) CC : -3.68363155E-01 -3.74417358E-11 1.09399258E-10
|
|
rho(1) VV1 : 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
rho(1) VV2 : 0.00000000E+00 0.00000000E+00 0.00000000E+00
|
|
rho(0) NL : -3.25269710E-02 -3.30834668E-12 9.66556713E-12
|
|
<L_R> terms : -4.06086549E-01 -4.12715185E-11 1.20589684E-10
|
|
<A0.An> terms : -3.29367631E-05 -8.36826372E-16 2.44510396E-15
|
|
|
|
|
|
Chern vector : -6.74228553E-05 -6.85203018E-15 2.00206698E-14
|
|
Ch CC : -4.86195537E-05 -4.94390000E-15 1.44451953E-14
|
|
Ch VV1 : -9.86019870E-07 -9.74536006E-17 2.84851313E-16
|
|
Ch VV2 : -1.78172818E-05 -1.81067658E-15 5.29062326E-15
|
|
|
|
====================================================
|
|
|
|
================================================================================
|
|
|
|
---- first-order wavefunction calculations are completed ----
|
|
|
|
|
|
respfn : d/dk was computed, but no 2DTE, so no DDB output.
|
|
|
|
== END DATASET(S) ==============================================================
|
|
================================================================================
|
|
|
|
-outvars: echo values of variables after computation --------
|
|
acell 2.0000000000E+01 2.0000000000E+01 2.0000000000E+01 Bohr
|
|
amu 2.01797000E+01
|
|
ecut 5.00000000E+00 Hartree
|
|
enunit 2
|
|
etotal1 -3.3612090135E+01
|
|
etotal2 -4.6226093584E-01
|
|
fcart1 9.9999999999E+99 9.9999999999E+99 9.9999999999E+99
|
|
- fftalg 512
|
|
getwfk1 0
|
|
getwfk2 -1
|
|
iscf1 17
|
|
iscf2 -3
|
|
istwfk 1
|
|
ixc -101130
|
|
jdtset 1 2
|
|
kptopt 0
|
|
lambsig1 0.00000000E+00
|
|
lambsig2 3.41860000E-04
|
|
P mkmem 1
|
|
P mkqmem 1
|
|
P mk1mem 1
|
|
natom 1
|
|
nband 4
|
|
ndtset 2
|
|
ngfft 45 45 45
|
|
ngfftdg 48 48 48
|
|
nkpt 1
|
|
nstep 500
|
|
nsym 1
|
|
ntypat 1
|
|
nucdipmom 1.0000000000E+00 0.0000000000E+00 0.0000000000E+00
|
|
occ 2.000000 2.000000 2.000000 2.000000
|
|
optdriver1 0
|
|
optdriver2 1
|
|
optforces 0
|
|
optstress 0
|
|
orbmag1 0
|
|
orbmag2 2
|
|
pawcpxocc 2
|
|
pawecutdg 6.00000000E+00 Hartree
|
|
prtpot1 0
|
|
prtpot2 1
|
|
rfddk1 0
|
|
rfddk2 1
|
|
spgroup 1
|
|
strten1 9.9999999999E+99 9.9999999999E+99 9.9999999999E+99
|
|
9.9999999999E+99 9.9999999999E+99 9.9999999999E+99
|
|
symmorphi 0
|
|
tolvrs1 1.00000000E-18
|
|
tolvrs2 0.00000000E+00
|
|
tolwfr1 0.00000000E+00
|
|
tolwfr2 1.00000000E-20
|
|
typat 1
|
|
usexcnhat 0
|
|
useylm 1
|
|
xangst 5.2917720859E+00 5.2917720859E+00 5.2917720859E+00
|
|
xcart 1.0000000000E+01 1.0000000000E+01 1.0000000000E+01
|
|
xred 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01
|
|
znucl 10.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] Orbital magnetism and chemical shielding in the projector augmented-wave formalism.
|
|
- J.W. Zwanziger, M. Torrent, and X. Gonze Phys. Rev. B 107, 165157 (2023).
|
|
- Comment: to be cited in case the computation of orbital magnetism is used, i.e. orbmag>0.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#zwanziger2023
|
|
-
|
|
- [2] Projector augmented-wave approach to density-functional perturbation theory.
|
|
- C. Audouze, F. Jollet, M. Torrent and X. Gonze, Phys. Rev. B 73, 235101 (2006).
|
|
- Comparison between projector augmented-wave and ultrasoft pseudopotential formalisms
|
|
- at the density-functional perturbation theory level.
|
|
- C. Audouze, F. Jollet, M. Torrent and X. Gonze, Phys. Rev. B 78, 035105 (2008).
|
|
- Comment: to be cited in case the computation of response function with PAW, i.e. (rfphon=1 or rfelfd=1) and usepaw=1.
|
|
- Strong suggestion to cite these papers.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#audouze2006,
|
|
- and https://docs.abinit.org/theory/bibliography/#audouze2008
|
|
-
|
|
- [3] Implementation of the Projector Augmented-Wave Method in the ABINIT code.
|
|
- M. Torrent, F. Jollet, F. Bottin, G. Zerah, and X. Gonze Comput. Mat. Science 42, 337, (2008).
|
|
- Comment: PAW calculations. Strong suggestion to cite this paper.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#torrent2008
|
|
-
|
|
- [4] Libxc: A library of exchange and correlation functionals for density functional theory.
|
|
- M.A.L. Marques, M.J.T. Oliveira, T. Burnus, Computer Physics Communications 183, 2227 (2012).
|
|
- Comment: to be cited when LibXC is used (negative value of ixc)
|
|
- Strong suggestion to cite this paper.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#marques2012
|
|
-
|
|
- [5] The Abinit project: Impact, environment and recent developments.
|
|
- Computer Phys. Comm. 248, 107042 (2020).
|
|
- X.Gonze, B. Amadon, G. Antonius, F.Arnardi, L.Baguet, J.-M.Beuken,
|
|
- J.Bieder, F.Bottin, J.Bouchet, E.Bousquet, N.Brouwer, F.Bruneval,
|
|
- G.Brunin, T.Cavignac, J.-B. Charraud, Wei Chen, M.Cote, S.Cottenier,
|
|
- J.Denier, G.Geneste, Ph.Ghosez, M.Giantomassi, Y.Gillet, O.Gingras,
|
|
- D.R.Hamann, G.Hautier, Xu He, N.Helbig, N.Holzwarth, Y.Jia, F.Jollet,
|
|
- W.Lafargue-Dit-Hauret, K.Lejaeghere, M.A.L.Marques, A.Martin, C.Martins,
|
|
- H.P.C. Miranda, F.Naccarato, K. Persson, G.Petretto, V.Planes, Y.Pouillon,
|
|
- S.Prokhorenko, F.Ricci, G.-M.Rignanese, A.H.Romero, M.M.Schmitt, M.Torrent,
|
|
- M.J.van Setten, B.Van Troeye, M.J.Verstraete, G.Zerah and J.W.Zwanzig
|
|
- Comment: the fifth generic paper describing the ABINIT project.
|
|
- Note that a version of this paper, that is not formatted for Computer Phys. Comm.
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT20.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze2020
|
|
-
|
|
- [6] ABINIT: Overview, and focus on selected capabilities
|
|
- J. Chem. Phys. 152, 124102 (2020).
|
|
- A. Romero, D.C. Allan, B. Amadon, G. Antonius, T. Applencourt, L.Baguet,
|
|
- J.Bieder, F.Bottin, J.Bouchet, E.Bousquet, F.Bruneval,
|
|
- G.Brunin, D.Caliste, M.Cote,
|
|
- J.Denier, C. Dreyer, Ph.Ghosez, M.Giantomassi, Y.Gillet, O.Gingras,
|
|
- D.R.Hamann, G.Hautier, F.Jollet, G. Jomard,
|
|
- A.Martin,
|
|
- H.P.C. Miranda, F.Naccarato, G.Petretto, N.A. Pike, V.Planes,
|
|
- S.Prokhorenko, T. Rangel, F.Ricci, G.-M.Rignanese, M.Royo, M.Stengel, M.Torrent,
|
|
- M.J.van Setten, B.Van Troeye, M.J.Verstraete, J.Wiktor, J.W.Zwanziger, and X.Gonze.
|
|
- Comment: a global overview of ABINIT, with focus on selected capabilities .
|
|
- Note that a version of this paper, that is not formatted for J. Chem. Phys
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT20_JPC.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#romero2020
|
|
-
|
|
- [7] Recent developments in the ABINIT software package.
|
|
- Computer Phys. Comm. 205, 106 (2016).
|
|
- X.Gonze, F.Jollet, F.Abreu Araujo, D.Adams, B.Amadon, T.Applencourt,
|
|
- C.Audouze, J.-M.Beuken, J.Bieder, A.Bokhanchuk, E.Bousquet, F.Bruneval
|
|
- D.Caliste, M.Cote, F.Dahm, F.Da Pieve, M.Delaveau, M.Di Gennaro,
|
|
- B.Dorado, C.Espejo, G.Geneste, L.Genovese, A.Gerossier, M.Giantomassi,
|
|
- Y.Gillet, D.R.Hamann, L.He, G.Jomard, J.Laflamme Janssen, S.Le Roux,
|
|
- A.Levitt, A.Lherbier, F.Liu, I.Lukacevic, A.Martin, C.Martins,
|
|
- M.J.T.Oliveira, S.Ponce, Y.Pouillon, T.Rangel, G.-M.Rignanese,
|
|
- A.H.Romero, B.Rousseau, O.Rubel, A.A.Shukri, M.Stankovski, M.Torrent,
|
|
- M.J.Van Setten, B.Van Troeye, M.J.Verstraete, D.Waroquier, J.Wiktor,
|
|
- B.Xu, A.Zhou, J.W.Zwanziger.
|
|
- Comment: the fourth generic paper describing the ABINIT project.
|
|
- Note that a version of this paper, that is not formatted for Computer Phys. Comm.
|
|
- is available at https://www.abinit.org/sites/default/files/ABINIT16.pdf .
|
|
- The licence allows the authors to put it on the Web.
|
|
- DOI and bibtex: see https://docs.abinit.org/theory/bibliography/#gonze2016
|
|
-
|
|
- Proc. 0 individual time (sec): cpu= 9.8 wall= 9.8
|
|
|
|
================================================================================
|
|
|
|
Calculation completed.
|
|
.Delivered 57 WARNINGs and 9 COMMENTs to log file.
|
|
+Overall time at end (sec) : cpu= 9.8 wall= 9.8
|