Modeling the polarization properties of synchrotron sources through simulations of relativistic MHD turbulence

The emission from the relativistically hot plasmas of high-energy astrophysical synchrotron sources, Pulsar Wind Nebulae (PWNe) in particular, depends on the level of magnetic fluctuations. Recent observations by the X-ray polarimeter IXPE support the presence of turbulence, with varying conditions even in different regions of a same source. We aim at modeling such synchrotron emission, and in particular the degree of linear polarization, by using for the first time 3D relativistic MHD turbulence simulations. Thanks to a novel accelerated version of the ECHO code, a series of 3D relativistic MHD simulations are performed assuming a relativistically hot plasma and various degrees of magnetization, mimicking different conditions encountered in synchrotron sources. Magnetic fluctuations at random directions with respect to a background field are initialized at large scales. After the turbulent cascade is fully developed, the statistical properties of the plasma and of the synthetic synchrotron emission maps are analyzed. Turbulence rapidly relaxes to a sort of Alfvénic equilibrium and a Kolmogorov cascade with slope -5/3 soon develops, with differences depending on the initial ratio η of magnetic fluctuations over the background field. Dissipation mostly occurs in thin current sheets, where (numerical) reconnection takes place, intermittency and deviation from isotropic Gaussian distributions are observed. Synthetic synchrotron maps and their statistical properties depend on η too, approaching analytical estimates for large η. The integrated degree of linear polarization is found to cover the whole range of observed values in PWNe, and its dependence on the relative amplitude of turbulent fluctuations shows a good agreement with analytical estimates, even in the presence of anisotropy.