Constraining Axion-Like Particle mediated Dark Matter with Observational Constraints: A Statistical and Machine Learning Approach
We present a comprehensive study of axion-like particle (ALP) mediated dark matter (DM) effects on neutron star (NS) structure within a relativistic mean-field framework with non-linear mesonic interactions constrained by nuclear and astrophysical data. We explore DM masses and Fermi momenta , generating over 30{,}000 equations of state using two representative hadronic models, a stiff EoS (EoS1) and a soft EoS (EoS18), including a consistent crust description. A multi-level statistical filtering scheme based on voting, likelihood, and kernel density estimation is applied using constraints from radio and X-ray pulsars, GW170817, and the low-mass compact object HESS~J1731347. We find that models satisfying the PSR~J06143329 radius constraint automatically comply with the HESS bound, allowing ALP-mediated DM to explain low-mass compact objects while remaining consistent with NSs. For the stiff EoS, we obtain a lower bound , with preferred values and , while the soft EoS yields no strict lower bound, though large and are disfavored. We also develop a supervised interpolation model using \texttt{AutoGluon} to infer DM parameters from NS mass--radius curves, achieving , and show that is mainly constrained by global radius ratios, whereas is driven by the tidal deformability .