Axionlike particles (ALPs) may undergo mixing with photons in the presence of astrophysical magnetic fields, leading to alterations in the observed high energy $\gamma$-ray spectra. In this study, we investigate the ALP-photon oscillation effect using the spectra of the blazar Mrk 421 over 15 observation periods measured by Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) and Fermi Large Area Telescope (Fermi-LAT). Compared with previous studies, we generate mock data under the ALP hypothesis and employ the ${\mathrm{CL}}_{\mathrm{s}}$ method to set constraints on the ALP parameters. This method is widely utilized in high energy experiments and avoids the exclusion of specific parameter regions where distinguishing between the null and ALP hypotheses is challenging. We find that the ALP-photon coupling $g_{a\gamma }$ is constrained to be smaller than $\sim 2\times {10}^{-11}{\mathrm{GeV}}^{-1}$ for ALP masses ranging from ${10}^{-9}\mathrm{eV}$ to ${10}^{-7}\mathrm{eV}$ at the 95% confidence level. We also present the constraints derived from the TS distribution under the null hypothesis, which is commonly utilized in previous astrophysical ALP studies. Our results reveal that the combined constraints of all the periods obtained from both methods are consistent. However, the ${\mathrm{CL}}_{\mathrm{s}}$ method remains effective in cases where the latter method fails to provide constraints for specific observation periods.

• Received 11 September 2023
• Accepted 12 February 2024

DOI:https://doi.org/10.1103/PhysRevD.109.063003