On the statistical characterization of the synchrotron multi-zone polarization of blazars

Multiwavelength polarimetric observations of blazars reveal complex, energy-dependent polarization behavior, including a decrease in polarization fraction from X-rays to millimeter bands and significant variability in the electric vector position angle (EVPA). These trends challenge simple single-zone synchrotron models and suggest a more intricate, turbulent jet structure with multiple emission zones. We develop a statistical framework to model the observed energy-dependent polarization patterns in blazars, focusing on the behavior captured by IXPE in the X-ray band and RoboPol in the optical. The goal is to statistically characterize multi-zone models in terms of the distributions of cell size and the physical parameters of the electron energy distribution (EED). A Monte Carlo approach, implemented with the JetSeT code, is used to generate synthetic multi-zone synchrotron emission from a spherical region filled with turbulent cells with randomly distributed physical properties. Simulations explore scenarios ranging from identical cells to power-law distributions of cell sizes and EED parameters with variable cutoff and low-energy slopes. The results show that a purely turbulent, multi-zone model can reproduce the observed energy-dependent polarization without requiring correlations between cell size and EED parameters. The polarization degree is primarily determined by the effective, flux-weighted, number of emitting cells, modulated by the dispersion in cell properties, particularly the EED cutoff energy at high frequencies and the low-energy spectral index at low frequencies. With a fractional dispersion in cutoff energy of about 90% and a low-energy spectral index dispersion of ~0.5-1.5, the model reproduces the chromatic mm-to-X-ray polarization trends seen by IXPE and the optical polarization limiting envelope observed in the RoboPol dataset.

💡 Research Summary

The paper addresses the puzzling energy‑dependent polarization behavior observed in blazars, namely the decrease of the polarization fraction from X‑rays to millimeter wavelengths and the large variability of the electric‑vector position angle (EVPA). Single‑zone synchrotron models cannot reproduce these trends, prompting the authors to develop a statistical multi‑zone framework that treats the jet as a collection of turbulent cells, each with its own magnetic field orientation, strength, Doppler factor, and electron energy distribution (EED).

Using the JetSeT code (v1.3.1) they compute the synchrotron spectral energy distribution (SED) and the Stokes parameters for each cell. The electron distribution is taken as a power‑law with an exponential cutoff, n(γ)=K γ⁻ᵖ exp(−γ/γ_cut). The cell radius follows a power‑law distribution (R_c∝r^q, q≤0) within a spherical region of radius R_S=10¹⁶ cm. Magnetic field strength B_c, Doppler factor δ_c, and the EED index p are either fixed or drawn from uniform (or log‑uniform) probability density functions; the magnetic field direction χ_r is random in