Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

Several of the many proposed Dark Matter candidate particles, already investigated with lower exposure and a higher software energy threshold, are further analyzed including the first DAMA/LIBRA--phase2 data release, with an exposure of 1.13 ton $\times$ yr and a lower software energy threshold (1 keV). The cumulative exposure above 2 keV considering also DAMA/NaI and DAMA/LIBRA--phase1 results is now 2.46 ton $\times$ yr. The analysis permits to constraint the parameters' space of the considered candidates restricting their values -- with respect to previous analyses -- thanks to the increase of the exposure and to the lower energy threshold.

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 first model independent results obtained by the DAMA/LIBRA-phase2 experiment are presented. The data have been collected over 6 annual cycles corresponding to a total exposure of 1.13 ton $\times$ yr, deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. The DAMA/LIBRA-phase2 apparatus, $\simeq$ 250 kg highly radio-pure NaI(Tl), profits from a second generation high quantum efficiency photomultipliers and of new electronics with respect to DAMA/LIBRA-phase1. The improved experimental configuration has also allowed to lower the software energy threshold. New data analysis strategies are presented. The DAMA/LIBRA-phase2 data confirm the evidence of a signal that meets all the requirements of the model independent Dark Matter (DM) annual modulation signature, at 9.5 $\sigma$ C.L. in the energy region (1-6) keV. In the energy region between 2 and 6 keV, where data are also available from DAMA/NaI and DAMA/LIBRA-phase1 (exposure $1.33$ ton $\times$ yr, collected over 14 annual cycles), the achieved C.L. for the full exposure (2.46 ton $\times$ yr) is 12.9 $\sigma$; the modulation amplitude of the single-hit scintillation events is: $(0.0103 \pm 0.0008)$ cpd/kg/keV, the measured phase is $(145 \pm 5)$ days and the measured period is $(0.999 \pm 0.001)$ yr, all these values are well in agreement with those expected for DM particles. No systematics or side reaction able to mimic the exploited DM signature (i.e. to account for the whole measured modulation amplitude and to simultaneously satisfy all the requirements of the signature), has been found or suggested by anyone throughout some decades thus far.

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.

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.

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.

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.

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.

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.

In the present paper the results obtained in the investigation of possible diurnal effects for low-energy single-hit scintillation events of DAMA/LIBRA-phase1 (1.04 ton $\times$ yr exposure) have been analysed in terms of an effect expected in case of Dark Matter (DM) candidates inducing nuclear recoils and having high cross-section with ordinary matter, which implies low DM local density in order to fulfill the DAMA/LIBRA DM annual modulation results. This effect is due to the different Earth depths crossed by those DM candidates during the sidereal day.

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

This paper summarizes in a simple and intuitive way why the neutrons, the muons and the solar neutrinos cannot give any significant contribution to the DAMA annual modulation results. A number of these elements have already been presented in individual papers; they are recalled here. Afterwards, few simple considerations are summarized which already demonstrate the incorrectness of the claim reported in PRL 113 (2014) 081302.

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.

MGRO J2019+37 within the Cygnus region is a bright and extended source revealed by Milagro at 12-35 TeV. This source is almost as bright as Crab Nebula in northern sky, while it is not confirmed by ARGO-YBJ around TeV. Up to now, no obvious counterpart at low energy wavelengths has been found. Hence, MGRO J2019+37 becomes mysterious and its VHE gamma-ray emission mechanism is attractive. In this paper, a brief summary of the multi-wavelength observations from Radio to gamma-ray is presented. All the available data from the XMM-Newton and INTEGRAL at X-ray, and the Fermi-LAT at gamma-ray bands were used to get constraint on its emission flux at low energy wavelengths. Then, its possible counterparts and the VHE emission mechanism are discussed.

The results obtained in the search for possible diurnal effect in the single-hit low energy data collected by DAMA/LIBRA-phase1 (total exposure: 1.04 ton x yr) deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. are presented. At the present level of sensitivity the presence of any significant diurnal variation and of diurnal time structures in the data can be excluded for both the cases of solar and sidereal time. In particular, the diurnal modulation amplitude expected, because of the Earth diurnal motion, on the basis of the DAMA Dark Matter annual modulation results is below the present sensitivity.

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.

The results obtained with the total exposure of 1.04 ton x yr collected by DAMA/LIBRA-phase1 deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. during 7 annual cycles (i.e. adding a further 0.17 ton x yr exposure) are presented. The DAMA/LIBRA-phase1 data give evidence for the presence of Dark Matter (DM) particles in the galactic halo, on the basis of the exploited model independent DM annual modulation signature by using highly radio-pure NaI(Tl) target, at 7.5 sigma C.L.. Including also the first generation DAMA/NaI experiment (cumulative exposure 1.33 ton x yr, corresponding to 14 annual cycles), the C.L. is 9.3 sigma and the modulation amplitude of the single-hit events in the (2-6) keV energy interval is: (0.0112 \pm 0.0012) cpd/kg/keV; the measured phase is (144 \pm 7) days and the measured period is (0.998 \pm 0.002) yr, values well in agreement with those expected for DM particles. No systematic or side reaction able to mimic the exploited DM signature has been found or suggested by anyone over more than a decade.

Experimental observations and theoretical arguments at Galaxy and larger scales have suggested that a large fraction of the Universe is composed by Dark Matter particles. This has motivated the DAMA experimental efforts to investigate the presence of such particles in the galactic halo by exploiting a model independent signature and very highly radiopure set-ups deep underground. Few introductory arguments are summarized before presenting a review of the present model independent positive results obtained by the DAMA/NaI and DAMA/LIBRA set-ups at the Gran Sasso National Laboratory of the INFN. Implications and model dependent comparisons with other different kinds of results will be shortly addressed. Some arguments put forward in literature will be confuted.

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 DAMA/LIBRA experiment, running at the Gran Sasso National Laboratory of the I.N.F.N. in Italy, has a sensitive mass of about 250 kg highly radiopure NaI(Tl). It is mainly devoted to the investigation of Dark Matter (DM) particles in the Galactic halo by exploiting the model independent DM annual modulation signature. The present DAMA/LIBRA experiment and the former DAMA/NaI one (the first generation experiment having an exposed mass of about 100 kg) have released so far results corresponding to a total exposure of 1.17 ton yr over 13 annual cycles. They provide a model independent evidence of the presence of DM particles in the galactic halo at 8.9 sigma C.L.. A short summary of the obtained results is presented and future perspectives of the experiment mentioned.

A few final comments on arXiv:1210.7548 are given to confute incorrect arguments claimed there.

We briefly remind references and arguments, already discussed in the past, which confute erroneous claims in arXiv:1210.5501.

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.

This paper gathers arguments and reasons why muons surviving the Gran Sasso mountain cannot mimic the Dark Matter annual modulation signature exploited by the DAMA/NaI and DAMA/LIBRA experiments. A number of these items have already been presented in individual papers. Further arguments have been addressed here in order to present a comprehensive collection and to enable a wider community to correctly approach this point.

Some theoretical and experimental aspects regarding the direct dark matter field are mentioned. In particular some arguments, which play a relevant role in the evaluation of model dependent interpretations of experimental results and in comparisons, are shortly addressed.

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.

A wide field of view Cherenkov/fluorescence telescope array is one of the main components of the Large High Altitude Air Shower Observatory project. To serve as Cherenkov and fluorescence detectors, a flexible and mobile design is adopted for easy reconfiguring of the telescope array. Two prototype telescopes have been constructed and successfully run at the site of the ARGO-YBJ experiment in Tibet. The features and performance of the telescopes are presented.

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.

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.

The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compositions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of unrecorded photons in the insensitive areas of the camera.

The present DAMA/LIBRA experiment and the former DAMA/NaI have cumulatively released so far the results obtained with the data collected over 13 annual cycles (total exposure: 1.17 ton $\times$ yr). They give a model independent evidence of the presence of DM particles in the galactic halo on the basis of the DM annual modulation signature at 8.9 $\sigma$ C.L. for the cumulative exposure.

DAMA/LIBRA is running at the Gran Sasso National Laboratory of the I.N.F.N.. Here the results obtained with a further exposure of 0.34 ton x yr are presented. They refer to two further annual cycles collected one before and one after the first DAMA/LIBRA upgrade occurred on September/October 2008. The cumulative exposure with those previously released by the former DAMA/NaI and by DAMA/LIBRA is now 1.17 ton x yr, corresponding to 13 annual cycles. The data further confirm the model independent evidence of the presence of Dark Matter (DM) particles in the galactic halo on the basis of the DM annual modulation signature (8.9 sigma C.L. for the cumulative exposure). In particular, with the cumulative exposure the modulation amplitude of the single-hit events in the (2 -- 6) keV energy interval measured in NaI(Tl) target is (0.0116 +- 0.0013) cpd/kg/keV; the measured phase is (146 +- 7) days and the measured period is (0.999 +- 0.002) yr, values well in agreement with those expected for the DM particles.

A variety of detectors has been proposed for dark matter direct detection, but most of them -- by the fact -- are still at R&D stage. In many cases, it is claimed that the lack of an adequate detectors' radio-purity might be compensated through heavy uses of MonteCarlo simulations, subtractions and handlings of the measured counting rates, in order to claim higher sensitivity (just for a particular scenario). The relevance of a correct evaluation of systematic effects in the use of MonteCarlo simulations at very low energy (which has always been safely discouraged in the field so far) and of multiple subtractions and handling procedures applied to the measured counting rate is shortly addressed here at some extent. Many other aspects would also deserve suitably deep investigations.

The DAMA/LIBRA set-up (about 250 kg highly radiopure NaI(Tl) sensitive mass) is running at the Gran Sasso National Laboratory of the I.N.F.N.. The first DAMA/LIBRA results confirm the evidence for the presence of a Dark Matter particle component in the galactic halo, as pointed out by the former DAMA/NaI set-up; cumulatively the data support such evidence at 8.2 sigma C.L. and satisfy all the many peculiarities of the Dark Matter annual modulation signature. The main aspects and prospects of this model independent experimental approach will be outlined.

The highly radiopure $\simeq$ 250 kg NaI(Tl) DAMA/LIBRA set-up is running at the Gran Sasso National Laboratory of the I.N.F.N.. In this paper the first result obtained by exploiting the model independent annual modulation signature for Dark Matter (DM) particles is presented. It refers to an exposure of 0.53 ton$\times$yr. The collected DAMA/LIBRA data satisfy all the many peculiarities of the DM annual modulation signature. Neither systematic effects nor side reactions can account for the observed modulation amplitude and contemporaneously satisfy all the several requirements of this DM signature. Thus, the presence of Dark Matter particles in the galactic halo is supported also by DAMA/LIBRA and, considering the former DAMA/NaI and the present DAMA/LIBRA data all together (total exposure 0.82 ton$\times$yr), the presence of Dark Matter particles in the galactic halo is supported at 8.2 $\sigma$ C.L..

The $\simeq$ 250 kg highly radiopure NaI(Tl) DAMA/LIBRA apparatus, running at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N., is described.

Some extensions of the Standard Model provide Dark Matter candidate particles with sub-GeV mass. These Light Dark Matter particles have been considered for example in Warm Dark Matter scenarios (e.g. the keV scale sterile neutrino, axino or gravitino). Moreover MeV scale DM candidates have been proposed in supersymmetric models and as source of the 511 keV line from the Galactic center. In this paper the possibility of direct detection of a Light Dark Matter candidate is investigated considering the inelastic scattering processes on the electron or on the nucleus targets. Some theoretical arguments are developed and related phenomenological aspects are discussed. Allowed volumes and regions for the characteristic phenomenological parameters of the considered scenarios are derived from the DAMA/NaI annual modulation data.

Some extensions of the Standard Model provide Dark Matter candidate particles which can have a dominant coupling with the lepton sector of the ordinary matter. Thus, such Dark Matter candidate particles ($\chi^{0}$) can be directly detected only through their interaction with electrons in the detectors of a suitable experiment, while they are lost by experiments based on the rejection of the electromagnetic component of the experimental counting rate. These candidates can also offer a possible source of the 511 keV photons observed from the galactic bulge. In this paper this scenario is investigated. Some theoretical arguments are developed and related phenomenological aspects are discussed. Allowed intervals and regions for the characteristic phenomenological parameters of the considered model and of the possible mediator of the interaction are also derived considering the DAMA/NaI data.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

The channeling effect of low energy ions along the crystallographic axes and planes of NaI(Tl) crystals is discussed in the framework of corollary investigations on WIMP Dark Matter candidates. In fact, the modeling of this existing effect implies a more complex evaluation of the luminosity yield for low energy recoiling Na and I ions. In the present paper related phenomenological arguments are developed and possible implications are discussed at some extent.

The effect pointed out by A. B. Migdal in the 40's (hereafter named Migdal effect) has so far been usually neglected in the direct searches for WIMP Dark Matter candidates. This effect consists in the ionization and the excitation of bound atomic electrons induced by the recoiling atomic nucleus. In the present paper the related theoretical arguments are developed and some consequences of the proper accounting for this effect are discussed by some examples of practical interest.

Galaxy hierarchical formation theories, numerical simulations, the discovery of the Sagittarius Dwarf Elliptical Galaxy (SagDEG) in 1994 and more recent investigations suggest that the dark halo of the Milky Way can have a rich phenomenology containing non thermalized substructures. In the present preliminary study, we investigate the case of the SagDEG (the best known satellite galaxy in the Milky Way crossing the solar neighbourhood) analyzing the consequences of its dark matter stream contribution to the galactic halo on the basis of the DAMA/NaI annual modulation data. The present analysis is restricted to some WIMP candidates and to some of the astrophysical, nuclear and particle Physics scenarios. Other candidates such as e.g. the light bosonic ones, we discussed elsewhere, and other non thermalized substructures are not yet addressed here.