In this short letter, we find that a magnetic transition dipole moment between tau and sterile neutrinos can account for the XENON1T excess events. Unlike the ordinary neutrino dipole moment, the introduction of the new sterile mass scale allows for astrophysical bounds to be suppressed. Interestingly, the best-fit regions that are compatible with the SN1987A imply either boron-8 or CNO neutrinos as the source flux. We find that sterile neutrinos of either $\sim$ 260 keV or in the $\sim$(500 - 800) keV mass range are capable of evading astrophysical constraints while being able to successfully explain the XENON1T event rate. The sterile neutrino in the best fit parameter space may have significant effects on big bang nucleosynthesis (BBN). We show the region in which a low reheating temperature of the Universe may allow the BBN constraints to be alleviated.

Indirect dark matter (DM) detection typically involves the observation of standard model (SM) particles emerging from DM annihilation/decay inside regions of high dark matter concentration. We consider an annihilation scenario in which this reaction has to be initiated by one of the DMs involved being boosted while the other is an ambient non-relativistic particle. This "trigger" DM must be created, for example, in a previous annihilation or decay of a heavier component of DM. Remarkably, boosted DM annihilating into gamma-rays at a specific point in a galaxy could actually have traveled from its source at another point in the same galaxy or even from another galaxy. Such a "non-local" behavior leads to a non-trivial dependence of the resulting photon signal on the galactic halo parameters, such as DM density and core size, encoded in the so-called "astrophysical" $J$-factor. These non-local $J$-factors are strikingly different than the usual scenario. A distinctive aspect of this model is that the signal from dwarf galaxies relative to the Milky Way tends to be suppressed from the typical value to various degrees depending on their characteristics. This feature can thus potentially alleviate the mild tension between the DM annihilation explanation of the observed excess of $\sim$ GeV photons from the Milky Way's galactic center vs. the apparent non-observation of the corresponding signal from dwarf galaxies.

We investigate the possibility of the gravitational-wave event GW170817 being a light, solar-mass black hole (BH) - neutron star (NS) merger. We explore two exotic scenarios involving primordial black holes (PBH) that could produce such an event, taking into account available observational information on NGC 4993. First, we entertain the possibility of dynamical NS-PBH binary formation where a solar-mass PBH and a NS form a binary through gravitational interaction. We find that while dynamical NS-PBH formation could account for the GW170817 event, the rate is highly dependent on unknown density contrast factors and could potentially be affected by galaxy mergers. We also find that PBH-PBH binaries would likely have a larger merger rate, assuming the density contrast boost factor of an order similar to the NS-PBH case. These exotic merger formations could provide new channels to account for the volumetric rate of compact-object mergers reported by LIGO/Virgo. Secondly, we consider the case where one of the NS's in a binary NS system is imploded by a microscopic PBH. We find that the predicted rate for NS implosion into a BH is very small, at least for the specific environment of NGC 4993. We point out that similar existing (e.g. GW190425 and GW190814) and future observations will shed additional light on these scenarios.

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.

At present, the strongest upper limit on $\sum m_{\nu}$, the sum of neutrino masses, is from cosmological measurements. However, this bound assumes that the neutrinos are stable on cosmological timescales, and is not valid if the neutrino lifetime is less than the age of the universe. In this paper, we explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation on timescales of order the age of the universe, and determine the bound on the sum of neutrino masses in this scenario. We focus on the case in which the neutrinos decay after becoming non-relativistic. We derive the Boltzmann equations that govern the cosmological evolution of density perturbations in the case of unstable neutrinos, and solve them numerically to determine the effects on the matter power spectrum and lensing of the cosmic microwave background. We find that the results admit a simple analytic understanding. We then use these results to perform a Monte Carlo analysis based on the current data to determine the limit on the sum of neutrino masses as a function of the neutrino lifetime. We show that in the case of decaying neutrinos, values of $\sum m_{\nu}$ as large as 0.9 eV are still allowed by the data. Our results have important implications for laboratory experiments that have been designed to detect neutrino masses, such as KATRIN and KamLAND-ZEN.

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.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation after becoming non-relativistic. We show that in this scenario, near-future large scale structure measurements from the Euclid satellite, when combined with cosmic microwave background data from Planck, may allow an independent determination of both the lifetime of the neutrinos and the sum of their masses. These parameters can be independently determined because the Euclid data will cover a range of redshifts, allowing the growth of structure over time to be tracked. If neutrinos are stable on cosmological timescales, these observations can improve the lower limit on the neutrino lifetime by seven orders of magnitude, from $\mathcal{O}(10)$ years to $2\times 10^8$ years ($95\%$ C.L.), without significantly affecting the measurement of neutrino mass. On the other hand, if neutrinos decay after becoming non-relativistic but on timescales less than $\mathcal{O}(100)$ million years, these observations may allow, not just the first measurement of the sum of neutrino masses, but also the determination of the neutrino lifetime from cosmology.

In this article we investigate the effects of single derivative mixing in massive bosonic fields. In the regime of large mixing, we show that this leads to striking changes of the field dynamics, delaying the onset of classical oscillations and decreasing, or even eliminating, the friction due to Hubble expansion. We highlight this phenomenon with a few examples. In the first example, we show how an axion like particle can have its number abundance parametrically enhanced. In the second example, we demonstrate that the QCD axion can have its number abundance enhanced allowing for misalignment driven axion dark matter all the way down to $f_a$ of order astrophysical bounds. In the third example, we show that the delayed oscillation of the scalar field can also sustain a period of inflation. In the last example, we present a situation where an oscillating scalar field is completely frictionless and does not dilute away in time.

Boosted dark matter (BDM) is a well-motivated class of dark matter candidates in which a small subcomponent is relativistic at the present time. In this work, we lay the foundation for BDM searches via its hadronic interactions at DUNE. We investigate BDM-nucleus scattering in detail by developing new dedicated event generation techniques with a parameterized detector simulation. We demonstrate that the excellent event reconstruction capability of DUNE enables a significantly broader discovery reach.

Phase transitions in the early universe can readily create an observable stochastic gravitational wave background. We show that such a background necessarily contains anisotropies analogous to those of the cosmic microwave background (CMB) of photons, and that these too may be within reach of proposed gravitational wave detectors. Correlations within the gravitational wave anisotropies and their cross-correlations with the CMB can provide new insights into the mechanism underlying primordial fluctuations, such as multi-field inflation, as well as reveal the existence of non-standard ``hidden sectors" of particle physics in earlier eras.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

During the next decade, gravitational waves will be observed from hundreds of binary inspiral events. When the redshifts of the host galaxies are known, these events can be used as `standard sirens', sensitive to the expansion rate of the Universe. Measurements of the Hubble constant $H_0$ from standard sirens can be done independently of other cosmological probes, and events occurring at $z<0.1$ will allow us to infer $H_0$ independently of cosmological models. The next generation of spectroscopic galaxy surveys will play a crucial role in reducing systematic uncertainties in $H_0$ from standard sirens, particularly for the numerous `dark sirens' which do not have an electromagnetic counterpart. In combination with large spectroscopic data sets, standard sirens with an EM counterpart are expected to constrain $H_0$ to $\sim 1-2\%$ precision within the next decade. This is competitive with the best estimates of $H_0$ obtained to date and will help illuminate the current tension between existing measurements. Information on the galaxies that host the gravitational wave events will also shed light on the origin and evolution of compact object binaries.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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.

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.

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.

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.

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.

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.

Primordial magnetic fields in the dark sector can be transferred to magnetic fields in the visible sector due to a gauge kinetic mixing term. We show that the transfer occurs when the evolution of magnetic fields is dominated by dissipation due to finite electric conductivity, and does not occur at later times if the magnetic fields evolve according to magnetohydrodynamics scaling laws. The efficiency of the transfer is suppressed by not only the gauge kinetic mixing coupling but also the ratio between the large electric conductivity and the typical momentum of the magnetic fields. We find that the transfer gives nonzero visible magnetic fields today. However, without possible dynamo amplifications, the field transfer is not efficient enough to obtain the intergalactic magnetic fields suggested by the gamma-ray observations, although there are plenty of possibilities for efficient dark magnetogenesis, which are experimentally unconstrained.

We identify new astrophysical signatures of dark matter that implodes neutron stars (NSs), which could decisively test whether NS-imploding dark matter is responsible for missing pulsars in the Milky Way galactic center, the source of some $r$-process elements, and the origin of fast-radio bursts. First, NS-imploding dark matter forms $\sim 10^{-10}$ solar mass or smaller black holes inside neutron stars, which proceed to convert neutron stars into $\sim$1.5 solar mass BHs. This decreases the number of neutron star mergers seen by LIGO/Virgo (LV) and associated merger kilonovae seen by telescopes like DES, BlackGEM, and ZTF, and instead, producing a population of "black mergers" containing $\sim$1.5 solar mass black holes. Second, dark matter-induced neutron star implosions may create a new kind of kilonovae that lacks a detectable, accompanying gravitational signal, which we call "quiet kilonovae." Using DES data and the Milky Way's r-process abundance, we constrain quiet kilonovae. Third, the spatial distribution of neutron star merger kilonovae and quiet kilonovae in galaxies can be used to detect dark matter. NS-imploding dark matter destroys most neutron stars at the centers of disc galaxies, so that neutron star merger kilonovae would appear mostly in a donut at large radii. We find that as few as ten neutron star merger kilonova events, located to $\sim$1 kpc precision could validate or exclude dark matter-induced neutron star implosions at $2 \sigma$ confidence, exploring dark matter-nucleon cross-sections 4-10 orders of magnitude below current direct detection experimental limits. Similarly, NS-imploding dark matter as the source of fast radio bursts can be tested at $2 \sigma$ confidence once 20 bursts are located in host galaxies by radio arrays like CHIME and HIRAX.

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.

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.

The standard paradigm of collisionless cold dark matter is in tension with measurements on large scales. In particular, the best fit values of the Hubble rate $H_0$ and the matter density perturbation $\sigma_8$ inferred from the cosmic microwave background seem inconsistent with the results from direct measurements. We show that both problems can be solved in a framework in which dark matter consists of two distinct components, a dominant component and a subdominant component. The primary component is cold and collisionless. The secondary component is also cold, but interacts strongly with dark radiation, which itself forms a tightly coupled fluid. The growth of density perturbations in the subdominant component is inhibited by dark acoustic oscillations due to its coupling to the dark radiation, solving the $\sigma_8$ problem, while the presence of tightly coupled dark radiation ameliorates the $H_0$ problem. The subdominant component of dark matter and dark radiation continue to remain in thermal equilibrium until late times, inhibiting the formation of a dark disk. We present an example of a simple model that naturally realizes this scenario in which both constituents of dark matter are thermal WIMPs. Our scenario can be tested by future stage-IV experiments designed to probe the CMB and large scale structure.

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.

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.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

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/

We consider a twin WIMP scenario whose twin sector contains a full dark copy of the SM hadrons, where the lightest twin particles are twin pions. By analogy to the standard WIMP paradigm, the dark matter (DM) freezes out through twin electroweak interactions, and annihilates into a dark shower of light twin hadrons. These are either stable or decay predominantly to standard model (SM) photons. We show that this 'hadrosymmetric' scenario can be consistent with all applicable astrophysical, cosmological and collider constraints. In order to decay the twin hadrons before the big-bang nucleosynthesis epoch, an additional portal between the SM and twin sector is required. In most cases we find this additional mediator is within reach of either the LHC or future intensity frontier experiments. Furthermore, we conduct simulations of the dark shower and consequent photon spectra. We find that fits of these spectra to the claimed galactic center gamma-ray excess seen by Fermi-LAT non-trivially coincide with regions of parameter space that both successfully generate the observed DM abundance and exhibit minimal fine-tuning.

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

We attempt to explain simultaneously the Galactic center gamma-ray excess and the 3.5 keV X-ray line from the Perseus cluster based on a class of non-abelian $SU(2)$ DM models, in which the dark matter and an excited state comprise a "dark" $SU(2)$ doublet. The non-abelian group kinetically mixes with the standard model gauge group via dimensions-5 operators. The dark matter particles annihilate into standard model fermions, followed by fragmentation and bremsstrahlung, and thus producing a continuous spectrum of gamma-rays. On the other hand, the dark matter particles can annihilate into a pair of excited states, each of which decays back into the dark matter particle and an X-ray photon, which has an energy equal to the mass difference between the dark matter and the excited state, which is set to be 3.5 keV. The large hierarchy between the required X-ray and $\gamma$-ray annihilation cross-sections can be achieved by a very small kinetic mixing between the SM and dark sector, which effectively suppresses the annihilation into the standard model fermions but not into the excited state.

In light of the recent discovery by the ATLAS and CMS experiments at the Large Hadron Collider (LHC) of a Higgs-like particle with a narrow mass range of 125-126 GeV, we perform an updated analysis on one of the popular scalar dark matter models, the Inert Higgs Doublet Model (IHDM). We take into account in our likelihood analysis of various experimental constraints, including recent relic density measurement, dark matter direct and indirect detection constraints as well as the latest collider constraints on the invisible decay width of the Higgs boson and monojet search at the LHC. It is shown that if the invisible decay of the standard model Higgs boson is open, LHC as well as direct detection experiments like LUX and XENON100 could put stringent limits on the Higgs boson couplings to dark matter. We find that the most favoured parameter space for IHDM corresponds to dark matter with a mass less than 100 GeV or so. In particular, the best-fit points are at the dark matter mass around 70 GeV where the invisible Higgs decay to dark matter is closed. Scalar dark matter in the higher mass range of 0.5-4 TeV is also explored in our study. Projected sensitivities for the future experiments of monojet at LHC-14, XENON1T and AMS-02 one year antiproton flux are shown to put further constraints on the existing parameter space of IHDM.

We explore the MSSM with 9 free parameters (p9MSSM) that have been selected as a minimum set that allows an investigation of neutralino dark matter and collider signatures while maintaining consistency with several constraints. These include measurement of the dark matter relic density from PLANCK, main properties of the discovered Higgs boson, LHC direct SUSY searches, recent evidence for a Standard Model-like BR(Bs -> \mu+ \mu-), and the measurement of \delta(g-2), plus a number of other electroweak and flavor physics constraints. We perform a simulation of two LHC direct SUSY searches at sqrt(s)=8 TeV: the CMS inclusive \alpha_T search for squarks and gluinos and the CMS electroweak production search with 3l+E_T^miss in the final state. We use the latter to identify the regions of the parameter space, consistent at 2\sigma\ with \delta(g-2), that are not excluded by the direct limits from the electroweak production. We find that they correspond to a neutralino mass in the window 200 GeV<m_\chi<500 GeV. We also implement the likelihood for the XENON100 exclusion bound, in which we consider for the first time the impact of a recent determination of the \Sigma_{\pi N} term from CHAOS data, \Sigma_{\pi N}=43\pm12 MeV. We show that in light of this measurement, the present statistical impact of the XENON100 bound is greatly reduced, although future sensitivities of the LUX and XENON1T experiments will have decisive impact on the mixed bino/higgsino composition of the neutralino. We point out some tension between the constraints from \delta(g-2) and XENON100. Finally, we present prospects for various indirect searches of dark matter, namely \gamma-ray fluxes from dSphs and the Galactic Center at Fermi-LAT, and the positron flux at AMS02. We also show the 5-year sensitivity on the spin-dependent neutralino-proton cross section due to neutrino fluxes from the Sun at IceCube.

Observation of gamma-rays from dwarf galaxies is an effective way to search for particle dark matter. Using 4-year data of Fermi-LAT observations on a series of Milky Way satellites, we develop a general way to search for the signals from dark matter annihilation in such objects. Instead of giving prior information about the energy spectrum of dark matter annihilation, we bin the Fermi-LAT data into several energy bins and build a likelihood map in the "energy bin - flux" plane. The final likelihood of any spectrum can be easily derived through combining the likelihood of all the energy bins. It gives consistent result with that directly calculated using the Fermi Scientific Tool. This method is very efficient for the study of any specific dark matter models with gamma-rays. We use the new likelihood map with Fermi-LAT 4 year data to fit the parameter space in three representative dark matter models: i) toy dark matter model, ii) effective dark matter operators, and iii) supersymmetric neutralino dark matter.

We discuss the event rate in DeepCore array due to neutrino flux produced by annihilations and decays of galactic dark matter. This event rate is calculated with a 10 GeV threshold energy, which is smaller than the threshold energy taken in previous works. Taking into account the background event rate due to the atmospheric neutrino flux, we evaluate the sensitivity of DeepCore array for probing dark matter annihilation cross section and decay time. The sensitivity studies include the annihilation modes $\chi\chi\to b\bar{b}, \ \tau^+ \tau^-$, \ $\mu^+\mu^-$, and $\nu\bar{\nu}$, and decay modes $\chi\to b\bar{b}, \ \tau^+ \tau^-$, \ $\mu^+\mu^-$, and $\nu\bar{\nu}$. We compare our results with corresponding constraints derived from observations of WMAP, ACT and Fermi-LAT.

We investigate the impact of recent limits from LHC searches for supersymmetry and from direct and indirect searches for dark matter on global Bayesian inferences of the parameter space of the Constrained Minimal Supersymmetric Standard Model (CMSSM). In particular we apply recent exclusion limits from the CMS \alpha_T analysis of 1.1/fb of integrated luminosity, the current direct detection dark matter limit from XENON100, as well as recent experimental constraints on \gamma-ray fluxes from dwarf spheroidal satellite galaxies of the Milky Way from the FermiLAT telescope, in addition to updating values for other non-LHC experimental constraints. We extend the range of scanned parameters to include a significant fraction of the focus point/hyperbolic branch region. While we confirm earlier conclusions that at present LHC limits provide the strongest constraints on the model's parameters, we also find that when the uncertainties are not treated in an excessively conservative way, the new bounds from dwarf spheroidal have the power to significantly constrain the focus point/hyperbolic branch region. Their effect is then comparable, if not stronger, to that from XENON100. We further analyze the effects of one-year projected sensitivities on the neutrino flux from the Sun in the 86-string IceCube+DeepCore configuration at the South Pole. We show that data on neutrinos from the Sun, expected for the next few months at IceCube and DeepCore, have the potential to further constrain the same region of parameter space independently of the LHC and can yield additional investigating power for the model.

We present a global study of the simplest scalar phantom dark matter model. The best fit parameters of the model are determined by simultaneously imposing (i) relic density constraint from WMAP, (ii) 225 live days data from direct experiment XENON100, (iii) upper limit of gamma-ray flux from Fermi-LAT indirect detection based on dwarf spheroidal satellite galaxies, and (iv) the Higgs boson candidate with a mass about 125 GeV and its invisible branching ratio no larger than 40% if the decay of the Higgs boson into a pair of dark matter is kinematically allowed. The allowed parameter space is then used to predict annihilation cross sections for gamma-ray lines, event rates for three processes mono-b jet, single charged lepton and two charged leptons plus missing energies at the Large Hadron Collider, as well as to evaluate the muon anomalous magnetic dipole moment for the model.

Elementary keV sterile Dirac neutrinos can be a natural ingredient of the composite neutrino scenario. For a certain class of composite neutrino theories, these sterile neutrinos naturally have the appropriate mixing angles to be resonantly produced warm dark matter (WDM). Alternatively, we show these sterile neutrinos can be WDM produced by an entropy-diluted thermal freeze-out, with the necessary entropy production arising not from an out-of-equilibrium decay, but rather from the confinement of the composite neutrino sector, provided there is sufficient supercooling.

An effective interaction approach is used to describe the interactions between the spin 0 or spin 1/2 dark matter particle and the degrees of freedom of the standard model. This approach is applicable to those models in which the dark matter particles do not experience the standard-model interactions, e.g., hidden-sector models. We explore the effects of these effective interaction operators on (i) dark matter relic density, (ii) spin-independent and spin-dependent dark matter-nucleon scattering cross sections, (iii) cosmic antiproton and gamma ray fluxes from the galactic halo due to dark matter annihilation, and (iv) monojet and monophoton production plus missing energy at the Tevatron and the Large Hadron Collider (LHC). We combine the experimental data of relic density from WMAP7, spin-independent cross section from XENON100, spin-dependent cross section from XENON10, ZEPLIN-III, and SIMPLE, cosmic antiproton flux from PAMELA, cosmic gamma-ray flux from ${\it Fermi}$-LAT, and the monojet and monophoton data from the Tevatron and the LHC, to put the most comprehensive limits on each effective operator.

Motivated by the recent re-confirmation by CoGENT of the low-energy excess of events observed last year, and the recent improved limits from the XENON-100 experiment that are in contention with the CoGENT data, we re-examine the low mass neutralino region of the Minimal Supersymmetric Standard Model and of the Next-to-Minimal Supersymmetric Standard Model, both without assuming gaugino mass unification. We make several focused scans for each model, determining conservative constraints on input parameters. We then determine how these constraints are made increasingly stringent as we re-invoke our experimental constraints involving the dark matter relic abundance, collider constraints from LEP and the Tevatron, and then from flavour physics, as a series of successive 2 sigma hard cuts. We find that for both models, when all relevant constraints are applied in this fashion, we do not generate neutralino LSPs that possess a spin-independent scattering cross section in excess of 10^-5 pb and a mass ~7 GeV < m_chi < ~9 GeV that is necessary in order to explain the CoGENT observations.

Dark matter pair production at high energy colliders may leave observable signatures in the energy and momentum spectra of the objects recoiling against the dark matter. We use LEP data on mono-photon events with large missing energy to constrain the coupling of dark matter to electrons. Within a large class of models, our limits are complementary to and competitive with limits on dark matter annihilation and on WIMP-nucleon scattering from indirect and direct searches. Our limits, however, do not suffer from systematic and astrophysical uncertainties associated with direct and indirect limits. For example, we are able to rule out light (< 10 GeV) thermal relic dark matter with universal couplings exclusively to charged leptons. In addition, for dark matter mass below about 80 GeV, LEP limits are stronger than Fermi constraints on annihilation into charged leptons in dwarf spheroidal galaxies. Within its kinematic reach, LEP also provides the strongest constraints on the spin-dependent direct detection cross section in models with universal couplings to both quarks and leptons. In such models the strongest limit is also set on spin independent scattering for dark matter masses below ~4 GeV. Throughout our discussion, we consider both low energy effective theories of dark matter, as well as several motivated renormalizable scenarios involving light mediators.

We investigate the extent to which the uncertainties associated with the propagation of Galactic cosmic rays impact upon estimates for the gamma-ray flux from the mid-latitude region. We consider contributions from both standard astrophysical background (SAB) processes as well as resolved point sources. We have found that the uncertainties in the total gamma-ray flux from the mid-latitude region relating to propagation parameter values consistent with local B/C and Be10/Be9 data dominate by 1-2 orders of magnitude. These uncertainties are reduced to less than an order of magnitude when the normalisations of the SAB spectral components are fitted to the corresponding Fermi LAT data. We have found that for many propagation parameter configurations (PPCs) our fits improve when an extragalactic background (EGB) component is simultaneously fitted to the data. We also investigate the improvement in our fits when a flux contribution from neutralino dark matter (DM), described by the Minimal Supersymmetric Standard Model, was simultaneously fitted to the data. We consider three representative cases of neutralino DM for both Burkert and Einasto DM density profiles, in each case simultaneously fitting a boost factor of the DM contribution together with the SAB and EGB components. We have found that for several PPCs there are significant improvements in our fits, yielding both substantial EGB and DM components, where for a few of these PPCs the best-fit EGB component is consistent with recent estimates by the Fermi Collaboration.

GAMER is a GPU-accelerated Adaptive-MEsh-Refinement code for astrophysical simulations. In this work, two further extensions of the code are reported. First, we have implemented the MUSCL-Hancock method with the Roe's Riemann solver for the hydrodynamic evolution, by which the accuracy, overall performance and the GPU versus CPU speed-up factor are improved. Second, we have implemented the out-of-core computation, which utilizes the large storage space of multiple hard disks as the additional run-time virtual memory and permits an extremely large problem to be solved in a relatively small-size GPU cluster. The communication overhead associated with the data transfer between the parallel hard disks and the main memory is carefully reduced by overlapping it with the CPU/GPU computations.

We present the newly developed code, GAMER (GPU-accelerated Adaptive MEsh Refinement code), which has adopted a novel approach to improve the performance of adaptive mesh refinement (AMR) astrophysical simulations by a large factor with the use of the graphic processing unit (GPU). The AMR implementation is based on a hierarchy of grid patches with an oct-tree data structure. We adopt a three-dimensional relaxing TVD scheme for the hydrodynamic solver, and a multi-level relaxation scheme for the Poisson solver. Both solvers have been implemented in GPU, by which hundreds of patches can be advanced in parallel. The computational overhead associated with the data transfer between CPU and GPU is carefully reduced by utilizing the capability of asynchronous memory copies in GPU, and the computing time of the ghost-zone values for each patch is made to diminish by overlapping it with the GPU computations. We demonstrate the accuracy of the code by performing several standard test problems in astrophysics. GAMER is a parallel code that can be run in a multi-GPU cluster system. We measure the performance of the code by performing purely-baryonic cosmological simulations in different hardware implementations, in which detailed timing analyses provide comparison between the computations with and without GPU(s) acceleration. Maximum speed-up factors of 12.19 and 10.47 are demonstrated using 1 GPU with 4096^3 effective resolution and 16 GPUs with 8192^3 effective resolution, respectively.

The diffuse Galactic gamma-ray data from the region of the Galactic Center has been collected by the LAT instrument on the Fermi Gamma-Ray Space Telescope. In this paper we argue that it may be able to provide an unambiguous evidence of originating, in addition to known astrophysical sources, from dark matter annihilations in the halo, independently of the mass and other properties of the dark matter particle. We also show that the recently released high precision data from mid-latitudes is already providing an upper bound, albeit still a weak one, on the cuspiness of the dark matter density profile as a function of the mass of the dark matter assumed to be a stable neutralino of minimal supersymmetry.

In this paper, we describe the architecture and performance of the GraCCA system, a Graphic-Card Cluster for Astrophysics simulations. It consists of 16 nodes, with each node equipped with 2 modern graphic cards, the NVIDIA GeForce 8800 GTX. This computing cluster provides a theoretical performance of 16.2 TFLOPS. To demonstrate its performance in astrophysics computation, we have implemented a parallel direct N-body simulation program with shared time-step algorithm in this system. Our system achieves a measured performance of 7.1 TFLOPS and a parallel efficiency of 90% for simulating a globular cluster of 1024K particles. In comparing with the GRAPE-6A cluster at RIT (Rochester Institute of Technology), the GraCCA system achieves a more than twice higher measured speed and an even higher performance-per-dollar ratio. Moreover, our system can handle up to 320M particles and can serve as a general-purpose computing cluster for a wide range of astrophysics problems.