------------------------------------------------------------------------
|                                                                      |
|    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                       |
|                                                                      |
------------------------------------------------------------------------
|                                                                      |
|   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: 14.00                                                  |
| sub  version: M-CPA                                                  |
| release     : 45                                                     |
|                                                                      |
------------------------------------------------------------------------
| date        : Sun Mar  2 15:45:55 2014                               |
| host        : unas                                                   |
------------------------------------------------------------------------
File =.sym exists!
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
    ={"14.00","M-CPA"};


    struct {char[*] hostname;char[*] date;char[*] by;} last_changes
    ={"unas","Sun Mar  2 15:45:46 2014","fedit (14.00-45-x86_64)"};

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


section structure_information{

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


    struct {int number;char[*] symbol;} spacegroup
    ={225,"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.
    };
    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,{40,40,40}};


    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"};

    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)"};

    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;
};


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 {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;} topological_insulator
    ={f};


    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}
    };

};


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

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

------------------------------------------------------------------------
Unit of length    : Bohr radii
Spacegroup        :  225 - FM3M
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 - FM3M
Point group      :  37 - 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,
 for RTG lattices the operations refer to the corresponding
 STG unit cell!)


Group Generators      : 3
Index   rotation          translation      : symbol
14:   (-Y  ,+X  ,+Z  ) + (0   ,0   ,0   )  :   C1/4(z)
 9:   (+Z  ,+X  ,+Y  ) + (0   ,0   ,0   )  :   C1/3(111)
32:   (-X  ,-Y  ,-Z  ) + (0   ,0   ,0   )  :   I

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   )  :   C1/3(-111)
 5:   (-Y  ,-Z  ,+X  ) + (0   ,0   ,0   )  :   C1/3(1-11)
 6:   (+Y  ,-Z  ,-X  ) + (0   ,0   ,0   )  :   C2/3(-1-11)
 7:   (+Y  ,+Z  ,+X  ) + (0   ,0   ,0   )  :   C2/3(111)
 8:   (-Z  ,-X  ,+Y  ) + (0   ,0   ,0   )  :   C2/3(-111)
 9:   (+Z  ,+X  ,+Y  ) + (0   ,0   ,0   )  :   C1/3(111)
10:   (+Z  ,-X  ,-Y  ) + (0   ,0   ,0   )  :   C2/3(1-11)
11:   (-Z  ,+X  ,-Y  ) + (0   ,0   ,0   )  :   C1/3(-1-11)
12:   (+Y  ,+X  ,-Z  ) + (0   ,0   ,0   )  :   C2(xy)
13:   (+Y  ,-X  ,+Z  ) + (0   ,0   ,0   )  :   C3/4(z)
14:   (-Y  ,+X  ,+Z  ) + (0   ,0   ,0   )  :   C1/4(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(x-z)
18:   (-Z  ,+Y  ,+X  ) + (0   ,0   ,0   )  :   C3/4(y)
19:   (+Z  ,+Y  ,-X  ) + (0   ,0   ,0   )  :   C1/4(y)
20:   (-X  ,-Z  ,-Y  ) + (0   ,0   ,0   )  :   C2(y-z)
21:   (+X  ,+Z  ,-Y  ) + (0   ,0   ,0   )  :   C3/4(x)
22:   (-X  ,+Z  ,+Y  ) + (0   ,0   ,0   )  :   C2(yz)
23:   (+X  ,-Z  ,+Y  ) + (0   ,0   ,0   )  :   C1/4(x)
32:   (-X  ,-Y  ,-Z  ) + (0   ,0   ,0   )  :   I
33:   (+X  ,-Y  ,+Z  ) + (0   ,0   ,0   )  :   s(y)
34:   (-X  ,+Y  ,+Z  ) + (0   ,0   ,0   )  :   s(x)
35:   (+X  ,+Y  ,-Z  ) + (0   ,0   ,0   )  :   s(z)
36:   (+Y  ,-Z  ,+X  ) + (0   ,0   ,0   )  :   S5/6(-111)
37:   (+Y  ,+Z  ,-X  ) + (0   ,0   ,0   )  :   S5/6(1-11)
38:   (-Y  ,+Z  ,+X  ) + (0   ,0   ,0   )  :   S1/6(-1-11)
39:   (-Y  ,-Z  ,-X  ) + (0   ,0   ,0   )  :   S1/6(111)
40:   (+Z  ,+X  ,-Y  ) + (0   ,0   ,0   )  :   S1/6(-111)
41:   (-Z  ,-X  ,-Y  ) + (0   ,0   ,0   )  :   S5/6(111)
42:   (-Z  ,+X  ,+Y  ) + (0   ,0   ,0   )  :   S1/6(1-11)
43:   (+Z  ,-X  ,+Y  ) + (0   ,0   ,0   )  :   S5/6(-1-11)
44:   (-Y  ,-X  ,+Z  ) + (0   ,0   ,0   )  :   s(xy)
45:   (-Y  ,+X  ,-Z  ) + (0   ,0   ,0   )  :   S1/4(z)
46:   (+Y  ,-X  ,-Z  ) + (0   ,0   ,0   )  :   S3/4(z)
47:   (+Y  ,+X  ,+Z  ) + (0   ,0   ,0   )  :   s(x-y)
48:   (-Z  ,+Y  ,-X  ) + (0   ,0   ,0   )  :   s(xz)
49:   (+Z  ,+Y  ,+X  ) + (0   ,0   ,0   )  :   s(x-z)
50:   (+Z  ,-Y  ,-X  ) + (0   ,0   ,0   )  :   S1/4(y)
51:   (-Z  ,-Y  ,+X  ) + (0   ,0   ,0   )  :   S3/4(y)
52:   (+X  ,+Z  ,+Y  ) + (0   ,0   ,0   )  :   s(y-z)
53:   (-X  ,-Z  ,+Y  ) + (0   ,0   ,0   )  :   S1/4(x)
54:   (+X  ,-Z  ,-Y  ) + (0   ,0   ,0   )  :   s(yz)
55:   (-X  ,+Z  ,-Y  ) + (0   ,0   ,0   )  :   S3/4(x)

                   TRANSLATION

lattice constants:  7.550000000000000      7.550000000000000      7.550000000000000
axis angles      :  90.000000000000000     90.000000000000000     90.000000000000000
bravais lattice  : FCC
primitive to bravais transformation
      b1  : 0    1/2  1/2
      b2  : 1/2  0    1/2
      b3  : 1/2  1/2  0
unit vectors of primitive lattice (SC)
      u1  :  1.000000000000000      0.000000000000000      0.000000000000000
      u2  :  0.000000000000000      1.000000000000000      0.000000000000000
      u3  :  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
reciprocial 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


------------------------------------------------------------------------
                              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
------------------------------------------------------------------------
Start: content of =.basdef
------------------------------------------------------------------------
    1 P=(0.638); N=(14); 1s: /  S2s Q=(0) P=(1);  S2p Q=(0) P=(1);  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



 BZMESH:  1661 irreducible k-points from  64000 ( 40 40 40 )

 TETIRR:     63869 inequivalent tetrahedra from    384000
------------------------------------------------------------------------

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 RECIPROCIAL 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                        1661
size of v-v coefficient matrix              960058
size of c-v coefficient matrix               56474
------------------------------------------------------------------------
                        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
------------------------------------------------------------------------

scalar relativistic calculation!


XC version :    9 - Perdew Wang 92               (LSDA)  


========================================================================
                     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:    541
Total number of skipped     points:    136
============================== 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.14 sec
 
 
 CPU   : Total One-center integrals: cpu time:           0.01 sec
 
 
 CPU   : Total Two-center integrals: cpu time:           0.05 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: 

------------------------------------------------------------------------
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.01 sec
 
calc_nv_erg: onsit env    :                  -376.39869073
calc_nv_erg: onsit exc0   :                   -17.44646464
calc_nv_erg: onsit envxc0 :                   -23.04645998
calc_nv_erg: onsit e0     :                  -176.29029926
calc_nv_erg: offsite      :                    -0.10175961
calc_nv_erg: n*vc-11      :                  -176.39205886
calc_nv_erg: n*vat        :                    -4.08609946
calc_nv_erg: offsit exc0  :                    -2.78180276
calc_nv_erg: offsit envxc0:                    -3.61415514
calc_nv_erg: n*excat      :                   -20.22826740
calc_nv_erg: n*vxcat      :                   -26.66061512
 
 CPU   : calc_nv_erg: cpu time:           0.06 sec
 


============================== START molpottogrid============================== 
molpottogrid: EXC(n)             :                   -17.70912651
molpottogrid: EXC(at)            :                    20.22834468
molpottogrid: EXC                :                   -17.70904924
molpottogrid: n*VXC(n)           :                   -23.38974555
molpottogrid: n*VXC(at)          :                    26.66071553
molpottogrid: ENVXC              :                   -23.38964513
molpottogrid: Int(n*v)           :                     4.30501781
molpottogrid: Int(n*VCoulomb)    :                    -0.00301253
molpottogrid: Int(Z*VCoulomb)    :                     0.00213432
 
 CPU   : MOLPOTpot -> grid: cpu time:           0.02 sec
 
============================== END   molpottogrid============================== 


erg_hart_0:secnd term  onsite  :          -289.992155
erg_hart_0: Int(Z*VCoulomb) tot:          -289.990021
 
 CPU   : make potential: cpu time:           0.09 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.01 sec
 
real site    1
real site    1
real site    1
real site    1
 
 CPU   : LOIThree-center integrals: cpu time:           0.32 sec
 
real site    1
real site    1
real site    1
real site    1
 
 CPU   : LOIThree-center vc-integrals: cpu time:           0.22 sec
 
============================== END   loi_int============================== 


 
 CPU   : LO integrals: cpu time:           0.56 sec
 


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


 
 CPU   : LOF integrals: cpu time:           0.23 sec
 
 
 CPU   : Total core transformation of S/Hcv: cpu time:           0.00 sec
 
DIAGONALIZATION:   0.06% done
 
K-vector     1 :     0.0000    0.0000    0.0000
Band energies at spin   1
    -3.8850707    -2.4954471    -2.4954471    -2.4954471    -0.5337943
     0.3659375     0.3659375     0.3659375     0.4218312     0.4218312
     0.4218312     0.6246157     0.6246157     0.8914051     0.9145126
     0.9145126     0.9145126
DIAGONALIZATION:   9.99% done
DIAGONALIZATION:  19.99% done
DIAGONALIZATION:  29.98% done
DIAGONALIZATION:  39.98% done
DIAGONALIZATION:  49.97% done
DIAGONALIZATION:  59.96% done
DIAGONALIZATION:  69.96% done
DIAGONALIZATION:  79.95% done
DIAGONALIZATION:  89.95% done
DIAGONALIZATION:  99.94% done
DIAGONALIZATION: 100.00% done
 EFRANG: ebmin=  -3.8850707,  ebmax=   2.4696113
 EFRANG: e1=  -0.2803904,  e2=   0.3659375
 BANDWT: ---- Tetrahedron Integration ---- 
              Estimating Fermi energy ...
      Estimate             Window            Tolerance     N(E_f)
     ------------------------------------------------------------------
    -0.1159784    -0.1188084  -0.1123451   0.64633D-02     10.754762
    -0.1159102    -0.1159646  -0.1158999   0.64633D-04     10.582840
    -0.1159102    -0.1159103  -0.1159096   0.64633D-06     10.585503
    -0.1159102    -0.1159102  -0.1159102   0.64633D-08     10.585482
    -0.1159102    -0.1159102  -0.1159102   0.64633D-10     10.583475
 BANDWT: Density of states at E_f:   N(E_f) =   10.583475 states/Ha
 TETWTS: Fermi energy: -0.115910; 11.000000 electrons
         Band energy:   -23.59226338, including Bloechl correction:-0.000384
------------------------------------------------------------------------
    Plasmon frequencies from Lindhard expression (intra band)

spin : 1
DOS(ef)      =     10.5855 [state/Ha]
<|vk|>       =     15.0284 [10^5 m/s]
sqrt(<vk^2>) =     15.5157 [10^5 m/s]
|vk| interval= [    1.4128,   19.0591] [10^5 m/s]

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

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   80.2453    80.2453    80.2453

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


------------------------------------------------------------------------
 BANDWT: spin:  1
 BANDWT: highest fully occupied band :     5
 BANDWT: lowest fully unoccupied band:     8
 BANDWT: WARNING: Gap information is insecure! Check DOS.
 BANDWT: total     gap :    0.00000 eV
 
 CPU   : det_bandweights: cpu time:           0.60 sec
 
START: optics matrix elements
OPTMATEL:   0.06% done
OPTMATEL:   9.99% done
OPTMATEL:  19.99% done
OPTMATEL:  29.98% done
OPTMATEL:  39.98% done
OPTMATEL:  49.97% done
OPTMATEL:  59.96% done
OPTMATEL:  69.96% done
OPTMATEL:  79.95% done
OPTMATEL:  89.95% done
OPTMATEL:  99.94% done
OPTMATEL: 100.00% done
 
 CPU   : optics matrix elements: cpu time:           0.88 sec
 
tetopt:   0.00% done
tetopt:  10.00% done
tetopt:  20.00% done
tetopt:  30.00% done
tetopt:  39.99% done
tetopt:  49.99% done
tetopt:  59.99% done
tetopt:  69.99% done
tetopt:  79.99% done
tetopt:  89.99% done
tetopt:  99.99% 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.02024    0.00207    0.67254    0.00310    0.04146     9.75308
N_gros Al   1    1.99993    6.00076    1.22027   -0.01898    1.54231   -0.04585    0.30155    11.00000    0.00000


total net   electron number:                                                                   9.75308
total gross electron number:                                                                  11.00000    0.00000


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

 
 CPU   : Kohn-Sham equation: cpu time:           4.35 sec
 
 
 CPU   : Net density: cpu time:           0.00 sec
 
 
 CPU   : Occupation numbers: cpu time:           0.48 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  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.61 sec
 
 
------------------------------------------------------------------------
==========                   TOTAL ENERGY                     ==========
------------------------------------------------------------------------
Ecore:     -110.41265739
EBand:     -134.00492077
T*S:          0.00000000
            total energy         kinetic energy    potential energy    ex.-corr. energy
EE:        -241.9192901388         242.17485162       -466.38509252        -17.70904924

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

estack     :        -241.91929014
eav  ,deav :        -241.91929014            0.00000000

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

 
 CPU   : fplo step: cpu time:           5.91 sec
 

------------------------------------------------------------------------
Writing =.dens ...!
 
 CPU   : fplo cycle: cpu time:           5.91 sec
 
SCF: 
SCF: iteration  1  dimension  0  last deviation u=  0.17E-07 CONVERGED

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