Recently the AMS-02 reported the precise measurements of the energy spectra of medium-mass compositions (Neon, Magnesium, Silicon) of primary cosmic rays, which reveal different properties from those of light compositions (Helium, Carbon, Oxygen). Here we propose a nearby source scenario, together with the background source contribution, to explain the newly measured spectra of cosmic ray Ne, Mg, Si, and particularly their differences from that of He, C, O. Their differences at high energies can be naturally accounted for by the element abundance of the nearby source. Specifically, the abundance ratio of the nearby source to the background of the Ne, Mg, Si elements is lower by a factor of $\sim1.7$ than that of the He, C, O elements. Such a difference could be due to the abundance difference of the stellar evolution of the progenitor star or the acceleration process/environment, of the nearby source. This scenario can simultaneously explain the high-energy spectral softening features of cosmic ray spectra revealed recently by CREAM/NUCLEON/DAMPE, as well as the energy-dependent behaviors of the large-scale anisotropies. It is predicted that the dipole anisotropy amplitudes below PeV energies of the Ne, Mg, Si group are smaller than that of the He, C, O group, which can be tested with future measurements.

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.

The observation of cosmic sub-PeV gamma-rays from the Crab Nebula opens up the possibility of testing cosmic ray photon transparency at the multi-hundred TeV scale. Assuming no deviation from a source gamma-ray emission due to accelerated electron inverse-Compton scattering, higher event energies can extend constraints on the effects of new physics; We consider oscillation between gamma-rays and axions/dark photons, plus attenuation effects from gamma-ray absorption in the case of dark photon dark matter. Combining the recent AS$\gamma$ and HAWC sub-PeV data with earlier MAGIC and HEGRA data, axion-like particles are most constrained in the $10^{-7}-10^{-6}$ eV mass range, where the coupling $g_{a\gamma}$ is constrained to be below ${1.8}\times 10^{-10}$ GeV$^{-1}$. In comparison, gamma ray flux attenuation due to oscillation with a dark photon leads to a very weak constraint on the mixing parameter; $\epsilon\lesssim$ 0.2 for dark photon mass between $10^{-7}$ and $10^{-6}$ eV. Direct scattering from dark photon dark matter limits $\epsilon\lesssim$ 0.01 for masses between $6$ and $400$ eV.

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

We study the reasonable requirements of two anomalous $U(1)$s in a flavored-axion framework for the anomaly cancellations of both $U(1)$-mixed gravity and $U(1)_Y\times[U(1)]^2$ which in turn determine the $U(1)_Y$ charges where $U(1)_Y$ is the hypercharge gauge symmetry of the standard model. We argue that, with a flavor symmetry group, axion-induced topology in symmetry-broken phases plays crucial roles in describing how quarks and leptons are organized at a fundamental level and make deep connections with each other. A unified model, as an example, is then proposed in a simple way to describe a whole spectrum of particles where both flavored-axion interactions with normal matter and the masses and mixings of fermions emerge from the spontaneous breaking of a given symmetry group. Once a scale of active neutrino mass defined at a seesaw scale is fixed by the commensurate $U(1)$ flavored-PQ charges of fermions, that of QCD axion decay constant $F_A$ is determined. In turn, fundamental physical parameters complementary to each other are predicted with the help of precision flavor experiments and astrophysical observations. Model predictions are extracted on the characteristics of neutrino and flavored-axion as a good candidate of dark matter: $F_A=3.05^{\,+1.01}_{\,-1.01}\times10^{10}$ GeV (consequently, QCD axion mass $m_a=1.78^{+0.88}_{-0.44}\times10^{-4}$ eV, and axion to electron coupling $|g_{a\gamma\gamma}|=2.51^{+0.62}_{-1.24}\times10^{-14}\,\text{GeV}^{-1}$, etc.); atmospheric mixing angle $\theta_{23}$, Dirac CP phase $\delta_{CP}$, and $0\nu\beta\beta$-decay rate for normal mass ordering and inverted one by taking quantum corrections into account. Future $0\nu\beta\beta$ decay searches and/or precise measurement of $\theta_{23}$ will exclude or favor the example model, giving predictions on $\delta_{CP}$ and mass orderings.

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.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

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 report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

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 hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

The precise observations of Galactic cosmic ray fluxes of the secondary family, such as Li, Be, B, are expected to have significant implications on our understanding of the cosmic ray origin and propagation. Here we employ the recent very precise measurements of those species by the Alpha Magnetic Spectrometer on the International Space Station, together with their parent species (C and O), as well as the data collected by the Voyager-1 spacecraft outside the heliosphere and the Advanced Composition Explorer, we investigate the propagation of cosmic rays in the Milky Way. We find that these new data strongly favor the existence of reacceleration of cosmic rays during the propagation rather than the convective transport. We further find that for the reacceleration model, the slope of rigidity dependence of the diffusion coefficient is close to the Kolmogorov-type of turbulence, and the spectral hardenings of both the primary and secondary particles can be well described by the injection hardening rather than the propagation hardening.

We build a comprehensive sample to statistically describe the properties of X-ray flashes (XRFs) and X-ray riches (XRRs) from the latest third Swift Burst Alert Telescope (BAT3) catalog of Gamma-ray bursts (GRBs). We obtain 81 XRFs, 540 XRRs, and 394 classical GRBs (C-GRBs). We statistically explore the different properties of the $\gamma$-ray prompt emission, the X-ray emission, the X-ray lightcurve type, the association with supernovae (SNe), and the host galaxy properties for these sources. We confirm that most XRFs/XRRs are long GRBs with low values of peak energy $E^{obs}_{peak}$ and they are low-luminosity GRBs. XRFs, XRRs, and C-GRBs follow the same $E_{X,iso}$-$E_{\gamma,iso}$-$E_{peak,z}$ correlations. Compared to the classical GRBs, XRFs are favorable to have the association with SN explosions. We do not find any significant differences of redshift distribution and host galaxy properties among XRFs, XRRs, and C-GRBs. We also discuss some observational biases and selection effects that may affect on our statistical results. The GRB detectors with wide energy range and low energy threshold are expected for the XRF/XRR research in the future.

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.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

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.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

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 remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

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

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

We propose an interesting framework for asymmetric scalar dark matter (ADM), which has novel collider phenomenology in terms of an unstable ADM bound state (ADMonium) produced via Higgs portals. ADMonium is a natural consequence of the basic features of ADM: the (complex scalar) ADM is charged under a dark local $U(1)_d$ symmetry which is broken at a low scale and provides a light gauge boson $X$. The dark gauge coupling is strong and then ADM can annihilate away into $X$-pair effectively. Therefore, the ADM can form bound state due to its large self-interaction via $X$ mediation. To explore the collider signature of ADMonium, we propose that ADM has a two-Higgs doublet portal. The ADMonium can have a sizable mixing with the heavier Higgs boson, which admits a large cross section of ADMonium production associated with $b\bar b$. The resulting signature at the LHC depends on the decays of $X$. In this paper we consider a case of particular interest: $pp\ra b\bar b+ {\rm ADMonium}$ followed by ${\rm ADMonium}\ra 2X\ra 2e^+e^-$ where the electrons are identified as (un)converted photons. It may provide a competitive explanation to heavy di-photon resonance searches at the LHC.

In this work we use the newly reported Boron-to-Carbon ratio (B/C) from AMS-02 and the time-dependent proton fluxes from PAMELA and AMS-02 to constrain the source and propagation parameters of cosmic rays in the Milky Way. A linear correlation of the solar modulation parameter with solar activities is assumed to account for the time-varying cosmic ray fluxes. A comprehensive set of propagation models, with/without reacceleration or convection, have been discussed and compared. We find that only the models with reacceleration can self-consistently fit both the proton and B/C data. The rigidity dependence slope of the diffusion coefficient, $\delta$, is found to be about $0.38-0.50$ for the diffusion-reacceleration models. The plain diffusion and diffusion-convection models fit the data poorly. We compare different model predictions of the positron and antiproton fluxes with the data. We find that the diffusion-reacceleration models over-produce low energy positrons, while non-reacceleration models give better fit to the data. As for antiprotons, reacceleration models tend to under-predict low energy antiproton fluxes, unless a phenomenological modification of the velocity-dependence of the diffusion coefficient is applied. Our results suggest that there could be important differences of the propagation for nuclei and leptons, in either the Milky Way or the solar heliosphere.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

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.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

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.

We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to gamma-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation time, from August 2008 to February 2013. In particular, thanks to the ARGO-YBJ and Fermi data, the whole energy range from 100 MeV to 10 TeV is covered without any gap. In the observation period, Mrk 421 showed both low and high activity states at all wavebands. The correlations among flux variations in different wavebands were analyzed. Seven large flares, including five X-ray flares and two GeV gamma-ray flares with variable durations (3-58 days), and one X-ray outburst phase were identified and used to investigate the variation of the spectral energy distribution with respect to a relative quiescent phase. During the outburst phase and the seven flaring episodes, the peak energy in X-rays is observed to increase from sub-keV to few keV. The TeV gamma-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV gamma-rays is found to vary depending on the flare, a feature that leads us to classify flares into three groups according to the GeV flux variation. Finally, the one-zone synchrotron self-Compton model was adopted to describe the emission spectra. Two out of three groups can be satisfactorily described using injected electrons with a power-law spectral index around 2.2, as expected from relativistic diffuse shock acceleration, whereas the remaining group requires a harder injected spectrum. The underlying physical mechanisms responsible for different groups may be related to the acceleration process or to the environment properties.

The events recorded by ARGO-YBJ in more than five years of data collection have been analyzed to determine the diffuse gamma-ray emission in the Galactic plane at Galactic longitudes 25{\deg} < l < 100{\deg} and Galactic latitudes . The energy range covered by this analysis, from ~350 GeV to ~2 TeV, allows the connection of the region explored by Fermi with the multi-TeV measurements carried out by Milagro. Our analysis has been focused on two selected regions of the Galactic plane, i.e., 40{\deg} < l < 100{\deg} and 65{\deg} < l < 85{\deg} (the Cygnus region), where Milagro observed an excess with respect to the predictions of current models. Great care has been taken in order to mask the most intense gamma-ray sources, including the TeV counterpart of the Cygnus cocoon recently identified by ARGO-YBJ, and to remove residual contributions. The ARGO-YBJ results do not show any excess at sub-TeV energies corresponding to the excess found by Milagro, and are consistent with the predictions of the Fermi model for the diffuse Galactic emission. From the measured energy distribution we derive spectral indices and the differential flux at 1 TeV of the diffuse gamma-ray emission in the sky regions investigated.

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 ARGO-YBJ experiment has been in stable data taking from November 2007 till February 2013 at the YangBaJing Cosmic Ray Observatory (4300 m a.s.l.). The detector consists of a single layer of Resistive Plate Chambers (RPCs) ( about 6700 m^2}) operated in streamer mode. The signal pick-up is obtained by means of strips facing one side of the gas volume. The digital readout of the signals, while allows a high space-time resolution in the shower front reconstruction, limits the measurable energy to a few hundred TeV. In order to fully investigate the 1-10 PeV region, an analog readout has been implemented by instrumenting each RPC with two large size electrodes facing the other side of the gas volume. Since December 2009 the RPC charge readout has been in operation on the entire central carpet (about 5800 m^2). In this configuration the detector is able to measure the particle density at the core position where it ranges from tens to many thousands of particles per m^2. Thus ARGO-YBJ provides a highly detailed image of the charge component at the core of air showers. In this paper we describe the analog readout of RPCs in ARGO-YBJ and discuss both the performance of the system and the physical impact on the EAS measurements.

The ARGO-YBJ experiment is a full-coverage air shower detector located at the Yangbajing Cosmic Ray Observatory (Tibet, People's Republic of China, 4300 m a.s.l.). The high altitude, combined with the full-coverage technique, allows the detection of extensive air showers in a wide energy range and offer the possibility of measuring the cosmic ray proton plus helium spectrum down to the TeV region, where direct balloon/space-borne measurements are available. The detector has been in stable data taking in its full configuration from November 2007 to February 2013. In this paper the measurement of the cosmic ray proton plus helium energy spectrum is presented in the region 3-300 TeV by analyzing the full collected data sample. The resulting spectral index is $\gamma = -2.64 \pm 0.01$. These results demonstrate the possibility of performing an accurate measurement of the spectrum of light elements with a ground based air shower detector.

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 operation of AMS-02 opens a new era for the study of cosmic ray physics with unprecedentedly precise data which are comparable with the laboratory measurements. The high precision data allow a quantitative study on the cosmic ray physics and give strict constraints on the nature of cosmic ray sources. However, the intrinsic errors from the theoretical models to interpret the data become dominant over the errors in the data. In the present work we try to give a systematic study on the uncertainties of the models to explain the AMS-02 positron fraction data, which shows the cosmic ray $e^+e^-$ excesses together with the PAMELA and Fermi-LAT measurements. The excesses can be attributed to contributions from the extra $e^+e^-$ sources, such as pulsars or the dark matter annihilation. The possible systematic uncertainties of the theoretical models considered include the cosmic ray propagation, the treatment of the low energy data, the solar modulation, the $pp$ interaction models, the nuclei injection spectrum and so on. We find that in general a spectral hardening of the primary electron injection spectrum above $\sim50-100$ GeV is favored by the data. Furthermore, the present model uncertainties may lead to a factor of $\sim2$ enlargement in the determination of the parameter regions of the extra source, such as the dark matter mass, annihilation rate and so on.

The search for Gamma Ray Burst (GRB) emission in the energy range 1-100 GeV in coincidence with the satellite detection has been carried out using the Astrophysical Radiation with Ground-based Observatory at YangBaJing (ARGO-YBJ) experiment. The high altitude location (4300 m a.s.l.), the large active surface ($\sim$ 6700 m$^2$ of Resistive Plate Chambers), the wide field of view ($\sim 2~$sr, limited only by the atmospheric absorption) and the high duty cycle ($>$ 86 %) make the ARGO-YBJ experiment particularly suitable to detect short and unexpected events like GRBs. With the scaler mode technique, i.e., counting all the particles hitting the detector with no measurement of the primary energy and arrival direction, the minimum threshold of $\sim$ 1 GeV can be reached, overlapping the direct measurements carried out by satellites. During the experiment lifetime, from December 17, 2004 to February 7, 2013, a total of 206 GRBs occurring within the ARGO-YBJ field of view (zenith angle $\theta$ $\le$ 45$^{\circ}$) have been analyzed. This is the largest sample of GRBs investigated with a ground-based detector. Two lightcurve models have been assumed and since in both cases no significant excess has been found, the corresponding fluence upper limits in the 1-100 GeV energy region have been derived, with values as low as 10$^{-5}~$erg cm$^{-2}$. The analysis of a subset of 24 GRBs with known redshift has been used to constrain the fluence extrapolation to the GeV region together with possible cutoffs under different assumptions on the spectrum.

The AMS-02 has just published the unprecedentedly precise measurement of the cosmic electron and positron spectra. In this paper we try to give a quantitative study on the AMS-02 results by a global fitting to the electron and positron spectra, together with the updated positron fraction data. The Markov Chain Monte Carlo algorithm is adopted to do the fitting. The primary electron spectrum and the parameters for pulsars or dark matter which contribute extra positrons are determined simultaneously. We find that there is a hardening of the primary electron spectrum at $\sim 60$ GeV. With such a new feature at the background spectrum both the pulsars and dark matter can explain the AMS-02 results very well. The dark matter scenario shows a drop at positron fraction at $\sim 300$ GeV, however, suffers very strong constraints from Fermi $\gamma$-ray observations. The fitting results also suggest that the propagation model with convection may be more favored by the lepton data than the reacceleration model.

The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding their origin. Locating the "knees" of the spectra is an important part of the approach and has yet to be achieved. Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640+/-87 TeV, with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration models.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

The AMS-02 collaboration has just released its first result of the cosmic positron fraction $e^+/(e^-+e^+)$ with high precision up to $\sim 350$ GeV. The AMS-02 result shows the same trend with the previous PAMELA result, which requires extra electron/positron sources on top of the conventional cosmic ray background, either from astrophysical sources or from dark matter annihilation/decay. In this paper we try to figure out the nature of the extra sources by fitting to the AMS-02 $e^+/(e^-+e^+)$ data, as well as the electron and proton spectra by PAMELA and the $(e^-+e^+)$ spectrum by Fermi and HESS. We adopt the GALPROP package to calculate the propagation of the Galactic cosmic rays and the Markov Chain Monte Carlo sampler to do the fit. We find that the AMS-02 data have implied essential difference from the PAMELA data. There is {\rm tension} between the AMS-02 $e^+/(e^-+e^+)$ data and the Fermi/HESS $(e^-+e^+)$ spectrum, that the AMS-02 data requires less contribution from the extra sources than Fermi/HESS. Then we redo the fit without including the Fermi/HESS data. In this case both the pulsars and dark matter annihilation/decay can explain the AMS-02 data. The pulsar scenario has a soft inject spectrum with the power-law index $\sim 2$, while the dark matter scenario needs $\tau^+\tau^-$ final state with mass $\sim 600$ GeV and a boost factor $\sim 200$.

The recently reported positron fraction up to $\sim 350$ GeV by AMS-02 seems to have tension with the total electron/positron spectra detected by Fermi and HESS, for either pulsar or dark matter annihilation/decay scenario as the primary positron sources. In this work we will show that the tension will be removed by an adjustment of the primary electron spectrum. If the primary electron spectrum becomes harder above $\sim50$ GeV, similar as the cosmic ray nuclei spectrum, the AMS-02 positron fraction and Fermi/HESS data can be well fitted by both the pulsar and dark matter models. This result indicates that there should be a common origin of the cosmic ray nuclei and the primary electrons. Furthermore, this study also implies that the properties of the extra sources derived from the fitting to the AMS-02 data should depend on the form of background.

The energy spectrum of cosmic Hydrogen and Helium nuclei has been measured, below the so-called "knee", by using a hybrid experiment with a wide field-of-view Cherenkov telescope and the Resistive Plate Chamber (RPC) array of the ARGO-YBJ experiment at 4300 m above sea level. The Hydrogen and Helium nuclei have been well separated from other cosmic ray components by using a multi-parameter technique. A highly uniform energy resolution of about 25% is achieved throughout the whole energy range (100 TeV - 700 TeV). The observed energy spectrum is compatible with a single power law with index gamma=-2.63+/-0.06.

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

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.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

The extragalactic background light (EBL) contains important information about stellar and galaxy evolution. It leaves imprint on the very high energy $\gamma$-ray spectra from sources at cosmological distances due to the process of pair production. In this work we propose to {\em measure} the EBL directly by extracting the collective attenuation effects in a number of $\gamma$-ray sources at different redshifts. Using a Markov Chain Monte Carlo fitting method, the EBL intensities and the intrinsic spectral parameters of $\gamma$-ray sources are derived simultaneously. No prior shape of EBL is assumed in the fit. With this method, we can for the first time to derive the spectral shape of the EBL model-independently. Our result shows the expected features predicted by the present EBL models and thus support the understanding of the EBL origin.

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.

Shocks of supernova remnants (SNRs) are important (and perhaps the dominant) agents for production of the Galactic cosmic rays. Recent $\gamma$-ray observations of several SNRs have made this case more compelling. However, these broadband high-energy measurements also reveal a variety of spectral shape demanding more comprehensive modeling of emissions from SNRs. According to the locally observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all kinds of accelerated particles, we propose a unified model that ascribes the distinct $\gamma$-ray spectra of different SNRs to variations of the medium density and the spectral difference between cosmic ray electrons and protons observed at Earth to transport effects. For low density environments, the $\gamma$-ray emission is inverse-Compton dominated. For high density environments like systems of high-energy particles interacting with molecular clouds, the $\gamma$-ray emission is $\pi^0$-decay dominated. The model predicts a hadronic origin of $\gamma$-ray emission from very old remnants interacting mostly with molecular clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density environment created by their progenitors via e.g., strong stellar winds. These results can be regarded as evidence in support of the SNR-origin of the Galactic cosmic rays.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

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

A self-consistent global fitting method based on the Markov Chain Monte Carlo technique to study the dark matter (DM) property associated with the cosmic ray electron/positron excesses was developed in our previous work. In this work we further improve the previous study to include the hadronic branching ratio of DM annihilation/decay. The PAMELA $\bar{p}/p$ data are employed to constrain the hadronic branching ratio. We find that the 95% ($2\sigma$) upper limits of the quark branching ratio allowed by the PAMELA $\bar{p}/p$ data is $\sim 0.032$ for DM annihilation and $\sim 0.044$ for DM decay respectively. This result shows that the DM coupling to pure leptons is indeed favored by the current data. Based on the global fitting results, we further study the neutrino emission from DM in the Galactic center. Our predicted neutrino flux is some smaller than previous works since the constraint from $\gamma$-rays is involved. However, it is still capable to be detected by the forth-coming neutrino detector such as IceCube. The improved points of the present study compared with previous works include: 1) the DM parameters, both the particle physical ones and astrophysical ones, are derived in a global fitting way, 2) constraints from various species of data sets, including $\gamma$-rays and antiprotons are included, and 3) the expectation of neutrino emission is fully self-consistent.

Recently PAMELA collaboration published the cosmic nuclei and electron spectra with high precision, together with the cosmic antiproton data updated, and the Fermi-LAT collaboration also updated the measurement of the total $e^+e^-$ spectrum to lower energies. In this paper we develop a Markov Chain Monte Carlo (MCMC) package {\it CosRayMC}, based on the GALPROP cosmic ray propagation model to study the implications of these new data. It is found that if only the background electrons and secondary positrons are considered, the fit is very bad with $\chi_{\rm red}^2 \approx 3.68$. Taking into account the extra $e^+e^-$ sources of pulsars or dark matter annihilation we can give much better fit to these data, with the minimum $\chi_{\rm red}^2 \approx 0.83$. This means the extra sources are necessary with a very high significance in order to fit the data. However, the data show little difference between pulsar and dark matter scenarios. Both the background and extra source parameters are well constrained with this MCMC method. Including the antiproton data, we further constrain the branching ratio of dark matter annihilation into quarks $B_q<0.5%$ at $2\sigma$ confidence level. The possible systematical uncertainties of the present study are discussed.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We investigate the constraints on the scalar, vector and spin-3/2 dark matter interaction with the standard model particles, from the observations of dark matter relic density, the direct detection experiments of CDMS and XENON, and the indirect detection of the antiproton-to-proton ratio by PAMELA. A model independent way is adopted by constructing general 4-particle operators up to dimension 6 for the effective interaction between dark matter and standard model particles. We find that the constraints from different experiments are complementary with each other. Comparison among these constraints may exclude some effective models of dark matter and limit some parameters of others. The spin-independent direct detection gives strong constraints for some operators, while the indirect detection of antiproton-to-proton data can be more sensitive than direct detection or relic density for light dark matter (whose mass less than 70 GeV) in some cases. The constraints on some operators for spin-3/2 dark matter are shown to be similar to those on their analogous operators for Dirac fermionic dark matter. There are still some operators not sensitive to the current dark matter direct and indirect search experiments.

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.

In this work we study the constraints on the dark matter interaction with the standard model particles, from the observations of dark matter relic density, the direct detection experiments of CDMS and XENON, and the indirect detection of the antiproton-to-proton ratio by PAMELA. A model independent way is adopted in the study by constructing the effective interaction operators between dark matter and standard model particles. The most general 4-fermion operators are investigated. We find that the constraints from different observations are complementary with each other. Especially the spin independent scattering gives very strong constraints for corresponding operators. In some cases the indirect detection of antiproton-to-proton data can actually be more sensitive than the direct detection or relic density for light dark matter (less than 70 GeV).

The cosmic ray (CR) energy spectra measured with ATIC, CREAM and PAMELA showed that there is remarkable hardening for rigidity of several hundred GV. We propose that this hardening is due to the superposition of spectra from a population of sources, e.g., supernova remnants (SNRs), whose injection spectral indices have a dispersion. Adopting proper model parameters the observational data can be well explained. It is interesting that the injection source parameters are similar with that derived from gamma-ray observations of SNRs, which may support the SNR-origin of CRs. Furthermore this mechanism provides an alternative explanation of the "ankle-cutoff" structure of the ultra high energy CR spectra.

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.

Recent cosmic ray (CR) experiments discovered that the CR spectra experience a remarkable hardening for rigidity above several hundred GV. We propose that this is caused by the superposition of the CR energy spectra of many sources that have a dispersion in the injection spectral indices. Adopting similar parameters as those of supernova remnants derived from the Fermi $\gamma$-ray observations, we can reproduce the observational CR spectra of different species well. This may be interpreted as evidence to support the supernova remnant origin of CRs below the knee. We further propose that the same mechanism may explain the "ankle" of the ultra high energy CR spectrum.

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.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

Emission mechanisms of the shell-type supernova remnant (SNR) RX J1713.7-3946 are studied with multi-wavelength observational data from radio, X-ray, GeV $\gamma$-ray to TeV $\gamma$-ray band. A Markov Chain Monte Carlo method is employed to explore the high-dimensional model parameter space systematically. Three scenarios for the $\gamma$-ray emission are investigated: the leptonic, the hadronic and a hybrid one. Thermal emission from the background plasma is also included to constrain the gas density, assuming ionization equilibrium, and a 2$\sigma$ upper limit of about 0.03 cm$^{-3}$ is obtained as far as thermal energies account for a significant fraction of the dissipated kinetic energy of the SNR shock. Although systematic errors dominate the $\chi^2$ of the spectral fit of all models, we find that 1) the leptonic model has the best constrained model parameters, whose values can be easily accommodated with a typical supernova, but gives relatively poor fit to the TeV $\gamma$-ray data; 2) The hybrid scenario has one more parameter than the leptonic one and improves the overall spectral fit significantly; 3) The hadronic one, which has three more parameters than the leptonic model, gives the best fit to the overall spectrum with relatively not-well-constrained model parameters and very hard spectra of accelerated particles. The uncertainties of the model parameters decrease significantly if the spectral indices of accelerated electrons and protons are the same. The hybrid and hadronic models also require an energy input into high-energy protons, which seems to be too high compared with typical values of a supernova explosion. Further investigations are required to reconcile these observations with SNR theories.

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 sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

Assuming an intrinsic `Band' shape spectrum and an intrinsic energy-independent emission profile we have investigated the connection between the evolution of the rest-frame spectral parameters and the spectral lags measured in gamma-ray burst (GRB) pulses by using a pulse model. We first focus our attention on the evolution of the peak energy, $E_{0,p}$, and neglect the effect of the curvature effect. It is found that the evolution of $E_{0,p}$ alone can produce the observed lags. When $E_{0,p}$ varies from hard to soft only the positive lags can be observed. The negative lags would occur in the case of $E_{0,p}$ varying from soft to hard. When the evolution of $E_{0,p}$ and the low-energy spectral index $\alpha_{0}$ varying from soft to hard then to soft we can find the aforesaid two sorts of lags. We then examine the combined case of the spectral evolution and the curvature effect of fireball and find the observed spectral lags would increase. A sample including 15 single pulses whose spectral evolution follows hard to soft has been investigated. All the lags of these pulses are positive, which is in good agreement with our theoretical predictions. Our analysis shows that only the intrinsic spectral evolution can produce the spectral lags and the observed lags should be contributed by the intrinsic spectral evolution and the curvature effect. But it is still unclear what cause the spectral evolution.

We calculate the CMB temperature distortion due to the energetic electrons and positrons produced by dark matter annihilation (Sunyaev-Zel'dovich effect), in dwarf spheroidal galaxies (dSphs). In the calculation we have included two important effects which were previously ignored. First we show that the electron-positron pairs with energy less than GeV, which were neglected in previous calculation, could contribute a significant fraction of the total signal. Secondly we also consider the full effects of diffusion loss, which could significantly reduce the density of electron-positron pairs at the center of cuspy halos. For neutralinos, we confirm that detecting such kind of SZ effect is beyond the capability of the current or even the next generation experiments. In the case of light dark matter (LDM) the signal is much larger, but even in this case it is only marginally detectable with the next generation of experiment such as ALMA. We conclude that similar to the case of galaxy clusters, in the dwarf galaxies the $SZ_2DM}$ effect is not a strong probe of DM annihilations.

The cosmic electron and positron excesses have been explained as possible dark matter (DM) annihilation products. In this work we investigate the possible effects of such a DM annihilation scenario during the evolution history of the Universe. We first calculate the extragalactic $\gamma$-ray background (EGRB), which is produced through the final state radiation of DM annihilation to charged leptons and the inverse Compton scattering between electrons/positrons and the cosmic microwave background. The DM halo profile and the minimal halo mass, which are not yet well determined from the current N-body simulations, are constrained by the EGRB data from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic DM models on cosmic evolution, such as the reionization and heating of intergalactic medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We show that the impact on the Hydrogen 21 cm signal might show interesting signatures of DM annihilation, but the influence on star formation is not remarkable. Future observations of the 21 cm signals could be used to place new constraints on the properties of DM.

Monochromatic gamma-rays are thought to be the smoking gun signal for identifying the dark matter annihilation. However, the flux of monochromatic gamma-rays is usually suppressed by the virtual quantum effects since dark matter should be neutral and does not couple with gamma-rays directly. In the work we study the detection strategy of the monochromatic gamma-rays in a future space-based detector. The monochromatic gamma-ray flux is calculated by assuming supersymmetric neutralino as a typical dark matter candidate. We discuss both the detection focusing on the Galactic center and in a scan mode which detects gamma-rays from the whole Galactic halo are compared. The detector performance for the purpose of monochromatic gamma-rays detection, with different energy and angular resolution, field of view, background rejection efficiencies, is carefully studied with both analytical and fast Monte-Carlo method.

In this work, we revisit the all-sky Galactic diffuse $\gamma$-ray emission taking into account the new measurements of cosmic ray electron/positron spectrum by PAMELA, ATIC and Fermi, which show excesses of cosmic electrons/positrons beyond the expected fluxes in the conventional model. Since the origins of the extra electrons/positrons are not clear, we consider three different scenarios to account for the excesses: the astrophysical sources such as the Galactic pulsars, dark matter decay and annihilation. Further, new results from Fermi-LAT of the (extra-)Galactic diffuse $\gamma$-ray are adopted. The background cosmic rays without the new sources give lower diffuse $\gamma$ rays compared to Fermi-LAT observation, which is consistent with previous analysis. The scenario with astrophysical sources predicts diffuse $\gamma$-rays with little difference with the background. The dark matter annihilation models with $\tau^{\pm}$ final state are disfavored by the Fermi diffuse $\gamma$-ray data, while there are only few constraints on the decaying dark matter scenario. Furthermore, these is always a bump at higher energies ($\sim$ TeV) of the diffuse $\gamma$-ray spectra for the dark matter scenarios due to final state radiation. Finally we find that the Fermi-LAT diffuse $\gamma$-ray data can be explained by simply enlarging the normalization of the electron spectrum without introduce any new sources, which may indicate that the current constraints on the dark matter models can be much stronger given a precise background estimate.

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.

In this paper we have analyzed the temporal and spectral behavior of 52 Fast Rise and Exponential Decay (FRED) pulses in 48 long-duration gamma-ray bursts (GRBs) observed by the CGRO/BATSE, using a pulse model with two shape parameters and the Band model with three shape parameters, respectively. It is found that these FRED pulses are distinguished both temporally and spectrally from those in long-lag pulses. Different from these long-lag pulses only one parameter pair indicates an evident correlation among the five parameters, which suggests that at least $\sim$4 parameters are needed to model burst temporal and spectral behavior. In addition, our studies reveal that these FRED pulses have correlated properties: (i) long-duration pulses have harder spectra and are less luminous than short-duration pulses; (ii) the more asymmetric the pulses are the steeper the evolutionary curves of the peak energy ($E_{p}$) in the $\nu f_{\nu}$ spectrum within pulse decay phase are. Our statistical results give some constrains on the current GRB models.

The paper investigates the overall and detailed features of cosmic ray (CR) spectra in the knee region using the scenario of nuclei-photon interactions around the acceleration sources. Young supernova remnants can be the physical realities of such kind of CR acceleration sites. The results show that the model can well explain the following problems simultaneously with one set of source parameters: the knee of CR spectra and the sharpness of the knee, the detailed irregular structures of CR spectra, the so-called "component B" of Galactic CRs, and the electron/positron excesses reported by recent observations. The coherent explanation serves as evidence that at least a portion of CRs might be accelerated at the sources similar to young supernova remnants, and one set of source parameters indicates that this portion mainly comes from standard sources or from a single source.

In this paper we develop a Markov Chain Monte Carlo code to study the dark matter properties in interpreting the recent observations of cosmic ray electron/positron excesses. We assume that the dark matter particles couple dominantly to leptons and consider two cases separately with annihilating or decaying into lepton pairs. The constraint on the central density profile from diffuse $\gamma$-rays around the Galactic center is also included in the Markov Chain Monte Carlo code self-consistently. In the numerical study, 7 parameters are introduced. We fit two data sets independently and find that for the Data set I (PAMELA+ATIC), dark matter with $m_{\chi}\approx0.7$ TeV for annihilation (or 1.4 TeV for decay) and pure $e^+e^-$ channel is favored, while for the Data set II (PAMELA+Fermi-LAT+H.E.S.S.) $m_{\chi}\approx 2$ TeV for annihilation (or 4 TeV for decay) and the combination of $\mu^+\mu^-$ and $\tau^+\tau^-$ final states can best fit the data. The H.E.S.S. observation of the Galactic center $\gamma$-rays puts a strong constraint on the central density profile of the dark matter halo for the annihilation dark matter scenario. In this case the NFW profile which is regarded as the typical predication from the cold dark matter scenario, is excluded with a high significance ($>5\sigma$). For the decaying dark matter scenario, the constraint is much weaker.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

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.

In this work we study how the cosmological parameter, the Hubble constant $H_0$, can be constrained by observation of very high energy (VHE) $\gamma$-rays at the TeV scale. The VHE $\gamma$-rays experience attenuation by background radiation field through $e^+e^-$ pair production during the propagation in the intergalactic space. This effect is proportional to the distance that the VHE $\gamma$-rays go through. Therefore the absorption of TeV $\gamma$-rays can be taken as cosmological distance indicator to constrain the cosmological parameters. Two blazars Mrk 501 and 1ES 1101-232, which have relatively good spectra measurements by the atmospheric Cerenkov telescope, are studied to constrain $H_0$. The mechanism constraining the Hubble constant adopted here is very different from the previous methods such as the observations of type Ia supernovae and the cosmic microwave background. However, at $2\sigma$ level, our result is consistent with other methods.

Dark matter annihilation is one of the leading explanations for the recently observed $e^\pm$ excesses in cosmic rays by PAMELA, ATIC, FERMI-LAT and HESS. Any dark matter annihilation model proposed to explain these data must also explain the fact that PAMELA data show excesses only in $e^\pm$ spectrum but not in anti-proton. It is interesting to ask whether the annihilation mode into anti-proton is completely disallowed or only suppressed at low energies. Most models proposed have negligible anti-protons in all energy ranges. We show that the leptocentric $U(1)_{B-3L_i}$ dark matter model can explain the $e^\pm$ excesses with suppressed anti-proton mode at low energies, but at higher energies there are sizable anti-proton excesses. Near future data from PAMELA and AMS can provide crucial test for this type of models. Cosmic $\gamma$ ray data can further rule out some of the models. We also show that this model has interesting cosmic neutrino signatures.

In this work we calculate the Sunyaev-Zel'dovich (SZ) effect due to the $e^+e^-$ from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, (1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) are investigated. For each case, we also consider several DM profiles with and without central cusp. We generally find smaller signals than previously reported. Moreover, the diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be negligible, is considered and found to have significant effect on the central electron/positron distribution for DM profile with large spatial gradient. We find that the SZ effect from WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the next generation SZ interferometer such as ALMA. Although the signal of the LDM is much larger than that of the WIMP, the final SZ effect is still very small due to the smoothing effect of diffusion. Only for the configuration with large central cusp and extremely small diffusion effect, the LDM induced SZ effect might have a bit chance of being detected.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

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.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants.

The perspective of the detectability of Galactic dark matter subhaloes on the Fermi satellite is investigated in this work. Under the assumptions that dark matter annihilation accounts for the "GeV excess" of the Galactic diffuse $\gamma$-rays discovered by EGRET and the $\gamma$-ray flux is dominated by the contribution from subhaloes of dark matter, we calculate the expected number of dark matter subhaloes that Fermi may detect. We show that Fermi may detect a few tens to several hundred subhaloes in 1-year all sky survey. Since EGRET observation is taken as a normalization, this prediction is independent of the particle physics property of dark matter. The uncertainties of the prediction are discussed in detail. We find that the major uncertainty comes from the mass function of subhaloes, i.e., whether the subhaloes are "point like" (high-mass rich) or "diffuse like" (low-mass rich). Other uncertainties like the background estimation and the observational errors will contribute a factor of $2\sim 3$.

In this paper we propose a dark matter model and study aspects of its phenomenology. Our model is based on a new dark matter sector with a U(1)' gauge symmetry plus a discrete symmetry added to the Standard Model of particle physics. The new fields of the dark matter sector have no hadronic charges and couple only to leptons. Our model can not only give rise to the observed neutrino mass hierarchy, but can also generate the baryon number asymmetry via non-thermal leptogenesis. The breaking of the new U(1)' symmetry produces cosmic strings. The dark matter particles are produced non-thermally from cosmic string loop decay which allows one to obtain sufficiently large annihilation cross sections to explain the observed cosmic ray positron and electron fluxes recently measured by the PAMELA, ATIC, PPB-BETS, Fermi-LAT, and HESS experiments while maintaining the required overall dark matter energy density. The high velocity of the dark matter particles from cosmic string loop decay leads to a low phase space density and thus to a dark matter profile with a constant density core in contrast to what happens in a scenario with thermally produced cold dark matter where the density keeps rising towards the center. As a result, the flux of gamma rays radiated from the final leptonic states of dark matter annihilation from the Galactic center is suppressed and satisfies the constraints from the HESS gamma-ray observations.

The GZK cutoff predicted at the Ultra High Energy Cosmic Ray (UHECR) spectrum as been observed by the HiRes and Auger experiments. The results put severe constraints on the effect of Lorentz Invariance Violation(LIV) which has been introduced to explain the absence of GZK cutoff indicated in the AGASA data. Assuming homogeneous source distribution with a single power law spectrum, we calculate the spectrum of UHECRs observed on Earth by taking the processes of photopion production, $e^+e^-$ pair production and adiabatic energy loss into account. The effect of LIV is also taken into account in the calculation. By fitting the HiRes monocular spectra and the Auger combined spectra, we show that the LIV parameter is constrained to $\xi=-0.8^{+3.2}_{-0.5}\times10^{-23}$ and $0.0^{+1.0}_{-0.4}\times10^{-23}$ respectively, which is well consistent with strict Lorentz Invariance up to the highest energy.

Recently the Milagro experiment observed diffuse multi-TeV gamma-ray emission in the Cygnus region, which is significantly stronger than what predicted by the Galactic cosmic ray model. However, the sub-GeV observation by EGRET shows no excess to the prediction based on the same model. This TeV excess implies possible high energy cosmic rays populated in the region with harder spectrum than that observed on the Earth. In the work we studied this theoretical speculation in detail. We find that, a diffuse proton source with power index $\alpha_p\lesssim 2.3$, or a diffuse electron source with power index $\alpha_e\lesssim2.6$ can reproduce the Milagro's observation without conflicting with the EGRET data. Further detections on neutrinos, a diagnostic of the hadronic model, and hard X-ray synchrontron radiation, a diagnostic of the lepton model, help to break this degeneracy. In combination with the gamma ray observations to several hundred GeV by Fermi, we will be able to understand the diffuse emission mechanisms in the Cygnus region better.

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.

The diffuse Galactic $\gamma$-rays from EGRET observation shows excesses above 1 GeV in comparison with the expectations from conventional Galactic cosmic ray (CR) propagation model. In the work we try to solve the "GeV excess" problem by dark matter (DM) annihilation in the frame of supersymmetry (SUSY). Compared with previous works, there are three aspects improved in this work: first, the direction-independent "boost factor" for diffuse $\gamma$-rays from dark matter annihilation (DMA) is naturally reproduced by taking the DM substructures into account; second, there is no need for renormalization of the diffuse $\gamma$-ray background produced by CRs; last but not the least, in this work our new propagation model can give consistent results of both diffuse $\gamma$-rays and antiprotons, by directly adding the signals from DMA to the diffuse $\gamma$-ray background. This is a self-consistent model among several possible scenarios at present, and can be tested or optimized by the forthcoming experiments such as GLAST, PAMELA and AMS02.

In this work we develop a new propagation model for the Galactic cosmic rays based on the GALPROP code, including contributions from dark matter annihilation. The model predicts compatible Galactic diffuse $\gamma$ ray spectra with EGRET data in all sky regions. It also gives consistent results of the diffuse $\gamma$ ray longitude and latitude distributions. Further the results for B/C, $^{10}$Be/$^9$Be, proton, electron and antiproton spectra are also consistent with cosmic ray measurements. In the model we have taken a universal proton spectrum throughout the Galaxy without introducing large fluctuation for the proton energy loss is negligible. The dark matter annihilation signals are `boosted' after taking the contribution from subhalos into account. Another interesting feature of the model is that it gives better description of the diffuse $\gamma$ rays when taking the source distribution compatible with supernova remnants data, which is different from previous studies.

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

Using the data taken from Tibet II High Density (HD) Array (1997 February-1999 September) and Tibet-III array (1999 November-2005 November), our previous northern sky survey for TeV $\gamma-$ray point sources has now been updated by a factor of 2.8 improved statistics. From $0.0^{\circ}$ to $60.0^{\circ}$ in declination (Dec) range, no new TeV $\gamma-$ray point sources with sufficiently high significance were identified while the well-known Crab Nebula and Mrk421 remain to be the brightest TeV $\gamma-$ray sources within the field of view of the Tibet air shower array. Based on the currently available data and at the 90% confidence level (C.L.), the flux upper limits for different power law index assumption are re-derived, which are approximately improved by 1.7 times as compared with our previous reported limits.

Anti-proton and positron Galactic cosmic ray spectra are among the key targets for indirect detection of dark matter. The boost factors, corresponding to an enhancement of the signal and linked to the clumpiness properties of the dark matter distribution, have been taken as high as thousands in the past. The dramatic impact of these boost factors for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks for their detailed calculation. We take into account the results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors linked to the cosmic ray propagation properties, for anti-protons and positrons. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are used to calculate the mean value and the variance of the boost factor for the primary fluxes. We take advantage of the statistical method introduced in Lavalle et al. (2007) and cross-check the results with Monte Carlo computations. By spanning some extreme configurations of sub-halo and propagation properties, we find that the average contribution of the clumps is negligible compared to that of the smooth dark matter component. Sub-halos do not lead to enhancement of the signals, unless they are taken with some extreme (unexpected) properties. This result is independent of the nature of the self-annihilating dark matter candidate considered, and provides precise estimates of the theoretical and the statistical uncertainties of the antimatter flux from dark matter substructures. Spectral distortions can still be expected in antimatter flux measurements, but scenarios invoking large and even mild clumpiness boost factors are strongly disfavoured by our analysis.

In this paper we study the detectability of $\gamma$-rays from dark matter annihilation in the subhalos of the Milky Way by the satellite-based experiments, EGRET and GLAST. We work in the frame of supersymmetric extension of the standard model and assume the lightest neutralino being the dark matter particles. Based on the N-body simulation of the evolution of dark matter subhalos we first calculate the average intensity distribution of this new class of $\gamma$-ray sources by neutralino annihilation. It is possible to detect these $\gamma$-ray sources by EGRET and GLAST. Conversely, if these sources are not detected the nature of the dark matter particls will be constrained by these experiments, which, however, depending on the uncertainties of the subhalo profile.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

We model the electrons/positrons produced by dark matter annihilations in the colliding galaxy cluster system 1E0657-56 (the bullet cluster). These charged particles, confined by the Magnetic field, clearly trace the path of the bullet, which passes through the main cluster with a speed of 3000-5000 km/s. Adding the effect of dark matter substructure (subhalos) in each cluster we find the annihilation rate is enhancedgreatly and the density of positrons in the trail is similar to that within the bullet cluster. This opens the door to unambiguous detection of the annihilation signal through the SZ effect, at significant separation from the thermal electrons.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

In this paper we calculate the Galactic positron flux from dark matter annihilation in the frame of supersymmetry, taking the enhancement of the flux by existence of dark matter substructures into account. The propagation of positrons in the Galactic magnetic field is solved in a realistic numerical model GALPROP. The secondary positron flux is recalculated in the GLAPROP model. The total positron flux from secondary products and dark matter annihilation can fit the HEAT data well when taking a cuspy density profile of the substructures.

In this paper we study the possibility of detecting gamma rays from dark matter annihilation in the subhalos of the Milky Way by the ground based EAS detectors within the frame of the minimal supersymmetric standard model. Based on the Monte Carlo simulation we also study the properties of two specific EAS detectors, the ARGO and HAWC, and the sensitivities of these detectors on the detection of dark matter annihilation. We find the ground EAS detectors have the possibility to observe such signals. Conversely if no signal observed we give the constraints on the supersymmetric parameter space, which however depends on the subhalos properties.

The CACTUS experiment recently observed a gamma ray excess above 50 GeV from the direction of the Draco dwarf spheroidal galaxy. Considering that Draco is dark matter dominated the gamma rays may be generated through dark matter annihilation in the Draco halo. In the framework of the minimal supersymmetric extension of the standard model we explore the parameter space to account for the gamma ray signals at CACTUS. We find that the neutralino mass is constrained to be approximately in the range between 100 GeV ~ 400 GeV and a sharp central cuspy of the dark halo profile in Draco is necessary to explain the CACTUS results. We then discuss further constraints on the supersymmetric parameter space by observations at the ground based ARGO detector. It is found that the parameter space can be strongly constrained by ARGO if no excess from Draco is observed above 100 GeV.

High resolution simulations reveal that in the cold dark matter scenario the structures form hierarchically and a large number of substructures survive in the galactic halos. The substructures can be probed if they emit gamma rays via dark matter annihilation. We calculated the gamma ray fluxes from the dark matter annihilations in the substructures of our Galaxy within the frame of the minimal supersymmetric extension of the standard model. The uncertainties of the prediction from both the low energy supersymmetry and especially from the density profiles of dark matter in the substructures are carefully investigated. The cumulative number of substructures emitting gamma rays above any given flux is calculated. Detectability of the gamma rays from the substructures is discussed. We propose the viability to detect these signals through the ground large field of view detectors.

In this paper we consider the cosmological implications of dark energy models with a coupled system of a dynamical scalar field (the quintessence) and the neutrinos. By detailed numerical calculations we study the various possibilities on the evolution and the fates of the universe in this class of models. Our results show that due to the interaction with quintessence, neutrinos could be dominant over the quintessence in the future universe, however would eventually decay away.

Adopting a recent calculation of the Galactic interstellar radiation field, we calculate the attenuation of the very high energy gamma rays from the Galactic sources. The infra-red radiation background near the Galactic Center is very intense due to the new calculation and our result shows that a cutoff of high energy gamma ray spectrum begins at about 20 TeV and reaches about 10% for 50 TeV gamma rays.

We study in this paper the possible influence of the dark energy on the detection of the dark matter particles. In models of dark energy described by a dynamical scalar field such as the Quintessence, its interaction with the dark matter will cause the dark matter particles such as the neutralino vary as a function of space and time. Given a specific model of the Quintessence and its interaction in this paper we calculate numerically the corrections to the neutralino masses and the induced spectrum of the neutrinos from the annihilation of the neutralinos pairs in the core of the Sun. This study gives rise to a possibility of probing for dark energy in the experiments of detecting the dark matter particles.

Quintessence is a possible candidate for dark energy. In this paper we study the phenomenologies of the fermionic partner of Quintessence, the Quintessino. Our results show that, for suitable choices of the model parameters, the Quintessino is a good candidate for cold or warm dark matter. In our scenario, dark energy and dark matter of the Universe are connected in one chiral superfield.

In this paper we consider the possibility of neutrino mass varying during the evolution of the Universe and study its implications on leptogenesis. Specifically, we take the minimal seesaw model of neutrino masses and introduce a coupling between the right-handed neutrinos and the dark energy scalar field, the Quintessence. In our model, the right-handed neutrino masses change as the Quintessence scalar evolves. We then examine in detail the parameter space of this model allowed by the observed baryon number asymmetry. Our results show that it is possible to lower the reheating temperature in this scenario in comparison with the case that the neutrino masses are unchanged, which helps solve the gravitino problem. Furthermore, a degenerate neutrino mass patten with $m_i$ larger than the upper limit given in the minimal leptogenesis scenario is permitted.

In this paper we study double scalar neutrino inflation in the minimal supersymmetric seesaw model in light of WMAP. Inflation in this model is firstly driven by the heavier sneutrino field $\tilde{N}_2$ and then the lighter field $\tilde{N}_1$. we will show that with the mass ratio $6\lesssim M_2/M_1 \lesssim 10$ the model predicts a suppressed primordial scalar spectrum around the largest scales and the predicted CMB TT quadrupole is much better suppressed than the single sneutrino model. So this model is more favored than the single sneutrino inflation model. We then consider the implications of the model on the reheating temperature, leptogenesis and lepton flavor violation. Our results show that the seesaw parameters are constrained strongly by the reheating temperature, together with the requirement by a successful inflation. The mixing between the first generation and the other two generations in the right-handed neutrino sector is tiny. The rates of lepton flavor violating processes in our scenario depend on only 4 unknown seesaw parameters through a 'reduced' seesaw formula, besides $U_{e3}$ and the supersymmetric parameters. We find that the branching ratio of $\mu \to e \gamma$ is generally near the present experimental limit, while ${\rm Br}(\tau \to \mu \gamma)$ is around $\mathcal{O}(10^{-10} - 10^{-9})$.