SYSTEM   = Al(111)-2x2-4Layer + O2 (fcc) single point MS-RPBEl

LMIXTAU  =.TRUE.       ! Include the kinetic energy density in the density mixing
LASPH    =.TRUE.       ! Include non-spherical contributions related to the gradient of the density in the PAW spheres

METAGGA  = SRMGGA      ! Parameters needed for the DF
SRVAL    = 1.0         ! The parameters employed here correspond to MS-RPBEl
FINTB    = 1.0
SRAC     = 1.0

SRTA0    = 1
SRA0MU   = 0.12345679012345678
SRA0KP   = 0.804
SRA0C    = 0.0767086

SRTA1    = 1
SRA1MU   = 0.12345679012345678
SRA1KP   = 0.804

SRTB0    = 1
SRB0MU   = 0.12345679012345678
SRB0KP   = 0.804
SRB0C    = 0.0767086

SRTB1    = 1
SRB1MU   = 0.12345679012345678
SRB1KP   = 0.804

PREC     = Accurate    ! OptAlgo precision

ALGO     = Conjugate   ! Better to use a more stable algorithm for spin-polarized calculations
ENCUT    = 600         ! Planewave cut-off

NELM     = 800         ! Max N of SCF cycles
NELMIN   = 4           !
NELMDL   =-12          !

ISTART   = 0           ! Initialization of wavefunction
ICHARG   = 1
INIWAV   = 1           !
ISPIN    = 2
ISYM     = 0
MAGMOM   = 16*1.0 2*4.0

IBRION   =-1           ! Single point

ISMEAR   = 1           ! Methfessel-Paxton smearing (1st order)
SIGMA    = 0.2         ! Smearing width (eV)

NSW      = 0           ! Maximum of ionic steps
EDIFF    = 1.E-7       ! Convergence w.r.t. energy (eV)
EDIFFG   =-0.005       ! Convergence w.r.t. forces
ADDGRID  =.TRUE.       ! Support grid for the forces

LCHARGE  =.FALSE.      ! Don't write large, useless files
LWAVE    =.FALSE.      ! Don't write large, useless files

NPAR     = 2           ! Parallelization option
KPAR     = 1           !
LPLANE   =.TRUE.       !
NSIM     = 4           !