2006ApJ...637..466N

We present Monte Carlo simulations of radiative transfer in magnetar atmospheres. We include the effects of vacuum polarization, electron and proton scattering, and free-free absorption. Simulations are performed for the atmosphere model with the magnetic field perpendicular and also tilted with respect to the neutron star surface, and we show that the average spectrum does not strongly depend on the orientation of the magnetic field. We investigate the region of the parameter space where the vacuum absorption-like feature appears in the spectrum, and we analyze the shape of the proton cyclotron line. Our results indicate that the existence of the vacuum polarization feature should be a general attribute of soft gamma-ray repeaters' burst spectra, provided that the energy release takes place at a sufficiently dense region and the atmosphere scale height is large enough. We discuss the existence of such a feature in recent observational data on these sources.