The ABINIT Users Mailing List ( CLOSED )

[i.galanakis@fz-juelich.de]

Dear All, 
I a new user of abinit and all my experience comes from all-electron methods. 
I had now problem with non-spinpolarized calculations but when I try to 
include spin a get always a non-magentic solution. I attach my input file. 
ANy idea why it does not work 

# SCF for NiMnSb 

# irdwfk 1 


#Definition of the unit cell
acell 3*11.2041        # 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

#Definition of the atom types
ntypat 3          # There is only one type of atom
znucl 28 25 51    # 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 3           # There is only one atom per cell
typat 1 2 3       # 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 Ni
   0.25  0.25  0.25     # Triplet giving the REDUCED coordinate of atom Mn
   0.75  0.75  0.75     # Triplet giving the REDUCED coordinate of atom Sb


#Exchange-correlation functional
ixc 1             # LDA Teter Pade parametrization

# Spin parameters
nsppol 2
nspden 2
nspinor 1
spinat 0.0 0.0 1.0
       0.0 0.0 3.0
       0.0 0.0 0.0
 
#Definition of the planewave basis set
ecut  12.0         # Maximal kinetic energy cut-off, in Hartree

# Calculation parameters
iscf 5            # Kind of mixing 
nstep 10          # Maximal number of SCF cycles
toldfe 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) 
nsym 0            # finds automatically the symmetry

#Definition of occupation numbers
occopt 4
tsmear 0.01
nband 20

#Definition of the k-point grid
kptopt 1          # Option for the automatic generation of k points, taking
                  # into account the symmetry
ngkpt 8 8 8       # 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