Search

HESS J1843–033 is a very high energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond 100 TeV from the HESS J1843–033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of 0.°34 ± 0.°12 is successfully detected above 25 TeV at (α, δ) = (281.°09 ± 0.°10, −3.°76 ± 0.°09) near HESS J1843–033 with a statistical significance of 6.2σ, and the source is named TASG J1844–038. The position of TASG J1844–038 is consistent with those of HESS J1843–033, eHWC J1842–035, and LHAASO J1843–0338. The measured gamma-ray energy spectrum in 25 TeV < E < 130 TeV is described with dN / dE = ( 9.70 ± 1.89 ) × 10 − 16 (E/40 TeV)−3.26±0.30 TeV−1 cm−2 s−1, and the spectral fit to the combined spectra of HESS J1843–033, LHAASO J1843–0338, and TASG J1844–038 implies the existence of a cutoff at 49.5 ± 9.0 TeV. Associations of TASG J1844–038 with SNR G28.6–0.1 and PSR J1844–0346 are also discussed in detail for the first time.

Cosmic rays (protons and other atomic nuclei) are believed to gain energies of petaelectronvolts (PeV) and beyond at astrophysical particle accelerators called 'PeVatrons' inside our Galaxy. Although a characteristic feature of a PeVatron is expected to be a hard gamma-ray energy spectrum that extends beyond 100 teraelectronvolts (TeV) without a cutoff, none of the currently known sources exhibits such a spectrum due to the low maximum energy of accelerated cosmic rays or insufficient detector sensitivity around 100 TeV. Here we report the observation of gamma-ray emission from the supernova remnant G106.3+2.7 above 10 TeV. This work provides flux data points up to and above 100 TeV and indicates that the very-high-energy gamma-ray emission above 10 TeV is well correlated with a molecular cloud rather than the pulsar PSR J2229+6114. Regarding the gamma-ray emission mechanism of G106.3+2.7, this morphological feature appears to favor a hadronic origin via the {\pi}0 decay caused by accelerated relativistic protons over a leptonic one via the inverse-Compton scattering by relativistic electrons. Furthermore, we point out that an X-ray flux upper limit on the synchrotron spectrum would provide important information to firmly establish the hadronic scenario as the mechanism of particle acceleration at the source., Comment: Published in Nature Astronomy

Using the Tibet-III air shower array, we search for steady TeV -rays from 18 pulsars in the Fermi Large Area Telescope pulsar catalog. Among them, we ob- serve 8 sources including the Crab instead of the expected

The observation of the extensive air showers has been carried out since 1990 by Tibet air-shower array to study the origin of the cosmic-rays. In order to investigate the chemical composition around the knee, air-shower core detector was used in two experimental phases aiming to separate light primary elements (proton and helium). In a first phase experiment made in 1996-1999, primary proton and helium spectra were measured using a hybrid detector consisted of Tibet II air-shower array, emulsion chambers, and burst detectors. The result indicated that the fraction of light elements among the all particle spectrum is less than 30%. We performed the second phase experiment in 2002-2005 to measure proton + helium together to confirm the previous result with higher statistics by improving the detection efficiency for the air-shower cores. The P + He spectrum was obtained in a good agreement with previous result by several times higher statistics. The flux of the light elements, together with the all-particle energy spectrum, indicates the dominance of nuclei heavier than helium in the energy range around the knee. We also present the outline of the next phase experiment, which is to measure the heavy elements around the knee, using low threshold air-shower core detector named YAC. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.

The Tibet-III air shower array, consisting of 533 scintillation detectors, has been operating successfully at Yangbajing in Tibet, China since 1999. Using the dataset collected by this array from 1999 November through 2005 November, we obtained the energy spectrum of $\gamma$-rays from the Crab Nebula, expressed by a power law as $(dJ/dE) = (2.09\pm0.32)\times10^{-12} (E/{\rm 3 TeV})^{-2.96\pm0.14} {\rm cm}^{-2} {\rm s}^{-1} {\rm TeV}^{-1}$ in the energy range of 1.7 to 40 TeV. This result is consistent with other independent $\gamma$-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 moon's shadow in comparison with a detailed Monte Carlo simulation. The shadow is shifted to the west of the moon's apparent position as an effect of the geomagnetic field, although the extent of this displacement depends on the primary energy 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 $\pm$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 $0\fdg011$. The deficit rate and position of the moon's shadow are shown to be very stable within a statistical error of $\pm$6% year by year. This guarantees the long-term stability of point-like source observation with the Tibet-III array. These systematic errors are adequately taken into account in our study of the Crab Nebula., Comment: 12 pages, 16 figures, accepted by ApJ

We report the detection of 10 new X-ray filaments using the data from the {\sl Chandra} X-ray satellite for the inner $6^{\prime}$ ($\sim 15$ parsec) around the Galactic center (GC). All these X-ray filaments are characterized by non-thermal energy spectra, and most of them have point-like features at their heads that point inward. Fitted with the simple absorbed power-law model, the measured X-ray flux from an individual filament in the 2-10 keV band is $\sim 2.8\times10^{-14}$ to $10^{-13}$ ergs cm$^{-2}$ s$^{-1}$ and the absorption-corrected X-ray luminosity is $\sim 10^{32}-10^{33}$ ergs s$^{-1}$ at a presumed distance of 8 kpc to the GC. We speculate the origin(s) of these filaments by morphologies and by comparing their X-ray images with the corresponding radio and infrared images. On the basis of combined information available, we suspect that these X-ray filaments might be pulsar wind nebulae (PWNe) associated with pulsars of age $10^3 \sim 3\times 10^5$ yr. The fact that most of the filament tails point outward may further suggest a high velocity wind blowing away form the GC.

Ulysses observed a stable strong CIR from early 1992 through 1994 during its first journey into the southern hemisphere. After the rapid latitude scan in early 1995, Ulysses observed a weaker CIR from early 1996 to mid-1997 in the northern hemisphere as it traveled back to the ecliptic at the orbit of Jupiter. These two CIRs are the observational basis of the investigation into the latitudinal structure of CIRs. The first CIR was caused by an extension of the northern coronal hole into the southern hemisphere during declining solar activity, whereas the second CIR near solar minimum activity was caused by small warps in the streamer belt. The latitudinal structure is described through the presentation of three 26-day periods during the southern CIR. The first at 24 Ss hows the full plasma interaction region including fast and slow wind streams, the compressed shocked flows with embedded stream interface and heliospheric current sheet (HCS), and the forward and reverse shocks with associated accelerated ions and electrons. The second at 40 S exhibits only the reverse shock, accelerated particles, and the 26-day modulation of cosmic rays. The third at 60 Ss hows only the accelerated particles and modulated cosmic rays. The possible mechanisms for the access of the accelerated particles and the CIR-modulated cosmic rays to high latitudes above the plasma

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., 18 pages, 2 figures, accepted by The Astrophysical Journal

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., Comment: This paper has been accepted by hepnp

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 $10^{14}$ eV and $10^{17}$ eV, in which the position of the knee is clearly seen at around 4 PeV. The spectral index is -2.68 $\pm$ 0.02(stat.) below 1PeV, while it is -3.12 $\pm$ 0.01(stat.) above 4 PeV in the case of QGSJET+HD model, and various systematic errors are under study now., Comment: 12 pages, 7 figures, accepted by Advances in space research

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., accepted to ApJL

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 m 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 x10^{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 $\sim 40 \sigma$ level. The center of the moon was detected in the direction away from the apparent center of the moon by 0.23$^\circ$ 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., Comment: 13pages,4figures; Accepted for publication in Astroparticle Physics

We consider standing wave solutions of the linearized magnetohydrodynamic (MHD) equations describing a perfectly conducting but viscous magnetofluid in a gravitationally-stratified cylindrical plasma. We assume a relatively strong homogeneous background magnetic field within the flux tube, which is oriented along the cylinder axis; the background atmosphere is in hydrostatic equilibrium, with gravity uniform and parallel to the cylinder axis; the vertical temperature profile (constrained by the stability conditions) is arbitrary. The MHD wave equations are separable, and appropriate boundary conditions are formulated. We employ a numerical shooting scheme to solve the coupled ordinary differential equations and to obtain the damping rates and eigenfrequency spectra for the standing MHD wave modes. Three series of modes are obtained for an isothermal flux tube: Alfvtnic, quasi-fast, and quasi-slow. The quasi-slow mode series is novel, having been missed previously. In the inviscid limit, with a st...

This paper considers two-dimensional nonlinear MHD waves of large horizontal spatial scales for a thin magnetofluid layer on the surface of a rotating sphere. The 'shallow fluid' hydrodynamic equations are generalized to include the effects of magnetic fields, and it is shown that the resulting MHD equations can be reduced to a single scalar equation for a stream function involving several free functions. For special choices of these free functions, two kinds of finite-amplitude MHD waves are obtained, propagating in the azimuthal direction relative to the uniformly rotating background atmosphere in the presence of a background zonal magnetic field and a steady differential zonal flow. These two kinds of MHD waves are fundamentally due to the joint effects of the uniform rotation of the background atmosphere and background magnetic field; the first is an inertial wave of the Rossby (1939) and Haurwitz (1940) type, modified by the presence of the background zonal magnetic field, while the second is a magnetic Alfven-like wave which is modified by the uniform rotation of the background atmosphere.

A computational code for solving two-dimensional, time-dependent MHD equations by the method of characteristics is presented. Its capabilities are demonstrated by solving two very different problems for which analytical solutions exist: linearized, standing MHD wave motions in a magnetized cylindrical plasma, and nonlinear self-similar expansion of a magnetized plasma ball. The nonlinear development of standing MHD wave solutions in a cylindrical plasma is also studied. The method can be naturally embedded in the computational architecture of massively parallel processors.