# Initial Parameters
    ISPIN = 2 # specifies spin polarization
    PREC = Accurate
    LREAL = Auto
    ISYM = -1
    ISTART = 2

    # Calculation
    ISMEAR = -5 # determines how partial occupancies f_nk are set for each orbital
    SIGMA = 0.05 # determines the width of the smearing in eV

    # Magnetic Moment
    #MAGMOM = 8*5.0 8*-5.0 3*-2.0 2*-4.0 2*-1.0 1*-4.0 32*0.0 
    IDIPOL = 0
    LDIPOL = .TRUE.

    # Electronic Relaxation
    ALGO = Fast # selects a fairly robust mixture of the Davidson and RMM-DIIS algorithms
    ENCUT = 550 # specifies the cutoff energy for the planewave basis set in eV
    NELM = 500 # sets the maximum number of electronic SC (selfconsistency) steps which may be performed
    EDIFF = 1e-06 # specifies the global break condition for the electronic SC-loop
    ICHARG = 11

    # Density Functional Theory
    LMAXMIX = 4 # controls up to which l-quantum number the one-center PAW charge densities are passed through the charge density mixer and written to the CHGCAR file

    # Hubbard Parameters
    LDAU = .TRUE.
    LDAUTYPE = 2 # LDAUTYPE = 1 Liechtenstein, LDAUTYPE = 2 Dudarev
    LDAUL = 2 2 2 2 2 0 
    LDAUU = 5.3 6.2 5.3 0 5.3 0 
    LDAUJ = 4.0 3.5 4.0 3.5 4.0 0  
    LDAUPRINT = 1

    # Ionic Relaxation
    IBRION = 2 # ionic relaxation (conjugate gradient algorithm)
    NSW = 0 # sets the maximum number of ionic steps
    EDIFFG = -1e-02 # force convergence, relaxation will stop if all forces are smaller than |EDIFFG|
    ISIF = 2 # ISIF = 3 for bulk and = 2 for surfaces

    # Output
    LORBIT = 11 # DOSCAR and lm-decomposed PROCAR written
    NWRITE = 2 # determines how much is written to the file OUTCAR
    LWAVE = .FALSE. # determines whether the wavefunctions are written to the WAVECAR file at the end of a run
    LCHARG = .FALSE. # determines whether the charge densities (files CHGCAR and CHG) are written

    # Mixing Parameters
    AMIX = 0.1              
    AMIX_MAG = 0.2 
    BMIX = 0.0001 
    BMIX_MAG = 0.0001  


    # Non Collinear Flags
    LNONCOLLINEAR =.TRUE. # specifies whether fully non-collinear magnetic calculations are performed.
    LSORBIT = .TRUE. # switches on SOC and automatically turns LNONCOLLINEAR = .TRUE.
    NBANDS = 640
    

    #Magnetic Moment
SAXIS = 1 0 0
    LORBMOM = .TRUE.

    #CONSTRAINMENT
    I_CONSTRAINED_M = 2
    LAMBDA = 0.01
M_CONSTR = 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 -5 0 0 -5 0 0 -5 0 0 -5 0 0 -5 0 0 -5 0 0 -5 0 0 -5 0 0 -2 0 0 -2 0 0 -2 0 0 -4 0 0 -4 0 0 -1 0 0 -1 0 0 -4 0 0 96*0.6