------------------------------------------------------------------------
|                                                                      |
|    FULL-POTENTIAL LOCAL-ORBITAL MINIMUM BASIS BANDSTRUCTURE CODE     |
|                      FULL RELATIVISTIC VERSION                       |
|                                                                      |
|             by K. Koepernik, A.Ernst and H.Eschrig (2003)            |
|                  relativistic version by Ingo Opahle                 |
|                       LSDA+U by Igor Chaplygin                       |
|                  dHvA module based on contributions by               |
|                             Vivien Petzold                           |
|                                                                      |
------------------------------------------------------------------------
|                                                                      |
|   Any publication of results obtained by this program                |
|   has to include the citation:                                       |
|                                                                      |
|        K.Koepernik and H.Eschrig, Phys. Rev. B 59, 1743 (1999)       |
|                                                                      |
|   Any publication of CPA results obtained by this program            |
|   additionally has to include the citation:                          |
|                                                                      |
|        K. Koepernik, B. Velicky, R. Hayn and H. Eschrig,             |
|        Phys. Rev. B 55, 5717 (1997)                                  |
|                                                                      |
|                                                                      |
------------------------------------------------------------------------
|                                                                      |
| main version: 21.00                                                  |
| sub  version: M-CPA                                                  |
| release     : 61                                                     |
|                                                                      |
------------------------------------------------------------------------
| compiled with ifort 1310                                             |
| compiled with own LAPACK copy                                        |
------------------------------------------------------------------------
|                                                                      |
------------------------------------------------------------------------
| date        : Fri Jun 24 11:42:29 2022                               |
| host        : unas                                                   |
------------------------------------------------------------------------

File =.in exists!
------------------------------------------------------------------------
Compound: Not specified
------------------------------------------------------------------------
------------------------------------------------------------------------
Start: content of =.in
------------------------------------------------------------------------
########################################################################
# this file is part of the FPLO bandstructure package
# it was created automatically
# please dont delete this file nor
# edit it directly, unless you really know what you are doing
# Use FEDIT to change its content
#
########################################################################
section header{

    struct {char[*] mainversion;char[*] subversion;} version
    ={"21.00","M-CPA"};


    struct {char[*] hostname;char[*] date;char[*] by;} last_changes
    ={"unas","Fri Jun 24 11:42:16 2022","fedit (21.00-61-x86_64)"};

    char[*] compound="Not specified";
};


section structure_information{

    struct {int type;char[*] description;} structure_type
    ={1,"Crystal   "};


    struct {char[*] number;char[*] setting;char[*] symbol;} spacegroup
    ={"221","default","P m -3 m"};

    int subgroupgenerators[*]={
    };

    struct {int type;char[*] description;} lengthunit
    ={1,"Bohr radii"};

    real lattice_constants[3]={5.4,5.4,5.4
    };
    real axis_angles[3]={90,90,90
    };

    struct {logical active;real newx[3];real newz[3];} cellrotation
    ={f,{1,0,0},{0,0,1}};

    int max_L=4;
    int nsort=2;

    struct {char[2] element;real tau[3];} wyckoff_positions[nsort]
    ={
        {"Fe",{0,0,0}}
        ,{"Fe",{1/2,1/2,1/2}}
    };

};


section structure_dependend{

    struct {char[2] atom;real concentration;} concentrations[nsort]
    ={
        {"Fe",1}
        ,{"Fe",1}
    };


    struct {char[2] element;real split;} initial_spin_split[nsort]
    ={
        {"Fe",2}
        ,{"Fe",-2}
    };


    struct {int sort;char[2] state;real spinup[7];real spindn[7];} core_occupation[*]
    ={
    };


    struct {char[20] label;real kpoint[3];} special_sympoints[*]
    ={
        {"$~G",{0,0,0}}
        ,{"X",{0,1/2,0}}
        ,{"M",{1/2,1/2,0}}
        ,{"$~G",{0,0,0}}
        ,{"R",{1/2,1/2,1/2}}
        ,{"X",{0,1/2,0}}
        ,{"<BREAK>",{0,0,0}}
        ,{"M",{1/2,1/2,0}}
        ,{"R",{1/2,1/2,1/2}}
    };

};


section mesh{

    struct {int nr;real rmin;real rmax;real rmax_factor;} radial_mesh
    ={200,1.0e-5,0.0,1.0};


    struct {int nr;real rmin;real rmax;} atom_radial_mesh
    ={2000,1.0e-7,0.0};

};


section brillouin{

    struct {logical metal;int nkxyz[3];} bzone_integration
    ={t,{12,12,12}};


    struct {int ne;real range;real width;int norder;} mp_bzone_integration
    ={1001,5.0,0.005,0};


    struct {int type;char[*] description;} bzone_method
    ={1,"tetrahedron method"};


    struct {int type;char[*] description;} bzintmethod
    ={1,"primitive cell"};

    real bzone_kT=100.0;
    int nband=-1;
};


section bandstructure_plot{

    struct {logical bandplot;logical read_sympoints;int ndivisions;real emin;real emax;int nptdos;logical plot_idos;logical plot_ndos;logical restrict_bands_to_window;logical coeffout;} bandplot_control
    ={f,t,50,-20.0,20.0,1000,f,f,f,f};


    struct {logical bandweights;logical bweight_rot;real z_axis[3];real x_axis[3];char[*] bwdfilename;
        struct {int type;char[*] description;} frelprojection;
    } bandweight_control
    ={f,f,{0.0,0.0,1.0},{1.0,0.0,0.0},""
        ,{1,"jmu"}
    };

};


section iteration_control{

    struct {int steps;real tolerance;real mixing;int maxdim;real progress;} iteration_control
    ={30,1.e-6,0.4,5,0.1};


    struct {int scheme;char[*] description;} iteration_version
    ={2,"Iterat   : Hyperspace interpolation"};


    struct {int type;char[*] description;} iteration_convergence_condition
    ={2,"Density"};

    real etot_tolerance=1.e-8;
    # ratio of things put into iteration vector compared
    # to density

    struct {real U;} iterationvector_ratio
    ={1.0};

};


section forces{

    struct {int steps;real tolerance;int maxdim;} force_iteration_control
    ={100,1.0e-3,0};


    struct {int scheme;char[*] description;} force_iteration_version
    ={6,"MinIterat: Find minimum with iterat"};

    logical force_in_each_step=f;
    logical approx_nodndL_ewald=f;

    struct {int mode;char[*] description;} force_mode
    ={1,"No forces"};

    int geomopt_sites[*]={
    };
};


section options{

    struct {int mspin;logical initial_spin_split;} spin
    ={2,f};


    struct {int type;char[*] description;} xc_potential
    ={9,"Perdew Wang 92               (LSDA)"};


    struct {logical overwrite_cmBJ;real cmBJ;} mBJcontrols
    ={f,1};

    real ex_field_strength=1.0;

    struct {int type;char[*] description;} relativistic
    ={2,"scalar relativistic"};


    struct {int type;char[*] description;} nuclear_model
    ={1,"Point charge"};


    struct {
        struct {int type;char[*] description;} chargemode;
        real ionicity;logical make_vca;
        struct {int sort;real Z;} vca[*];
    } charges
    ={
         {1,"None"}
        ,0.0,f
        ,{
        }
    };


    struct {int type;char[*] description;} spin_pop_type
    ={3,"exact expression (expensive)"};

    real quantization_axis[3]={0,0,1
    };
    flag options[*]={CALC_DOS(-),NOT_USED(-),FULLBZ(-),CALC_PLASMON_FREQ(-)
    ,EMPTY_LATTICE_TEST(-),NO_DOS(-),PLOT_REALFUNC(-),PLOT_BASIS(-)
    ,TEST_LOI(-),TEST_DIAGNO(-),TEST_SYMMETRIZATION(-),TEST_HS_SYM(-)
    ,PROT_PRINT_BASIS(-),PROT_MAKELATTICE(-),PROT_STRUCTURE_PRNT(-),PROT_TCI(-)
    ,NOT_USED(-),NOT_USED(-),NOT_USED(-),NO_SYMMETRYTEST(+)
    ,NO_POTENTIAL(-),NOT_USED(-),NO_CORE(-),NOT_USED(-)
    ,NO_POPANALYSIS(-),NO_LOI(-),NO_BASIS(-),NO_EFCH(-)
    ,NOT_USED(-),NOT_USED(-),NOT_USED(-)
    };

    struct {int type;char[*] description;} verbosity_level
    ={200,"more information"};


    struct {real moment;logical fsm;} fsmoment
    ={1.0,f};

    logical calc_etot=t;
    logical usedatadirectories=f;
};


section cpa{

    struct {real emin0;real emax0;int nec;int max_se_loops;real eimag;} cpa_control
    ={-1.60,0.4,30,10,1.0e-5};


    struct {real efermi0;real slop_0;real slop_min;int max_ef_loops;} cpa_fermi_control
    ={0.08,30.0,0.3,10};


    struct {int type;char[*] description;} cpa_fermi_method
    ={1,"secant method"};


    struct {int ne;real eimag;} cpa_dos
    ={1000,1.e-3};


    struct {int ne;real eimag;} cpa_bloch_spectral_density
    ={100,1.e-3};

};


section numerics{

    struct {int degree;} interpolation
    ={8};


    struct {int nr;} one_center_integrals
    ={200};


    struct {int nxi;int neta;} coulomb_energy_calc
    ={26,46};


    struct {int nr;int neta;} overlap_density_calc
    ={50,20};


    struct {real dens;real pot;real orbital;} cutoff_tolerances
    ={1.e-9,1.e-12,1.e-8};

    #WARNING: cutoff is obsolete now
    #WARNING: cutoff is obsolete now

    struct {real cutoff;real eps_r;real eps_g;} ewald_parameter
    ={2.0,1.0e-5,1.0e-5};


    struct {int nr;int angular_mesh_type_min;int angular_mesh_type_max;logical use_symmetry;} three_center_integrals
    ={80,3,15,t};


    struct {real tol;} symtest
    ={5.0e-13};

    # WARNING: nfft_ewp_per_site is obsolete now
    # WARNING: nfft_ewp_per_site is obsolete now
    int nfft_ewp_per_site=-500;
};


section LSDA_U{
    logical make_lsdau=f;

    struct {int type;char[*] description;} lsdau_projection
    ={1,"orthogonal"};


    struct {int type;char[*] description;} lsdau_functional
    ={1,"LSDA+U Around mean field (AMF/OP)"};


    struct {int type;char[*] description;} lsdau_dmat_mode
    ={1,"self consistent"};

    real lsdau_dmat_deviation=1.e-2;

    struct {int sort;char[2] state;real f_slater[4];} lsdau_defs[*]
    ={
    };

};


section OPC{
    logical make_opc=f;

    struct {int type;char[*] description;} opc_projection
    ={1,"orthogonal"};


    struct {int type;char[*] description;} opc_functional
    ={1,"Spin dependend (Nordstroem, Brooks, Johansson, JPhysC 4, 3261 (1992))"};


    struct {int sort;char[2] state;} opc_defs[*]
    ={
    };

};


section Advanced_output{
    int ldossites[*]={
    };
    logical grids_active=f;
    logical grids_stop_after=f;

    struct {
        struct {int type;char[*] description;} basis;
        real direction1[3];real direction2[3];real direction3[3];real origin[3];int subdivision[3];char[*] file;logical InclPeriodicPoints;logical createOpenDX;
        struct {int type;char[*] description;} quantities[*];
        struct {int type;char[*] description;} outputdata[*];
        real scfwavefnctwindow[2];
        struct {char[*] name;
            struct {real emin;real emax;real de;
                struct {int type;char[*] description;} spin;
            } sections[*];
        } energywindows[*];
        struct {char[*] name;real kpoint[3];int bandindices[*];real energywindow[2];} kresolved[*];
    } grids[*]
    ={
    };


    struct {logical active;logical forceusewanniercenters;int homooffset;
        struct {int nint;int nky;} wanniercenters;
    } topological_insulator
    ={f,f,4
        ,{40,100}
    };


    struct {logical active;logical jointdos;logical bandoutput;logical stopafter;
        struct {int ne;real emin;real emax;} energy;
    } optics
    ={f,f,f,f
        ,{1000,0,10}
    };

};


section BasisSettings{

    struct {
        struct {int type;char[*] description;} version;
        int extensionlevel;logical add3d;logical addf;char[*] core4f_by_el[*];int core4f_by_sort[*];char[*] core4fNovalenceF_by_el[*];int core4fNovalenceF_by_sort[*];
        struct {real Q;real S;} multicore[*];
        struct {real Q;real S;} multisemicore[*];
    } basis
    ={
         {1,"default FPLO9 basis"}
        ,1,f,f,{},{},{},{}
        ,{
        }
        ,{
        }
    };

};


section dHvA{

    struct {logical active;int bandsup[*];int bandsdown[*];real isovalue;int subdiv[3];int bisections;char[*] wandata;logical deleteoldfiles;logical needbandweights;char[*] statefile;
        struct {int weightindex;real coef;} energyshift[*];
    } isosurface
    ={f,{},{},0.,{12,12,12},0,"",t,f,""
        ,{
        }
    };


    struct {
        struct {char[*] label;real field[3];} fields[*];
        int bands[*];int parts[*];
        struct {int type;char[*] description;} spin;
        int nphi;int phifrom;int phito;int nplane;int planefrom;int planeto;real areachainthreshold;real orbitcomparisonsearchoffset;real orbitsampleportion;real arearadiusfactor;real smoothnessthreshold;flag outputoptions[*];flag debugoptions[*];
    } dhva
    ={
         {
            {"[010]",{0.,1.,0.}}
            ,{"[001]",{0.,0.,1.}}
        }
        ,{},{}
        ,{0,"both"}
        ,30,-1,-2,100,-1,-2,0.1,0.2,0.1,1.,1.,{show_non_closed(-),output_orbits(-),output_extremal_orbits(+)},{determine_area_chain_extrema(-),sort_area_chains(-),sort_extrema_chains(-),extend_periodic_area_chains(-)
        ,collect_orbits(-),calc_orbit(-),run(+)}
    };

};


------------------------------------------------------------------------
End  : content of =.in
------------------------------------------------------------------------
========================================================================
                      CALCULATION OF CRYSTALL STRUCTURE

------------------------------------------------------------------------
                                INPUT DATA

------------------------------------------------------------------------
Unit of length    : Bohr radii
Spacegroup        : 221 - P m -3 m
Lattice constants : 5.4        5.4        5.4
Axis angles       : 90         90         90

               Wyckoff positions

Number of Wyckoff positions :    2
No.   Element   X          Y          Z
   1    Fe      0          0          0
   2    Fe      1/2        1/2        1/2

------------------------------------------------------------------------
------------------------------------------------------------------------
                            SYMMETRY CREATION

------------------------------------------------------------------------

                   SYMMETRYGROUP

Space group      : 221, HM: P m -3 m, Hall: -P 4 2 3
Point group      : Oh
Inversion        : yes
Symmorphic       : yes

                   OPERATIONS


(The XYZ-operation symbols "(A,B,C)+(t1,t2,t3)" mean:
 if a vector is given relative to the lattice basis vectors
 r=X*a1+Y*a2+Z*a3 then it is transformed into
 r'=(A+t1)*a1+(B+t2)*a2+(C+t3)*a2.


Full Group Operations : 48
Index   rotation          translation      : symbol
 0:   (+X  ,+Y  ,+Z  ) + (0   ,0   ,0   )  :   E
 1:   (-X  ,+Y  ,-Z  ) + (0   ,0   ,0   )  :   C2(y)
 2:   (+X  ,-Y  ,-Z  ) + (0   ,0   ,0   )  :   C2(x)
 3:   (-X  ,-Y  ,+Z  ) + (0   ,0   ,0   )  :   C2(z)
 4:   (-Y  ,+Z  ,-X  ) + (0   ,0   ,0   )  :   C3-(1-1-1)
 5:   (-Y  ,-Z  ,+X  ) + (0   ,0   ,0   )  :   C3-(-11-1)
 6:   (+Y  ,-Z  ,-X  ) + (0   ,0   ,0   )  :   C3-(-1-11)
 7:   (+Y  ,+Z  ,+X  ) + (0   ,0   ,0   )  :   C3-(111)
 8:   (-Z  ,-X  ,+Y  ) + (0   ,0   ,0   )  :   C3(1-1-1)
 9:   (+Z  ,+X  ,+Y  ) + (0   ,0   ,0   )  :   C3(111)
10:   (+Z  ,-X  ,-Y  ) + (0   ,0   ,0   )  :   C3(-11-1)
11:   (-Z  ,+X  ,-Y  ) + (0   ,0   ,0   )  :   C3(-1-11)
12:   (+Y  ,+X  ,-Z  ) + (0   ,0   ,0   )  :   C2(xy)
13:   (+Y  ,-X  ,+Z  ) + (0   ,0   ,0   )  :   C4-(z)
14:   (-Y  ,+X  ,+Z  ) + (0   ,0   ,0   )  :   C4(z)
15:   (-Y  ,-X  ,-Z  ) + (0   ,0   ,0   )  :   C2(x-y)
16:   (+Z  ,-Y  ,+X  ) + (0   ,0   ,0   )  :   C2(xz)
17:   (-Z  ,-Y  ,-X  ) + (0   ,0   ,0   )  :   C2(-xz)
18:   (-Z  ,+Y  ,+X  ) + (0   ,0   ,0   )  :   C4-(y)
19:   (+Z  ,+Y  ,-X  ) + (0   ,0   ,0   )  :   C4(y)
20:   (-X  ,-Z  ,-Y  ) + (0   ,0   ,0   )  :   C2(y-z)
21:   (+X  ,+Z  ,-Y  ) + (0   ,0   ,0   )  :   C4-(x)
22:   (-X  ,+Z  ,+Y  ) + (0   ,0   ,0   )  :   C2(yz)
23:   (+X  ,-Z  ,+Y  ) + (0   ,0   ,0   )  :   C4(x)
32:   (-X  ,-Y  ,-Z  ) + (0   ,0   ,0   )  :   I
33:   (+X  ,-Y  ,+Z  ) + (0   ,0   ,0   )  :   m(y)
34:   (-X  ,+Y  ,+Z  ) + (0   ,0   ,0   )  :   m(x)
35:   (+X  ,+Y  ,-Z  ) + (0   ,0   ,0   )  :   m(z)
36:   (+Y  ,-Z  ,+X  ) + (0   ,0   ,0   )  :   -3-(1-1-1)
37:   (+Y  ,+Z  ,-X  ) + (0   ,0   ,0   )  :   -3-(-11-1)
38:   (-Y  ,+Z  ,+X  ) + (0   ,0   ,0   )  :   -3-(-1-11)
39:   (-Y  ,-Z  ,-X  ) + (0   ,0   ,0   )  :   -3-(111)
40:   (+Z  ,+X  ,-Y  ) + (0   ,0   ,0   )  :   -3(1-1-1)
41:   (-Z  ,-X  ,-Y  ) + (0   ,0   ,0   )  :   -3(111)
42:   (-Z  ,+X  ,+Y  ) + (0   ,0   ,0   )  :   -3(-11-1)
43:   (+Z  ,-X  ,+Y  ) + (0   ,0   ,0   )  :   -3(-1-11)
44:   (-Y  ,-X  ,+Z  ) + (0   ,0   ,0   )  :   m(xy)
45:   (-Y  ,+X  ,-Z  ) + (0   ,0   ,0   )  :   -4-(z)
46:   (+Y  ,-X  ,-Z  ) + (0   ,0   ,0   )  :   -4(z)
47:   (+Y  ,+X  ,+Z  ) + (0   ,0   ,0   )  :   m(x-y)
48:   (-Z  ,+Y  ,-X  ) + (0   ,0   ,0   )  :   m(xz)
49:   (+Z  ,+Y  ,+X  ) + (0   ,0   ,0   )  :   m(-xz)
50:   (+Z  ,-Y  ,-X  ) + (0   ,0   ,0   )  :   -4-(y)
51:   (-Z  ,-Y  ,+X  ) + (0   ,0   ,0   )  :   -4(y)
52:   (+X  ,+Z  ,+Y  ) + (0   ,0   ,0   )  :   m(y-z)
53:   (-X  ,-Z  ,+Y  ) + (0   ,0   ,0   )  :   -4-(x)
54:   (+X  ,-Z  ,-Y  ) + (0   ,0   ,0   )  :   m(yz)
55:   (-X  ,+Z  ,-Y  ) + (0   ,0   ,0   )  :   -4(x)

                   TRANSLATION

Length units throughout calculations are Bohr radii.

lattice constants [aB] :  5.400000000000000      5.400000000000000      5.400000000000000
lattice constants [Ang]:  2.857556938320000      2.857556938320000      2.857556938320000
axis angles            :  90.000000000000000     90.000000000000000     90.000000000000000
bravais lattice        : Simple Cubic
conventional to primitive transformation
      b1  : 1    0    0
      b2  : 0    1    0
      b3  : 0    0    1
coordinate rotation: cell=(a1,a2,a3)
                         =R*default_cell
      R=       1.000000000000000      0.000000000000000      0.000000000000000
               0.000000000000000      1.000000000000000      0.000000000000000
               0.000000000000000      0.000000000000000      1.000000000000000
lattice vectors
      a1  :  5.400000000000000      0.000000000000000      0.000000000000000
      a2  :  0.000000000000000      5.400000000000000      0.000000000000000
      a3  :  0.000000000000000      0.000000000000000      5.400000000000000
reciprocal lattice vectors / 2*Pi
      g1  :  0.185185185185185      0.000000000000000      0.000000000000000
      g2  :  0.000000000000000      0.185185185185185      0.000000000000000
      g3  :  0.000000000000000      0.000000000000000      0.185185185185185

Volume (real cell) [aB^3] :  157.464000000000027
Volume (real cell) [Ang^3]:  23.333757393627479

------------------------------------------------------------------------
                              UNIT CELL CREATION

------------------------------------------------------------------------

               Wyckoff positions
Number of Wyckoff positions :    2
No.   Element    X                      Y                      Z                      Concentration
   1    Fe       0.000000000000000      0.000000000000000      0.000000000000000      1.000000000000000
   2    Fe       0.500000000000000      0.500000000000000      0.500000000000000      1.000000000000000


                 Atom sites

Number of sites :    2
Number of real sites :    2
Number of nonempty real sites :    2
No.  Element WPS CPA-Block    X                      Y                      Z
   1   Fe     1      1        0.000000000000000      0.000000000000000      0.000000000000000
   2   Fe     2      2       -2.700000000000000     -2.700000000000000     -2.700000000000000

Structure type: Crystal


Averaged WS radius :    2.65881 Bohr radii
Vuc                :            157.464000000000027 Bohr radii


------------------------------------------------------------------------

========================================================================
------------------------------------------------------------------------
Number of diferent sites (Naps) :     2
maximum atoms at one site       :     1
pointer ---    tau(1..3,1..nsite)    --- atoms
   1           0.0000  0.0000  0.0000        1      isrt :    1
                                                    isa  :    1
   2          -2.7000 -2.7000 -2.7000        1      isrt :    2
                                                    isa  :    2
------------------------------------------------------------------------

 ------- mapping of sort onto site ------------
sort    1 site    1
sort    2 site    2
 ----------------------------------------------
 ------- classification of sorts   ------------
number of inequivalent sites    2
icis     1 representing sort    1
icis     2 representing sort    2

number of inequivalent atoms    2
sort    1 representing sort    1
sort    2 representing sort    2
 ----------------------------------------------

========================================================================
                     RADIAL BASIS FUNCTION MESH
------------------------------------------------------------------------
  Mesh type        =Root
------------------------------------------------------------
 El. sort        nr      rmin     default_rmax   factor     rmax 
 Fe     1       200    0.100E-05    7.531871   1.000000   7.531871
 Fe     2       200    0.100E-05    7.531871   1.000000   7.531871

========================================================================

Nuclear charges
site sort real-site El  nuclear-charge averaged-charge
   1    1    1      Fe      26.00000    26.00000
   2    2    2      Fe      26.00000    26.00000
BASIS DEFINITION:

INFORMATION: Next output is =.basdef as it result from fedit basis definitions.
------------------------------------------------------------------------
Start: content of =.basdef
------------------------------------------------------------------------
    1 P=(0.62); N=(14);  1s2s2p : 3s3p / D4s Q=(0,4.8) P=(0.95,0.9); D3d Q=(0,3.27) P=(0.94,0.95); S4p Q=(2.9) P=(1);
    2 P=(0.62); N=(14);  1s2s2p : 3s3p / D4s Q=(0,4.8) P=(0.95,0.9); D3d Q=(0,3.27) P=(0.94,0.95); S4p Q=(2.9) P=(1);
------------------------------------------------------------------------
End  : content of =.basdef
------------------------------------------------------------------------
--------------------------------------------------------------------------
--------------------------------------------------------------------------
         Vatom-classes and ionicities
--------------------------------------------------------------------------
   sort     state           vat-class   Q      S     N       p        px
--------------------------------------------------------------------------
    1 Fe     1 1s      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    1 Fe     2 2s      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    1 Fe     3 2p      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    1 Fe     1 3s      val      1   0.0000   0.0000 14.00   0.6200   1.0000
    1 Fe     2 3p      val      1   0.0000   0.0000 14.00   0.6200   1.0000
    1 Fe     3 4s      val      1   0.0000   0.0000 14.00   0.6200   0.9500
    1 Fe     4 5s      val      2   4.8000   0.0000 14.00   0.6200   0.9000
    1 Fe     5 3d      val      1   0.0000   0.0000 14.00   0.6200   0.9400
    1 Fe     6 4d      val      3   3.2700   0.0000 14.00   0.6200   0.9500
    1 Fe     7 4p      val      4   2.9000   0.0000 14.00   0.6200   1.0000
    2 Fe     1 1s      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    2 Fe     2 2s      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    2 Fe     3 2p      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    2 Fe     1 3s      val      1   0.0000   0.0000 14.00   0.6200   1.0000
    2 Fe     2 3p      val      1   0.0000   0.0000 14.00   0.6200   1.0000
    2 Fe     3 4s      val      1   0.0000   0.0000 14.00   0.6200   0.9500
    2 Fe     4 5s      val      2   4.8000   0.0000 14.00   0.6200   0.9000
    2 Fe     5 3d      val      1   0.0000   0.0000 14.00   0.6200   0.9400
    2 Fe     6 4d      val      3   3.2700   0.0000 14.00   0.6200   0.9500
    2 Fe     7 4p      val      4   2.9000   0.0000 14.00   0.6200   1.0000
--------------------------------------------------------------------------

------------------------------------------------------------------------
|   OCCUPATION NUMBERS OF CORE STATES FOR CRYSTAL CALCULATION          |
------------------------------------------------------------------------
Sort:   1 Element: Fe
------------------------------------------------------------------------
 State  Spin   Occu    M =  -3     -2     -1      0      1      2      3
  1:1s  up      1.00   :                        1.00
        down    1.00   :                        1.00
  2:2s  up      1.00   :                        1.00
        down    1.00   :                        1.00
  3:2p  up      3.00   :                 1.00   1.00   1.00
        down    3.00   :                 1.00   1.00   1.00
Sort:   2 Element: Fe
------------------------------------------------------------------------
 State  Spin   Occu    M =  -3     -2     -1      0      1      2      3
  1:1s  up      1.00   :                        1.00
        down    1.00   :                        1.00
  2:2s  up      1.00   :                        1.00
        down    1.00   :                        1.00
  3:2p  up      3.00   :                 1.00   1.00   1.00
        down    3.00   :                 1.00   1.00   1.00
------------------------------------------------------------------------

Number of occupied valence bands (averaged):     32.00000


AUTOSUBDIVISION: |ga|,|gb|,|gc|  =     0.18518519      0.18518519      0.18518519
AUTOSUBDIVISION: input nkxyz     =             12              12              12
AUTOSUBDIVISION: centering factor=              1               1               1
AUTOSUBDIVISION: NO AUTO SUBDIVISION REQUESTED (if some nkxyz=0)
AUTOSUBDIVISION: 
AUTOSUBDIVISION: SUBDIVISIONS (NBZ = prod(nkxyz)*  1~=      1728):
AUTOSUBDIVISION: input nkxyz    =        12        12        12
AUTOSUBDIVISION: actual nkxyz   =        12        12        12 NBZ=      1728
AUTOSUBDIVISION: k-subdivisions =      0.0154      0.0154      0.0154

 BZMESH:    84 irreducible k-points from   1728 (   12   12   12 )
PREPGRADIENT: Numerical gradient formula, order=         1
PREPGRADIENT: 
PREPGRADIENT:  grad f(k)=sum_n ~b_n*(f(k+b_n)-f(k))
PREPGRADIENT: 
PREPGRADIENT:                    b                                      ijk                        ~b
PREPGRADIENT:          1 :       0.0000      0.0000      0.0833         0         0         1      0.0000      0.0000      6.0000
PREPGRADIENT:          2 :      -0.0833      0.0000      0.0000        -1         0         0     -6.0000      0.0000      0.0000
PREPGRADIENT:          3 :       0.0833      0.0000      0.0000         1         0         0      6.0000      0.0000      0.0000
PREPGRADIENT:          4 :       0.0000      0.0833      0.0000         0         1         0      0.0000      6.0000      0.0000
PREPGRADIENT:          5 :       0.0000     -0.0833      0.0000         0        -1         0      0.0000     -6.0000      0.0000
PREPGRADIENT:          6 :       0.0000      0.0000     -0.0833         0         0        -1      0.0000      0.0000     -6.0000

 TETIRR:      1728 inequivalent tetrahedra from     10368
------------------------------------------------------------------------

make_lattice: N1         6
make_lattice: N2         6
make_lattice: N3         6

 LOCAL SITE POTENTIAL IS NONSPHERICAL

AUTOMATIC MODE: TOLERANCE DETERMINED NUMBER OF FOURIER COMPONENTS
TOLERANCE FOR REAL    SPACE CUTOFF:  0.10E-004
TOLERANCE FOR FOURIER SPACE CUTOFF:  0.10E-004
EWALD PARAMETER                   :    0.66805
MAXIMUM RECIPROCAL VECTOR NORM    :    5.63663
 NUMBER OF FT. COMPONENTS FOR EWALD POTENIAL:    461

DIMENSIONS OF MATRICES
valence matrix dimension                        38
core    matrix dimension                        10
max number of valence states per atom           19
max number of core    states per atom            5
max number of bands                             38
    number of k-points                          84
size of v-v coefficient matrix              485184
size of c-v coefficient matrix              127680
------------------------------------------------------------------------
                        matrix index structure
------------------------------------------------------------------------
               core block
 site  El start_index  end_index  length 
    1  Fe      1           5         5
    2  Fe      6          10         5
              valence block
 site  El start_index  end_index  length 
    1  Fe      1          19        19
    2  Fe     20          38        19
------------------------------------------------------------------------

The following file has been created:
./+points
scalar relativistic calculation!


XC version :    9 - Perdew Wang 92               (LSDA)


ISO: active: off
ISO: deleteoldfiles: on
ISO: bandsup  : 
ISO: bandsdown: 
ISO: isoE: 0
ISO: initial subdiv: 12 12 12
ISO: bisections: 0
ISO: wandata: 
ISO: needbandweights: off
ISO: statefile: 
ISO: energyshift: 
ISO:
ISO: Reciprocal iso cell vectors/(2*pi):
ISO: g1:     -0.00000     -0.18519     -0.00000,  |g1|:      0.18519
ISO: g2:      0.18519      0.00000      0.00000,  |g2|:      0.18519
ISO: g3:      0.00000      0.00000      0.18519,  |g3|:      0.18519
ISO:
ISO: Most homogenous initial subdivision would be: 12 12 12 or multiples of these.
ISO:

========================================================================
                     RADIAL BASIS FUNCTION MESH
------------------------------------------------------------------------
  Mesh type        =Root
------------------------------------------------------------
 El. sort        nr      rmin     default_rmax   factor     rmax 
 Fe     1       200    0.100E-05    7.531871   1.000000   7.531871
 Fe     2       200    0.100E-05    7.531871   1.000000   7.531871

========================================================================

========================================================================
LOADING DENSITY from file =.dens

 
 CPU   : densin: cpu time:           0.01 sec
 

END DENSITY LOADING
========================================================================
Maximum radial mesh radius :              7.531871484566357
Maximum LOI Cut-off radius :             15.063742969132715
Maximum of       pairs:  88

-----------------------------------------------------------------------------
==========                      STRUCTURE DATA                     ==========
-----------------------------------------------------------------------------
| ATOM SITE  R_NN/2    nonempty          SHELLS             NEIGHBOURS      |
-----------------------------------------------------------------------------
|  Fe    1      2.3383     2.3383         10                168             |
-----------------------------------------------------------------------------
|  Fe    2      2.3383     2.3383         10                168             |
-----------------------------------------------------------------------------
-------------------------------------------------------------------------------
==========             TABLE OF NEIGHBOURS                           ==========
-------------------------------------------------------------------------------
| EL SITE  EL SITE       PAIRS                  SHELLS         RMIN     RMAX  |
-------------------------------------------------------------------------------
| Fe   1   Fe   1        80                      6             5.4000  13.2272|
-------------------------------------------------------------------------------
| Fe   1   Fe   2        88                      4             4.6765  14.0296|
-------------------------------------------------------------------------------
| Fe   2   Fe   2        80                      6             5.4000  13.2272|
-------------------------------------------------------------------------------


============================== START molinit============================== 
isa=   1
isa=   2
Total number of irreducible points:   1022
Total number of skipped     points:    312
============================== END   molinit============================== 


Preparing basis potentials ...
(Wanna see details ... set verbosity level to second highest.)
... Done

------------------------------------------------------------------------
==========         RELATIVISTIC ATOMIC ENERGIES               ==========
------------------------------------------------------------------------
| ATOM  ELECTRON           ENERGY      COMPRESSION     TYPE    SORT    |
------------------------------------------------------------------------
|  Fe     1s1 1/2        -256.57124911      none        core    1      |
------------------------------------------------------------------------
|  Fe     2s1 1/2         -30.04502293      none        core    1      |
------------------------------------------------------------------------
|  Fe     2p1 1/2         -25.94523337      none        core    1      |
------------------------------------------------------------------------
|  Fe     2p1 3/2         -25.48771739      none        core    1      |
------------------------------------------------------------------------
|  Fe     3s 1/2           -3.43436995  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3p 1/2           -2.23874794  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3p 3/2           -2.18167726  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4s 1/2           -0.13670452  4.436  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     5s 1/2           -0.77704156  4.203  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3d 3/2           -0.28486504  4.390  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3d 5/2           -0.27915229  4.390  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4d 3/2           -0.59912462  4.436  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4d 5/2           -0.59772687  4.436  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4p 1/2           -0.90805444  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4p 3/2           -0.89982858  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     1s1 1/2        -256.57124911      none        core    2      |
------------------------------------------------------------------------
|  Fe     2s1 1/2         -30.04502293      none        core    2      |
------------------------------------------------------------------------
|  Fe     2p1 1/2         -25.94523337      none        core    2      |
------------------------------------------------------------------------
|  Fe     2p1 3/2         -25.48771739      none        core    2      |
------------------------------------------------------------------------
|  Fe     3s 1/2           -3.43436995  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3p 1/2           -2.23874794  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3p 3/2           -2.18167726  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4s 1/2           -0.13670452  4.436  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     5s 1/2           -0.77704156  4.203  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3d 3/2           -0.28486504  4.390  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3d 5/2           -0.27915229  4.390  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4d 3/2           -0.59912462  4.436  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4d 5/2           -0.59772687  4.436  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4p 1/2           -0.90805444  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4p 3/2           -0.89982858  4.670  14.000 valence   2      |
------------------------------------------------------------------------
==========       SCALARRELATIVISTIC ATOMIC ENERGIES           ==========
------------------------------------------------------------------------
|  Fe     1s1            -256.57124911      none        core    1      |
------------------------------------------------------------------------
|  Fe     2s1             -30.04502293      none        core    1      |
------------------------------------------------------------------------
|  Fe     2p1             -25.64022272      none        core    1      |
------------------------------------------------------------------------
|  Fe     3s               -3.43436995  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3p               -2.20070082  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4s               -0.13670452  4.436  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     5s               -0.77704156  4.203  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     3d               -0.28143739  4.390  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4d               -0.59828597  4.436  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     4p               -0.90257054  4.670  14.000 valence   1      |
------------------------------------------------------------------------
|  Fe     1s1            -256.57124911      none        core    2      |
------------------------------------------------------------------------
|  Fe     2s1             -30.04502293      none        core    2      |
------------------------------------------------------------------------
|  Fe     2p1             -25.64022272      none        core    2      |
------------------------------------------------------------------------
|  Fe     3s               -3.43436995  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3p               -2.20070082  4.670  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4s               -0.13670452  4.436  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     5s               -0.77704156  4.203  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     3d               -0.28143739  4.390  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4d               -0.59828597  4.436  14.000 valence   2      |
------------------------------------------------------------------------
|  Fe     4p               -0.90257054  4.670  14.000 valence   2      |
------------------------------------------------------------------------
--------------------------------------------------
--------------------------------------------------
bascheck: core    states are o.k.
bascheck: valence states are o.k.

------------------------------------------------------------------------
|                  MINIMAL RADIUS FOR THE LOCAL BASIS FUNCTIONS        |
------------------------------------------------------------------------
| ATOM   ELECTRON    RMIN_LA  RMIN_SA   RMIN_LD  RMIN_SD  TYPE    SORT |
------------------------------------------------------------------------
|  Fe  1s1           0.9928   0.8899    0.0000   0.0000     Core     1 |
|  Fe  2s1           2.7851   2.3487    0.0000   0.0000     Core     1 |
|  Fe  2p1           2.9200   2.4435    2.5721   2.1832     Core     1 |
------------------------------------------------------------------------
|  Fe                2.9200   2.4435    2.5721   2.1832     Core     1 |
------------------------------------------------------------------------
|  Fe  3s            6.0534   5.5658    0.0000   0.0000   Valence    1 |
|  Fe  3p            6.2460   5.8863    5.9593   5.4795   Valence    1 |
|  Fe  4s            6.3787   6.2117    0.0000   0.0000   Valence    1 |
|  Fe  5s            6.3365   6.1407    0.0000   0.0000   Valence    1 |
|  Fe  3d            6.2309   6.0284    5.8949   5.7453   Valence    1 |
|  Fe  4d            6.3959   6.2249    5.9399   5.9496   Valence    1 |
|  Fe  4p            6.6444   6.4450    6.3277   6.1346   Valence    1 |
------------------------------------------------------------------------
|  Fe                6.6444   6.4450    6.3277   6.1346   Valence    1 |
------------------------------------------------------------------------
|  Fe  1s1           0.9928   0.8899    0.0000   0.0000     Core     2 |
|  Fe  2s1           2.7851   2.3487    0.0000   0.0000     Core     2 |
|  Fe  2p1           2.9200   2.4435    2.5721   2.1832     Core     2 |
------------------------------------------------------------------------
|  Fe                2.9200   2.4435    2.5721   2.1832     Core     2 |
------------------------------------------------------------------------
|  Fe  3s            6.0534   5.5658    0.0000   0.0000   Valence    2 |
|  Fe  3p            6.2460   5.8863    5.9593   5.4795   Valence    2 |
|  Fe  4s            6.3787   6.2117    0.0000   0.0000   Valence    2 |
|  Fe  5s            6.3365   6.1407    0.0000   0.0000   Valence    2 |
|  Fe  3d            6.2309   6.0284    5.8949   5.7453   Valence    2 |
|  Fe  4d            6.3959   6.2249    5.9399   5.9496   Valence    2 |
|  Fe  4p            6.6444   6.4450    6.3277   6.1346   Valence    2 |
------------------------------------------------------------------------
|  Fe                6.6444   6.4450    6.3277   6.1346   Valence    2 |
------------------------------------------------------------------------
------------------------------------------------------------------------
| max extension of core    states=     2.91996589                      |
| max extension of valence states=     6.64435521                      |
| max extension of state         =     6.64435521                      |
------------------------------------------------------------------------
 
 CPU   : radial equation: cpu time:           0.10 sec
 
 
 CPU   : Total One-center integrals: cpu time:           0.01 sec
 
 
 CPU   : Total Two-center integrals: cpu time:           0.02 sec
 
SCF: Iteration version (   2) : Iterat   : Hyperspace interpolation             
SCF: Tolerance         =    0.10E-05
SCF: Vector dimension  =        2400
SCF: Max hyperspace dim=     5
SCF: Max mixing        =    0.40E+00
SCF: Max steps         =    30
SCF: 
SCF: iteration  0  dimension  0  last deviation u=  0.00E+00
SCF:               interpolated  new deviation  u=  0.00E+00  step p= 0.400
SCF: 
best Fname=.dens

------------------------------------------------------------------------
Writing =.dens ...!
------------------------------------------------------------------------------
                                Density Analysis
site elem. Spin up       Spin down    Sum           Nuclear        Excess
           (electrons)   (electrons)  (electrons)   (protons)      electrons
------------------------------------------------------------------------------
  1  Fe   13.75426994   12.24573010   26.00000004   26.00000000    0.00000004
  2  Fe   12.24572470   13.75427489   25.99999959   26.00000000   -0.00000041
------------------------------------------------------------------------------
sum  :    25.99999464   26.00000499   51.99999962   52.00000000   -0.00000038
------------------------------------------------------------------------------
 
------------------------------------------------------------------------
==========                      CHARGE                        ==========
------------------------------------------------------------------------
  nu. charge    el. charge    mag. moment     spin up       spin down
       52.00    51.9999996     -0.0000103    25.9999946     26.0000050
 
-------------------------------------------------
|  ATOM   |  SITE  |  MAG.MOMENT  |  NU.CHARGE  |
-------------------------------------------------
|   Fe    |   1    |  1.5085398   |     26.00   |
|   Fe    |   2    | -1.5085502   |     26.00   |
-------------------------------------------------
(Scalar)Relativistic and no finite nucleus: no r=0 densities.
 
 CPU   : Coulomb-potential: cpu time:           0.00 sec
 
calc_nv_erg: onsit env    :                 -4072.19828586
calc_nv_erg: onsit exc0   :                  -106.47474688
calc_nv_erg: onsit envxc0 :                  -140.94544606
calc_nv_erg: onsit e0     :                 -1962.70160614
calc_nv_erg: offsite      :                    -0.26919208
calc_nv_erg: n*vc-11      :                 -1962.97079821
calc_nv_erg: n*vat        :                   -17.63508399
calc_nv_erg: offsit exc0  :                   -12.25424556
calc_nv_erg: offsit envxc0:                   -15.91926862
calc_nv_erg: n*excat      :                  -118.72899244
calc_nv_erg: n*vxcat      :                  -156.86471468
 
 CPU   : calc_nv_erg: cpu time:           0.03 sec
 


============================== START molpottogrid============================== 
molpottogrid: EXC(n)             :                  -107.26135864
molpottogrid: EXC(at)            :                   118.72882877
molpottogrid: EXC                :                  -107.26152230
molpottogrid: n*VXC(n)           :                  -141.97792875
molpottogrid: n*VXC(at)          :                   156.86450197
molpottogrid: ENVXC              :                  -141.97814146
molpottogrid: Int(n*v)           :                    21.81104622
molpottogrid: Int(n*VCoulomb)    :                    -0.05127604
molpottogrid: Int(Z*VCoulomb)    :                     0.06542014
 
 CPU   : MOLPOTpot -> grid: cpu time:           0.02 sec
 
============================== END   molpottogrid============================== 


erg_hart_0:secnd term  onsite  :         -3045.312128
erg_hart_0: Int(Z*VCoulomb) tot:         -3045.246707
 
 CPU   : make potential: cpu time:           0.05 sec
 


============================== START loi_int============================== 
Allocated p_gdiff(    89,    48,     2,     2) =    0.06835 MByte
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
 
 CPU   : LOICore one-center integrals: cpu time:           0.00 sec
 
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
 
 CPU   : LOIThree-center integrals: cpu time:           0.21 sec
 
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
real site    1
real site    2
 
 CPU   : LOIThree-center vc-integrals: cpu time:           0.09 sec
 
============================== END   loi_int============================== 


 
 CPU   : LO integrals: cpu time:           0.30 sec
 
 
 CPU   : Total core transformation of S/Hcv: cpu time:           0.00 sec
 
DIAGONALIZATION:   1.19% done
 
K-vector     1 :     0.0000    0.0000    0.0000
Band energies at spin   1
    -3.3360556    -3.2795257    -2.1025548    -2.1025548    -2.1025548
    -2.0468177    -2.0468177    -2.0468177    -0.4501802    -0.2927229
    -0.2927229    -0.2022767    -0.2022767    -0.2022767    -0.1372874
    -0.1372874    -0.1005749    -0.1005749    -0.1005749     0.2351369
     0.2351369     0.2351369     0.4427112     0.4427112     0.6070370
     1.0292664     1.0292664     1.0362858     1.0362858     1.0362858
     1.0947760     1.0947760     1.0947760     1.7705817     1.7705817
     1.7705817     1.8701366     2.2568602
DIAGONALIZATION:   9.52% done
DIAGONALIZATION:  19.05% done
DIAGONALIZATION:  28.57% done
DIAGONALIZATION:  38.10% done
DIAGONALIZATION:  47.62% done
DIAGONALIZATION:  57.14% done
DIAGONALIZATION:  66.67% done
DIAGONALIZATION:  76.19% done
DIAGONALIZATION:  85.71% done
DIAGONALIZATION:  95.24% done
DIAGONALIZATION: 100.00% done
DIAGONALIZATION:   1.19% done
 
K-vector     1 :     0.0000    0.0000    0.0000
Band energies at spin   2
    -3.3360559    -3.2795259    -2.1025551    -2.1025551    -2.1025551
    -2.0468179    -2.0468179    -2.0468179    -0.4501802    -0.2927230
    -0.2927230    -0.2022769    -0.2022769    -0.2022769    -0.1372875
    -0.1372875    -0.1005750    -0.1005750    -0.1005750     0.2351368
     0.2351368     0.2351368     0.4427111     0.4427111     0.6070369
     1.0292664     1.0292664     1.0362857     1.0362857     1.0362857
     1.0947759     1.0947759     1.0947759     1.7705816     1.7705816
     1.7705816     1.8701365     2.2568601
DIAGONALIZATION:   9.52% done
DIAGONALIZATION:  19.05% done
DIAGONALIZATION:  28.57% done
DIAGONALIZATION:  38.10% done
DIAGONALIZATION:  47.62% done
DIAGONALIZATION:  57.14% done
DIAGONALIZATION:  66.67% done
DIAGONALIZATION:  76.19% done
DIAGONALIZATION:  85.71% done
DIAGONALIZATION:  95.24% done
DIAGONALIZATION: 100.00% done
 EFRANG: ebmin=  -3.3360559,  ebmax=   2.9132883
 EFRANG: e1=  -0.1619811,  e2=  -0.1198408
 BANDWT: ---- Tetrahedron Integration ---- 
              Estimating Fermi energy ...
      Estimate             Window            Tolerance     N(E_f)
     ------------------------------------------------------------------
    -0.1387971    -0.1388039  -0.1383825   0.42140D-03    173.681225
    -0.1387971    -0.1387997  -0.1387955   0.42140D-05    173.720069
    -0.1387971    -0.1387971  -0.1387971   0.42140D-07    173.720874
    -0.1387971    -0.1387971  -0.1387971   0.42140D-09    173.721013
 BANDWT: Density of states at E_f:   N(E_f) =  173.721013 states/Ha
 TETWTS: Fermi energy: -0.138797; 32.000000 electrons
         Band energy:   -41.716861881070, including Bloechl correction:-0.000963
 BANDWT: spin:  1
 BANDWT: highest fully occupied band :    14
 BANDWT: lowest fully unoccupied band:    18
 BANDWT: spin:  2
 BANDWT: highest fully occupied band :    14
 BANDWT: lowest fully unoccupied band:    18
 BANDWT: WARNING: Gap information is not reliable! Check DOS.
 BANDWT: spin up   gap :    0.00000 eV
 BANDWT: spin down gap :    0.00000 eV
 BANDWT: total     gap :    0.00000 eV
 
 CPU   : det_bandweights: cpu time:           0.01 sec
 
 
 CPU   : Total KS-DOS: cpu time:           0.00 sec
 
------------------------------------------------------------------------
                              Population Analysis 
                           (valence electron numbers) 


       El Site                                                                                  SUM       Excess electrons
                    3s         3p         4s         5s         3d         4d         4p    
N_net  Fe   1    2.01163    6.03823    0.39656    0.00033    6.13401    0.01784    0.15321    14.75181
N_gros Fe   1    2.00092    6.00001    0.74075   -0.00414    6.55356    0.11661    0.59229    16.00000    0.00000

S_net  Fe   1   -0.00027   -0.00356    0.02760    0.00003    1.51016   -0.00480    0.01692     1.54609
S_gros Fe   1    0.00005   -0.00159    0.02833   -0.00069    1.45481   -0.02086    0.03973     1.49979

                    3s         3p         4s         5s         3d         4d         4p    
N_net  Fe   2    2.01163    6.03823    0.39656    0.00033    6.13401    0.01784    0.15321    14.75181
N_gros Fe   2    2.00092    6.00001    0.74075   -0.00414    6.55356    0.11661    0.59229    16.00000   -0.00000

S_net  Fe   2    0.00027    0.00356   -0.02760   -0.00003   -1.51018    0.00480   -0.01692    -1.54610
S_gros Fe   2   -0.00005    0.00159   -0.02833    0.00069   -1.45482    0.02086   -0.03973    -1.49980


total net   electron number:                                                                  29.50362
total gross electron number:                                                                  32.00000    0.00000

total net   spin  :                                                                          -0.00002
total gross spin  :                                                                          -0.00001

------------------------------------------------------------------------

 
 CPU   : Kohn-Sham equation: cpu time:           0.12 sec
 
 
 CPU   : Net Occupation Numbers: cpu time:           0.02 sec
 
 
 CPU   : Net density: cpu time:           0.01 sec
 
 
 CPU   : Occupation numbers: cpu time:           0.03 sec
 
 
 CPU   : Overlap density: cpu time:           0.13 sec
 
------------------------------------------------------------------------------
                                Density Analysis
site elem. Spin up       Spin down    Sum           Nuclear        Excess
           (electrons)   (electrons)  (electrons)   (protons)      electrons
------------------------------------------------------------------------------
  1  Fe   13.75427136   12.24572889   26.00000025   26.00000000    0.00000025
  2  Fe   12.24572127   13.75427811   25.99999938   26.00000000   -0.00000062
------------------------------------------------------------------------------
sum  :    25.99999262   26.00000700   51.99999962   52.00000000   -0.00000038
------------------------------------------------------------------------------
 
------------------------------------------------------------------------
==========                      CHARGE                        ==========
------------------------------------------------------------------------
  nu. charge    el. charge    mag. moment     spin up       spin down
       52.00    51.9999996     -0.0000144    25.9999926     26.0000070
 
-------------------------------------------------
|  ATOM   |  SITE  |  MAG.MOMENT  |  NU.CHARGE  |
-------------------------------------------------
|   Fe    |   1    |  1.5085425   |     26.00   |
|   Fe    |   2    | -1.5085568   |     26.00   |
-------------------------------------------------
(Scalar)Relativistic and no finite nucleus: no r=0 densities.
 
 CPU   : Density calculation: cpu time:           0.17 sec
 
 
------------------------------------------------------------------------
==========                   TOTAL ENERGY                     ==========
------------------------------------------------------------------------
Ecore:    -1451.90798853
EBand:    -1493.62485041
-T*S:          0.00000000
this is actually E_entropy= -T S

            total energy         kinetic energy    potential energy    ex.-corr. energy
EE:     -2541.132830797662        2574.39747322      -5008.26878171       -107.26152230
this is the total energy without -T S

free energy E-T S:     -2541.132830797662

corrected energy (E+E-T S)/2:     -2541.132830797662

free energy should be force consistent,
except for Bloechl correction (if it applies)


------------------------------------------------------------------------

estack     :       -2541.13283080
eav  ,deav :       -2541.13283080            0.00000000

------------------------------------------------------------------------

 
 CPU   : fplo step: cpu time:           0.67 sec
 
best Fname=.dens

------------------------------------------------------------------------
Writing =.dens ...!
 
 CPU   : fplo cycle: cpu time:           0.67 sec
 
best Tname=.dens_best

------------------------------------------------------------------------
Writing =.dens_best ...!
SCF: 
SCF: iteration  1  dimension  0  last deviation u=  0.90E-06 CONVERGED
best Fname=.dens

------------------------------------------------------------------------
Writing =.dens ...!
 
 CPU   : total fplo calculation: cpu time:           0.67 sec
 
TERMINATION: Finished : SCF calculation
