Dwarf spheroidal galaxies (dSphs) are promising targets of indirect detection experiments searching for dark matter (DM) at present universe. Toward robust prediction for the amount of signal flux originating in DM annihilation inside dSphs, a precise determination of DM distributions as well as $J$-factors of the dSphs is particularly important. In this work, we estimate those of Draco, Sculptor, and Ursa Minor dSphs by an improved statistical method in which both foreground stars and dSph member stars are simultaneously taken into account. We define the likelihood function of the method as the so-called conditional one to remove sampling bias of observed stellar data. This improved method enables us to estimate DM distributions and $J$-factors of the dSphs directly from observed stellar data contaminated by foreground stars without imposing stringent membership criteria on the measured quantities.

Dwarf spheroidal galaxies (dSphs) are promising targets for the gamma-ray dark matter (DM) search. In particular, DM annihilation signal is expected to be strong in some of the recently discovered nearby ultra-faint dSphs, which potentially give stringent constraints on the $\mathcal{O}(1)$ TeV WIMP DM. However, various non-negligible systematic uncertainties complicate the estimation of the astrophysical factors relevant for the DM search in these objects. Among them, the effects of foreground stars particularly attract attention because the contamination is unavoidable even for the future kinematical survey. In this article, we assess the effects of the foreground contamination on the astrophysical $J$-factor estimation by generating mock samples of stars in the four ultra-faint dSphs and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and apply a likelihood analysis which takes member and foreground stellar distributions into account. We show that the foreground star contaminations in the signal region (the region of interest) and their statistical uncertainty can be estimated by interpolating the foreground star distribution in the control region where the foreground stars dominate the member stars. Such regions can be secured at future spectroscopic observations utilizing a multiple object spectrograph with a large field of view; e.g. the Prime Focus Spectrograph mounted on Subaru Telescope. The above estimation has several advantages: The data-driven estimation of the contamination makes the analysis of the astrophysical factor stable against the complicated foreground distribution. Besides, foreground contamination effect is considered in the likelihood analysis.

A novel idea of the direct detection to search for a ultralight dark matter based on the interaction between the dark matter and a nucleon is proposed. Solar system bodies feel the dark matter wind and it acts as a resistant force opposing their motions. The astronomical ephemeris of solar system bodies is so precise that it has a strong capability to detect a dark matter whose mass is much lighter than O(1) eV. We have estimated the resistant force based on the calculation of the elastic scattering cross section between the dark matter and the bodies beyond the Born approximation, and show that the astronomical ephemeris indeed put a very strong constraint on the interaction between the dark matter and a nucleon, depending on how smoothly the ultralight dark matter is distributed at the scale smaller than the celestial bodies in our solar system.

The gamma-ray observation of the dwarf spheroidal galaxies (dSphs) is a promising approach to search for the dark matter annihilation (or decay) signal. The dSphs are the nearby satellite galaxies with a clean environment and dense dark matter halo so that they give stringent constraints on the 1 TeV dark matter. However, recent studies have revealed that current estimation of astrophysical factors relevant for the dark matter searches are not conservative, where the various non-negligible systematic uncertainties are not taken into account. Among them, the effect of foreground stars on the astrophysical factors has not been paid much attention, which becomes more important for deeper and wider stellar surveys in the future. In this article, we assess the effects of the foreground contamination by generating the mock samples of stars and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and find that the cuts on the velocity and surface gravity can efficiently eliminate the contamination. We also propose a new likelihood function which includes the foreground distribution function. We apply this likelihood function to the fit of the three types of the mock data (Ursa Minor, Draco with large dark matter halo, and Draco with small halo) and three cases of the observation. The likelihood successfully reproduces the input $J$-factor value while the fit without considering the foreground distribution gives large deviation from the input value by a factor of three.

The dwarf spheroidal galaxies (dSphs) in the Milky Way are the primary targets in the indirect searches for particle dark matter. To set robust constraints on candidate of dark matter particles, understanding the dark halo structure of these systems is of substantial importance. In this paper, we first evaluate the astrophysical factors for dark matter annihilation and decay for 24 dSphs, with taking into account a non-spherical dark halo, using generalized axisymmetric mass models based on axisymmetric Jeans equations. First, from a fitting analysis of the most recent kinematic data available, our axisymmetric mass models are a much better fit than previous spherical ones, thus, our work should be the most realistic and reliable estimator for astrophysical factors. Secondly, we find that among analysed dSphs, the ultra-faint dwarf galaxies Triangulum II and Ursa Major II are the most promising but large uncertain targets for dark matter annihilation while the classical dSph Draco is the most robust and detectable target for dark matter decay. It is also found that the non-sphericity of luminous and dark components influences the estimate of astrophysical factors, even though these factors largely depend on the sample size, the prior range of parameters and the spatial extent of the dark halo. Moreover, owing to these effects, the constraints on the dark matter annihilation cross-section are more conservative than those of previous spherical works. These results are important for optimizing and designing dark matter searches in current and future multi-messenger observations by space and ground-based telescopes.

The AMS-02 collaboration has recently reported the antiproton to proton ratio with improved accuracy. In view of uncertainties of the production and the propagation of the cosmic rays, the observed ratio is still consistent with the secondary astrophysical antiproton to proton ratio. However, it is nonetheless enticing to examine whether the observed spectrum can be explained by a strongly motivated dark matter, the wino dark matter. As we will show, we find that the antiproton flux from the wino annihilation can explain the observed spectrum well for its mass range 2.5-3 TeV. The fit to data becomes particularly well compared to the case without the annihilation for the thermal wino dark matter case with a mass about 3 TeV. The ratio is predicted to be quickly decreased at the energy several hundreds of GeV, if this possibility is true, and it will be confirmed or ruled out in near future when the AMS-02 experiment accumulates enough data at this higher energy region.

We discuss the indirect detection of the wino dark matter utilizing gamma-ray observations of dwarf spheroidal galaxies (dSphs). After carefully reviewing current limits with particular attention to astrophysical uncertainties, we show prospects of the wino mass limit in future gamma-ray observation by the Fermi-LAT and the GAMMA-400 telescopes. We find that the improvement of the so-called $J$-factor of both the classical and the ultra-faint dSphs will play a crucial role to cover whole mass range of the wino dark matter. For example, with $\delta (\log_{10}J) = 0.1$ for both the classical and the ultra-faint dSphs, whole wino dark matter mass range can be covered by 15 years and 10 years data at the Fermi-LAT and GAMMA-400 telescopes, respectively.

The large value of the tensor-to-scalar ratio in the cosmic microwave background radiation reported by the BICEP2 collaboration gives strong impact on models of supersymmetry (SUSY). The large ratio indicates inflation with a high-energy scale and thus a high reheating temperature in general, and various SUSY models suffer from the serious gravitino and Polonyi problems. In this article, we discuss a class of the high-scale SUSY breaking models which are completely free from those problems. With especially focusing on the dark matter relic abundance, we examine how the BICEP2 result narrows down the parameter space of the models, assuming the simplest chaotic inflation model. We find that the mass of the dark matter is predicted to be less than about 1 TeV thanks to the non-thermal production in the early universe through the decay of abundant gravitinos produced after the reheating process. We also discuss implications in some details to dark matter searches at collider and indirect dark matter detection experiments.

We study phenomenological aspects of the bino-wino co-annihilation scenario in high-scale supersymmetry breaking models. High-scale SUSY breaking scenarios are considered to be promising possibility after the discovery of the Higgs boson with a mass around 126 GeV. In this paper, we discuss the bino lightest supersymmetric particle (LSP) accompanied by the at most around 30 GeV heavier wino. With the suitable mass splitting between the bino and the wino, the bino LSP has the correct relic abundance of dark matter. For the smaller mass splitting, the late-time decay of the gravitino can provide the correct abundance of the bino dark matter. It is extremely challenging to find signals from the bino dark matter in direct and indirect detections. By utilizing multi-jets plus missing transverse momentum events at the LHC, we can constraint the gluino mass and thus probe the bino-wino co-annihilation scenario indirectly. The collider experiment, however, can not search the bino dark matter directly. In this paper, we suggest the direct probe of the bino dark matter. We show that the bino dark matter leaves imprints on the small-scale matter power spectrum when the bino dark matter is produced by the decay of the gravitino. The non-thermal bino dark matter behaves as mixed (cold+warm) dark matter.

A monochromatic line in the cosmic neutrino spectrum would be a smoking gun signature of dark matter. It is intriguing that the IceCube experiment has recently reported two PeV neutrino events with energies that may be equal up to experimental uncertainties, and which have a probability of being a background fluctuation estimated to be less than a percent. Here we explore prospects for these events to be the first indication of a monochromatic line signal from dark matter. While measurable annihilation signatures would seem to be impossible at such energies, we discuss the dark matter quantum numbers, effective operators, and lifetimes which could lead to an appropriate signal from dark matter decays. We will show that the set of possible decay operators is rather constrained, and will focus on several viable candidates which could explain the IceCube events; R-parity violating gravitinos, hidden sector gauge bosons, and singlet fermions in an extra dimension. In essentially all cases we find that a PeV neutrino line signal from dark matter would be accompanied by a potentially observable continuum spectrum of neutrinos rising towards lower energies.

We study how "warm" the wino dark matter is when it is non-thermally produced by the decays of the gravitino in the early Universe. We clarify the energy distribution of the wino at the decay of the gravitino and the energy loss process after their production. By solving the Boltzmann equation, we show that a sizable fraction of the wino dark matter can be "warm" for the wino mass m_{\tilde w} \sim 100-500 GeV. The "warmness" of the wino dark matter leaves imprints on the matter power spectra and may provide further insights on the origin of dark matter via the future 21 cm line survey. Our calculations can be applied to other non-thermal wino production scenarios such as the wino dark matter produced by the decay of the moduli fields.

One of important properties of dark matter is its stability. The U(1)$_{\rm B-L}$ gauge symmetry is the most attractive symmetry to guarantee the stability. Though the symmetry is expected to be broken at very high energy scale to account for tiny neutrino masses through the seesaw mechanism, the residual discrete symmetry of U(1)$_{\rm B-L}$ can stabilize the dark matter naturally. We prove that, when there is new physics connecting B$-$L charges of dark matter and standard model particles at the scale between the electroweak and the U(1)$_{\rm B-L}$ breaking, the mass of dark matter is definitely predicted to be (5--7)/$Q_{\rm DM}$ GeV ($Q_{DM}$ is the B$-$L charge of dark matter) independent not only of details of the new physics but also of its energy scale. We also show two attractive examples. First one is the scalar dark matter with the B$-$L charge of one, which turns out to be very consistent with current CoGeNT results. Another one is the fermionic dark matter having the B$-$L charge of one third, which is also attractive from the viewpoint of model building.

We study cosmic-rays in decaying dark matter scenario, assuming that the dark matter is the lightest superparticle and it decays through a R-parity violating operator. We calculate the fluxes of cosmic-rays from the decay of the dark matter and those from the standard astrophysical phenomena in the same propagation model using the GALPROP package. We reevaluate the preferred parameters characterizing standard astrophysical cosmic-ray sources with taking account of the effects of dark matter decay. We show that, if energetic leptons are produced by the decay of the dark matter, the fluxes of cosmic-ray positron and electron can be in good agreements with both PAMELA and Fermi-LAT data in wide parameter region. It is also discussed that, in the case where sizable number of hadrons are also produced by the decay of the dark matter, the mass of the dark matter is constrained to be less than 200-300 GeV in order to avoid the overproduction of anti-proton. We also show that the cosmic gamma-ray flux can be consistent with the results of Fermi-LAT observation if the mass of the dark matter is smaller than nearly 4 TeV.

The positron anomaly recently reported by the cosmic-ray measurements suggests that, if explained by the decay of dark matter particle, the decay source is closely linked up with the leptonic sector of the standard model. It is observed that, with a simple dimensional analysis, the lifetime of dark matter for the anomaly is expressed by the energy scale of neutrino masses. We present two scenarios in which these two matter at issue (the dark matter width and the tiny neutrino masses) stem from a single operator involving a gauge-singlet scalar field.

The positron anomaly recently reported by the cosmic-ray measurements can be explained by the decaying dark matter scenario, where it decays mainly into leptons with the lifetime of O(10^26) second. When the dark matter is a fermionic particle, the lifetime of this order is known to be obtained by a dimension 6 operator suppressed by the unification scale 10^16 GeV, while such decay operators do not necessarily involve only leptons. In addition, the scenario would be spoiled if there exist lower-dimensional operators inducing the dark matter decay. We show in this letter that a single non-Abelian discrete symmetry such as A_4 is possible to prohibit all such harmful (non-leptonically coupled and lower-dimensional) operators. Moreover, the dark matter decays into charged leptons in a flavor-blind fashion due to the non-Abelian flavor symmetry, which results in perfect agreements not only with the PAMELA data but also with the latest Fermi-LAT data reported very recently. We also discuss some relevance between the discrete symmetry and neutrino physics.

The CDMS Collaboration has reported two candidate events for dark matter. If the events are due to the elastic scattering of dark matter, the dark matter would be a WIMP dark matter with its mass of the order of 10-100GeV and its scattering cross section with a nucleon is about 10^-43cm^2. We show that such a dark matter is properly realized as a neutralino dark matter in the light higgs boson scenario of the MSSM. The lightest higgs boson mass can be lighter than 114.4GeV in the scenario because of a suppressed interaction between higgs boson and Z bosons. As a result, a large scattering cross section between the dark matter and ordinary matter is obtained.

Motivated by the PAMELA anomaly in the fluxes of cosmic-ray positron and electron, we study the cosmic gamma-ray induced by the inverse Compton (IC) scattering process in unstable dark matter scenario assuming that the anomaly is due to the positron and electron emission by the decay of dark matter. We calculate the fluxes of IC-induced gamma-ray produced in our Galaxy and that from cosmological distance, and show that both of them are significant. We discuss a possibility that large dark matter mass over TeV scale might be constrained by the gamma-ray observation by Fermi Gamma-ray Space Telescope.

Motivated by the PAMELA anomaly in the fluxes of cosmic-ray electron and positron, we study the cosmic gamma-ray induced by the inverse Compton (IC) scattering process in unstable dark matter scenario assuming that the anomaly is due to the electron and positron emission by the decay of dark matter. We calculate the fluxes of IC-induced gamma-ray produced in our Galaxy and that from cosmological distance, and show that both of them are significant. If the gamma-ray flux is precisely determined by Fermi Gamma-ray Space Telescope for various line-of-sight directions, it will provide an important test of the decaying dark matter scenario.

Motivated by the recent PAMELA and ATIC results, we calculate the electron and positron fluxes from the decay of lightest-superparticle (LSP) dark matter. We assume that the LSP is the dominant component of dark matter, and consider the case that the R-parity is very weakly violated so that the lifetime of the LSP becomes of the order of 10^26 sec. We will see that, with such a choice of the lifetime, the cosmic-ray electron and positron from the decay can be the source of the anomalous electron and positron fluxes observed by PAMELA and ATIC. We consider the possibilities that the LSP is the gravitino, the lightest neutralino, and scalar neutrino, and discuss how the resultant fluxes depend on the dark-matter model. We also discuss the fluxes of gamma-ray and anti-proton, and show that those fluxes can be consistent with the observed value in the parameter region where the PAMELA and ATIC anomalies are explained.

Motivated by the anomalous positron flux recently reported by the PAMELA collaboration, we study the cosmic-ray positron produced by the pair annihilation and the decay of superparticle dark matter. We calculate the cosmic-ray positron flux and discuss implications of the PAMELA data. We show that the positron excess observed by the PAMELA can be explained with some of the superparticle dark matter.

We investigate a possibility of precision measurements for parameters of the Littlest Higgs model with T-parity at the International Linear Collider (ILC). The model predicts new gauge bosons (AH, ZH, and WH), among which the heavy photon (AH) is a candidate for dark matter. The masses of these new gauge bosons strongly depend on the vacuum expectation value that breaks a global symmetry of the model. Through Monte Carlo simulations of the processes: e+ e- ->AH ZH and e+ e- -> WH+ WH-, we show how precisely the masses can be determined at the ILC for a representative parameter point of the model. We also discuss the determination of the Little Higgs parameters and its impact on the future measurement of the thermal abundance of the dark matter relics in our universe.

We discuss the synchrotron radiation flux from the Galactic center in unstable dark matter scenario. Motivated by the anomalous excess of the positron fraction recently reported by the PAMELA collaboration, we consider the case that the dark matter particle is unstable (and long-lived), and that energetic electron and positron are produced by the decay of dark matter. Then, the emitted electron and positron becomes the source of the synchrotron radiation. We calculate the synchrotron radiation flux for models of decaying dark matter, which can explain the PAMELA positron excess. Taking the lifetime of the dark matter of O(10^26 sec), which is the suggested value to explain the PAMELA anomaly, the synchrotron radiation flux is found to be O(1 kJy/str) or smaller, depending on the particle-physics and cosmological parameters.

We study gamma ray and positron in high energy cosmic ray from the decay of the gravitino dark matter in the framework of supersymmetric model with R-parity violation. Even though R-parity is violated, the lifetime of the gravitino, which is assumed to be the lightest superparticle, can be longer than the present age of the universe if R-parity violating interactions are weak enough. In such a case, gravitino can be dark matter of the universe and its decay produces high energy cosmic rays. We calculate the fluxes of gamma ray and positron from the decay of the gravitino dark matter and discuss implications of such a scenario to present and future observations. In particular, we show that excesses of the fluxes of gamma ray and positron observed by EGRET and HEAT experiments, respectively, can be simultaneously explained as the cosmic rays from the decay of the gravitino dark matter.

If the Dark Matter (DM) in the Universe has interactions with the standard-model particle, the pair annihilation may give the imprints in the cosmic ray. In this paper we study the pair annihilation processes of the DM, which is neutral, however has the electroweak (EW) gauge non-singlet. In this estimation the non-relativistic (NR) effective theory in the EW sector is a suitable technique. We find that if the DM mass is larger than about 1 TeV, the attractive Yukawa potentials induced by the EW gauge bosons have significant effects on the DM annihilation processes, and the cross sections may be enhanced by several orders of magnitude, due to the zero energy resonance under the potentials. Especially, the annihilation to two $\gamma$'s might have a comparable cross section to other tree-level processes, while the cross section under the conventional calculation is suppressed by a loop factor. We also discuss future sensitivities to the $\gamma$ ray from the galactic center by the GLAST satellite detector and the Air Cerenkov Telescope (ACT) arrays.

The neutralino pair annihilation into two photons in our galactic halo gives a robust dark matter signal, since it would give a quasi-monotonic gamma ray. This process is radiatively-induced, and the full-one loop calculation was done previously. However, for the heavy wino-like or Higgsino-like neutralino, the one-loop cross section violates unitarity, therefore the higher-order corrections may be important. We construct a non-relativistic theory for chargino and neutralino two-body states, and estimate all-order QED corrections and two-loop corrections by $Z$ and/or $W$ exchange. We find that the critical mass, above that the two-loop contribution is larger than one-loop one, is about 8 TeV (O(10) TeV) in the limit where neutralino is wino (Higgsino)-like, respectively. Around and above the critical mass, the all-order Z and/or W exchange must be included to estimate the cross section. On the other hand, the QED corrections depend on the mass difference between the neutralino and chargino. In the wino-like limit where neutralino is highly degenerate with chargino in mass, we find that QED corrections enhance the pair annihilation cross section by 1.5-2.