2019ApJ...876..109Z

The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS 0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multiwavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and neutrino flux of hadronic models, and predict their X-ray to γ-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multiwavelength polarization variability of the proton synchrotron (PS) model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and γ-ray bands, and to infer the magnetic field strength in the emission region. In addition, the temporal multiwavelength polarization signatures in the PS shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multiwavelength polarimetry facilities such as the Imaging X-ray Polarimetry Explorer and All-sky Medium Energy Gamma-ray Observatory together with neutrino telescopes such as IceCube can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.