# By setting 'Unrestricted Yes' we do a spin polarized calculation. Normally
# this would converge to the ferromagnetic solution.
# With the SpinFlipRegion keyword we make sure that we start with an antiferromagnetic
# For antiferromagnetic iron we need a larger unit cell of two atoms. Since
# these atoms appear to the program as symmetry equivalent we have to specify
# them as separate types.
Fe 0.0 0.0 0.0
Fe -1.435 -1.435 1.435 region=flip
-1.435 1.435 1.435
1.435 -1.435 1.435
2.87 2.87 -2.87
Title Beta iron
SpinFlipRegion flip # Flip (startup) spin density at second atom