------------------------------------------------------------------------
|                                                                      |
|    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: 22.00                                                  |
| sub  version: M-CPA                                                  |
| release     : 62                                                     |
|                                                                      |
------------------------------------------------------------------------
| compiled with ifort 1310                                             |
| compiled with own LAPACK copy                                        |
------------------------------------------------------------------------
|                                                                      |
------------------------------------------------------------------------
| date        : Sun Jul  3 18:19:25 2022                               |
| host        : unas                                                   |
------------------------------------------------------------------------

File =.in exists!
------------------------------------------------------------------------
Compound: Al example
------------------------------------------------------------------------
------------------------------------------------------------------------
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
    ={"22.00","M-CPA"};


    struct {char[*] hostname;char[*] date;char[*] by;} last_changes
    ={"unas","Wed Mar 31 15:05:13 2021","fplo (22.00-62)"};

    char[*] compound="Al example";
};


section structure_information{

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


    struct {char[*] number;char[*] setting;char[*] symbol;} spacegroup
    ={"225","default","FM3M"};

    int subgroupgenerators[*]={
    };

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

    real lattice_constants[3]={7.55,7.55,7.55
    };
    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=1;

    struct {char[2] element;real tau[3];} wyckoff_positions[nsort]
    ={
        {"Al",{0.,0.,0.}}
    };

};


section structure_dependend{

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


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


    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",{1,0,0}}
        ,{"W",{1,1/2,0}}
        ,{"K",{3/4,3/4,0}}
        ,{"$~G",{0,0,0}}
        ,{"L",{1/2,1/2,1/2}}
        ,{"W",{1,1/2,0}}
        ,{"U",{1,1/4,1/4}}
        ,{"X",{1,0,0}}
    };

};


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,{60,60,60}};


    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
    ={1,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

    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
    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
    ={t,f,f,f
        ,{1000,0,20}
    };

};


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        : 225 - F m -3 m
Lattice constants : 7.55       7.55       7.55
Axis angles       : 90.        90.        90.

               Wyckoff positions

Number of Wyckoff positions :    1
No.   Element   X          Y          Z
   1    Al      0.         0.         0.

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

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

                   SYMMETRYGROUP

Space group      : 225, HM: F m -3 m, Hall: -F 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] :  7.550000000000000      7.550000000000000      7.550000000000000
lattice constants [Ang]:  3.995287941540000      3.995287941540000      3.995287941540000
axis angles            :  90.000000000000000     90.000000000000000     90.000000000000000
bravais lattice        : Face Centered Cubic
conventional to primitive transformation
      b1  : 0    1/2  1/2
      b2  : 1/2  0    1/2
      b3  : 1/2  1/2  0
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  :  0.000000000000000      3.775000000000000      3.775000000000000
      a2  :  3.775000000000000      0.000000000000000      3.775000000000000
      a3  :  3.775000000000000      3.775000000000000      0.000000000000000
reciprocal lattice vectors / 2*Pi
      g1  : -0.132450331125828      0.132450331125828      0.132450331125828
      g2  :  0.132450331125828     -0.132450331125828      0.132450331125828
      g3  :  0.132450331125828      0.132450331125828     -0.132450331125828

Volume (real cell) [aB^3] :  107.592218749999986
Volume (real cell) [Ang^3]:  15.943521882808749

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

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

               Wyckoff positions
Number of Wyckoff positions :    1
No.   Element    X                      Y                      Z                      Concentration
   1    Al       0.000000000000000      0.000000000000000      0.000000000000000      1.000000000000000


                 Atom sites

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

Structure type: Crystal


Averaged WS radius :    2.95051 Bohr radii
Vuc                :            107.592218749999986 Bohr radii


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

========================================================================
------------------------------------------------------------------------
Number of diferent sites (Naps) :     1
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
------------------------------------------------------------------------

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

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

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

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

Nuclear charges
site sort real-site El  nuclear-charge averaged-charge
   1    1    1      Al      13.00000    13.00000
BASIS DEFINITION:

INFORMATION: Next output is =.basdef as it result from fedit basis definitions.
------------------------------------------------------------------------
Start: content of =.basdef
------------------------------------------------------------------------
    1 P=(0.638); N=(14);  1s : 2s2p / D3s Q=(0,3) P=(0.95,1.05); D3p Q=(0,2.5) P=(1,0.97); S3d Q=(3) P=(1.05);
------------------------------------------------------------------------
End  : content of =.basdef
------------------------------------------------------------------------
--------------------------------------------------------------------------
--------------------------------------------------------------------------
         Vatom-classes and ionicities
--------------------------------------------------------------------------
   sort     state           vat-class   Q      S     N       p        px
--------------------------------------------------------------------------
    1 Al     1 1s      core     1   0.0000   0.0000  0.00   0.0000   0.0000
    1 Al     1 2s      val      1   0.0000   0.0000 14.00   0.6380   1.0000
    1 Al     2 2p      val      1   0.0000   0.0000 14.00   0.6380   1.0000
    1 Al     3 3s      val      1   0.0000   0.0000 14.00   0.6380   0.9500
    1 Al     4 4s      val      2   3.0000   0.0000 14.00   0.6380   1.0500
    1 Al     5 3p      val      1   0.0000   0.0000 14.00   0.6380   1.0000
    1 Al     6 4p      val      3   2.5000   0.0000 14.00   0.6380   0.9700
    1 Al     7 3d      val      2   3.0000   0.0000 14.00   0.6380   1.0500
--------------------------------------------------------------------------

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

Number of occupied valence bands (averaged):     11.00000


AUTOSUBDIVISION: |ga|,|gb|,|gc|  =     0.22941070      0.22941070      0.22941070
AUTOSUBDIVISION: input nkxyz     =             60              60              60
AUTOSUBDIVISION: centering factor=              1               1               1
AUTOSUBDIVISION: NO AUTO SUBDIVISION REQUESTED (if some nkxyz=0)
AUTOSUBDIVISION: 
AUTOSUBDIVISION: SUBDIVISIONS (NBZ = prod(nkxyz)*  1~=    216000):
AUTOSUBDIVISION: input nkxyz    =        60        60        60
AUTOSUBDIVISION: actual nkxyz   =        60        60        60 NBZ=    216000
AUTOSUBDIVISION: k-subdivisions =      0.0038      0.0038      0.0038

 BZMESH:  5216 irreducible k-points from 216000 (   60   60   60 )
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.0167      0.0167     -0.0167         0         0         1      7.5000      7.5000     -7.5000
PREPGRADIENT:          2 :      -0.0167     -0.0167     -0.0167        -1        -1        -1     -7.5000     -7.5000     -7.5000
PREPGRADIENT:          3 :       0.0167     -0.0167     -0.0167        -1         0         0      7.5000     -7.5000     -7.5000
PREPGRADIENT:          4 :       0.0167     -0.0167      0.0167         0         1         0      7.5000     -7.5000      7.5000
PREPGRADIENT:          5 :      -0.0167     -0.0167      0.0167         0         0        -1     -7.5000     -7.5000      7.5000
PREPGRADIENT:          6 :      -0.0167      0.0167      0.0167         1         0         0     -7.5000      7.5000      7.5000
PREPGRADIENT:          7 :       0.0167      0.0167      0.0167         1         1         1      7.5000      7.5000      7.5000
PREPGRADIENT:          8 :      -0.0167      0.0167     -0.0167         0        -1         0     -7.5000      7.5000     -7.5000

 TETIRR:    215809 inequivalent tetrahedra from   1296000
------------------------------------------------------------------------

OPTICS:
OPTICS: structure of tensor
OPTICS:
OPTICS: X  0  0 
OPTICS: 0  X  0 
OPTICS: 0  0  X 
OPTICS:
OPTICS: P(1,1)= P(2,2)= P(3,3)=  0.333*P(1,1)   0.333*P(2,2)   0.333*P(3,3) 
OPTICS:
OPTICS:
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.60212
MAXIMUM RECIPROCAL VECTOR NORM    :    5.08034
 NUMBER OF FT. COMPONENTS FOR EWALD POTENIAL:    259

DIMENSIONS OF MATRICES
valence matrix dimension                        17
core    matrix dimension                         1
max number of valence states per atom           17
max number of core    states per atom            1
max number of bands                             17
    number of k-points                        5216
size of v-v coefficient matrix             3014848
size of c-v coefficient matrix              177344
------------------------------------------------------------------------
                        matrix index structure
------------------------------------------------------------------------
               core block
 site  El start_index  end_index  length 
    1  Al      1           1         1
              valence block
 site  El start_index  end_index  length 
    1  Al      1          17        17
------------------------------------------------------------------------

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.13245      0.13245      0.13245,  |g1|:      0.22941
ISO: g2:      0.13245     -0.13245      0.13245,  |g2|:      0.22941
ISO: g3:      0.13245      0.13245     -0.13245,  |g3|:      0.22941
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 
 Al     1       200    0.100E-05    8.356605   1.000000   8.356605

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

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

 
 CPU   : densin: cpu time:           0.00 sec
 

END DENSITY LOADING
========================================================================
Maximum radial mesh radius :              8.356604954364068
Maximum LOI Cut-off radius :             16.713209908728135
Maximum of       pairs: 176

-----------------------------------------------------------------------------
==========                      STRUCTURE DATA                     ==========
-----------------------------------------------------------------------------
| ATOM SITE  R_NN/2    nonempty          SHELLS             NEIGHBOURS      |
-----------------------------------------------------------------------------
|  Al    1      2.6693     2.6693          9                176             |
-----------------------------------------------------------------------------
-------------------------------------------------------------------------------
==========             TABLE OF NEIGHBOURS                           ==========
-------------------------------------------------------------------------------
| EL SITE  EL SITE       PAIRS                  SHELLS         RMIN     RMAX  |
-------------------------------------------------------------------------------
| Al   1   Al   1       176                      9             5.3387  16.0160|
-------------------------------------------------------------------------------


============================== START molinit============================== 
isa=   1
Total number of irreducible points:    531
Total number of skipped     points:    146
============================== END   molinit============================== 


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

------------------------------------------------------------------------
==========         RELATIVISTIC ATOMIC ENERGIES               ==========
------------------------------------------------------------------------
| ATOM  ELECTRON           ENERGY      COMPRESSION     TYPE    SORT    |
------------------------------------------------------------------------
|  Al     1s1 1/2         -55.27785854      none        core    1      |
------------------------------------------------------------------------
|  Al     2s 1/2           -3.95013389  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     2p 1/2           -2.57235621  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     2p 3/2           -2.55614407  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3s 1/2           -0.26980869  5.065  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     4s 1/2           -0.37428095  5.598  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3p 1/2           -0.05952439  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3p 3/2           -0.05854706  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     4p 1/2           -0.08000254  5.172  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     4p 3/2           -0.07847538  5.172  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3d 3/2           -0.57658215  5.598  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3d 5/2           -0.57650774  5.598  14.000 valence   1      |
------------------------------------------------------------------------
==========       SCALARRELATIVISTIC ATOMIC ENERGIES           ==========
------------------------------------------------------------------------
|  Al     1s1             -55.27785854      none        core    1      |
------------------------------------------------------------------------
|  Al     2s               -3.95013389  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     2p               -2.56154812  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3s               -0.26980869  5.065  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     4s               -0.37428095  5.598  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3p               -0.05887284  5.332  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     4p               -0.07898443  5.172  14.000 valence   1      |
------------------------------------------------------------------------
|  Al     3d               -0.57653751  5.598  14.000 valence   1      |
------------------------------------------------------------------------
--------------------------------------------------
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 |
------------------------------------------------------------------------
|  Al  1s1           2.0016   1.7091    0.0000   0.0000     Core     1 |
------------------------------------------------------------------------
|  Al                2.0016   1.7091    0.0000   0.0000     Core     1 |
------------------------------------------------------------------------
|  Al  2s            6.3678   5.4253    0.0000   0.0000   Valence    1 |
|  Al  2p            6.7127   6.0528    6.0649   5.3254   Valence    1 |
|  Al  3s            7.1273   6.9054    0.0000   0.0000   Valence    1 |
|  Al  4s            7.8462   7.6153    0.0000   0.0000   Valence    1 |
|  Al  3p            7.4776   7.2491    6.9508   6.8684   Valence    1 |
|  Al  4p            7.2413   7.1007    6.7564   6.7682   Valence    1 |
|  Al  3d            7.7589   7.5084    6.8866   7.1540   Valence    1 |
------------------------------------------------------------------------
|  Al                7.8462   7.6153    6.9508   7.1540   Valence    1 |
------------------------------------------------------------------------
------------------------------------------------------------------------
| max extension of core    states=     2.00164122                      |
| max extension of valence states=     7.84621620                      |
| max extension of state         =     7.84621620                      |
------------------------------------------------------------------------
 
 CPU   : radial equation: cpu time:           0.04 sec
 
 
 CPU   : Total One-center integrals: cpu time:           0.01 sec
 
 
 CPU   : Total Two-center integrals: cpu time:           0.01 sec
 
SCF: Iteration version (   2) : Iterat   : Hyperspace interpolation             
SCF: Tolerance         =    0.10E-05
SCF: Vector dimension  =         600
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  Al    6.50000286    6.50000286   13.00000573   13.00000000    0.00000573
------------------------------------------------------------------------------
sum  :     6.50000286    6.50000286   13.00000573   13.00000000    0.00000573
------------------------------------------------------------------------------
 
------------------------------------------------------------------------
==========                      CHARGE                        ==========
------------------------------------------------------------------------
  nu. charge    el. charge    mag. moment     spin up       spin down
       13.00    13.0000057      0.0000000     6.5000029      6.5000029
 
-------------------------------------------------
|  ATOM   |  SITE  |  MAG.MOMENT  |  NU.CHARGE  |
-------------------------------------------------
|   Al    |   1    |  0.0000000   |     13.00   |
-------------------------------------------------
(Scalar)Relativistic and no finite nucleus: no r=0 densities.
 
 CPU   : Coulomb-potential: cpu time:           0.00 sec
 
calc_nv_erg: onsit env    :                  -376.39914797
calc_nv_erg: onsit exc0   :                   -17.44649537
calc_nv_erg: onsit envxc0 :                   -23.04650035
calc_nv_erg: onsit e0     :                  -176.29022485
calc_nv_erg: offsite      :                    -0.10175629
calc_nv_erg: n*vc-11      :                  -176.39198113
calc_nv_erg: n*vat        :                    -4.08605266
calc_nv_erg: offsit exc0  :                    -2.78178317
calc_nv_erg: offsit envxc0:                    -3.61412962
calc_nv_erg: n*excat      :                   -20.22827855
calc_nv_erg: n*vxcat      :                   -26.66062997
 
 CPU   : calc_nv_erg: cpu time:           0.01 sec
 


============================== START molpottogrid============================== 
molpottogrid: EXC(n)             :                   -17.70914565
molpottogrid: EXC(at)            :                    20.22835579
molpottogrid: EXC                :                   -17.70906841
molpottogrid: n*VXC(n)           :                   -23.38977087
molpottogrid: n*VXC(at)          :                    26.66073034
molpottogrid: ENVXC              :                   -23.38967050
molpottogrid: Int(n*v)           :                     4.30548446
molpottogrid: Int(n*VCoulomb)    :                    -0.00304946
molpottogrid: Int(Z*VCoulomb)    :                     0.00218319
 
 CPU   : MOLPOTpot -> grid: cpu time:           0.01 sec
 
============================== END   molpottogrid============================== 


erg_hart_0:secnd term  onsite  :          -289.992644
erg_hart_0: Int(Z*VCoulomb) tot:          -289.990461
 
 CPU   : make potential: cpu time:           0.02 sec
 


============================== START loi_int============================== 
Allocated p_gdiff(   177,    48,     1,     1) =    0.03398 MByte
real site    1
real site    1
real site    1
real site    1
 
 CPU   : LOICore one-center integrals: cpu time:           0.00 sec
 
real site    1
real site    1
real site    1
real site    1
 
 CPU   : LOIThree-center integrals: cpu time:           0.08 sec
 
real site    1
real site    1
real site    1
real site    1
 
 CPU   : LOIThree-center vc-integrals: cpu time:           0.06 sec
 
============================== END   loi_int============================== 


 
 CPU   : LO integrals: cpu time:           0.14 sec
 


============================== START loif_int============================== 
 
 CPU   : One-center force integrals: cpu time:           0.00 sec
 
 
 CPU   : Two-center force integrals: cpu time:           0.05 sec
 
 
 CPU   : Total One-center force integrals: cpu time:           0.00 sec
 
 
 CPU   : Total Two-center force integrals: cpu time:           0.05 sec
 
============================== END   loif_int============================== 


 
 CPU   : LOF integrals: cpu time:           0.05 sec
 
 
 CPU   : Total core transformation of S/Hcv: cpu time:           0.00 sec
 
DIAGONALIZATION:   0.02% done
 
K-vector     1 :     0.0000    0.0000    0.0000
Band energies at spin   1
    -3.8850279    -2.4954049    -2.4954049    -2.4954049    -0.5337801
     0.3659528     0.3659528     0.3659528     0.4218489     0.4218489
     0.4218489     0.6246334     0.6246334     0.8914247     0.9145261
     0.9145261     0.9145261
DIAGONALIZATION:   9.99% done
DIAGONALIZATION:  19.98% done
DIAGONALIZATION:  29.97% done
DIAGONALIZATION:  39.95% done
DIAGONALIZATION:  49.94% done
DIAGONALIZATION:  59.93% done
DIAGONALIZATION:  69.92% done
DIAGONALIZATION:  79.91% done
DIAGONALIZATION:  89.90% done
DIAGONALIZATION:  99.88% done
DIAGONALIZATION: 100.00% done
 EFRANG: ebmin=  -3.8850279,  ebmax=   2.4696250
 EFRANG: e1=  -0.2803757,  e2=   0.3659528
 BANDWT: ---- Tetrahedron Integration ---- 
              Estimating Fermi energy ...
      Estimate             Window            Tolerance     N(E_f)
     ------------------------------------------------------------------
    -0.1160995    -0.1187935  -0.1123303   0.64633D-02     10.842031
    -0.1160664    -0.1160790  -0.1160143   0.64633D-04     10.882948
    -0.1160664    -0.1160667  -0.1160660   0.64633D-06     10.881741
    -0.1160664    -0.1160664  -0.1160664   0.64633D-08     10.881741
    -0.1160664    -0.1160664  -0.1160664   0.64633D-10     10.881032
 BANDWT: Density of states at E_f:   N(E_f) =   10.881032 states/Ha
 TETWTS: Fermi energy: -0.116066; 11.000000 electrons
         Band energy:   -23.591882319178, including Bloechl correction:-0.000174
------------------------------------------------------------------------
    Plasmon frequencies from Lindhard expression (intra band)

spin : 1
DOS(ef)      =     10.8817 [state/Ha]
<|vk|>       =     14.8868 [10^5 m/s]
sqrt(<vk^2>) =     15.4290 [10^5 m/s]
|vk| interval= [    2.2204,   18.9655] [10^5 m/s]

Fermi velocity tensor:
<v o v> [10^5 m/s]^2    79.3508     0.0000     0.0000
                         0.0000    79.3508     0.0000
                         0.0000     0.0000    79.3508

1st main axis            1.0000     0.0000     0.0000
2nd main axis            0.0000     1.0000     0.0000
3rd main axis            0.0000     0.0000     1.0000

eigenval.[10^5 m/s]^2   79.3508    79.3508    79.3508

Plasma frequencies along main axis from Lindhard (intra band)
w_p              12.4912     12.4912     12.4912   eV


------------------------------------------------------------------------
 BANDWT: spin:  1
 BANDWT: highest fully occupied band :     5
 BANDWT: lowest fully unoccupied band:     8
 BANDWT: WARNING: Gap information is not reliable! Check DOS.
 BANDWT: total     gap :    0.00000 eV
 
 CPU   : det_bandweights: cpu time:           0.41 sec
 
START: optics matrix elements
OPTMATEL:   0.02% done
OPTMATEL:   9.99% done
OPTMATEL:  19.98% done
OPTMATEL:  29.97% done
OPTMATEL:  39.95% done
OPTMATEL:  49.94% done
OPTMATEL:  59.93% done
OPTMATEL:  69.92% done
OPTMATEL:  79.91% done
OPTMATEL:  89.90% done
OPTMATEL:  99.88% done
OPTMATEL: 100.00% done
 
 CPU   : optics matrix elements: cpu time:           0.72 sec
 
tetopt:   0.00% done
tetopt:  10.00% done
tetopt:  20.00% done
tetopt:  30.00% done
tetopt:  40.00% done
tetopt:  50.00% done
tetopt:  60.00% done
tetopt:  70.00% done
tetopt:  80.00% done
tetopt:  90.00% done
tetopt: 100.00% done
tetopt: 100.00% done
 
 CPU   : Total KS-DOS: cpu time:           0.00 sec
 
------------------------------------------------------------------------
                              Population Analysis 
                           (valence electron numbers) 


       El Site                                                                                  SUM       Excess electrons
                    2s         2p         3s         4s         3p         4p         3d    
N_net  Al   1    2.00299    6.01067    1.02195    0.00208    0.67197    0.00310    0.04146     9.75422
N_gros Al   1    1.99993    6.00076    1.22116   -0.01905    1.54153   -0.04585    0.30152    11.00000   -0.00000


total net   electron number:                                                                   9.75422
total gross electron number:                                                                  11.00000   -0.00000


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

 
 CPU   : Kohn-Sham equation: cpu time:           3.26 sec
 
 
 CPU   : Net Occupation Numbers: cpu time:           0.07 sec
 
 
 CPU   : Net density: cpu time:           0.01 sec
 
 
 CPU   : Occupation numbers: cpu time:           0.32 sec
 
 
 CPU   : Overlap density: cpu time:           0.02 sec
 
------------------------------------------------------------------------------
                                Density Analysis
site elem. Spin up       Spin down    Sum           Nuclear        Excess
           (electrons)   (electrons)  (electrons)   (protons)      electrons
------------------------------------------------------------------------------
  1  Al    6.50000286    6.50000286   13.00000573   13.00000000    0.00000573
------------------------------------------------------------------------------
sum  :     6.50000286    6.50000286   13.00000573   13.00000000    0.00000573
------------------------------------------------------------------------------
 
------------------------------------------------------------------------
==========                      CHARGE                        ==========
------------------------------------------------------------------------
  nu. charge    el. charge    mag. moment     spin up       spin down
       13.00    13.0000057      0.0000000     6.5000029      6.5000029
 
-------------------------------------------------
|  ATOM   |  SITE  |  MAG.MOMENT  |  NU.CHARGE  |
-------------------------------------------------
|   Al    |   1    |  0.0000000   |     13.00   |
-------------------------------------------------
(Scalar)Relativistic and no finite nucleus: no r=0 densities.
 
 CPU   : Density calculation: cpu time:           0.35 sec
 
 
------------------------------------------------------------------------
==========                   TOTAL ENERGY                     ==========
------------------------------------------------------------------------
Ecore:     -110.41256499
EBand:     -134.00444731
-T*S:          0.00000000
this is actually E_entropy= -T S

            total energy         kinetic energy    potential energy    ex.-corr. energy
EE:      -241.919290784995         242.17526886       -466.38549124        -17.70906841
this is the total energy without -T S

free energy E-T S:      -241.919290784995

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

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


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

estack     :        -241.91929078
eav  ,deav :        -241.91929078            0.00000000

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

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

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

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

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