Version 4.4.3 of ABINIT (sequential version, prepared for a P6/Linux computer) Copyright (C) 1998-2004 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/Infos/copyright 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/Infos/contributors. Please read ~ABINIT/Infos/acknowledgments.htm for suggested acknowledgments of the ABINIT effort. For more information, see http://www.abinit.org . Starting date : Sat 21 May 2005. - input file -> gs.097 - output file -> gs_097.out - root for input files -> fexi - root for output files -> fexo Symmetries : space group Im -3 m (#229); Bravais cI (body-center cubic) ================================================================================ Values of the parameters that define the memory need of the present run intxc = 0 ionmov = 0 iscf = 5 ixc = 11 lmnmax = 4 lnmax = 4 mband = 12 mffmem = 1 P mgfft = 24 mkmem = 6 mpssoang= 4 mpw = 471 mqgrid = 1201 natom = 1 nfft = 13824 nkpt = 6 nloalg = 4 nspden = 1 nspinor = 1 nsppol = 1 nsym = 48 n1xccc = 0 ntypat = 1 occopt = 3 ================================================================================ P This job should need less than 5.228 Mbytes of memory. Rough estimation (10% accuracy) of disk space for files : WF disk file : 0.519 Mbytes ; DEN or POT disk file : 0.107 Mbytes. ================================================================================ -outvars: echo values of preprocessed input variables -------- acell 5.3687119211E+00 5.3687119211E+00 5.3687119211E+00 Bohr amu 5.58470000E+01 ecut 2.50000000E+01 Hartree ixc 11 kpt 1.25000000E-01 1.25000000E-01 1.25000000E-01 3.75000000E-01 1.25000000E-01 1.25000000E-01 -1.25000000E-01 1.25000000E-01 1.25000000E-01 -3.75000000E-01 3.75000000E-01 1.25000000E-01 -3.75000000E-01 -3.75000000E-01 1.25000000E-01 -3.75000000E-01 3.75000000E-01 3.75000000E-01 kptrlen 1.85977636E+01 kptopt 1 kptrlatt 4 0 0 0 4 0 0 0 4 P mkmem 6 natom 1 nband 12 ngfft 24 24 24 nkpt 6 nstep 35 nsym 48 ntypat 1 occ 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 occopt 3 rprim -5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 -5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 -5.0000000000E-01 shiftk 5.00000000E-01 5.00000000E-01 5.00000000E-01 spgroup 229 symrel 1 0 0 0 1 0 0 0 1 -1 0 0 0 -1 0 0 0 -1 -1 -1 -1 0 0 1 0 1 0 1 1 1 0 0 -1 0 -1 0 0 1 0 1 0 0 -1 -1 -1 0 -1 0 -1 0 0 1 1 1 0 0 1 -1 -1 -1 1 0 0 0 0 -1 1 1 1 -1 0 0 0 1 0 1 0 0 0 0 1 0 -1 0 -1 0 0 0 0 -1 -1 -1 -1 0 0 1 1 0 0 1 1 1 0 0 -1 -1 0 0 1 0 0 0 1 0 -1 -1 -1 -1 0 0 0 -1 0 1 1 1 0 0 1 -1 -1 -1 0 1 0 0 0 -1 1 1 1 0 -1 0 -1 -1 -1 1 0 0 0 0 1 1 1 1 -1 0 0 0 0 -1 0 1 0 0 0 1 1 0 0 0 -1 0 0 0 -1 -1 0 0 1 0 0 -1 -1 -1 0 1 0 -1 0 0 1 1 1 0 -1 0 0 0 1 0 1 0 -1 -1 -1 0 0 -1 0 -1 0 1 1 1 1 0 0 -1 -1 -1 0 0 1 -1 0 0 1 1 1 0 0 -1 0 1 0 0 0 1 -1 -1 -1 0 -1 0 0 0 -1 1 1 1 -1 -1 -1 1 0 0 0 1 0 1 1 1 -1 0 0 0 -1 0 0 0 1 0 1 0 1 0 0 0 0 -1 0 -1 0 -1 0 0 0 1 0 -1 -1 -1 0 0 1 0 -1 0 1 1 1 0 0 -1 1 0 0 0 0 1 -1 -1 -1 -1 0 0 0 0 -1 1 1 1 -1 -1 -1 0 1 0 1 0 0 1 1 1 0 -1 0 -1 0 0 0 0 1 1 0 0 0 1 0 0 0 -1 -1 0 0 0 -1 0 -1 -1 -1 0 1 0 0 0 1 1 1 1 0 -1 0 0 0 -1 1 0 0 0 0 1 0 1 0 -1 0 0 0 0 -1 0 -1 0 0 1 0 -1 -1 -1 1 0 0 0 -1 0 1 1 1 -1 0 0 0 0 1 1 0 0 -1 -1 -1 0 0 -1 -1 0 0 1 1 1 toldfe 1.00000000E-06 Hartree tsmear 1.83746623E-03 Hartree typat 1 wtk 0.12500 0.18750 0.09375 0.37500 0.12500 0.09375 znucl 26.00000 ================================================================================ chkinp: Checking input parameters for consistency. ================================================================================ == DATASET 1 ================================================================== Real(R)+Recip(G) space primitive vectors, cartesian coordinates (Bohr,Bohr^-1): R(1)= -2.6843560 2.6843560 2.6843560 G(1)= 0.0000000 0.1862644 0.1862644 R(2)= 2.6843560 -2.6843560 2.6843560 G(2)= 0.1862644 0.0000000 0.1862644 R(3)= 2.6843560 2.6843560 -2.6843560 G(3)= 0.1862644 0.1862644 0.0000000 Unit cell volume ucvol= 7.7371374E+01 bohr^3 Angles (23,13,12)= 1.09471221E+02 1.09471221E+02 1.09471221E+02 degrees getcut: wavevector= 0.0000 0.0000 0.0000 ngfft= 24 24 24 ecut(hartree)= 25.000 => boxcut(ratio)= 2.29338 getcut : COMMENT - Note that boxcut > 2.2 ; recall that boxcut=Gcut(box)/Gcut(sphere) = 2 is sufficient for exact treatment of convolution. Such a large boxcut is a waste : you could raise ecut e.g. ecut= 32.872352 Hartrees makes boxcut=2 --- Pseudopotential description ------------------------------------------------ - pspini: atom type 1 psp file is fe.fhi - pspatm: opening atomic psp file fe.fhi OPIUM generated Fe potential 26.00000 8.00000 50520 znucl, zion, pspdat 6 11 3 3 1129 0.00000 pspcod,pspxc,lmax,lloc,mmax,r2well 0.00000000000000 0.00000000000000 0.00000000000000 rchrg,fchrg,qchrg 1.013085 amesh (Hamman grid) pspatm: epsatm= 57.34358424 --- l ekb(1:nproj) --> 0 0.275000 1 191.755585 2 -7.766120 pspatm: atomic psp has been read and splines computed 4.58748674E+02 ecore*ucvol(ha*bohr**3) -------------------------------------------------------------------------------- P newkpt: treating 12 bands with npw= 471 for ikpt= 1 P newkpt: treating 12 bands with npw= 460 for ikpt= 2 P newkpt: treating 12 bands with npw= 454 for ikpt= 3 P newkpt: treating 12 bands with npw= 464 for ikpt= 4 P newkpt: treating 12 bands with npw= 454 for ikpt= 5 P newkpt: treating 12 bands with npw= 454 for ikpt= 6 setup2: Arith. and geom. avg. npw (full set) are 461.000 460.966 ================================================================================ iter Etot(hartree) deltaE(h) residm vres2 diffor maxfor ETOT 1 -18.427053113671 -1.843E+01 1.611E+01 1.905E+03 0.000E+00 0.000E+00 ETOT 2 -19.137292809749 -7.102E-01 1.436E+00 5.905E+02 0.000E+00 0.000E+00 ETOT 3 -19.423571896407 -2.863E-01 7.591E-01 3.971E+01 0.000E+00 0.000E+00 ETOT 4 -19.427712152486 -4.140E-03 4.692E-02 1.355E+00 0.000E+00 0.000E+00 ETOT 5 -19.428165557056 -4.534E-04 8.931E-01 1.584E+00 0.000E+00 0.000E+00 ETOT 6 -19.428136659070 2.890E-05 4.747E-01 4.564E-01 0.000E+00 0.000E+00 ETOT 7 -19.428369422775 -2.328E-04 2.871E-01 1.794E-02 0.000E+00 0.000E+00 ETOT 8 -19.428372559689 -3.137E-06 1.342E-01 4.256E-04 0.000E+00 0.000E+00 ETOT 9 -19.428372985115 -4.254E-07 3.915E-02 4.363E-05 0.000E+00 0.000E+00 ETOT 10 -19.428372994661 -9.546E-09 1.119E-02 2.758E-04 0.000E+00 0.000E+00 At SCF step 10, etot is converged : for the second time, diff in etot= 9.546E-09 < toldfe= 1.000E-06 Cartesian components of stress tensor (hartree/bohr^3) sigma(1 1)= 6.17750752E-05 sigma(3 2)= 0.00000000E+00 sigma(2 2)= 6.17750752E-05 sigma(3 1)= 0.00000000E+00 sigma(3 3)= 6.17750752E-05 sigma(2 1)= 0.00000000E+00 ================================================================================ ----iterations are completed or convergence reached---- Mean square residual over all n,k,spin= 2.0118E-04; max= 1.1191E-02 0.1250 0.1250 0.1250 1 2.02858E-03 kpt; spin; max resid(k); each band: 1.99E-09 4.78E-10 4.47E-10 3.66E-10 3.43E-10 3.96E-10 8.75E-10 7.30E-10 1.95E-10 5.82E-09 8.02E-05 2.03E-03 0.3750 0.1250 0.1250 1 1.52747E-05 kpt; spin; max resid(k); each band: 5.56E-10 6.43E-10 6.07E-10 5.21E-10 3.20E-10 6.48E-10 8.81E-10 6.03E-10 1.44E-10 6.80E-09 7.82E-06 1.53E-05 -0.1250 0.1250 0.1250 1 1.00683E-08 kpt; spin; max resid(k); each band: 4.74E-10 2.58E-10 2.65E-10 3.31E-10 1.74E-10 1.15E-10 1.19E-10 8.30E-10 1.13E-09 1.47E-09 1.40E-09 1.01E-08 -0.3750 0.3750 0.1250 1 1.11915E-02 kpt; spin; max resid(k); each band: 4.84E-10 4.48E-10 2.44E-10 6.74E-10 1.09E-10 2.27E-10 6.58E-10 3.16E-10 1.28E-09 4.04E-07 2.84E-05 1.12E-02 -0.3750 -0.3750 0.1250 1 1.95529E-09 kpt; spin; max resid(k); each band: 2.56E-10 2.53E-10 5.61E-10 5.61E-10 6.45E-10 4.05E-10 4.15E-10 5.10E-10 3.00E-10 2.43E-10 2.27E-10 1.96E-09 -0.3750 0.3750 0.3750 1 1.13293E-03 kpt; spin; max resid(k); each band: 2.80E-11 1.28E-11 2.64E-10 2.47E-10 1.24E-10 1.06E-10 2.62E-10 1.92E-10 4.12E-11 8.00E-10 6.50E-10 1.13E-03 reduced coordinates (array xred) for 1 atoms 0.000000000000 0.000000000000 0.000000000000 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 0.00000000000000 0.00000000000000 0.00000000000000 cartesian forces (hartree/bohr) at end: 1 0.00000000000000 0.00000000000000 0.00000000000000 frms,max,avg= 0.0000000E+00 0.0000000E+00 0.000E+00 0.000E+00 0.000E+00 h/b cartesian forces (eV/Angstrom) at end: 1 0.00000000000000 0.00000000000000 0.00000000000000 frms,max,avg= 0.0000000E+00 0.0000000E+00 0.000E+00 0.000E+00 0.000E+00 e/A length scales= 5.368711921107 5.368711921107 5.368711921107 bohr = 2.841000000000 2.841000000000 2.841000000000 angstroms Fermi (or HOMO) energy (hartree) = -0.20779 Average Vxc (hartree)= -0.46223 Eigenvalues (hartree) for nkpt= 6 k points: kpt# 1, nband= 12, wtk= 0.12500, kpt= 0.1250 0.1250 0.1250 (reduced coord) -0.40134 -0.29950 -0.29950 -0.21218 -0.21218 -0.18265 0.42801 0.42801 0.47972 0.87269 1.11023 1.11348 occupation numbers for kpt# 1 2.00000 2.00000 2.00000 1.83207 1.83207 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 kpt# 2, nband= 12, wtk= 0.18750, kpt= 0.3750 0.1250 0.1250 (reduced coord) -0.34349 -0.31557 -0.27279 -0.20804 -0.20309 -0.04512 0.10276 0.33780 0.35551 0.95184 1.11707 1.16269 occupation numbers for kpt# 2 2.00000 2.00000 2.00000 1.06723 0.14389 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 kpt# 3, nband= 12, wtk= 0.09375, kpt= -0.1250 0.1250 0.1250 (reduced coord) -0.47797 -0.26962 -0.26962 -0.25716 -0.23457 -0.20334 0.62082 0.69253 0.69253 0.70170 0.89495 0.89495 occupation numbers for kpt# 3 2.00000 2.00000 2.00000 2.00000 2.00000 0.16271 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 kpt# 4, nband= 12, wtk= 0.37500, kpt= -0.3750 0.3750 0.1250 (reduced coord) -0.36663 -0.29172 -0.28324 -0.22734 -0.18668 -0.16629 0.21858 0.48896 0.50338 0.81730 0.90106 0.96963 occupation numbers for kpt# 4 2.00000 2.00000 2.00000 1.99995 0.00002 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 kpt# 5, nband= 12, wtk= 0.12500, kpt= -0.3750 -0.3750 0.1250 (reduced coord) -0.31690 -0.31690 -0.29992 -0.20191 -0.20191 0.07180 0.07180 0.16580 0.69734 0.74008 0.74008 1.04256 occupation numbers for kpt# 5 2.00000 2.00000 2.00000 0.07829 0.07829 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 kpt# 6, nband= 12, wtk= 0.09375, kpt= -0.3750 0.3750 0.3750 (reduced coord) -0.34706 -0.33910 -0.21212 -0.21212 -0.19121 -0.01801 0.24858 0.24858 0.47393 0.50128 0.76058 1.31979 occupation numbers for kpt# 6 2.00000 2.00000 1.82698 1.82698 0.00024 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ,Min el dens= 1.9714E-04 el/bohr^3 at reduced coord. 0.0000 0.0000 0.0000 , next min= 1.3170E-02 el/bohr^3 at reduced coord. 0.0417 0.0000 0.0000 ,Max el dens= 6.2632E-01 el/bohr^3 at reduced coord. 0.7917 0.9583 0.9583 , next max= 6.2632E-01 el/bohr^3 at reduced coord. 0.9583 0.7917 0.9583 -------------------------------------------------------------------------------- Components of total free energy (in Hartree) : Kinetic energy = 2.45574332529892E+01 Hartree energy = 2.52725228091096E+00 XC energy =-4.06838665073547E+00 Ewald energy =-2.16915103837943E+01 PspCore energy = 5.92917836166837E+00 Loc. psp. energy= 6.73509062087529E-01 NL psp energy=-2.73544769215961E+01 >>>>> Internal E=-1.94270009984699E+01 -kT*entropy =-1.37199619132389E-03 >>>>>>>>> Etotal=-1.94283729946612E+01 Other information on the energy : Total energy(eV)=-5.28672927767142E+02 ; Band energy (Ha)= -2.3452413636E+00 -------------------------------------------------------------------------------- Cartesian components of stress tensor (hartree/bohr^3) sigma(1 1)= 6.17750752E-05 sigma(3 2)= 0.00000000E+00 sigma(2 2)= 6.17750752E-05 sigma(3 1)= 0.00000000E+00 sigma(3 3)= 6.17750752E-05 sigma(2 1)= 0.00000000E+00 -Cartesian components of stress tensor (GPa) [Pressure= -1.8175E+00 GPa] - sigma(1 1)= 1.81748522E+00 sigma(3 2)= 0.00000000E+00 - sigma(2 2)= 1.81748522E+00 sigma(3 1)= 0.00000000E+00 - sigma(3 3)= 1.81748522E+00 sigma(2 1)= 0.00000000E+00 == END DATASET(S) ============================================================== ================================================================================ -outvars: echo values of variables after computation -------- acell 5.3687119211E+00 5.3687119211E+00 5.3687119211E+00 Bohr amu 5.58470000E+01 ecut 2.50000000E+01 Hartree etotal -1.9428372995E+01 fcart 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 ixc 11 kpt 1.25000000E-01 1.25000000E-01 1.25000000E-01 3.75000000E-01 1.25000000E-01 1.25000000E-01 -1.25000000E-01 1.25000000E-01 1.25000000E-01 -3.75000000E-01 3.75000000E-01 1.25000000E-01 -3.75000000E-01 -3.75000000E-01 1.25000000E-01 -3.75000000E-01 3.75000000E-01 3.75000000E-01 kptrlen 1.85977636E+01 kptopt 1 kptrlatt 4 0 0 0 4 0 0 0 4 P mkmem 6 natom 1 nband 12 ngfft 24 24 24 nkpt 6 nstep 35 nsym 48 ntypat 1 occ 2.000000 2.000000 2.000000 1.832070 1.832070 0.000002 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 1.067226 0.143888 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 2.000000 1.999999 0.162708 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 1.999952 0.000021 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 2.000000 0.078294 0.078294 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 2.000000 2.000000 1.826979 1.826979 0.000242 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 occopt 3 rprim -5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 -5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 5.0000000000E-01 -5.0000000000E-01 shiftk 5.00000000E-01 5.00000000E-01 5.00000000E-01 spgroup 229 strten 6.1775075152E-05 6.1775075152E-05 6.1775075152E-05 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 -1 -1 0 0 1 0 1 0 1 1 1 0 0 -1 0 -1 0 0 1 0 1 0 0 -1 -1 -1 0 -1 0 -1 0 0 1 1 1 0 0 1 -1 -1 -1 1 0 0 0 0 -1 1 1 1 -1 0 0 0 1 0 1 0 0 0 0 1 0 -1 0 -1 0 0 0 0 -1 -1 -1 -1 0 0 1 1 0 0 1 1 1 0 0 -1 -1 0 0 1 0 0 0 1 0 -1 -1 -1 -1 0 0 0 -1 0 1 1 1 0 0 1 -1 -1 -1 0 1 0 0 0 -1 1 1 1 0 -1 0 -1 -1 -1 1 0 0 0 0 1 1 1 1 -1 0 0 0 0 -1 0 1 0 0 0 1 1 0 0 0 -1 0 0 0 -1 -1 0 0 1 0 0 -1 -1 -1 0 1 0 -1 0 0 1 1 1 0 -1 0 0 0 1 0 1 0 -1 -1 -1 0 0 -1 0 -1 0 1 1 1 1 0 0 -1 -1 -1 0 0 1 -1 0 0 1 1 1 0 0 -1 0 1 0 0 0 1 -1 -1 -1 0 -1 0 0 0 -1 1 1 1 -1 -1 -1 1 0 0 0 1 0 1 1 1 -1 0 0 0 -1 0 0 0 1 0 1 0 1 0 0 0 0 -1 0 -1 0 -1 0 0 0 1 0 -1 -1 -1 0 0 1 0 -1 0 1 1 1 0 0 -1 1 0 0 0 0 1 -1 -1 -1 -1 0 0 0 0 -1 1 1 1 -1 -1 -1 0 1 0 1 0 0 1 1 1 0 -1 0 -1 0 0 0 0 1 1 0 0 0 1 0 0 0 -1 -1 0 0 0 -1 0 -1 -1 -1 0 1 0 0 0 1 1 1 1 0 -1 0 0 0 -1 1 0 0 0 0 1 0 1 0 -1 0 0 0 0 -1 0 -1 0 0 1 0 -1 -1 -1 1 0 0 0 -1 0 1 1 1 -1 0 0 0 0 1 1 0 0 -1 -1 -1 0 0 -1 -1 0 0 1 1 1 toldfe 1.00000000E-06 Hartree tsmear 1.83746623E-03 Hartree typat 1 wtk 0.12500 0.18750 0.09375 0.37500 0.12500 0.09375 znucl 26.00000 ================================================================================ - Total cpu time (s,m,h): 15.6 0.26 0.004 - Total wall clock time (s,m,h): 15.6 0.26 0.004 - For major independent code sections, cpu and wall times (sec), - as well as % of the total time and number of calls - routine cpu % wall % number of calls - (-1=no count) - fourwf(pot) 9.789 62.8 9.831 63.1 4314 - nonlop(apply) 1.268 8.1 1.217 7.8 4314 - projbd 0.830 5.3 0.837 5.4 6900 - pspini 0.525 3.4 0.526 3.4 1 - cgwf-O(npw) 0.480 3.1 0.473 3.0 -1 - timing timab 0.420 2.7 0.420 2.7 10 - fourwf(den) 0.412 2.6 0.404 2.6 326 - fourdp 0.268 1.7 0.283 1.8 139 - vtowfk(ssdiag) 0.254 1.6 0.248 1.6 -1 - getghc-other 0.240 1.5 0.244 1.6 -1 - xc:pot/=fourdp 0.223 1.4 0.216 1.4 11 - nonlop(forces) 0.184 1.2 0.189 1.2 720 - invars2 0.117 0.8 0.117 0.8 1 - forces 0.109 0.7 0.108 0.7 10 - vtorho(4)-mkrho 0.094 0.6 0.097 0.6 20 - 42 others 0.291 1.9 0.289 1.9 - subtotal 15.504 99.5 15.499 99.4 ================================================================================ Calculation completed. Delivered 0 WARNINGs and 2 COMMENTs to log file. +Overall time at end (sec) : cpu= 15.6 wall= 15.6