A fundamental understanding of the radiation damage effects in solids is of great importance in assisting the development of improved materials with ultra-high strength, toughness, and radiation resistance for nuclear energy applications. In this presentation, we show our recent theoretical investigation on the magnetic structure evolution of bulk iron in the region surrounding the radiation defects. We applied the locally self-consistent multiple scattering method (LSMS), a linear scaling ab-initio method based on density functional theory with local spin density approximation, to the study of the magnetic structure in a low energy cascade in a 10,000-atom sample for a series of time steps for the evolution of the defects. The primary damage state and the evolution of all defects in the sample were simulated using molecular dynamics with empirical, embedded-atom inter-atomic potentials. We also discuss the importance of thermal effect on the magnetic structure evolution.