qpt52 qpt53 qpt54 qpt55 qpt56 qpt57 qpt58 qpt59 qpt60 qpt61 qpt62 qpt63 qpt64 qpt65 qpt66 qpt67 qpt68 qpt69 qpt70

qpt71 qpt72 qpt73 kptopt2 2 kptopt 3 tolvrs 1.0d-15 iscf 5 getwfk 1

0.50000000000000 0.25000000000000 0.83333333333333 -0.41666666666667 0.25000000000000 0.83333333333333 -0.33333333333333 0.25000000000000 0.83333333333333 -0.25000000000000 0.25000000000000 0.83333333333333 -0.16666666666667 0.25000000000000 0.83333333333333 0.41666666666667 0.33333333333333 0.83333333333333 0.50000000000000 0.33333333333333 0.83333333333333 -0.41666666666667 0.33333333333333 0.83333333333333 -0.33333333333333 0.33333333333333 0.83333333333333 -0.25000000000000 0.33333333333333 0.83333333333333 0.50000000000000 0.41666666666667 0.83333333333333 -0.41666666666667 0.41666666666667 0.83333333333333 -0.33333333333333 0.41666666666667 0.83333333333333 -0.41666666666667 0.50000000000000 0.83333333333333 0.50000000000000 0.33333333333333 0.16666666666667 -0.41666666666667 0.33333333333333 0.16666666666667 -0.33333333333333 0.33333333333333 0.16666666666667 -0.41666666666667 0.41666666666667 0.16666666666667 -0.33333333333333 0.41666666666667 0.16666666666667 -0.25000000000000 0.41666666666667 0.16666666666667 -0.33333333333333 0.50000000000000 0.16666666666667 -0.25000000000000 0.50000000000000 0.25000000000000 # Automatic generation of k points, taking # into account the time-reversal symmetry only # SCF stopping criterion # Self-consistent calculation, using algorithm 5 # Read the ground-state wavefunctions

#Definition of occupation numbers occopt 4 tsmear 0.01 #Definition of the unit cell #acell 3*7.5041 # This is equivalent to 7.60 7.60 7.60 acell 3*7.6497601142 # This is equivalent to 7.60 7.60 7.60 rprim 0.0 0.5 0.5 # FCC primitive vectors (to be scaled by acell) 0.5 0.0 0.5 0.5 0.5 0.0

##Optimization of the lattice parameters #optcell 1 #ionmov 3 #ntime 10 #dilatmx 1.05 #ecutsm 0.5 #Definition of the atom types ntypat 1 # There is only one type of atom znucl 13 # The keyword "znucl" refers to the atomic number of the # possible type(s) of atom. The pseudopotential(s) # mentioned in the "files" file must correspond # to the type(s) of atom. Here, the only type is Aluminum

#Definition of the atoms natom 1 # There is only one atom per cell typat 1 # This atom is of type 1, that is, Aluminum xred # This keyword indicate that the location of the atoms # will follow, one triplet of number for each atom 0.0 0.0 0.0 # Triplet giving the REDUCED coordinate of atom 1. #Definition of the planewave basis set ecut 50.0 # Maximal kinetic energy cut-off, in Hartree #Exchange-correlation functional ixc 11 # LDA Teter Pade parametrization nband 6 #Definition of the k-point grid ngkpt 12 12 12 # This is a 2x2x2 FCC grid, based on the primitive vectors nshiftk 4 # of the reciprocal space. For a FCC real space lattice, # like the present one, it actually corresponds to the # so-called 4x4x4 Monkhorst-Pack grid, if the following shifts # are used : shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 #Definition of the SCF procedure nstep 100 # Maximal number of SCF cycles #to

ldfe 1.0d-6 # Will stop when, twice in a row, the difference # between two consecutive evaluations of total energy # differ by less than toldfe (in Hartree) # This value is way too large for most realistic studies of materials