Pulsar Polarization Array Limits on Ultralight Axionlike Dark Matter

We conduct the first-ever pulsar polarization array (PPA) analysis to detect the ultralight axionlike dark matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes pulsar timing array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of ρ_{0}=0.4  GeV/cm^{3}, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely ≲10^{-13.5}-10^{-12.2}  GeV^{-1}, for the mass range of 10^{-22}-10^{-21}  eV. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted.