Search

We investigate the production of axion-like particles (ALPs) in stellar cores, where they interact with electromagnetic fields and electrons, with typical masses between $\mathcal O(0.1)$ and $\mathcal O(10)$ keV. These low-energy ALPs are gravitationally trapped in the orbits of stars and subsequently decay into two photons that we detect as monochromatic X-ray lines. We propose to search for these gravitationally trapped ALPs in the Alpha Centauri binary system, our closest stellar neighbor, using sensitive X-ray detectors like Chandra and eROSITA. Our search for ALP decay signals in the energy range of 0.2 keV to 10 keV yielded null results, thus establishing the most stringent limits on ALP interactions to date. Specifically, if ALPs are mainly produced by Compton or bremsstrahlung processes (ALP-electron coupling $g_{aee}$ being significant), we have improved the limits on the ALP-photon coupling $g_{a\gamma\gamma}$ by two to three orders of magnitude, in ALP mass range between 0.2 keV to 5 keV, compared to previous measurements, including those from GW170817, SN 2023ixf, and other sources., Comment: 8 pages, 3 figures

The precision measurements of the muon magnetic moment and the $W$ boson mass have sparked interest in the potential deviations from standard model (SM) predictions. While it may be premature to attribute any excesses in these precision measurements to new physics, they do offer a valuable indication of potential directions for physics beyond the SM. Additionally, the particle nature of dark matter (DM) remains a crucial enigma. Despite the absence of any definitive DM signal in direct detection and collider experiments, the Galactic Center GeV $\gamma$-ray excess and the AMS-02 antiproton ($\overline{p}$) excess could potentially offer hints related to the evidence of DM. Motivated by these observations, we propose a simple DM model that addresses all these issues. This model extends the SM by incorporating singlet and doublet Dirac fermion fields, along with a doublet complex scalar field. For the viable parameter regions in this model, we find that future upgrades of the Large Hadron Collider and DM direct detection experiments can only partially probe them, while future high-energy muon colliders hold promise for exploring the unexplored parameter space., Comment: 33 pages, 10 figures, 2 tables. Comments are welcome

This study reanalyzes the detection prospects of dark matter (DM) annihilation signals in the Galactic Center, focusing on velocity-dependent dynamics within a spike density near the supermassive black hole (Sgr~A$^{\star}$). We investigate three annihilation processes -- $p$-wave, resonance, and forbidden annihilation -- under semi-relativistic velocities, leveraging gamma-ray data from Fermi and DAMPE telescopes. Our analysis integrates a fermionic DM model with an electroweak axion-like particle (ALP) portal, exploring annihilation into two or four photons. Employing a comprehensive six-dimensional integration, we precisely calculate DM-induced gamma-ray fluxes near Sgr~A$^{\star}$, incorporating velocity and positional dependencies in the annihilation cross-section and photon yield spectra. Our findings highlight scenarios of resonance and forbidden annihilation, where the larger ALP-DM-DM coupling constant $C_{a\chi\chi}$ can affect spike density, potentially yielding detectable gamma-ray line spectra within Fermi and DAMPE energy resolution. We set upper limits for $C_{a\chi\chi}$ across these scenarios, offering insights into the detectability and spectral characteristics of DM annihilation signals from the Galactic Center., Comment: 24 pages, 13 figures

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 dark matter (DM) and the standard model 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: 40 pages, 10 figures. Matches PRD accepted version

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: 35 pages, 4 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: 24 pages, 11 figures, published in JCAP

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 their cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. In a specific example with a coupling strength of $b_{\chi}=0.1$ and a DM mass of $m_{\chi}=0.1$ GeV, considering a simplified scenario where DM interacts exclusively with Oxygen, a notable modification in the Boron-to-Carbon spectrum due to the DM-CR interaction is observed. Particularly, the peak within the spectrum, spanning from $0.1$ GeV to $10$ GeV, experiences an enhancement of approximately 1.5 times. However, in a more realistic scenario where DM particles interact with all CRs, this peak can be amplified to twice its original value.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$ MeV, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$., Comment: 19 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

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 paper, 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 nanohertz gravitational waves reported by pulsar timing array (PTA) 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 extragalactic magnetic field inferred with GRB 221009A., Comment: 7 pages, 3 figures, 1 table

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 published

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. 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

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

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

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

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

We perform a comprehensive analysis for the light scalar dark matter (DM)in the Inert two Higgs doublet model (i2HDM) with compressed mass spectra, small mass splittings among three $\mathbb{Z}_2$ odd particles---scalar $S$, pseudo-scalar $A$, and charged Higgs $H^\pm$. In such a case, the co-annihilation processes play a significant role to reduce DM relic density. As long as a co-annihilation governs the total interaction rate in the early universe, a small annihilation rate is expected to reach a correct DM relic density and its coupling $\lambda_S$ between DM pair and Higgs boson shall be tiny. Consequently, a negligible DM-nucleon elastic scattering cross section is predicted at the tree-level. In this work, we include the one-loop quantum corrections of the DM-nucleon elastic scattering cross section. We found that the quartic self-coupling $\lambda_2$ between $\mathbb{Z}_2$ odd particles indeed contributes to the one-loop quantum correction and behaves non-trivially for the co-annihilation scenario. Interestingly, the parameter space, which is allowed by the current constraints considered in this study, can predict the DM mass and annihilation cross section at the present compatible with the AMS-02 antiproton excess. The parameter space can be further probed at the future high luminosity LHC., Comment: 41 pages, JHEP accepted version

The complex scalar dark matter (DM) candidate in the gauged two Higgs doublet model (G2HDM), stabilized by a peculiar hidden parity ($h$-parity), is studied in detail. We explore the parameter space for the DM candidate by taking into account the most recent DM constraints from various experiments, in particular, the PLANCK relic density measurement and the current DM direct detection limit from XENON1T. We separate our analysis in three possible compositions for the mixing of the complex scalar. We first constrain our parameter space with the vacuum stability and perturbative unitarity conditions for the scalar potential, LHC Higgs measurements, plus Drell-Yan and electroweak precision test constraints on the gauge sector. We find that DM dominated by composition of the inert doublet scalar is completely excluded by further combining the previous constraints with both the latest results from PLANCK and XENON1T. We also demonstrate that the remaining parameter space with two other DM compositions can be further tested by indirect detection like the future CTA gamma-ray telescope., Comment: 64 pages, 10 figures, 3 tables

We study the constraints on the new parameters in the gauge sector of gauged two Higgs doublet model using the electroweak precision test data collected from the Large Electron Positron Collider (LEP) at and off the Z-pole as well as the current Drell-Yan and high-mass dilepton resonance data from the Large Hadron Collider (LHC). Impacts on the new parameters by the projected sensitivities of various electroweak observables at the Circular Electron Positron Collider (CEPC) proposed to be built in China are also discussed. We also clarify why the St\"{u}eckelberg mass $M_Y$ for the hypercharge $U(1)_Y$ is set to be zero in the model by showing that it would otherwise lead to the violation of the standard charge assignments for the elementary quarks and leptons., Comment: 51 pages, 11 figures, update references and correct typos

A light fermionic weakly interacting massive particle (WIMP) dark matter is investigated by studying its minimal renormalizable model, where it requires a scalar mediator to have an interaction between the WIMP and standard model particles. We perform a comprehensive likelihood analysis of the model involving the latest but robust constraints and those will be obtained in the near future. In addition, we pay particular attention to properly take the kinematically equilibrium condition into account. It is shown that near-future experiments and observations such as low-mass direct dark matter detections, flavor experiments and CMB observations play important roles to test the model. Still, a wide parameter region will remain even if no WIMP and mediator signals are detected there. We also show that precise Higgs boson measurements at future lepton colliders will significantly test this remaining region., Comment: 50 pages, 52 figures

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta., Comment: 32 pages, 8 figures, 331 references

We study the theoretical and phenomenological constraints imposed on the scalar sector of the gauged two Higgs doublet model proposed recently as a variant of the popular inert Higgs doublet model of dark matter. The requirements of tree-level vacuum stability and perturbative unitarity in the scalar sector are analyzed in detail. Furthermore, taking into account the constraints from the 125 GeV Higgs boson measurements at the Large Hadron Collider, we map out the allowed ranges for the fundamental parameters of the scalar potential in the model., Comment: 41 pages, 31 figures. V2: Content added to section III.A, references added, content matches version published in PRD

In this paper we study sneutrino dark matter in a recently proposed supersymmetric electroweak-scale inverse seesaw model, in which the majority of the sneutrino dark matter particle is a mixture of the right-handed sneutrino $\tilde{N}^c$ and the singlet field $\tilde{S}$. The scalar field $X$ responsible for the generation of neutrino masses can simultaneously play a crucial role for sneutrino annihilation in the early Universe via the pseudoscalar mediator $A_X$ into neutrinos. We focus here on the dominant annihilation channels and provide all the formulas together with analytic estimates in order to identify the relevant parameters. Furthermore, we show that the direct detection scattering cross section is many orders of magnitude below the current limits, and estimate the indirect detection annihilation rate, which is only a few orders of magnitude below the current limits., Comment: 24 pages, 8 figures, and 2 tables. https://sarah.hepforge.org/trac/wiki/inverse-Seesaw-SUSY-EW, v2: correspond to published version

The potential antiproton excess in the AMS-02 data is of much interest and can probably come from dark matter annihilations. Based on the effective field theory approach, in this work we investigate the compatibility of the DM interpretation of the AMS-02 antiproton excess and the null results from direct detection experiments, LUX, PandaX-II, and XENON1T. We focus on dimension-five and -six operators with fermion DM. Only one of dimension-five and one of dimension-six operators can successfully account for the antiproton excess, while the rest either are excluded by direct detection or require very small cut-off scales which invalidate the effective field theory approach.

The recently announced results on the 21-cm absorption spectrum by the EDGES experiment can place very stringent limits on dark matter annihilation cross-sections. We properly take into account the heating energy released from dark matter annihilation from the radiation epoch to the 21-cm observation redshifts in the radiative transfer to compute the evolution of the gas temperature. Our results show that the global 21-cm absorption profile is a powerful cosmological probe of the dark matter interactions. For dark matter annihilating into electron-positron pairs, the EDGES results give a more stringent upper limit than the PLANCK result on the annihilation cross section at the lower dark matter mass region., Comment: v3, accepted version

We propose a flavored $U(1)_{e\mu}$ neutrino mass and dark matter~(DM) model to explain the recent DArk Matter Particle Explorer (DAMPE) data, which feature an excess on the cosmic ray electron plus positron flux around 1.4 TeV. Only the first two lepton generations of the Standard Model are charged under the new $U(1)_{e\mu}$ gauge symmetry. A vector-like fermion $\psi$, which is our DM candidate, annihilates into $e^{\pm}$ and $\mu^{\pm}$ via the new gauge boson $Z'$ exchange and accounts for the DAMPE excess. We have found that the data favors a $\psi$ mass around 1.5~TeV and a $Z'$ mass around 2.6~TeV, which can potentially be probed by the next generation lepton colliders and DM direct detection experiments., Comment: 7 pages, 3 figures. V2: version accepted by Physics Letters B

Weakly interacting massive particles are a widely well-probed dark matter candidate by the dark matter direct detection experiments. Theoretically, there are a large number of ultraviolet completed models that consist of a weakly interacting massive particle dark matter. The variety of models makes the comparison with the direct detection data complicated and often non-trivial. To overcome this, in the non-relativistic limit, the effective theory was developed in the literature which works very well to significantly reduce the complexity of dark matter-nucleon interactions and to better study the nuclear response functions. In the effective theory framework for a spin-1/2 dark matter, we combine three independent likelihood functions from the latest PandaX, LUX, and XENON1T data, and give a joint limit on each effective coupling. The astrophysical uncertainties of the dark matter distribution are also included in the likelihood. We further discuss the isospin violating cases of the interactions. Finally, for both dimension-five and dimension-six effective theories above the electroweak scale, we give updated limits of the new physics mass scales., Comment: 33 pages, 11 figures, PandaX run10 data included and version accepted in JHEP, "code is available at the LikeDM website, https://likedm.hepforge.org/"

With recent Lyman-alpha forest data from BOSS and XQ-100, some studies suggested that the lower mass limit on the fuzzy dark matter (FDM) particles is lifted up to $10^{-21}\,\mathrm{eV}$. However, such a limit was obtained by $\Lambda$CDM simulations with the FDM initial condition and the quantum pressure (QP) was not taken into account which could have generated non-trivial effects in large scales structures. We investigate the QP effects in cosmological simulations systematically, and find that the QP leads to further suppression of the matter power spectrum at small scales, as well as the halo mass function in the low mass end. We estimate the suppressing effect of QP in the 1D flux power spectrum of Lyman-alpha forest and compare it with data from BOSS and XQ-100. The rough uncertainties of thermal gas properties in the flux power spectrum model calculation were discussed. We conclude that more systematic studies, especially with QP taken into account, are necessary to constrain FDM particle mass using Lyman-alpha forest., Comment: 22 pages, 8 figures, "for simulation code, see https://github.com/liambx/Axion-Gadget", Accepted by ApJ, there are minor changes in the published version, main content is the same

Recently a gauged two Higgs doublet model, in which the two Higgs doublets are embedded into the fundamental representation of an extra local $SU(2)_H$ group, is constructed. Both the new gauge bosons $Z^\prime$ and $W^{\prime (p,m)}$ are electrically neutral. While $Z^\prime$ can be singly produced at colliders, $W^{\prime (p,m)}$, which is heavier, must be pair produced. We explore the constraints of $Z^\prime$ using the current Drell-Yan type data from the Large Hadron Collider. Anticipating optimistically that $Z^\prime$ can be discovered via the clean Drell-Yan type signals at high luminosity upgrade of the collider, we explore the detectability of extra heavy fermions in the model via the two leptons/jets plus missing transverse energy signals from the exotic decay modes of $Z^\prime$. For the $W^{\prime (p,m)}$ pair production in a future 100 TeV proton-proton collider, we demonstrate certain kinematical distributions for the two/four leptons plus missing energy signals have distinguishable features from the Standard Model background. In addition, comparisons of these kinematical distributions between the gauged two Higgs doublet model and the littlest Higgs model with T-parity, the latter of which can give rise to the same signals with competitive if not larger cross sections, are also presented., Comment: 39 pages, 23 figures, 7 tables and two new appendixes, to appear in EPJC

In this paper we propose a novel supersymmetric inverse seesaw model which has only one additional $Z_6$ symmetry. The field content is minimal to get a viable neutrino spectrum at tree-level. Interestingly, the inverse seesaw scale in our model is related to the scale of electroweak symmetry breaking. Due to that origin we are less biased about hierarchies and discuss three different types of the inverse seesaw mechanism with different phenomenologies. We can successfully reproduce neutrino masses and mixing and our model is consistent with current bounds on neutrinoless double beta decay, non-unitarity of the PMNS matrix and charged lepton flavor violation., Comment: 20 pages, 1 figure; version published in JHEP

The Ultra-Light Axion (ULA) is a dark matter candidate with mass $\mathcal{O}(10^{-22})$eV and de Broglie wavelength of order kpc. Such an axion, also called the Fuzzy Dark Matter (FDM), thermalizes via the gravitational force and forms a Bose-Einstein condensate. Recent studies suggested that the quantum pressure from the FDM can significantly affect the structure formation in small scales, thus alleviating the so-called "small-scale crisis". In this paper, we develop a new technique to discretize the quantum pressure and illustrate the interactions among FDM particles in the $N$-body simulation, which accurately simulates the formation of the dark-matter halo and its inner structure in the region outside the softening length. In a self-gravitationally-bound virialized halo, we find a constant density, solitonic core, which is consistent with the theoretical prediction. The existence of the solitonic core reveals the non-linear effect of quantum pressure and impacts the structure formation in the FDM model., Comment: 17 pages, 9 figures, Accepted by ApJ, there are minor changes in the published version, main content is the same

Using the latest AMS-02 cosmic ray antiproton flux data, we search for potential dark matter annihilation signal. The background parameters about the propagation, source injection, and solar modulation are not assumed {\it a priori}, but based on the results inferred from the recent B/C ratio and proton data measurements instead. The possible dark matter signal is incorporated into the model self-consistently under a Bayesian framework. Compared with the astrophysical background only hypothesis, we find that a dark matter signal is favored. The rest mass of the dark matter particles is $\sim 20-80$ GeV and the velocity-averaged hadronic annihilation cross section is about $(0.2-5)\times 10^{-26}$ cm$^{3}$s$^{-1}$, in agreement with that needed to account for the Galactic center GeV excess and/or the weak GeV emission from dwarf spheroidal galaxies Reticulum 2 and Tucana III. Tight constraints on the dark matter annihilation models are also set in a wide mass region., Comment: 4 pages + 3 pages (supplemental material); published in Phys. Rev. Lett

We consider indirect detection of meta-stable dark matter particles decaying into a stable neutral particle and a pair of standard model fermions. Due to the softer energy spectra from the three-body decay, such models could potentially explain the AMS-02 positron excess without being constrained by the Fermi-LAT gamma-ray data and the cosmic ray anti-proton measurements. We scrutinize over different final state fermions, paying special attention to handling of the cosmic ray background and including various contributions from cosmic ray propagation with the help of the \textsc{LikeDM} package. It is found that primary decays into an electron-positron pair and a stable neutral particle could give rise to the AMS-02 positron excess and, at the same time, stay unscathed against the gamma-ray and anti-proton constraints. Decays to a muon pair or a mixed flavor electron-muon pair may also be viable depending on the propagation models. Decays to all other standard model fermions are severely disfavored., Comment: 22 pages, 7 figures and 2 tables; updated to AMS-02 antiproton data; added the comparison with the newest AMS-02 positron fluxes; accepted for publication in JCAP

We enumerate the set of simplified models which match onto the complete set of gauge invariant effective operators up to dimension six describing interactions of a singlet-like Majorana fermion dark matter with the standard model. Tree level matching conditions for each case are worked out in the large mediator mass limit, defining a one to one correspondence between the effective operator coefficients and the simplified model parameters for weakly interacting models. Utilizing such a mapping, we compute the dark matter annihilation rate in the early universe, as well as other low-energy observables like nuclear recoil rates using the effective operators, while the simplified models are used to compute the dark matter production rates at high energy colliders like LEP, LHC and future lepton colliders. Combining all relevant constraints with a profile likelihood analysis, we then discuss the currently allowed parameter regions and prospects for future searches in terms of the effective operator parameters, reducing the model dependence to a minimal level. In the parameter region where such a model-independent analysis is applicable, and leaving aside the special dark matter mass regions where the annihilation proceeds through an s-channel Z or Higgs boson pole, the current constraints allow effective operator suppression scales ($\Lambda$) of the order of a few hundred GeV for dark matter masses $m_\chi >$ 20 GeV at $95\%$ C.L., while the maximum allowed scale is around $3$ TeV for $m_\chi \sim \mathcal{O}(1\,{\rm TeV})$. An estimate of the future reach of ton-scale direct detection experiments and planned electron-positron colliders show that most of the remaining regions can be probed, apart from dark matter masses near half of the Z-boson mass (with $500\,{\rm GeV} < \Lambda < 2\,{\rm TeV} $) and those beyond the kinematic reach of the future lepton colliders., Comment: 30 pages, 6 figures

With the large progress in searches for dark matter (DM) particles with indirect and direct methods, we develop a numerical tool that enables fast calculations of the likelihoods of specified DM particle models given a number of observational data, such as charged cosmic rays from space-borne experiments (e.g., PAMELA, AMS-02), gamma-rays from the Fermi space telescope, and underground direct detection experiments. The purpose of this tool --- LikeDM, likelihood calculator for dark matter detection --- is to bridge the gap between a particle model of DM and the observational data. The intermediate steps between these two, including the astrophysical backgrounds, the propagation of charged particles, the analysis of Fermi gamma-ray data, as well as the DM velocity distribution and the nuclear form factor, have been dealt with in the code. We release the first version (v1.0) focusing on the constraints from indirect detection of DM with charged cosmic and gamma rays. Direct detection will be implemented in the next version. This manual describes the framework, usage, and related physics of the code. The code LikeDM can be download from https://likedm.hepforge.org/, Comment: v2, published in Computer Physics Communications, "LikeDM can be download from https://likedm.hepforge.org/"