2021MNRAS.500.5369B

The study of polarized radiation transfer in the highly magnetized surface locales of neutron stars is of great interest to the understanding of accreting X-ray pulsars, rotation-powered pulsars, and magnetars. This paper explores scattering transport in the classical magnetic Thomson domain that is of broad applicability to these neutron star classes. The development of a Monte Carlo simulation for the polarized radiative transfer is detailed: it employs an electric field vector formalism to enable a breadth of utility in relating linear, circular, and elliptical polarizations. The simulation can be applied to any neutron star surface locale, and is adaptable to accretion column and magnetospheric problems. Validation of the code for both intensity and Stokes parameter determination is illustrated in a variety of ways. Representative results for emergent polarization signals from surface layers are presented for both polar and equatorial magnetic locales, exhibiting contrasting signatures between the two regions. There is also a strong dependence of these characteristics on the ratio of the frequency ωof a photon to the cyclotron frequency ω_B=eB/mc. Polarization signatures for high-opacity domains are presented, highlighting compact analytical approximations for the Stokes parameters and anisotropy relative to the local field direction for an extended range of frequencies. These are very useful in defining injection conditions deep in the simulation slab geometries, expediting the generation of emission signals from highly opaque stellar atmospheres. The results are interpreted throughout using the polarization characteristics of the magnetic Thomson differential cross-section.