Search

We have developed a set of four fully coupled Boltzmann equations to precisely determine the relic density and temperature of dark matter by including three distinct sectors: dark matter, light scalar, and standard model sectors. The intricacies of heat transfer between DM and the SM sector through a light scalar particle are explored, inspired by stringent experimental constraints on the scalar-Higgs mixing angle and the DM-scalar coupling. Three distinct sectors emerge prior to DM freeze-out, requiring fully coupled Boltzmann equations to accurately compute relic density. Investigation of forbidden, resonance, and secluded DM scenarios demonstrates significant deviations between established methods and the novel approach with fully coupled Boltzmann equations. Despite increased computational demands, this emphasizes the need for improved precision in relic density calculations, underlining the importance of incorporating these equations in comprehensive analyses., Comment: 38 pages, 10 figures

We investigate the elastic scattering cross section between dark matter and protons using the DES Year 3 weak lensing data. This scattering induces a dark acoustic oscillation structure in the matter power spectra. To address non-linear effects at low redshift, we utilize principal component analysis alongside a limited set of $N$-body simulations, improving the reliability of our matter power spectrum prediction. We further perform a robust Markov Chain Monte Carlo analysis to derive the upper bounds on the DM-proton elastic scattering cross-section, assuming different velocity dependencies. Our results, presented as the first Frequentist upper limits, are compared with the ones obtained by Bayesian approach. Compared with the upper limits derived from the Planck cosmic microwave background data, our findings from DES Year 3 data exhibit improvements of up to a factor of five. In addition, we forecast the future sensitivities of the China Space Station Telescope, the upcoming capabilities of this telescope could improve the current limits by approximately one order of magnitude., Comment: 10 pages, 5 figures

This study explores a minimal renormalizable dark matter (DM) model, incorporating a sub-GeV Majorana DM and a singlet scalar particle $\phi$. Using scalar and pseudo-scalar interactions (couplings $c_s$ and $c_p$), we investigate implications for DM detection, considering $s$-wave, $p$-wave, and combined ($s$+$p$ wave) contributions in DM annihilation cross-section, as well as loop-correction contributions to DM-nucleon elastic scattering. Identifying a broad parameter space ($10 \,\rm{MeV} < m_\chi \lesssim m_\phi$) within the $2\sigma$ allowed region, we explore scenarios ($\left|c_s\right|\gg \left|c_p\right|$, $\left|c_s\right|\ll \left|c_p\right|$, and $\left|c_s\right|\approx \left|c_p\right|$). We find that (i) a non-zero pseudo-scalar coupling alleviates direct detection constraints as a comparison with the previous pure scalar coupling case; (ii) CMB observations set stringent limits on pseudo-scalar interaction dominant cases, making $s$-wave annihilation viable only for $m_\chi>1\,\rm{GeV}$; (iii) the preferred $\phi$-resonance region can be tested in the future indirect detection experiments, such as e-ASTROGAM., Comment: 34 pages, 4 figures

We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey and the Planck CMB and BAO data. In the early Universe, the mirror baryon and mirror radiation behave as dark matter and dark radiation, and their presence modifies the Universe's expansion history. Additionally, the scattering between mirror baryon and photon generates the dark acoustic oscillation process, suppressing the matter power spectrum from the cosmic shear measurement. We demonstrate how current data constrain these corrections to the $\Lambda$CDM cosmology and find that for a viable solution to the little hierarchy problem, the proportion of MTH dark matter cannot exceed about $30\%$ of the total dark matter density, unless the temperature of twin photon is less than $30\%$ of that of the standard model photon. While the MTH model is presently not a superior solution to the observed $H_0$ tension compared to the $\Lambda$CDM+$\Delta N_{\rm eff}$ model, we demonstrate that it has the potential to alleviate both the $H_0$ and $S_8$ tensions, especially if the $S_8$ tension persists in the future and approaches the result reported by the Planck SZ (2013) analysis. In this case, the MTH model can relax the tensions while satisfying the DES power spectrum constraint up to $k \lesssim 10~h\rm {Mpc}^{-1}$. If the MTH model is indeed accountable for the $S_8$ and $H_0$ tensions, we show that the future China Space Station Telescope (CSST) can determine the twin baryon abundance with a $10\%$ level precision., Comment: 32 pages, 12 figures, 4 tables. Published in JCAP

Several Pulsar Timing Array (PTA) collaborations have recently provided strong evidence for a nHz Stochastic Gravitational-Wave Background (SGWB). Here we investigate the implications of a first-order phase transition occurring within the early universe's dark quantum chromodynamics (dQCD) epoch, specifically within the framework of the mirror twin Higgs dark sector model. Our analysis indicates a distinguishable SGWB signal originating from this phase transition, which can explain the measurements obtained by PTAs. Remarkably, a significant portion of the parameter space for the SGWB signal also effectively resolves the existing tensions in both the $H_0$ and $S_8$ measurements in Cosmology. This intriguing correlation suggests a possible common origin of these three phenomena for $0.2 < \Delta N_{\rm eff} < 0.5$, where the mirror dark matter component constitutes about $30\%$ of the total dark matter abundance. Next generation CMB experiment such as CMB-S4 is able to test this parameter region., Comment: 7 pages, 4 figures

The origin of interstellar and intergalactic magnetic fields remains largely unknown. One possibility is that they are related to the primordial magnetic fields (PMFs) produced by, for instance, the phase transitions of the early Universe. In this letter, we show that the PMF-induced turbulence generated at around the QCD phase transition epoch -- the characteristic magnetic field strength $B_{\rm ch}^* \sim \mathcal{O}(1)~\rm{\mu G}$ and coherent length scale $\ell_{\rm ch}^* \sim \mathcal{O}(1)~\rm{pc}$, can naturally accommodate with nanohertz gravitational waves reported by pulsar timing array collaborations. Moreover, the evolution of the PMFs to the recombination era with the form of $B_{\rm ch}\sim \ell_{\rm ch}^{-\alpha}$ can induce baryon density inhomogeneities, alter the recombination history, and alleviate the tension of the Hubble parameter $H_0$ and the matter clumpiness parameter $S_8$ between early and late-time measurements for $0.88\leq \alpha \leq 1.17$ (approximate 95\% credible region based on three PTA likelihoods). The further evolved PMFs may account for the $\sim {\cal O}(10^{-16})$ Gauss extra-galactic magnetic field inferred with GRB 221009A., Comment: 7 pages, 3 figures, 1 table

We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey and the Planck CMB and BAO data. In the early Universe, the mirror baryon and mirror radiation behave as dark matter and dark radiation, and their presence modifies the Universe's expansion history. Additionally, the scattering between mirror baryon and photon generates the dark acoustic oscillation process, suppressing the matter power spectrum from the cosmic shear measurement. We demonstrate how current data constrain these corrections to the $\Lambda$CDM cosmology and find that for a viable solution to the little hierarchy problem, the proportion of MTH dark matter cannot exceed about $30\%$ of the total dark matter density, unless the temperature of twin photon is less than $30\%$ of that of the standard model photon. While the MTH model is presently not a superior solution to the observed $H_0$ tension compared to the $\Lambda$CDM+$\Delta N_{\rm eff}$ model, we demonstrate that it has the potential to alleviate both the $H_0$ and $S_8$ tensions, especially if the $S_8$ tension persists in the future and approaches the result reported by the Planck SZ (2013) analysis. In this case, the MTH model can relax the tensions while satisfying the DES power spectrum constraint up to $k \lesssim 10~h\rm {Mpc}^{-1}$. If the MTH model is indeed accountable for the $S_8$ and $H_0$ tensions, we show that the future China Space Station Telescope (CSST) can determine the twin baryon abundance with a $10\%$ level precision., Comment: 32 pages, 12 figures, 4 tables

New long-lived particles produced at the colliders may escape from conventional particle detectors. Using satellites or ground telescopes, we can detect the photons generated from the annihilation of the star-captured dark matter into a pair of long-lived particles. When the propagation length of these long-lived particles surpasses the interplanetary distance between the Sun and Jupiter, it becomes unfeasible to detect such dark matter signals originating from the Sun or Jupiter on Earth. Our analysis of the dark matter-induced photons produced by prompt radiation, inverse Compton scattering, and synchrotron radiation mechanisms reveals that a decay length of about $10^{-3}$ pc for long-lived particles is required for maximum detectability. We investigate the parameters that allow the long-lived particle's lifetime to be consistent with Big Bang nucleosynthesis while also allowing it to escape the confines of our solar system. The Sirius system is proposed as a promising target for the indirect detection of such long-lived particles. Utilizing the prompt, inverse Compton scattering, and synchrotron radiation, upper limits on the dark matter-proton spin-independent and spin-dependent cross section are estimated with the Fermi-LAT null-signal observation and the capabilities of the upcoming Square Kilometre Array radio telescope., Comment: 30 pages, 6 figures, JCAP accepted

We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR Boron and Carbon. We derive new CR spectra from DM-CR collisions by incorporating these DM-CR cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. Utilizing the latest data from AMS-02 and DAMPE on the Boron-to-Carbon ratio, we estimate a 95\% upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass. Our findings reveal that at $m_\chi \simeq 2 \,\rm{MeV}$, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$., Comment: 25 pages, 8 figures

Several Pulsar Timing Array (PTA) collaborations have recently provided strong evidence for a nHz Stochastic Gravitational-Wave Background (SGWB). Here we investigate the implications of a first-order phase transition occurring within the early universe's dark quantum chromodynamics (dQCD) epoch, specifically within the framework of the mirror twin Higgs dark sector model. Our analysis indicates a distinguishable SGWB signal originating from this phase transition, which can explain the measurements obtained by PTAs. Remarkably, a significant portion of the parameter space for the SGWB signal also effectively resolves the existing tensions in both the $H_0$ and $S_8$ measurements in Cosmology. This intriguing correlation suggests a possible common origin of these three phenomena for $0.2 < \Delta N_{\rm eff} < 0.5$, where the mirror dark matter component constitutes about $30\%$ of the total dark matter abundance. Next generation CMB experiment such as CMB-S4 is able to test this parameter region., Comment: 7 pages, 4 figures

This study investigates the antineutrinos production by $\beta$-decay of $r$-process nuclei in two astrophysical sites that are capable of producing gamma-ray bursts (GRBs): binary neutron star mergers (BNSMs) and collapsars, which are promising sites for heavy element nucleosynthesis. We employ a simplified method to compute the $\beta$-decay $\bar\nu_e$ energy spectrum and consider a number of different representative thermodynamic trajectories for $r$-process simulations, each with four sets of $Y_e$ distribution. The time evolution of the $\bar\nu_e$ spectrum is derived for both the dynamical ejecta and the disk wind for BNSMs and collapsar outflow, based on approximated mass outflow rates. Our results show that the $\bar\nu_e$ has an average energy of approximately 3 to 9~MeV, with a high energy tail of up to 20 MeV. The $\bar\nu_e$ flux evolution is primarily determined by the outflow duration, and can thus remain large for $\mathcal{O}(10)$~s and $\mathcal{O}(100)$~s for BNSMs and collapsars, respectively. For a single merger or collapsar at 40~Mpc, the $\bar\nu_e$ flux is $\mathcal{O}(10-100)$~cm$^{-2}$~s$^{-1}$, indicating a possible detection horizon up to $0.1-1$~Mpc for Hyper-Kamiokande. We also estimate their contributions to the diffuse $\bar\nu_e$ background, and find that both sources should only contribute subdominantly to the diffuse background when compared to that expected from core-collapse supernovae., Comment: 14 pages, 8 figures

We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey and the Planck CMB and BAO data. In the early Universe, the mirror baryon and mirror radiation behave as dark matter and dark radiation, and their presence modifies the Universe's expansion history. Additionally, the scattering between mirror baryon and photon generates the dark acoustic oscillation process, suppressing the matter power spectrum from the cosmic shear measurement. We demonstrate how current data constrain these corrections to the $\Lambda$CDM cosmology and find that for a viable solution to the little hierarchy problem, the proportion of MTH dark matter cannot exceed about $30\%$ of the total dark matter density, unless the temperature of twin photon is less than $30\%$ of that of the standard model photon. While the MTH model is presently not a superior solution to the observed $H_0$ tension compared to the $\Lambda$CDM+$\Delta N_{\rm eff}$ model, we demonstrate that it has the potential to alleviate both the $H_0$ and $S_8$ tensions, especially if the $S_8$ tension persists in the future and approaches the result reported by the Planck SZ (2013) analysis. In this case, the MTH model can relax the tensions while satisfying the DES power spectrum constraint up to $k \lesssim 10~h\rm {Mpc}^{-1}$. If the MTH model is indeed accountable for the $S_8$ and $H_0$ tensions, we show that the future China Space Station Telescope (CSST) can determine the twin baryon abundance with a $10\%$ level precision., Comment: 32 pages, 12 figures, 4 tables. Published in JCAP

Precise measurements of the stellar orbits around Sagittarius A* have established the existence of a supermassive black hole (SMBH) at the Galactic centre (GC). Due to the interplay between the SMBH and dark matter (DM), the DM density profile in the innermost region of the Galaxy, which is crucial for the DM indirect detection, is still an open question. Among the most popular models in the literature, the theoretical spike profile proposed by Gondolo and Silk (1999; GS hereafter) is well adopted. In this work, we investigate the DM spike profile using updated data from the Keck and VLT telescopes considering that the presence of such an extended mass component may affect the orbits of the S-stars in the Galactic center. We examine the radius and slope of the generalized NFW spike profile, analyze the Einasto spike, and discuss the influence of DM annihilation on the results. Our findings indicate that an initial slope of $\gamma \gtrsim 0.92$ for the generalized NFW spike profile is ruled out at a 95% confidence level. Additionally, the spike radius $R_{\rm sp}$ larger than 21.5 pc is rejected at 95% probability for the Einasto spike with $\alpha=0.17$, which also contradicts the GS spike model. The constraints with the VLT/GRAVITY upper limits are also projected. Although the GS NFW spike is well constrained by the Keck and VLT observation of S2, an NFW spike with a weak annihilation cusp may still be viable, as long as the DM annihilation cross section satisfies $\left< \sigma v \right> \gtrsim 7.7\times 10^{-27}~{\rm cm^3\,s^{-1}} (m_{\rm DM}/100~{\rm GeV})$ at 95% level., Comment: 12 pages, 10 figures and 1 table. Accepted for publication in MNRAS

New long-lived particles produced at the colliders may escape from conventional particle detectors. Using satellites or ground telescopes, we can detect the photons generated from the annihilation of the star-captured dark matter into a pair of long-lived particles. When the propagation length of these long-lived particles surpasses the interplanetary distance between the Sun and Jupiter, it becomes unfeasible to detect such dark matter signals originating from the Sun or Jupiter on Earth. Our analysis of the dark matter-induced photons produced by prompt radiation, inverse Compton scattering, and synchrotron radiation mechanisms reveals that a decay length of about $10^{-3}$ pc for long-lived particles is required for maximum detectability. We investigate the parameters that allow the long-lived particle's lifetime to be consistent with Big Bang nucleosynthesis while also allowing it to escape the confines of our solar system. The Sirius system is proposed as a promising target for the indirect detection of such long-lived particles. Utilizing the prompt, inverse Compton scattering, and synchrotron radiation, upper limits on the dark matter-proton spin-independent and spin-dependent cross section are estimated with the Fermi-LAT null-signal observation and the capabilities of the upcoming Square Kilometre Array radio telescope., Comment: 30 pages, 6 figures, JCAP accepted

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at Zenodo.

The $W$-boson mass, which was recently measured at FermiLab, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02., Comment: 8 pages, 5 figures

Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon anomalous magnetic moment ($g$-2) and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider data, B-physics, dark matter relic density and direct detection experiments, we find that a neutralino in the mass range of $\sim 160-270$ GeV is a viable solution. Moreover, the favored parameter region can be effectively probed by the ongoing direct detection experiments like LZ, PandaX-4T and XENON-nT. The velocity averaged annihilation cross section of the dark matter particles, however, is suppressed., Comment: 14 pages, 6 figures

Ultralight scalar dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with ordinary matter may be detectable. In this work, we search for possible evidence of the scalar cloud using accurate orbital measurements of the S2 star around Sgr A*. We solve the first order Post-Newtonian equation, considering simultaneously the extended mass distribution of the scalar cloud and the additional coupling via Higgs-portal or photon-portal interaction. We find that the astrometric and spectroscopic data of the S2 star are well consistent with the scenario of point-like mass of Sgr A*. We thus derive upper limits of the extended mass and the coupling of the new interaction. The limits of the Higgs coupling and extended mass of the scalar cloud are the most stringent ones for the scalar mass window between $3.2\times 10^{-19}$ eV and $1.6\times 10^{-18}$ eV., Comment: 11 pages, 6 figures, comments welcome

Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon anomalous magnetic moment ($g-2$) and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider data, B-physics, dark matter relic density and direct detection experiments, we find that a neutralino in the mass range of $\sim 150-290$ GeV is a viable solution. Moreover, the favored parameter region can be effectively probed by the ongoing direct detection experiments like LZ, PandaX-4T and XENON-nT. The velocity averaged annihilation cross section of the dark matter particles, however, is suppressed., Comment: 17 pages, 4 figures

The $W$-boson mass, which was recently measured at FermiLab, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02., Comment: 8 pages, 5 figures

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at Zenodo.

Ultralight bosonic dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with the ordinary matter may be detectable. In this work, we search for possible evidence of the scalar cloud using accurate orbital measurements of the S2 star around Sgr~A*. We solve the first order Post-Newtonian equation, considering simultaneously the extended mass distribution of the scalar cloud and the frequency shift induced by the additional coupling via Higgs portal or photon portal interaction. Furthermore, we also investigate the impact of an astrophysical power-law component from the gas and stellar remnants. We find that the astrometric and spectroscopic data of the S2 star are well consistent with the scenario of a point-like mass of Sgr~A*. We thus derive upper limits of the coupling of the new interaction and the extended mass, with and without the contribution from the astrophysical component. The limits of the Higgs/photon coupling and the extended mass of the scalar cloud are the most stringent ones for the scalar mass window between $3.2\times 10^{-19}$~eV and $1.6\times 10^{-18}$~eV., Comment: 12 pages, 6 figures, published at PRD

Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon $g-2$ and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider data, flavor physics, dark matter relic density, and direct detection experiments, we find that lighter electroweakinos and sleptons can generate sufficient contributions to muon $g-2$ and $m_W$. Moreover, the corresponding bino-like neutralino dark matter mass is in the $\sim 180-280$ GeV range. Interestingly, the favored DM mass region can soon be entirely probed by ongoing direct detection experiments like PandaX-4T, XENONnT, LUX-ZEPLIN, and DARWIN., Comment: 21 pages, 7 figures

The $W$-boson mass, which was recently measured at FermiLab with an unprecedented precision, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02., Comment: 9 pages, 7 figures, Replaced to match published version in PRL

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma ray data, and so on. It has been claimed that only the DM annihilation to a muon-pair, namely muonphilic DM, is compatible with the null detection of all the corresponding astrophysical measurements [Di Mauro and Winkle (2021)]. On the other hand, a muonphilic DM model may also lead to a signal in the recent muon g-2 measurement or the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including different DM and mediator spins. In agreement with GCE (not only $2\mu$ but also $4\mu$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the mediator resonance region can explain GCE and relic density simultaneously. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon g-2 excess on top of GCE, relic density, and other DM and mediator constraints., Comment: 30 pages, 8 figures

The null signal from collider and dark matter (DM) direct detector experiments makes the interaction between DM and visible matter too small to reproduce the correct relic density for many thermal DM models. The remaining parameter space indicates that two almost degenerated states in the dark sector, the inelastic DM scenario, can co-annihilate in the early universe to produce the correct relic density. Regarding the direct detection of the inelastic DM scenario, the virialized DM component from the nearby halo is nonrelativistic and not able to excite the DM ground state, even if the relevant couplings can be considerable. Thus, a DM with a large mass splitting can evade traditional virialized DM direct detection. In this study, we connect the concept of cosmic-ray accelerated DM in our Milky Way and the direct detection of inelastic scattering in underground detectors to explore spectra that result from several interaction types of the inelastic DM. We find that the mass splitting $\delta<\mathcal{O}(1~{\rm MeV})$ can still be reachable for cosmic ray accelerated DM with mass range $1~{\rm MeV}

Conventionally, one can constrain the dark matter (DM) interaction with DM mass heavier than GeV by searching for DM induced synchrotron emission in the radio frequency band. However, an MeV DM can also generate detectable radio emission if electrons and positrons produced by DM annihilation or decay undergoes inverse Compton scattering (ICS) with the cosmic microwave background. The upcoming radio telescope Square Kilometre Array (SKA) is designed to operate with extremely high sensitivity. We investigate the capability of the SKA to detect DM particles in a board mass range from MeV to TeV, for both annihilation and decay scenarios. In this paper, we consider the sensitivities of the future SKA first and second phase (SKA1 and SKA2). As a comprehensive study, we systematically study the impacts on the DM-induced signal computation from the magnetic field strengths and particle diffusion coefficients. We compare the detection potential of four very different sources: two dwarf spheroidal galaxies (Draco and Segue 1), one radio-poor cluster (A2199), and one DM-rich ultra-diffuse galaxy (Dragonfly 44). We project the SKA1 and SKA2 sensitivities with the exposure of 100 hours on the annihilation cross section and decay time for DM mass from MeV to TeV by considering two different leptonic final states $e^+ e^-$ and $\mu^+\mu^-$., Comment: 24 pages, 6 figures. Replaced to match version accepted for publication in JCAP, the SKA sensitivity updated

Conventionally, one can constrain the dark matter (DM) interaction with DM mass heavier than GeV by searching for DM induced synchrotron emission in the radio frequency band. However, an MeV DM can also generate detectable radio emission if electrons and positrons produced by DM annihilation or decay undergoes inverse Compton scattering (ICS) with the cosmic microwave background. The upcoming radio telescope Square Kilometre Array (SKA) is designed to operate with extremely high sensitivity. We investigate the capability of the SKA to detect DM particles in a board mass range from MeV to TeV, for both annihilation and decay scenarios. In this paper, we consider the sensitivities of the future SKA first and second phase (SKA1 and SKA2). As a comprehensive study, we systematically study the impacts on the DM-induced signal computation from the magnetic field strengths and particle diffusion coefficients. We compare the detection potential of four very different sources: two dwarf spheroidal galaxies (Draco and Segue 1), one radio-poor cluster (A2199), and one DM-rich ultra-diffuse galaxy (Dragonfly 44). We project the SKA1 and SKA2 sensitivities with the exposure of 100 hours on the annihilation cross section and decay time for DM mass from MeV to TeV by considering two different leptonic final states $e^+ e^-$ and $\mu^+\mu^-$., Comment: 24 pages, 6 figures. Replaced to match version accepted for publication in JCAP, the SKA sensitivity updated

The supersymmetric model is one of the most attractive extensions of the Standard Model of particle physics. In light of the most recently reported anomaly of the muon g-2 measurement by the FermiLab E989 experiment, and the excesses of gamma rays at the Galactic center observed by Fermi-LAT space telescope, as well as the antiproton excess observed by the Alpha Magnetic Spectrometer, we propose to account for all these anomalies or excesses in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental constraints including the Higgs mass, B-physics, collider data, dark matter relic density and direct detections, we find that a $\sim 60$ GeV bino-like neutralino is able to successfully explain all these observations. Our scenario can be sensitively probed by future direct detection experiments., Comment: 5 pages, 3 figures, accepted version

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2\mu$ but also $4\mu$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints., Comment: 34 pages, 8 figures, Matches JHEP accepted version. Major revisions to address referee's comments

The null signal from collider and dark matter (DM) direct detector experiments makes the interaction between DM and visible matter too small to reproduce the correct relic density for many thermal DM models. The remaining parameter space indicates that two almost degenerated states in the dark sector, the inelastic DM scenario, can co-annihilate in the early universe to produce the correct relic density. Regarding the direct detection of the inelastic DM scenario, the virialized DM component from the nearby halo is nonrelativistic and not able to excite the DM ground state, even if the relevant couplings can be considerable. Thus, a DM with a large mass splitting can evade traditional virialized DM direct detection. In this study, we connect the concept of cosmic-ray accelerated DM in our Milky Way and the direct detection of inelastic scattering in underground detectors to explore spectra that result from several interaction types of the inelastic DM. We find that the mass splitting $\delta<\mathcal{O}(1~{\rm MeV})$ can still be reachable for cosmic ray accelerated DM with mass range $1~{\rm MeV}

The null signal from collider and dark matter (DM) direct detector experiments makes the interaction between DM and visible matter too small to reproduce the correct relic density for many thermal DM models. The remaining parameter space indicates that two almost degenerated states in the dark sector, the inelastic DM scenario, can co-annihilate in the early universe to produce the correct relic density. Regarding the direct detection of the inelastic DM scenario, the virialized DM component from the nearby halo is nonrelativistic and not able to excite the DM ground state, even if the relevant couplings can be considerable. Thus, a DM with a large mass splitting can evade traditional virialized DM direct detection. In this study, we connect the concept of cosmic-ray accelerated DM in our Milky Way and the direct detection of inelastic scattering in underground detectors to explore spectra that result from several interaction types of the inelastic DM. We find that the mass splitting $\delta<\mathcal{O}(1~{\rm MeV})$ can still be reachable for cosmic ray accelerated DM with mass range $1~{\rm MeV}

Recently, we have witnessed two hints of physics beyond the standard model: A 3.3$\sigma$ local excess ($M_{A_0} = 52$ GeV) in the search for $H_0\to A_0 A_0\to b\bar{b}\mu^{+}\mu^{-}$ and a 4.2$\sigma$ deviation from the SM prediction in the $(g-2)_\mu$ measurement. The first excess was found by the ATLAS Collaboration using 139 fb$^{-1}$ data at $\sqrt{s}=13$ TeV. The second deviation is a combination of the results from the Brookhaven E821 and the recently reported Fermilab E989 experiment. We attempt to explain these deviations in terms of a renormalizable simplified dark matter model. Inspired by the null signal result from dark matter (DM) direct detection, we interpret the possible new particle, $A_0$, as a pseudoscalar mediator connecting DM and the standard model. On the other hand, a new vector-like muon lepton can explain these two excesses at the same time while contributiong to the DM phenomenology., Comment: 28 pages, 8 figures. Recasting of SR5L added and figures updated accordingly. More references added. v3: Matches published version

We show that the impact of energy injection by dark matter annihilation on the cosmic microwave background power spectra can be apprehended via a residual likelihood map. By resorting to convolutional neural networks that can fully discover the underlying pattern of the map, we propose a novel way of constraining dark matter annihilation based on the Planck 2018 data. We demonstrate that the trained neural network can efficiently predict the likelihood and accurately place bounds on the annihilation cross-section in a $\textit{model-independent}$ fashion. The machinery will be made public in the near future., Comment: 25 pages, 8 figures; to match the published version

The null signal from collider and dark matter (DM) direct detector experiments makes the interaction between DM and visible matter too small to reproduce the correct relic density for many thermal DM models. The remaining parameter space indicates that two almost degenerated states in the dark sector, the inelastic DM scenario, can co-annihilate in the early universe to produce the correct relic density. Regarding the direct detection of the inelastic DM scenario, the virialized DM component from the nearby halo is nonrelativistic and not able to excite the DM ground state, even if the relevant couplings can be considerable. Thus, a DM with a large mass splitting can evade traditional virialized DM direct detection. In this study, we connect the concept of cosmic-ray accelerated DM in our Milky Way and the direct detection of inelastic scattering in underground detectors to explore spectra that result from several interaction types of the inelastic DM. We find that the mass splitting $\delta<\mathcal{O}(1~{\rm MeV})$ can still be reachable for cosmic ray accelerated DM with mass range $1~{\rm MeV}

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2\mu$ but also $4\mu$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints., Comment: 34 pages, 8 figures, Matches JHEP accepted version. Major revisions to address referee's comments

The null signal from collider and dark matter (DM) direct detector experiments makes the interaction between DM and visible matter too small to reproduce the correct relic density for many thermal DM models. The remaining parameter space indicates that two almost degenerated states in the dark sector, the inelastic DM scenario, can co-annihilate in the early universe to produce the correct relic density. Regarding the direct detection of the inelastic DM scenario, the virialized DM component from the nearby halo is nonrelativistic and not able to excite the DM ground state, even if the relevant couplings can be considerable. Thus, a DM with a large mass splitting can evade traditional virialized DM direct detection. In this study, we connect the concept of cosmic-ray accelerated DM in our Milky Way and the direct detection of inelastic scattering in underground detectors to explore spectra that result from several interaction types of the inelastic DM. We find that the mass splitting $\delta<\mathcal{O}(1~{\rm MeV})$ can still be reachable for cosmic ray accelerated DM with mass range $1~{\rm MeV}

Conventionally, one can constrain the dark matter (DM) interaction with DM mass heavier than GeV by searching for DM induced synchrotron emission in the radio frequency band. However, an MeV DM can also generate detectable radio emission if electrons and positrons produced by DM annihilation or decay undergoes inverse Compton scattering (ICS) with the cosmic microwave background. The upcoming radio telescope Square Kilometre Array (SKA) is designed to operate with extremely high sensitivity. We investigate the capability of the SKA to detect DM particles in a board mass range from MeV to TeV, for both annihilation and decay scenarios. In this paper, we consider the sensitivities of the future SKA first and second phase (SKA1 and SKA2). As a comprehensive study, we systematically study the impacts on the DM-induced signal computation from the magnetic field strengths and particle diffusion coefficients. We compare the detection potential of four very different sources: two dwarf spheroidal galaxies (Draco and Segue 1), one radio-poor cluster (A2199), and one DM-rich ultra-diffuse galaxy (Dragonfly 44). We project the SKA1 and SKA2 sensitivities with the exposure of 100 hours on the annihilation cross section and decay time for DM mass from MeV to TeV by considering two different leptonic final states $e^+ e^-$ and $\mu^+\mu^-$., Comment: 24 pages, 6 figures. Replaced to match version accepted for publication in JCAP, the SKA sensitivity updated

The supersymmetric model is one of the most attractive extensions of the Standard Model of particle physics. In light of the most recently reported anomaly of the muon g-2 measurement by the FermiLab E989 experiment, and the excesses of gamma rays at the Galactic center observed by Fermi-LAT space telescope, as well as the antiproton excess observed by the Alpha Magnetic Spectrometer, we propose to account for all these anomalies or excesses in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental constraints including the Higgs mass, B-physics, collider data, dark matter relic density and direct detections, we find that a $\sim 60$ GeV bino-like neutralino is able to successfully explain all these observations. Our scenario can be sensitively probed by future direct detection experiments., Comment: 5 pages, 3 figures, accepted version

We show that the impact of energy injection by dark matter annihilation on the cosmic microwave background power spectra can be apprehended via a residual likelihood map. By resorting to convolutional neural networks that can fully discover the underlying pattern of the map, we propose a novel way of constraining dark matter annihilation based on the Planck 2018 data. We demonstrate that the trained neural network can efficiently predict the likelihood and accurately place bounds on the annihilation cross-section in a $\textit{model-independent}$ fashion. The machinery will be made public in the near future., Comment: 25 pages, 8 figures; to match the published version

Conventionally, one can constrain the dark matter (DM) interaction with DM mass heavier than GeV by searching for DM induced synchrotron emission in the radio frequency band. However, an MeV DM can also generate detectable radio emission if electrons and positrons produced by DM annihilation or decay undergoes inverse Compton scattering (ICS) with the cosmic microwave background. The upcoming radio telescope Square Kilometre Array (SKA) is designed to operate with extremely high sensitivity. We investigate the capability of the SKA to detect DM particles in a board mass range from MeV to TeV, for both annihilation and decay scenarios. In this paper, we consider the sensitivities of the future SKA first and second phase (SKA1 and SKA2). As a comprehensive study, we systematically study the impacts on the DM-induced signal computation from the magnetic field strengths and particle diffusion coefficients. We compare the detection potential of four very different sources: two dwarf spheroidal galaxies (Draco and Segue 1), one radio-poor cluster (A2199), and one DM-rich ultra-diffuse galaxy (Dragonfly 44). We project the SKA1 and SKA2 sensitivities with the exposure of 100 hours on the annihilation cross section and decay time for DM mass from MeV to TeV by considering two different leptonic final states $e^+ e^-$ and $\mu^+\mu^-$., Comment: 24 pages, 6 figures. Replaced to match version accepted for publication in JCAP, the SKA sensitivity updated

The cold dark matter (CDM) candidate with weakly interacting massive particles can successfully explain the observed dark matter relic density in cosmic scale and the large-scale structure of the Universe. However, a number of observations at the satellite galaxy scale seem to be inconsistent with CDM simulation. This is known as the small-scale problem of CDM. In recent years, it has been demonstrated that self-interacting dark matter (SIDM) with a light mediator offers a reasonable explanation for the small-scale problem. We adopt a simple model with SIDM and focus on the effects of Sommerfeld enhancement. In this model, the dark matter candidate is a leptonic scalar particle with a light mediator. We have found several regions of the parameter space with proper masses and coupling strength generating a relic density that is consistent with the observed CDM relic density. Furthermore, this model satisfies the constraints of recent direct searches and indirect detection for dark matter as well as the effective number of neutrinos and the observed small-scale structure of the Universe. In addition, this model with the favored parameters can resolve the discrepancies between astrophysical observations and $N$-body simulations., Comment: REVTEX, 32 pages, 11 figures. To be published in Phys. Rev. D

The direct detection of particle dark matter through its scattering with nucleons is of fundamental importance to understand the nature of DM. In this work, we propose that the high-energy neutrino detectors like IceCube can be used to uniquely probe the DM-nucleon cross-section for high-energy DM of $\sim$ PeV, up-scattered by the high-energy cosmic rays. We derive for the first time strong constraints on the DM-nucleon cross-section down to $\sim 10^{-32}$ cm$^2$ at this energy scale for sub-GeV DM candidates. Such independent probe at energy scale far exceeding other existing direct detection experiments can therefore provide useful insights complementary to other searches., Comment: 16 pages, 6 figures. Published version

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283., Comment: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

We propose a machine learning method to investigate the propagation of cosmic rays based on the precisely measured spectra of the primary and secondary cosmic ray nuclei of Li, Be, B, C, and O from AMS-02, ACE, and Voyager-1. We train two convolutional neural networks. One network learns how to infer propagation and source parameters from the energy spectra of cosmic rays, and the other network, which is similar to the former, has the flexibility to learn from the data with added artificial fluctuations. Together with the simulated data generated by GALPROP, we find that both networks can properly invert the propagation process and infer the propagation and source parameters reasonably well. This approach can be much more efficient than the traditional Markov chain Monte Carlo fitting method for deriving the propagation parameters if users choose to update confidence intervals with new experimental data. Both of the trained networks are available at (https://github.com/alan200276/CR_ML)., Comment: 33 pages, 11 figures. Matches JCAP accepted version

We revisit the signatures from collisions of cosmic-rays on sub-GeV dark matter (DM) in the Milky Way. In addition to the upscattered DM component that can be probed by existing DM and neutrino experiments widely discussed, we examine the associated signals in $\gamma$-rays and neutrinos that span a wide energy range due to the inelastic scatterings. Assuming a simple vector portal DM model for illustration, we compute both the upscattered DM flux by cosmic-ray protons, and the resulting emission of secondary $\gamma$-rays and high-energy neutrinos from proton excitation, hadronization, and the subsequent meson decay. We derive limits on coupling constants in the vector portal model using data from the $\gamma$-ray and high-energy neutrino telescopes including Fermi, H.E.S.S. and IceCube. These limits are compared to those obtained by considering the upscattered DM signals at the low-energy DM/neutrino detectors XENON1T/MiniBooNE and the IceCube. For this particular model, the limits are set predominantly by non-detection of the upscattered DM events in XENON1T, for most of the DM mass range due to the large scattering cross section at low energies. Nevertheless, our study demonstrates that the $\gamma$-ray and neutrino signals, traditionally considered as indirect probes for DM annihilation and decay, can also be directly used to constrain the DM--nucleon interaction in complementary to the direct search experiments., Comment: 17 pages, 7 figures; Published version

The cold dark matter (CDM) candidate with weakly interacting massive particles can successfully explain the observed dark matter relic density in cosmic scale and the large-scale structure of the Universe. However, a number of observations at the satellite galaxy scale seem to be inconsistent with CDM simulation. This is known as the small-scale problem of CDM. In recent years, it has been demonstrated that self-interacting dark matter (SIDM) with a light mediator offers a reasonable explanation for the small-scale problem. We adopt a simple model with SIDM and focus on the effects of Sommerfeld enhancement. In this model, the dark matter candidate is a leptonic scalar particle with a light mediator. We have found several regions of the parameter space with proper masses and coupling strength generating a relic density that is consistent with the observed CDM relic density. Furthermore, this model satisfies the constraints of recent direct searches and indirect detection for dark matter as well as the effective number of neutrinos and the observed small-scale structure of the Universe. In addition, this model with the favored parameters can resolve the discrepancies between astrophysical observations and $N$-body simulations., Comment: REVTEX, 32 pages, 11 figures. To be published in Phys. Rev. D

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283., Comment: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

We propose a machine learning method to investigate the propagation of cosmic rays based on the precisely measured spectra of the primary and secondary cosmic ray nuclei of Li, Be, B, C, and O from AMS-02, ACE, and Voyager-1. We train two convolutional neural networks. One network learns how to infer propagation and source parameters from the energy spectra of cosmic rays, and the other network, which is similar to the former, has the flexibility to learn from the data with added artificial fluctuations. Together with the simulated data generated by GALPROP, we find that both networks can properly invert the propagation process and infer the propagation and source parameters reasonably well. This approach can be much more efficient than the traditional Markov chain Monte Carlo fitting method for deriving the propagation parameters if users choose to update confidence intervals with new experimental data. Both of the trained networks are available at (https://github.com/alan200276/CR_ML)., Comment: 33 pages, 11 figures. Matches JCAP accepted version

We revisit the signatures from collisions of cosmic-rays on sub-GeV dark matter (DM) in the Milky Way. In addition to the upscattered DM component that can be probed by existing DM and neutrino experiments widely discussed, we examine the associated signals in $\gamma$-rays and neutrinos that span a wide energy range due to the inelastic scatterings. Assuming a simple vector portal DM model for illustration, we compute both the upscattered DM flux by cosmic-ray protons, and the resulting emission of secondary $\gamma$-rays and high-energy neutrinos from proton excitation, hadronization, and the subsequent meson decay. We derive limits on coupling constants in the vector portal model using data from the $\gamma$-ray and high-energy neutrino telescopes including Fermi, H.E.S.S. and IceCube. These limits are compared to those obtained by considering the upscattered DM signals at the low-energy DM/neutrino detectors XENON1T/MiniBooNE and the IceCube. For this particular model, the limits are set predominantly by non-detection of the upscattered DM events in XENON1T, for most of the DM mass range due to the large scattering cross section at low energies. Nevertheless, our study demonstrates that the $\gamma$-ray and neutrino signals, traditionally considered as indirect probes for DM annihilation and decay, can also be directly used to constrain the DM--nucleon interaction in complementary to the direct search experiments., Comment: 17 pages, 7 figures; Published version