Dear users

I have a problem which can be (i think) solved by fixing the spin value at each atom. The problem is as follows;

I am trying to calculate the ground state of an anti-ferromagnetic (AFM) system. Although i know the value of final spins approximately, i do not know the lattice constant and position of the atoms for that system. When i start from the spin unpolarized geometry or from AFM geometry obtained with some other pseudopotential i get 0 magnetic moment at each atom. Therefore, i thought, to be sure that i cannot have AFM state with this pseudopotential, i need to fix the spin of each atom to the value that i expect and relax the system for the geometry. Then, after obtaining the optimized geometry i can relax the geometry and spin values together to see if it is going to go zero again or give me some spin value.

I could not find a variable that tells abinit to keep the spin at each atom fixed. Does anyone know how to do it or how to solve my problem woth some other method?

Cheers,

sener