ABSTRACT This paper reports on the measurement of the large-scale anisotropy in the distribution of cosmic-ray arrival directions using the data collected by the air shower detector ARGO-YBJ from ...2008 January to 2009 December, during the minimum of solar activity between cycles 23 and 24. In this period, more than 2 × 1011 showers were recorded with energies between ∼1 and 30 TeV. The observed two-dimensional distribution of cosmic rays is characterized by two wide regions of excess and deficit, respectively, both of relative intensity ∼10−3 with respect to a uniform flux, superimposed on smaller size structures. The harmonic analysis shows that the large-scale cosmic-ray relative intensity as a function of R.A. can be described by the first and second terms of a Fouries series. The high event statistics allow the study of the energy dependence of the anistropy, showing that the amplitude increases with energy, with a maximum intensity at ∼10 TeV, and then decreases while the phase slowly shifts toward lower values of R.A. with increasing energy. The ARGO-YBJ data provide accurate observations over more than a decade of energy around this feature of the anisotropy spectrum.
Context: 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. Aims: We take into account the state-of-the-art results of high resolution N-body dark matter simulations to calculate the most likely energy dependent boost factors, which are linked to the cosmic ray propagation properties, for anti-protons and positrons. The results from extreme, but still possible, configurations of the clumpy dark matter component are also discussed. Methods: 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. Results: 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. Dark matter clumps 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 sub-halos. Conclusions: 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. Some very extreme configurations could still lead to large enhancements, e.g. (i) very small clumps with masses ⪉10-6 M following a M-α mass distribution with α⪆ 2, highly concentrated with internal r-β profiles with β⪆ 1.5, and spatially distributed according to the smooth component; or (ii) a big sub-halo of mass ⪆ 07 M within a distance of ⪉ 1 kpc from the Earth. However, they are very unlikely from either theoretical or statistical arguments. Appendices A and B are only available in electronic form at www.aanda.org
We report on the analysis of the 10-1000 TeV large-scale sidereal anisotropy of Galactic cosmic rays (GCRs) with the data collected by the Tibet Air Shower Array from 1995 October to 2010 February. ...In this analysis, we improve the energy estimate and extend the decl. range down to −30°. We find that the anisotropy maps above 100 TeV are distinct from that at a multi-TeV band. The so-called tail-in and loss-cone features identified at low energies get less significant, and a new component appears at ∼100 TeV. The spatial distribution of the GCR intensity with an excess (7.2 pre-trial, 5.2 post-trial) and a deficit (−5.8 pre-trial) are observed in the 300 TeV anisotropy map, in close agreement with IceCube's results at 400 TeV. Combining the Tibet results in the northern sky with IceCube's results in the southern sky, we establish a full-sky picture of the anisotropy in hundreds of TeV band. We further find that the amplitude of the first order anisotropy increases sharply above ∼100 TeV, indicating a new component of the anisotropy. All these results may shed new light on understanding the origin and propagation of GCRs.
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.
ABSTRACT 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° < l ...< 100° 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° < l < 100° and 65° < l < 85° (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.
We present an updated all-particle energy spectrum of primary cosmic rays in a wide range from image to image eV using image events collected from 2000 November through 2004 October by the Tibet-III ...air-shower array located 4300 m above sea level (an atmospheric depth of 606 g cm(-2)). The size spectrum exhibits a sharp knee at a corresponding primary energy around 4 PeV. This work uses increased statistics and new simulation calculations for the analysis. We discuss our extensive Monte Carlo calculations and the model dependencies involved in the final result, assuming interaction models QGSJET01c and SIBYLL2.1, and heavy dominant (HD) and proton dominant (PD) primary composition models. Pure proton and pure iron primary models are also examined as extreme cases. A detector simulation was also performed to improve our accuracy in determining the size of the air showers and the energy of the primary particle. We confirmed that the all-particle energy spectra obtained under various plausible model parameters are not significantly different from each other, which was the expected result given the characteristics of the experiment at high altitude, where the air showers of the primary energy around the knee reach near-maximum development, with their features dominated by electromagnetic components, leading to a weak dependence on the interaction model or the primary mass. This is the highest statistical and the best systematics-controlled measurement covering the widest energy range around the knee energy region.
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 2007 November to 2013 January, 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 10degrees x 10degrees, finding a source extension sigma sub(ext)= 1degrees.8 + or - 0degrees.5. The observed differential energy spectrum is dN/dE = (2.5 + or - 0.4) x 10 super(-11) (E/1 TeV) super(-2.6+ or -0.3) photons cm super(-2) s super(-1) TeV super(-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 (SNRs) indicates that the particle acceleration inside a superbubble is similar to that in an 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.
Summary
Human brucellosis is a re‐emerging bacterial anthropozoonotic disease, which remains a public health concern in China with the growing number of cases and more widespread natural foci. The ...purpose of this study was to short‐term forecast the incidence of human brucellosis with a prediction model. We collected the annual and monthly laboratory data of confirmed cases from January 2004 to December 2013 in Shandong Diseases Reporting Information System (SDRIS). Autoregressive integrated moving average (ARIMA) model was fitted based on the monthly human brucellosis incidence from 2004 to 2013. Finally, monthly brucellosis incidences in 2014 were short‐term forecasted by the obtained model. The incidence of brucellosis was increasing from 2004 to 2013. For the ARIMA (0, 2, 1) model, the white noise diagnostic check (x2 = 5.58 P = 0.35) for residuals obtained was revealed by the optimum goodness‐of‐fit test. The monthly incidences that fitted by ARIMA (0, 2, 1) model were closely consistent with the real incidence from 2004 to 2013. And forecasting incidences from January 2014 to December 2014 were, respectively, 0.101, 0.118, 0.143, 0.166, 0.160, 0.172, 0.169, 0.133, 0.122, 0.105, 0.103 and 0.079 per100 000 population, with standard error 0.011–0.019 and mean absolute percentage error (MAPE) of 58.79%.
ABSTRACT We report on the extensive multi-wavelength observations of the blazar Markarian 421 (Mrk 421) covering radio to γ-rays, during the 4.5 year period of ARGO-YBJ and Fermi common operation ...time, from 2008 August to 2013 February. These long-term observations, extending over an energy range of 18 orders of magnitude, provide a unique chance to study the variable emission of Mrk 421. In particular, due to the ARGO-YBJ and Fermi data, the entire 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. The X-ray flux is clearly correlated with the TeV γ-ray flux, while the GeV γ-rays only show a partial correlation with the TeV γ-rays. Radio and UV fluxes seem to be weakly or not correlated with the X-ray and γ-ray fluxes. Seven large flares, including five X-ray flares and two GeV γ-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 a few keV. The TeV γ-ray flux increases up to 0.9-7.2 times the flux of the Crab Nebula. The behavior of GeV γ-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.