Probing Ultralight Axion-like Dark Matter $-$ A PTA-PPA Synergy

Abstract: Ultralight axion-like dark matter (ALDM) is a leading candidate in the dark matter realm, characterized by its prominent wave properties on astronomical scales. Pulsar Timing Arrays (PTAs) and Pulsar Polarization Arrays (PPAs) aim to detect this dark matter through timing and polarization measurements, respectively, of pulsars. The PTA relies on gravitational effects, as the ALDM halo perturbs the spacetime metric within the Milky Way, while the PPA detects non-gravitational effects, namely cosmological birefringence induced by the ALDM Chern-Simons coupling with photons. These two methods complement each other, synergistically enhancing the pulsar array's capability to identify the ALDM signals in the data. In this article, we provide a foundational development of this synergy. We begin by revisiting previously derived two-point correlation functions for both PTA and PPA, and expand our analysis to include correlations between timing and polarization signals. We then construct likelihood functions for PTA and combined PTA-PPA analyses within a Bayesian framework, aimed at detecting the characteristic correlations of ALDM signals. We emphasize the non-Gaussianity of the ALDM timing signals, which arises from their non-linear dependence on the field, in contrast to the Gaussian nature of its polarization signals. To address the complexities introduced, we approach this investigation in two ways: one involves a Gaussian approximation with proper justifications, while the other derives the formalism from the generic Gaussian characteristics of the ALDM field. We anticipate that these efforts will lead to further developments in PTA and PTA-PPA analysis methods.