As a next-generation complex extensive air shower array with a large field of view, the large high altitude air shower observatory (LHAASO) is very sensitive to the very high energy gamma-rays from $\sim$ 300 GeV to 1 PeV, and may thus serve as an important probe for the heavy dark matter (DM) particles. In this study, we make a forecast for the LHAASO sensitivities to the gamma-ray signatures resulting from DM decay in dwarf spheroidal satellite galaxies (dSphs) within the LHAASO field of view. Both individual and combined limits for 19 dSphs incorporating the uncertainties of the DM density profile are explored. Owing to the large effective area and strong capability of the photon-proton discrimination, we find that LHASSSO is sensitive to the signatures from decaying DM particles above $\mathcal{O}(1)$ TeV. The LHAASO sensitivity to the DM decay lifetime reaches $\mathcal{O} (10^{26}) \sim \mathcal{O} (10^{28})$ s for several decay channels at the DM mass scale from 1 TeV to 100 TeV.

The hardening and softening features in the DAMPE proton spectrum are very likely to be originated from a nearby supernova remnant (SNR). The proton spectrum from the nearby SNR is required to be very hard below $\approx10$ TeV. To reproduce this feature, we illustrate that anomalously slow-diffusion zone for cosmic rays (CRs) must be existed in the local interstellar medium (ISM) after also taking the dipole anisotropy of CRs into account. Assuming that the diffusion coefficient is homogeneous in the nearby ISM, we show that the diffusion coefficient is constrained to the magnitude of $10^{26}$ cm$^2$ s$^{-1}$ when normalized to 1 GeV, which is about 100 times smaller than the average value in the Galaxy. We further discuss the spatial distribution of the slow diffusion and find two distinct possibilities. In one case, the SNR is several hundred of parsecs away from the solar system, meanwhile both the SNR and the solar system are required to be included in a large slow-diffusion zone. The homogeneous diffusion belongs to this case. In the other case, the SNR is very close with a distance of $\sim50$ pc and the slow-diffusion zone is only limited around the SNR. The required diffusion coefficient is further smaller in the latter case. This work provides a new way of studying the CR diffusion in the local ISM.

We investigate the axion like particle (ALP)-photon oscillation effect in the high energy $\gamma$-ray spectra of PG 1553+113 and PKS 2155-304 measured by Fermi-LAT and H.E.S.S.. The choice of extragalactic background light (EBL) model, which induces the attenuate effect in observed $\gamma$-ray spectra, would affect the ALP implication. For the ordinary EBL model that prefers a null hypothesis, we set constraint on the ALP-photon coupling constant at 95% C.L. as $g_{a\gamma}\lesssim 5\times 10^{-11} ~\rm{GeV}^{-1}$ for the ALP mass $\sim 10$ neV. We also consider the CIBER observation of the cosmic infrared radiation, which shows an excess at the wave wavelength of $\sim 1~\mu$m after the substraction of foregrounds. The high energy gamma-rays from extragalactic sources at high redshifts would suffer from a more significant attenuate effect caused by this excess. In this case, we find that the ALP-photon oscillation would improve the fit to the observed spectra of PKS 2155-304 and PG 1553+113 and find a favored parameter region at 95% C.L..

The Large High Altitude Air Shower Observatory (LHAASO) is a next-generation observatory for high energy gamma rays and cosmic rays with wide field of view. It will detect gamma rays with high sensitivity in the energy range from 300 GeV to 1 PeV. Therefore, it is promising for LHAASO to search for the high-energy gamma rays induced by dark matter (DM) self-annihilation in dwarf spheroidal satellite galaxies (dSphs), which are ideal objects for the DM indirect detection. In this work, we investigate the LHAASO sensitivity to DM self-annihilation signatures for 19 dSphs and take the uncertainties on the spatial DM distribution of dSphs into account. We perform a joint likelihood analysis for the 19 dSphs and find that the LHAASO sensitivity to the DM annihilation cross section will reach $\mathcal{O}(10^{-24})\sim \mathcal{O}(10^{-25})$ cm$^3$ s$^{-1}$ at the mass scale above TeV for several annihilation modes, which is larger than the canonical thermal relic cross section by a factor of 10 to 100.

Geminga pulsar is surrounded by a multi-TeV $\gamma$-ray halo radiated by the high energy electrons and positrons accelerated by the central pulsar wind nebula (PWN). The angular profile of the $\gamma$-ray emission reported by HAWC indicates an anomalously slow diffusion for the cosmic-ray electrons and positrons in the halo region around Geminga. In the paper we study the possible mechanism for the origin of the slow diffusion. At first, we consider the self-generated Alfv\'en waves due to the streaming instability of the electrons and positrons released by Geminga. However, even considering a very optimistic scenario for the wave growth, we find this mechanism DOES NOT work to account for the extremely slow diffusion at the present day if taking the proper motion of Geminga pulsar into account. The reason is straightforward as the PWN is too weak to generate enough high energy electrons and positrons to stimulate strong turbulence at the late time. We then propose an assumption that the strong turbulence is generated by the shock wave of the parent supernova remnant (SNR) of Geminga. Geminga may still be inside the SNR, and we find that the SNR can provide enough energy to generate the slow-diffusion circumstance. The TeV halos around PSR B0656+14, Vela X, and PSR J1826-1334 may also be explained under this assumption.

As the TeV halos around Geminga and PSR B0656+14 have been confirmed by HAWC, slow diffusion of cosmic rays could be general around pulsars, and the cosmic positron spectrum from pulsars could be significantly changed. As a consequence, the most likely pulsar source of the positron excess, Geminga, is no more a viable candidate under the additional constraint from Fermi-LAT. Moreover, the latest measurement by AMS-02 shows a clear cutoff in the positron spectrum, which sets a strict constraint on the age of the pulsar source. Considering these new developments we reanalyze the scenario in this work. By checking all the observed pulsars under the two-zone diffusion scenario, we propose for the first time that PSR B1055-52 is a very promising source of the positron excess. B1055-52 can well reproduce both the intensity and the high-energy cut of the AMS-02 positron spectrum, and may also explain the H.E.S.S $e^-+e^+$ spectrum around 10 TeV. Moreover, if the slow diffusion is universal in the local interstellar medium, B1055-52 will be the unique reasonable source of the AMS-02 positron spectrum among the observed pulsars.

We study the time-dependent modulation effect and derive the local interstellar spectra (LIS) for the cosmic ray (CR) proton, helium, boron and carbon. A two-dimensional modulation model including the variation of the interplanetary environment with time is adopted to describe modulation process. The propagation equation of CRs in the heliosphere is numerically solved by the package Solarprop. We derive the LIS by fitting the latest results of several experiments, including Voyager 1, PAMELA, BESS-POLARII and ACE, during low solar activity periods. We further study the modulation in the polarity reversal periods with the PAMELA proton data. We find that the rigidity dependence of the diffusion coefficient is critical to explain the modulation effect during reversal periods. Our results also indicate a power law relation between the diffusion coefficient and the magnitude of the heliospheric magnetic field (HMF) at the Earth.

We investigate the implications on the dark matter (DM) signal from the AMS-02 cosmic antiproton flux. Global fits to the data are performed under different propagation and hadronic interaction models. The uncertainties from the injection spectrum, propagation effects and solar modulation of the cosmic rays are taken into account comprehensively. Since we need to investigate extended parameter regions with multiple free parameters in the fit, the machine learning method is adopted to maintain a realistic time cost. We find all the effects considered in the fitting process interplay with each other, among which the hadronic interaction model is the most important factor affecting the result. In most hadronic interaction and CR propagation models no DM signal is found with significance larger than $2\sigma$ except that the EPOS-LHC interaction model requires a more than $3\sigma$ DM signal with DM mass around $1\,\mathrm{TeV}$. For the diffusive reacceleration propagation model there is a highly significant DM signal with mass around $100\,\mathrm{GeV}$. However, the signal becomes less than $1\sigma$ if we take a charge dependent solar modulation potential in the analysis.

The recent HAWC observations of very-high-energy $\gamma$-ray halo around Geminga and Monogem indicate a very slow diffusion of cosmic rays, which results in tiny contribution of positrons from these two pulsars to the local flux. This makes the cosmic positron excess anomaly observed by PAMELA and AMS-02 even more puzzling. However, from the Boron-to-Carbon ratio data one can infer that the average diffusion coefficient in the Galaxy should be much larger. In this work we propose a two-zone diffusion model that the diffusion is slow only in a small region around the source, outside of which the propagation is as fast as usual. We find that such a scenario can naturally explain the positron excess data with parameters even more reasonable than that in the conventional one-zone diffusion model. The reason is that during the life time of Geminga ($\sim 300$ kyr) the electrons/positrons have propagated too far away with a fast diffusion and lead to a low local flux. The slow diffusion region in the two-zone model helps to confine the electrons/positrons for a long time and lead to an enhancement of the local flux. So under the constraint of the HAWC observations, pulsars are still the probable origin of the cosmic-ray positron excess.

The Fermi-LAT detects no significant anisotropy of the cosmic-ray (CR) electrons and positrons ($e^-+e^+$) with seven years of data, which provides the strongest restriction to the $e^-+e^+$ anisotropy up to now. As next generation CR observatory, HERD is expected to have a better capability of anisotropy detection than Fermi-LAT. In this paper, we discuss several models aimed to explain the AMS-02 data by the present and future anisotropy measurements. We find that the upper limits of Fermi-LAT disfavor Vela SNR as the dominant source in sub-TeV, while other cases that remain safe under the constraint of Fermi-LAT are expected to be distinguished from each other by HERD. We then discuss the possibilities of remarkable TeV spectral features, and test the corresponding anisotropies. We find the conditions under which the TeV model can have a prominent spectral feature and avoid the constraint of Fermi-LAT at the same time. Furthermore, the expected performance of HERD is sensitive enough to detect the anisotropies of all these TeV models, and even for the case of a featureless TeV spectrum. Thus HERD may play a crucial part in the study of the origin of cosmic electrons and positrons.

Searching for gamma rays from dwarf spheroidal galaxies (dSphs) is a promising approach to detect dark matter (DM) due to the high DM densities and low baryon components in dSphs. The Fermi-LAT observations from dSphs have set stringent constraints on the velocity independent annihilation cross section. However, the constraints from dSphs may change in velocity dependent annihilation scenarios because of the different velocity dispersions in galaxies. In this work, we study how to set constraints on the velocity dependent annihilation cross section from the combined Fermi-LAT observations of dSphs with the kinematic data. In order to calculate the gamma ray flux from the dSph, the correlation between the DM density profile and velocity dispersion at each position should be taken into account. We study such correlation and the relevant uncertainty from kinematic observations by performing a Jeans analysis. Using the observational results of three ultrafaint dSphs with large J-factors, including Willman 1, Reticulum II, and Triangulum II, we set constraints on the p-wave annihilation cross section in the Galaxy as an example.

Recently, many new dwarf spheroidal satellites (dSphs) have been discovered by the Dark Energy Survey (DES). These dSphs are ideal candidates for probing for gamma-ray emissions from dark matter (DM) annihilation. However, no significant signature has been found by the Fermi-LAT dSph observations. In this work, we reanalyze the Fermi-LAT Pass 8 data from the direction of Reticulum II, where a slight excess has been reported by some previous studies. We treat Reticulum II (DES J0335.6-5403) as a spatially extended source, and and that no significant gamma-ray signature is observed. Based on this result, we set upper-limits on the DM annihilation cross section.

The cosmic antideuteron is a promising probe for the dark matter annihilation signature. In order to determine the DM signature, the background astrophysical antideuteron flux should be carefully studied. In this work we provide a new calculation of the secondary antideuteron flux, and pay special attention to the uncertainties from hadronic interaction models by using several Monte Carlo generators. The uncertainties from propagation effects are also carefully investigated for both the astrophysical background and DM annihilation signature in several scenarios, which are constrained by the latest B/C ratio measured by AMS-02. Considering these uncertainties, we find that the secondary antideuteron flux is hard to detect in the near future detectors. However, the antideuteron signature from dark matter annihilation will be detectable even considering the constraint from the AMS-02 observation of the $\bar{p}/p$ ratio.

Many studies have shown that either the nearby astrophysical source or dark matter (DM) annihilation/decay is required to explain the origin of high energy cosmic ray (CR) $e^\pm$, which are measured by many experiments, such as PAMELA and AMS-02. Recently, the Dark Matter Particle Explorer (DAMPE) collaboration has reported its first result of the total CR $e^\pm$ spectrum from $25 \,\mathrm{GeV}$ to $4.6 \,\mathrm{TeV}$ with high precision. In this work, we study the DM annihilation and pulsar interpretations of the DAMPE high energy $e^\pm$ spectrum. In the DM scenario, the leptonic annihilation channels to $\tau^+\tau^-$, $4\mu$, $4\tau$, and mixed charged lepton final states can well fit the DAMPE result, while the $\mu^+\mu^-$ channel has been excluded. In addition, we find that the mixed charged leptons channel would lead to a sharp drop at $\sim$ $\mathrm{TeV}$. However, these DM explanations are almost excluded by the observations of gamma-ray and CMB, unless some complicated DM models are introduced. In the pulsar scenario, we analyze 21 nearby known pulsars and assume that one of them is the primary source of high energy CR $e^\pm$.Considering the constraint from the Fermi-LAT observation of the $e^\pm$ anisotropy, we find that two pulsars are possible to explain the DAMPE data. Our results show that it is difficult to distinguish between the DM annihilation and single pulsar explanations of high energy $e^\pm$ with the current DAMPE result.

The DArk Matter Particle Explorer (DAMPE), a space-based high precision cosmic ray detector, has just reported the new measurement of the total electron plus positron energy spectrum up to 4.6 TeV. A notable feature in the spectrum is the spectral break at $\sim 0.9$ TeV, with the spectral index softening from $-3.1$ to $-3.9$. Such a feature is very similar to the knee at the cosmic nuclei energy spectrum. In this work we propose that the knee-like feature can be explained naturally by assuming that the electrons are accelerated at the supernova remnants (SNRs) and released when the SNRs die out with lifetime around $10^5$ years. The cut-off energy of those electrons have already decreased to several TeV due to radiative cooling, which may induce the observed TeV spectral break. Another possibility is that the break is induced by a single nearby old SNR. Such a scenario may bring a large electron flux anisotropy that is observable by the future detectors. We also show that a minor part of electrons escaping during the acceleration in young and nearby SNRs are able to contribute to several TeV or higher energy region of the spectrum.

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed.

The cosmic-ray (CR) $e^\pm$ excess observed by AMS-02 can be explained by dark matter (DM) annihilation. However, the DM explanation requires a large annihilation cross section which is strongly disfavored by other observations, such as the Fermi-LAT gamma-ray observation of dwarf galaxies and the Planck observation of the cosmic microwave background (CMB). Moreover, the DM annihilation cross section required by the CR $e^\pm$ excess is also too large to generate the correct DM relic density with thermal production. In this work we use the Breit-Wigner mechanism with a velocity dependent DM annihilation cross section to reconcile these tensions. If DM particles accounting for the CR $e^\pm$ excess with $v\sim \mathcal{O}(10^{-3})$ are very close to a resonance in the physical pole case, their annihilation cross section in the Galaxy reaches a maximal value. On the other hand, the annihilation cross section would be suppressed for DM particles with smaller relative velocities in dwarf galaxies and at recombination, which may affect the gamma-ray and CMB observations, respectively. We find a proper parameter region that can simultaneously explain the AMS-02 results and the thermal relic density, while satisfying the Fermi-LAT and Planck constraints.

Recently the AMS-02 collaboration has published the measurement of the cosmic antiproton to proton ratio $\bar{p}/p$ and the $\bar{p}$ flux with a high precision up to $\sim 450\,\mathrm{GeV}$. In this work, we perform a systematic analysis of the secondary antiproton flux generated by the cosmic ray interaction with the interstellar gas. The uncertainty of the prediction originates from the cosmic ray propagation process and the hadronic interaction models. Although the cosmic ray propagation parameters have been well controlled by the AMS-02 $B/C$ ratio data for a specified model, different propagation models can not be discriminated by the $B/C$ data. The $\bar{p}$ flux is also calculated for several hadronic interaction models, which are generally adopted by the cosmic ray community. However, the results for different hadronic models do not converge. We find the EPOS LHC model, which seems to fit the collider data very well, predicts a slightly lower $\bar{p}/p$ ratio than the AMS-02 data at the high energy end. Finally we derive the constraints on the dark matter annihilation cross section from the AMS-02 $\bar{p}/p$ ratio for different propagation and hadronic interaction models.

The AMS-02 has measured the cosmic ray electron (plus positron) spectrum up to ~TeV with an unprecedent precision. The spectrum can be well described by a power law without any obvious features above 10 GeV. The satellite instrument Dark Matter Particle Explorer (DAMPE), which was launched a year ago, will measure the electron spectrum up to 10 TeV with a high energy resolution. The cosmic electrons beyond TeV may be attributed to few local cosmic ray sources, such as supernova remnants. Therefore, spectral features, such as cutoff and bumps, can be expected at high energies. In this work we give a careful study on the perspective of the electron spectrum beyond TeV. We first examine our astrophysical source models on the latest leptonic data of AMS-02 to give a self-consistent picture. Then we focus on the discussion about the candidate sources which could be electron contributors above TeV. Depending on the properties of the local sources (especially on the nature of Vela), DAMPE may detect interesting features in the electron spectrum above TeV in the future.

Growing evidence reveals universal hardening on various cosmic ray spectra, e.g. proton, positron, as well as antiproton fraction. Such universality may indicate they have a common origin. In this paper, we argue that these widespread excesses can be accounted for by a nearby supernova remnant surrounded by a giant molecular cloud. Secondary cosmic rays ($\rm p$, $\rm e^+$) are produced through the collisions between the primary cosmic ray nuclei from this supernova remnant and the molecular gas. Different from the background, which is produced by the ensemble of large amount of sources in the Milky Way, the local injected spectrum can be harder. The time-dependent transport of particles would make the propagated spectrum even harder. Under this scenario, the anomalies of both primary ($\rm p$, $\rm e^-$) and secondary ($\rm e^+$, $\rm \bar{p}/p$) cosmic rays can be properly interpreted. We further show that the TeV to sub-PeV anisotropy of proton is consistent with the observations if the local source is relatively young and lying at the anti-Galactic center direction.

The gamma-ray observation of dwarf spheroidal satellites (dSphs) is an ideal approach for probing the dark matter (DM) annihilation signature. The latest Fermi-LAT dSph searches have set stringent constraints on the velocity independent annihilation cross section in the small DM mass range, which gives very strong constraints on the scenario to explain the the AMS-02 positron excess by DM annihilation. However, the dSph constraints would change in the velocity dependent annihilation scenarios, because the velocity dispersion in the dSphs varies from that in the Milky Way. In this work, we use a likelihood map method to set constraints on the velocity dependent annihilation cross section from the Fermi-LAT observation of six dSphs. We consider three typical forms of the annihilation cross section, i.e. p-wave annihilation, Sommerfeld enhancement, and Breit-Wigner resonance. For the p-wave annihilation and Sommerfeld-enhancement, the dSph limits would become much weaker and stronger compared with those for the velocity independent annihilation, respectively. For the Breit-Wigner annihilation, the dSph limits would vary depending on the model parameters. We show that the scenario to explain the AMS-02 positron excess by DM annihilation is still viable in the velocity dependent cases.

In this work we study the constraints on dark matter (DM) annihilation/decay from the Fermi-LAT Isotropic Gamma-Ray Background (IGRB) observation. We consider the contributions from both extragalactic and galactic DM components. For DM annihilation, the evolutions of extragalactic DM halos are taken into account. We find that the IGRB constraints under some DM subhalo models can be comparable to those derived from the observations of dwarf spheroidal galaxies. We also use the IGRB results to constrain the parameter regions accounting for the latest AMS-02 electron-positron anomaly. We find that the majority of DM annihilation/decay channels are strongly disfavored by the latest Fermi-LAT IGRB observation; only DM annihilation/decay to $\mu^+\mu^-$ may be valid.

The AMS-02 collaboration has just released the cosmic antiproton to proton ratio $\bar{p}/p$ with a high precision up to $\sim 450$ GeV. In this work, we calculate the secondary antiprotons generated by cosmic ray interactions with the interstellar medium taking into account the uncertainties from the cosmic ray propagation. The $\bar{p}/p$ ratio predicted by these processes shows some tension with the AMS-02 data in some regions of propagation parameters, but the excess is not significant. We then try to derive upper bounds on the dark matter annihilation cross section from the $\bar{p}/p$ data or signal regions favored by the data. It is shown that the constraint derived by the AMS-02 data is similar to that from Fermi-LAT observations of dwarf galaxies. The signal region for dark matter is usually required $m_\chi \sim O(10)$ TeV and $\left<\sigma v\right>\sim\mathcal{O}(10^{-23})~\cm^3~\sec^{-1}$.

The gamma-ray excess observed by the Fermi-LAT in the Galactic Center can be interpreted by the dark matter annihilation to $b\bar{b}$ via a light pseudoscalar in the NMSSM. It is interesting to note that the corresponding singlet scalar is useful to achieve a strongly first order phase transition required by the electroweak baryogenesis. In this paper, we investigate the possibility that the NMSSM model can simultaneously accommodate these two issues. The phase transition strength can be characterized by the vacua energy gap at zero temperature and be sufficiently enhanced by the tree-level effect in the NMSSM. We find that the annihilation of Singlino/Higgsino DM particles occurring close to the light pseudoscalar resonance is favored by the galactic center excess and the observed DM relic density, and a resulting small $\kappa/\lambda$ and a negative $A_\kappa$ can also lead to a successful strongly first order electroweak phase transition.

The three-body force (TBF) effect on the neutrino emissivity in neutron star matter and the total neutrino emissivity of neutron stars have been investigated within the framework of the Brueckner-Hartree-Fock approach by adopting the AV18 two-body interaction plus a microscopic TBF. The neutrino emissivity from the direct Urca process turns out to be much larger than that from the modified Urca process. Inclusion of the TBF reduces strongly the density thresholds of the direct Urca processes involving electrons and muons. The TBF effect on the total neutrino emissivity of neutron stars is shown to be negligibly weak for neutron stars with small masses. For neutron stars with large masses, the TBF effect becomes visible and inclusion of the TBF may enhance the total neutrino emissivity by about 50% for neutron stars with a given mass of $M=1.6M_{\odot}$.

The AMS-02 collaboration has just published a high precision measurement of the cosmic positron fraction $e^+/(e^- + e^+)$, which rises with energy from $\sim 5$ GeV to $\sim 350$ GeV. The result indicates the existence of primary electron/positron sources to account for the positron excess. In this work, we investigate the possibility that the nearby mature pulsars are the primary positron sources. By fitting the data we find that the positrons from a single nearby pulsar, such as Geminga or Monogem, with the spectral index $\alpha \sim 2$ can interpret the AMS-02 result. We also investigate the possibility that high energy positrons are generated by multiple known pulsars in the ATNF catalogue. Such a scenario can also fit the AMS-02 data well. Future precise measurements of fine structures in the positron spectrum would be a support to the pulsar scenario.

We have investigated the three-body force (TBF) effect on the neutron and proton momentum distributions in asymmetric nuclear matter within the framework of the extended Brueckner-Hartree-Fock approach by adopting the $AV18$ two-body interaction plus a microscopic TBF. In asymmetric nuclear matter, it is shown that the neutron and proton momentum distributions become different from their common distribution in symmetric nuclear matter. The predicted depletion of the proton hole states increases while the neutron one decreases as a function of isospin-asymmetry. The TBF effect on the neutron and proton momentum distributions turns out to be negligibly weak at low densities around and below the normal nuclear density. The TBF effect is found to become sizable only at high densities well above the saturation density, and inclusion of the TBF leads to an overall enhancement of the depletion of the neutron and proton Fermi seas.

The weakly interacting massive particle (WIMP) often serves as a candidate for the cold dark matter, however when produced non-thermally it could behave like warm dark matter. In this paper we study the properties of the $\gamma$-ray emission from annihilation of WIMP dark matter in the halo of our own Milky-Way Galaxy with high resolution $N$-body simulations of a Milky-Way like dark matter halo, assuming different nature of WIMPs. Due to the large free-streaming length in the scenario of warm WIMPs, the substructure contend of the dark matter halo is significantly different from that of the cold WIMP counterpart, resulting in distinct predictions of the $\gamma$-ray signals from the dark matter annihilation. We illustrate these by comparing the predicted $\gamma$-ray signals from the warm WIMP annihilation to that of cold WIMPs. Pronounced differences from the subhalo skymap and statistical properties between two WIMP models are demonstrated. Due to the potentially enhanced cross section of the non-thermal production mechanism in warm WIMP scenario, the Galactic center might be prior for the indirect detection of warm WIMPs to dwarf galaxies, which might be different from the cold dark matter scenario. As a specific example we consider the non-thermally produced neutralino of supersymmetric model and discuss the detectability of warm WIMPs with Fermi $\gamma$-ray telescope.

It was recently reported that there may exist monochromatic $\gamma$-ray emission at $\sim 130$ GeV from the Galactic center in the Fermi Large Area Telescope data, which might be related with dark matter (DM) annihilation. In this work we carry out a comprehensive check of consistency of the results with the DM annihilation scenario, using the 3.7 yrs Fermi observation of the inner Galaxy, Galactic halo, clusters of galaxies and dwarf galaxies. The results found are as follows. 1) Very strong constraints on the DM annihilation into continuous $\gamma$-rays from the Galactic center are set, which are as stringent as the "natural" scale assuming thermal freeze-out of DM. Such limit sets strong constraint on the DM models to explain the line emission. 2) No line emission from the Galactic halo is found in the Fermi data, and the constraints on line emission is marginally consistent with the DM annihilation interpretation of the $\sim 130$ GeV line emission from the inner Galaxy. 3) No line emission from galaxy clusters and dwarf galaxies is detected, although possible concentration of photons from clusters in 120-140 GeV is revealed. The constraints from clusters and dwarf galaxies are weak and consistent with the DM annihilation scenario to explain the $\sim 130$ GeV line emission.

The globular clusters are probably good targets for dark matter (DM) searches in $\gamma$-rays due to the possible adiabatic contraction of DM by baryons. In this work we analyse the three-year data collected by {\it Fermi} Large Area Telescope of globular clusters NGC 6388 and M 15 to search for possible DM signals. For NGC 6388 the detection of $\gamma$-ray emission was reported by {\it Fermi} collaboration, which is consistent with the emission of a population of millisecond pulsars. The spectral shape of NGC 6388 is also shown to be consistent with a DM contribution if assuming the annihilation final state is $b\bar{b}$. No significant $\gamma$-ray emission from M 15 is observed. We give the upper limits of DM contribution to $\gamma$-ray emission in both NGC 6388 and M 15, for annihilation final states $b\bar{b}$, $W^+W^-$, $\mu^+\mu^-$, $\tau^+\tau^-$ and monochromatic line. The constraints are stronger than that derived from observation of dwarf galaxies by {\it Fermi}.

The superluminal neutrinos detected by OPERA indicates Lorentz invariance violation (LIV) of the neutrino sector at the order of $10^{-5}$. We study the implications of the result in this work. We find that such a large LIV implied by OPERA data will make the neutrino production process $\pi \to \mu + \nu_\mu$ kinematically forbidden for neutrino energy greater than about 5 GeV. The OPERA detection of neutrinos at 40 GeV can constrain the LIV parameter to be smaller than $3\times 10^{-7}$. Furthermore the neutrino decay in the LIV framework will modify the neutrino spectrum greatly. The atmospheric neutrino spectrum measured by IceCube can constrain the LIV parameter to the level of $10^{-12}$. The future detection of astrophysical neutrinos of Galactic sources is expected to be able to give even stronger constraint on the LIV parameter of neutrinos.

Recently the AGILE and Fermi/LAT detectors uncovered giant $\gamma$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $\gamma$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

In the paper we study the Breit-Wigner enhancement of dark matter (DM) annihilation considering the kinetic decoupling in the evolution of DM freeze-out at the early universe. Since the DM temperature decreases much faster (as $1/R^2$) after kinetic decoupling than that in kinetic equilibrium (as 1/R) we find the Breit-Wigner enhancement of DM annihilation rate after the kinetic decoupling will affect the DM relic density significantly. Focusing on the model parameters that trying to explain the anomalous cosmic positron/electron excesses observed by PAMELA/Fermi/ATIC we find the elastic scattering $Xf\to Xf$ is not efficient to keep dark matter in kinetic equilibrium, and the kinetic decoupling temperature $T_{kd}$ is comparable to the chemical decoupling temperature $T_f\sim O(10) GeV$. The reduction of the relic density after $T_{kd}$ is significant and leads to a limited enhancement factor $\sim O(10^2)$. Therefore it is difficult to explain the anomalous positron/electron excesses in cosmic rays by DM annihilation and give the correct DM relic density simultaneously in the minimal Breit-Wigner enhancement model.

The Fermi $\gamma$-ray space telescope reported the observation of several Galactic supernova remnants recently, with the $\gamma$-ray spectra well described by hadronic $pp$ collisions. The possible neutrino emissions from these Fermi detected supernova remnants are discussed in this work, assuming the hadronic origin of the $\gamma$-ray emission. The muon event rates induced by the neutrinos from these supernova remnants on typical km$^3$ neutrino telescopes, such as the IceCube and the KM3NeT, are calculated. The results show that for most of these supernova remnants the neutrino signals are too weak to be detected by the on-going or up-coming neutrino experiment. Only for the TeV bright sources RX J1713.7-3946 and possibly W28 the neutrino signals can be comparable with the atmospheric background in the TeV region, if the protons can be accelerated to very high energies. The northern hemisphere based neutrino telescope might detect the neutrinos from these two sources.

A statistical scenario is proposed to explain the $\gamma$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $\gamma$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $\gamma$-ray flares over a period of $\sim$ year and the monthly scale $\gamma$-ray flux fluctuations as observed by the Fermi/LAT. The observed $\gamma$-ray spectra in both the steady and flaring states are also well reproduced.

The $\gamma$-ray and neutrino emissions from dark matter (DM) annihilation in galaxy clusters are studied. After about one year operation of Fermi-LAT, several nearby clusters are reported with stringent upper limits of GeV $\gamma$-ray emission. We use the Fermi-LAT upper limits of these clusters to constrain the DM model parameters. We find that the DM model distributed with substructures predicted in cold DM (CDM) scenario is strongly constrained by Fermi-LAT $\gamma$-ray data. Especially for the leptonic annihilation scenario which may account for the $e^{\pm}$ excesses discovered by PAMELA/Fermi-LAT/HESS, the constraint on the minimum mass of substructures is of the level $10^2-10^3$ M$_{\odot}$, which is much larger than that expected in CDM picture, but is consistent with a warm DM scenario. We further investigate the sensitivity of neutrino detections of the clusters by IceCube. It is found that neutrino detection is much more difficult than $\gamma$-rays. Only for very heavy DM ($\sim 10$ TeV) together with a considerable branching ratio to line neutrinos the neutrino sensitivity is comparable with that of $\gamma$-rays.

We investigate the possibility of detecting light long-lived particle (LLP) produced by high energy cosmic ray colliding with atmosphere. The LLP may penetrate the atmosphere and decay into a pair of muons near/in the neutrino telescope. Such muons can be treated as the detectable signal for neutrino telescope. This study is motivated by recent cosmic electron/positron observations which suggest the existence of $O(TeV)$ dark matter and new light $O(GeV)$ particle. It indicates that dark sector may be complicated, and there may exist more than one light particles, for example the dark gauge boson $A'$ and associated dark Higgs boson $h'$. In this work, we discuss the scenario with $A'$ heavier than $h'$ and $h'$ is treated as LLP. Based on our numerical estimation, we find that the large volume neutrino telescope IceCube has the capacity to observe several tens of di-muon events for favorable parameters if the decay length of LLP can be comparable with the depth of atmosphere. The challenge here is how to suppress the muon backgrounds induced by cosmic rays and atmospheric neutrinos.

Boost factors of dark matter annihilation into antiprotons and electrons/positrons due to the clumpiness of dark matter distribution are studied in detail in this work, taking the Sommerfeld effect into account. It has been thought that the Sommerfeld effect, if exists, will be more remarkable in substructures because they are colder than the host halo, and may result in a larger boost factor. We give a full calculation of the boost factors based on the recent N-body simulations. Three typical cases of Sommerfeld effects, the non-resonant, moderately resonant and strongly resonant cases are considered. We find that for the non-resonant and moderately resonant cases the enhancement effects of substructures due to the Sommerfeld effect are very small ($\sim \mathcal{O}(1)$) because of the saturation behavior of the Sommerfeld effect. For the strongly resonant case the boost factor is typically smaller than $\sim \mathcal{O}(10)$. However, it is possible in some very extreme cases that DM distribution is adopted to give the maximal annihilation the boost factor can reach up to $\sim 1000$. The variances of the boost factors due to different realizations of substructures distribution are also discussed in the work.

The excesses of the cosmic positron fraction recently measured by PAMELA and the electron spectra by ATIC, PPB-BETS, Fermi and H.E.S.S. indicate the existence of primary electron and positron sources. The possible explanations include dark matter annihilation, decay, and astrophysical origin, like pulsars. In this work we show that these three scenarios can all explain the experimental results of the cosmic $e^\pm$ excess. However, it may be difficult to discriminate these different scenarios by the local measurements of electrons and positrons. We propose possible discriminations among these scenarios through the synchrotron and inverse Compton radiation of the primary electrons/positrons from the region close to the Galactic center. Taking typical configurations, we find the three scenarios predict quite different spectra and skymaps of the synchrotron and inverse Compton radiation, though there are relatively large uncertainties. The most prominent differences come from the energy band $10^4\sim 10^9$ MHz for synchrotron emission and $\gtrsim 10$ GeV for inverse Compton emission. It might be able to discriminate at least the annihilating dark matter scenario from the other two given the high precision synchrotron and diffuse $\gamma$-ray skymaps in the future.

Recent PAMELA data show that positron fraction has an excess above several GeV while anti-proton one is not. Moreover ATIC data indicates that electron/positron flux have a bump from 300 GeV to 800 GeV. Both annihilating dark matter (DM) with large boost factor and decaying DM with the life around $ 10^{26} s$ can account for the PAMELA and ATIC observations if their main final products are charged leptons ($e$, $\mu$ and $\tau$). In this work, we calculated the neutrino flux arising from $\mu$ and $\tau$ which originate from annihilating/decaying DM, and estimated the final muon rate in the neutrino telescopes, namely Antares and IceCube. Given the excellent angular resolution, Antares and IceCube are promising to discover the neutrino signals from Galactic center and/or large DM subhalo in annihilating DM scenario, but very challenging in decaying DM scenario.

Recently PAMELA released their first results on the positron and antiproton ratios. Stimulated by the new data, we studied the cosmic ray propagation models and calculated the secondary positron and antiproton spectra. The low energy positron ratio can be consistent with data in the convection propagation model. Above $\sim 10$ GeV PAMELA data shows a clear excess on the positron ratio. However, the secondary antiproton is roughly consistent with data. The positron excess may be a direct evidence of dark matter annihilation or decay. We compare the positron and anti-proton spectra with data by assuming dark matter annihilates or decays into different final states. The PAMELA data actually excludes quark pairs being the main final states, disfavors gauge boson final states. Only in the case of leptonic final states the positron and anti-proton spectra can be explained simultaneously. We also compare the decaying and annihilating dark matter scenarios to account for the PAMELA results and prefer to the decaying dark matter. Finally we consider a decaying neutralino dark matter model in the frame of supersymmetry with R-parity violation. The PAMELA data is well fitted with neutralino mass $600\sim 2000$ GeV and life time $\sim 10^{26}$ seconds. We also demonstrate that neutralino with mass around 2TeV can fit PAMELA and ATIC data simultaneously.

Dark matter (DM) annihilations in the Galaxy may produce high energy neutrinos, which can be detected by the neutrino telescopes, for example IceCube, ANTARES and Super-Kamiokande. The neutrinos can also arise from hadronic interaction between cosmic ray and atmosphere around the Earth, known as atmospheric neutrino. Current measurements on neutrino flux is consistent with theoretical prediction of atmospheric neutrino within the uncertainties. In this paper, by requiring that the DM annihilation neutrino flux is less than the current measurements, we obtain an upper bound on the cross section of dark matter annihilation $ < {\sigma v} >$. Compared with previous investigations, we improve the bound by including DM substructure contributions. In our paper, two kinds of substructure effects are scrutinized. One is the substructure average contribution over all directions. The other is point source effect by single massive sub-halo. We found that the former can improve the bound by several times, while the latter can improve the bound by $ 10^1 \sim 10^4$ utilizing the excellent angular resolution of neutrino telescope IceCube. The exact improvement depends on the DM profile and the sub-halo concentration model. In some model, IceCube can achieve the sensitivity of $ < {\sigma v} > \sim 10^{- 26} cm^3 s^{- 1} $.

The lightest neutralino, as the dark matter candidate, can be gravitationally captured by the Sun. In this paper, we studied the high energy neutrino signals from solar neutralino annihilations in the core of the Sun in the anomaly mediated supersymmetry (SUSY) breaking (AMSB) model. Based on the event-by-event monte carlo simulation code WimpSim, we studied the detailed energy and angular spectrum of the final muons at large neutrino telescope IceCube. More precisely we simulated the processes since the production of neutrino via neutralino annihilation in the core of the Sun, neutrino propagation from the Sun to the Earth, as well as the converting processes from neutrino to muon. Our results showed that in the AMSB model it is possible to observe the energetic muons at IceCube, provided that the lightest neutralio has relatively large higgsino component, as a rule of thumb $ N_{13}^2 + N_{14}^2 > 4%$ or equivalently $ \sigma_{SD} > 10^{- 5} pb$. Especially, for our favorable parameters the signal annual events can reach 102 and the statistical significance can reach more than 20. We pointed out that the energy spectrum of muons may be used to distinguish among the AMSB model and other SUSY breaking scenarios.