Unveiling the graviton mass bounds through the analysis of 2023 pulsar timing array data releases

Strong evidence for the Helling-Downs correlation curves have been reported by multiple pulsar timing array (PTA) collaborations in the middle of 2023. In this work, we investigate the graviton mass bounds by analyzing the observational data of the overlap reduction functions from the NANOGrav 15-year data release and CPTA first data release. The results from our data analysis display the state-of-the-art upper limits on the graviton mass at 90% confidence level, namely, mg8.6×10^-24 eV from NANOGrav and mg3.8×10^-23 eV from CPTA. We also study the cosmic-variance limit on the graviton mass bounds, i.e., ${{σ}}_{{m}g}CV=4.8×{}{10}^{{-}24} eV×{}f/(10 year{)}^{{-}1}$, with $f$ being a typical frequency band of PTA observations. This is equivalent to the cosmic-variance limit on the speed of gravitational waves, i.e., ${{σ}}_{{v}g}CV=0.07c$, with $c$ being the speed of light. Moreover, we discuss potential implications of these results for scenarios of ultralight tensor dark matter.