Z. H. Ye

The All-Particle Spectrum of Primary Cosmic Rays in the Wide Energy Range from 1014 to 1017 eV Observed with the Tibet-III Air-Shower Array

We present an updated all-particle energy spectrum of primary cosmic rays in a wide range from 10 14 to 10 17 eVusing 5:5 ; 10 7 events collected from 2000 November through 2004 October by the Tibet-III air-shower array located 4300 m abovesealevel(anatmosphericdepthof 606gcm � 2 ).Thesizespectrumexhibitsasharpkneeatacorrespondingprimary 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 thefinal result, assuming interaction models QGSJET01c and SIBYLL2.1, and heavy dominant (HD) and proton dominant (PD) primary composition models. Pure protonandpureironprimarymodelsarealsoexaminedasextremecases.A detector simulationwasalsoperformedtoimproveouraccuracyindeterminingthesizeof theairshowersandtheenergyof theprimaryparticle.Weconfirmedthatthe 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,leadingtoaweakdependenceontheinteractionmodel or theprimarymass.Thisisthehighest statistical and the best systematics-controlled measurement covering the widest energy range around the knee energy region.

Multi-TeV Gamma-Ray Observation from the Crab Nebula Using the Tibet-III Air Shower Array Finely Tuned by the Cosmic Ray Moon's Shadow

The Tibet-III air shower array, consisting of 533 scintillation detectors, has been operating successfully at Yangbajing in Tibet, China since 1999. Using the data set collected by this array from 1999 November through 2005 November, we obtained the energy spectrum of γ-rays from the Crab Nebula, expressed by a power law as (dJ/dE) = (2.09 ± 0.32) × 10–12(E/3 TeV)–2.96±0.14 cm–2 s–1 TeV–1 in the energy range of 1.7-40 TeV. This result is consistent with other independent γ-ray observations by imaging air Cherenkov telescopes. In this paper, we carefully checked and tuned the performance of the Tibet-III array using data on the Moons shadow in comparison with a detailed Monte Carlo (MC) simulation. The shadow is shifted to the west of the Moons apparent position as an effect of the geomagnetic field, although the extent of this displacement depends on the primary energy of positively charged cosmic rays. This finding enables us to estimate the systematic error in determining the primary energy from its shower size. This error is estimated to be less than ±12% in our experiment. This energy scale estimation is the first attempt among cosmic ray experiments at ground level. The systematic pointing error is also estimated to be smaller than 0011. The deficit rate and the position of the Moons shadow are shown to be very stable within a statistical error of ±6% year by year. This guarantees the long-term stability of pointlike source observation with the Tibet-III array. These systematic errors are adequately taken into account in our study of the Crab Nebula.

On temporal variations of the multi-tev cosmic ray anisotropy using the tibet iii air shower array

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays (CRs) by the 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 with 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 CRs remains insensitive to solar activities since the observation period covers more than half of the 23rd solar cycle.

Observation of TeV gamma rays from the fermi bright galactic sources with the tibet air shower array

Using the Tibet-III air shower array, we search for TeV γ-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 seven sources instead of the expected 0.61 sources at a significance of 2σ or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 × 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 × 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 γ-ray excesses. We also find that all seven sources are associated with pulsars, and six of them are coincident with sources detected by the Milagro experiment at a significance of 3σ 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 ≥3σ 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.

Implication of the sidereal anisotropy of ∼5 TeV cosmic ray intensity observed with the Tibet III air shower array

We show that the large‐scale anisotropy of ∼5 TeV galactic cosmic ray (GCR) intensity observed by Tibet Air Shower experiment can be reproduced by the superposition of a bi‐directional and uni‐directional flows (UDF and BDF) of GCRs. The heliosphere is located inside the local interstellar cloud (LIC) very close to the inner edge of the LIC. If the GCR population is lower inside the LIC than outside, the BDF flow is expected from the parallel diffusion of GCRs into LIC along the local interstellar magnetic field (LISMF) connecting the heliosphere with the region outside the LIC, where the GCR population is higher. A type of the UDF, on the other hand, is expected from the B×∇n drift flux driven by a gradient of GCR density (n) in the LISMF (B). The LISMF orientation deduced from the best‐fit direction of the BDF is almost parallel to the galactic plane and more consistent with the suggestion of Frisch (1996) than that of Lallement et al. (2005). We note that the model, if holds, yields the LISMF polarity ...

Moon shadow by cosmic rays under the influence of geomagnetic field and search for antiprotons at multi-TeV energies

We have observed the shadowing of galactic cosmic ray flux in the direction of the moon, the so-called moon shadow, using the Tibet-III air shower array operating at Yangbajing (4300 in a.s.l.) in Tibet since 1999. Almost all cosmic rays are positively charged; for that reason, they are bent by the geomagnetic field, thereby shifting the moon shadow westward. The cosmic rays will also produce an additional shadow in the eastward direction of the moon if cosmic rays contain negatively charged particles, such as antiprotons, with some fraction. We selected 1.5 x 10(10) air shower events with energy beyond about 3 TeV from the dataset observed by the Tibet-III air shower array and detected the moon shadow at similar to 40 sigma level. The center of the moon was detected in the direction away from the apparent center of the moon by 0.23 degrees to the west. Based on these data and a full Monte Carlo simulation, we searched for the existence of the shadow produced by antiprotons at the multi-TeV energy region. No evidence of the existence of antiprotons was found in this energy region. We obtained the 90% confidence level upper limit of the flux ratio of antiprotons to protons as 7% at multi-TeV energies.

Tibet air shower array: results and future plan

The Tibet air shower array, which has an effective area of 36,900 m2, has been in operation at Yangbajing in Tibet, China at an altitude of 4,300 m above sea level. In this paper, we will briefly introduce the recent gamma-ray observation with the present Tibet air shower array and our future plan which is called the Tibet muon detector (MD) project.

The energy spectrum of all-particle cosmic rays around the knee region observed with the Tibet-III air-shower array

We have already reported the first result on the all-particle spectrum around the knee region based on data from 2000 November to 2001 October observed by the Tibet-III air-shower array. In this paper, we present an updated result using data set collected in the period from 2000 November through 2004 October in a wide range over 3 decades between 1014 and 1017 eV, in which the position of the knee is clearly seen at around 4 PeV. The spectral index is 2.68 ± 0.02(stat.) below 1 PeV, while it is 3.12 ± 0.01(stat.) above 4 PeV in the case of QGSJET+HD model, and various systematic errors should be carefully done in the very near future.

Comparison of methods for a 3-D density inversion from airborne gravity gradiometry

In this study, we apply Tikhonov’s regularization algorithm for a 3-D density inversion from the gravity-gradiometry data. To reduce the non-uniqueness of the inverse solution (carried out without additional information from geological evidence), we implement the depth-weighting empirical function. However, the application of an empirical function in the inversion equation brings the bias problem of the regularization factor when a traditional Tikhonov’s algorithm is applied. To solve the bias problem of regularization factor selection, we present a standardized solution that comprises two parts for solving a 3-D constrained inversion equation, specifically the linear matrix transformation and Tikhonov’s regularization algorithm. Since traditional regularization techniques become numerically inefficient when dealing with large number of data, we further apply methods which include the Simultaneous Iterative Reconstruction Technique (SIRT) and the wavelet compression combined with Least Squares QR-decomposition (LSQR). In our simulation study, we demonstrate that SIRT as well as the wavelet compression plus LSQR algorithm improve the computation efficiency, while provide results which closely agree with that obtained from applying Tikhonov’s regularization. In particular, the algorithm of wavelet compression plus LSQR shows the best computing efficiency, because it combines the advantages of coefficients compression of big matrix and fast solution of sparse matrix. Similar findings are confirmed from the vertical gravity gradient data inversion for detecting potential deposits at the Kauring (near Perth, Western Australia) testing site.

TeV Gamma Ray Survey on the Direction of Fermi-LAT Pulsars with the Tibet Air Shower Array

We search for steady TeV γ-rays from 18 pulsars in the Fermi-LAT pulsar catalog using the Tibet air shower array. We observe 8 sources including the Crab against expected 0.41 sources at a significance of 2σ or more among the 18 pulsars. The chance probability from Poisson statistics is estimated to be 1.4×10-8.With the Crab excluded, it is estimated to be 1.8×10-7. These low chance probabilities clearly show that the Fermi pulsars have a statistically significant correlation with TeV γ-ray excesses observed by the Tibet air shower array.