The crystal structure of Fe16N2 was optimized by calculating atomic forces on the basis of an augmented-plane-wave force method. We found interesting cooperative structural relaxations of the iron atoms around the nitrogen, though the optimized structure is not far from the experimental proposal. Calculated hyperfine fields for each iron site for the optimized Fe16N2 show an interesting systematic feature when compared with experimental results.