宇宙线直接探测进展概述

宇宙线从发现起至今已超过百年。在20世纪上半叶,大型粒子加速器技术成熟以前,对宇宙线的研究引领着基本粒子物理的发展,从宇宙线研究中取得的多项成果斩获诺贝尔奖。21世纪,宇宙线因其与极端高能的物理规律和暗物质等新物理现象联系密切而绽放出新的活力,宇宙线起源、加速、传播等相关的天文学及物理学问题也备受关注。简述了近年来在空间直接观测宇宙线实验方面取得的进展,以及其对理解宇宙线物理问题的推动。最后概述了中国在相关领域的研究历程和现状。     

Overview of Direct Measurements of Cosmic Rays

The history of cosmic ray studies can be traced back to the 1910s when Hess and other scientists first discovered them. Cosmic rays are very important laboratories of particle physics, and have led to many important discoveries of fundamental particles, such as the positrons, muons, pions, and a series of strange particles. Cosmic rays are nowadays the key probes of the extremely high-energy physics and dark matter particles. A brief review about the history and recent progresses of direct observations of cosmic rays is presented. In recent years, the new space-borne experiments such as PAMELA and AMS-02, as well as a few of balloon-borne experiments, have measured the energy spectra of cosmic rays very precisely, and revealed several new features/anomalies. Remarkable excesses of positron fraction in the total electron plus positron fluxes have been observed, which may be caused by the annihilation/decay of dark matter particles or by astrophysical pulsars. The cosmic ray antiprotons, which are expected to have the same secondary origin as that of positrons, do not show significant excesses compared with the background prediction. This result also constrains the modeling of the positron excesses. In addition, the spectral hardening above several hundred GeV of cosmic ray nuclei has been revealed. These results have important and interesting implications on our understandings of the origin, acceleration, and propagation of cosmic rays. In particular, China has launched the Dark Matter Particle Explorer (DAMPE) to indirectly search for the dark matter and explore the high-energy universe in the TeV window. Most recently, the DAMPE collaborators reported the new measurements of the cosmic ray electron plus positron fluxes up to about 5 TeV with a very high precision. The DAMPE data revealed clearly a deflection around 0.9 TeV in the electron energy spectrum. Possible fine structures of the electron plus positron spectra can be critically addressed with the accumulation of data in the coming years.     

Inverse Compton e± pair cascade model for the γ-ray production in massive binary LSI +61° 303

We apply an inverse Compton   e ^±  pair cascade model for γ-ray production in the massive binary system LSI +61° 303 assuming that electrons are accelerated already inside the inner part of the jet launched by the compact object. γ-ray spectra, affected by the cascade process, and lower energy spectra, from the synchrotron cooling of the highest energy electrons in the jet, are calculated as a function of the phase of this binary system. γ-ray spectra expected in such a model have different shape than those ones produced by electrons in the jet directly to observer. Moreover, the model predicts clear anticorrelation between γ-ray fluxes in the GeV (1–10 GeV) and TeV (>200 GeV) energy ranges with the peak of the TeV emission at the phase ∼0.5 (the peak half-width ranges between the phases ∼0.4–0.9, for the inclination of the binary system equal to 60°, and ∼0.4–0.1 for 30°). The fine features of TeV γ-ray emission (fluxes and spectral shapes) as a function of the phase of the binary system are consistent with recent observations reported by the MAGIC collaboration. Future simultaneous observations in the GeV energies (by the GLAST and AGILE telescopes) and in the TeV energies (by the MAGIC and VERITAS telescopes) should test other predictions of the considered model supporting or disproving the hypothesis of acceleration of electrons already in the inner part of the microquasar jets.     

TeV gamma-rays from accreting magnetars in massive binaries

This paper focuses on neutron stars (NS) of the magnetar type inside massive binary systems. We determine the conditions under which the matter from the stellar wind can penetrate the inner magnetosphere of the magnetar. At a certain distance from the NS surface, the magnetic pressure can balance the gravitational pressure of the accreting matter, creating a very turbulent, magnetized transition region. It is suggested that this region provides good conditions for the acceleration of electrons to relativistic energies. These electrons lose energy due to the synchrotron process and inverse Compton (IC) scattering of the radiation from the nearby massive stellar companion, producing high-energy radiation from X-rays up to ∼TeV γ-rays. The primary γ-rays can be further absorbed in the stellar radiation field, developing an IC  e^±  pair cascade. We calculate the synchrotron X-ray emission from primary electrons and secondary  e^±  pairs and the IC γ-ray emission from the cascade process. It is shown that quasi-simultaneous observations of the TeV γ-ray binary system LSI +61 303 in the X-ray and TeV γ-ray energy ranges can be explained with such an accreting magnetar model.     

Hard X-ray observations of PSR J1833−1034 and its associated pulsar wind nebula

PSR J1833−1034 and its associated pulsar wind nebula (PWN) have been investigated in depth through X-ray observations ranging from 0.1 to 200 keV. The low-energy X-ray data from Chandra reveal a complex morphology that is characterized by a bright central plerion, no thermal shell and an extended diffuse halo. The spectral emission from the central plerion softens with radial distance from the pulsar, with the spectral index ranging from  Γ= 1.61  in the central region to  Γ= 2.36  at the edge of the PWN. At higher energy, INTEGRAL detected the source in the 17–200 keV range. The data analysis clearly shows that the main contribution to the spectral emission in the hard X-ray energy range is originated from the PWN, while the pulsar is dominant above 200 keV. Recent High Energy Stereoscopic System (HESS) observations in the high-energy gamma-ray domain show that PSR J1833−1034 is a bright TeV emitter, with a flux corresponding to ∼2 per cent of the Crab in 1–10 TeV range. In addition, the spectral shape in the TeV energy region matches well with that in the hard X-rays observed by INTEGRAL . Based on these findings, we conclude that the emission from the pulsar and its associated PWN can be described in a scenario where hard X-rays are produced through synchrotron light of electrons with Lorentz factor  γ∼ 10⁹  in a magnetic field of ∼10 μG. In this hypothesis, the TeV emission is due to inverse-Compton interaction of the cooled electrons off the cosmic microwave background photons. Search for PSR J1833−1034 X-ray pulsed emission, via RXTE and Swift X-ray observations, resulted in an upper limit that is about 50 per cent.     

Modulation of the high energy γ-ray flux from PSR B1259-63/SS2883 due to the orbital variation of the maximum energy of accelerated electrons

The inverse Compton (IC) scattering of ultrarelativistic electrons accelerated at the pulsar wind termination shock is generally believed to be responsible for TeV gamma-ray signal recently reported from the binary system PSR B1259-63/SS2883. In such a system the acceleration takes place in the presence of a dense radiation field provided by a companion Be2-type star. Thus it is natural to expect an orbital phase dependence of the acceleration efficiency in the system. The HESS collaboration reported the tendency of reduction of TeV γ-rays around the periastron. In this paper we study a possible explanation of this effect by the “early” (sub-TeV) cutoffs in the energy spectrum of accelerated electrons due to the enhanced rate of Compton losses close to the periastron.     

On the formation of TeV radiation in LS 5039

The recent detections of TeV gamma-rays from compact binary systems show that relativistic outflows (jets or winds) are sites of effective acceleration of particles up to multi-TeV energies. In this paper, we discuss the conditions of acceleration and radiation of ultrarelativistic electrons in LS 5039, the gamma-ray emitting binary system for which the highest quality TeV data are available. Assuming that the gamma-ray emitter is a jet-like structure, we performed detailed numerical calculations of the energy spectrum and light curves accounting for the acceleration efficiency, the location of the accelerator, the speed of the emitting flow, the inclination angle of the system, as well as specific features related to anisotropic inverse Compton (IC) scattering and pair production. We conclude that the accelerator should not be deep inside the binary system unless we assume a very efficient acceleration rate. We show that within the IC scenario both the gamma-ray spectrum and flux are strongly orbital phase dependent. Formally, our model can reproduce, for specific sets of parameter values, the energy spectrum of gamma-rays reported by HESS for wide orbital phase intervals. However, the physical properties of the source can be constrained only by observations capable of providing detailed energy spectra for narrow orbital phase intervals (Δφ≪ 0.1).     

Acceleration of particles in the vicinity of a massive black hole

Recent results of the gamma-ray Cherenkov astronomy definitely prove the existence of fast variability in the very high energy (V.H.E.) gamma-ray flux of some active galactic nuclei. The BL Lac PKS 2155-304 for instance showed variations down to a few minutes time scale. From standard light travel time argument, these variations put extremely strong constraints on the size of the TeV emitting zone, which has to be of the order of a few Schwarzschild radius, even for high values of the relativistic Doppler factor of the emitting jets. Such discovery is a challenge for particle acceleration scenarios, which have to imagine efficient acceleration processes at work in a very compact zone. Eventually, the immediate vicinity of the central black hole appears as the most conservative choice for the location of the TeV emission region of active galactic nuclei. In this paper, we propose a two-step mechanism for charged particle acceleration in the magnetosphere of a massive black hole surrounded by an accretion disk. Particles first gain energy by a stochastic process during the accretion phase. It is shown that effective proton acceleration up to energies 10¹⁷–10¹⁹ eV is possible in a low-luminosity magnetized accretion disk with 2D turbulent motion. The distribution function of energetic protons over energies is a power law function with typical index ≃−1. Here electrons are not very efficiently accelerated because of their drastic losses by synchrotron radiation. In a second time, part of the fast particles escape from the disk and are then entrained by the magnetic structure above the disk, in the rotating black hole magnetosphere. They thus gain additional energy by direct centrifugal mechanism, up to about 10²⁰ eV for the protons and to 10–100 TeV for the electrons when they cross the light cylinder surface. Such energetic particles can further radiate in the TeV spectral range observed by Cherenkov experiments as HESS, MAGIC and VERITAS. Energetic protons can produce γ-radiation in the energy band 1 GeV–100 TeV and above mainly by nuclei collisions with the disk matter, clouds, or ambient low energy photons. Energetic electrons can also reach the required spectral range by inverse Compton emission. However their acceleration is less efficient due to heavy radiation losses, and only gained by centrifugal process during the second phase of the whole mechanism we describe. Our present analysis would therefore favor hadronic scenarios for TeV emission of active galactic nuclei. It is tempting to relate long term variability over years of TeV active galactic nuclei to the first stochastic acceleration phase, which also provides the needed power law particle distributions, while short term variability over minutes is more likely due to perturbations of the second fast direct acceleration phase.     

The gamma ray background from large scale structure formation

Hierarchical clustering of dark matter halos is thought to describe well the large scale structure of the universe. The baryonic component of the halos is shock heated to the virial temperature while a small fraction of the energy flux through the shocks may be energized through the first order Fermi process to relativistic energy per particle. It has been proposed that the electrons accelerated in this way may upscatter the photons of the universal microwave background to gamma ray energies and indeed generate a diffuse background of gamma rays that compares well to the observations. In this paper, we calculate the spectra of the particles accelerated at the merger shocks and re-evaluate the contribution of structure formation to the extragalactic diffuse gamma ray background (EDGRB), concluding that this contribution adds up to at most 10% of the observed EDGRB.     

In-flight measurement of the absolute energy scale of the Fermi Large Area Telescope

The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope is a pair-conversion telescope designed to survey the gamma-ray sky from 20 MeV to several hundreds of GeV. In this energy band there are no astronomical sources with sufficiently well known and sharp spectral features to allow an absolute calibration of the LAT energy scale. However, the geomagnetic cutoff in the cosmic ray electron-plus-positron (CRE) spectrum in low Earth orbit does provide such a spectral feature. The energy and spectral shape of this cutoff can be calculated with the aid of a numerical code tracing charged particles in the Earth’s magnetic field. By comparing the cutoff value with that measured by the LAT in different geomagnetic positions, we have obtained several calibration points between ∼6 and ∼13 GeV with an estimated uncertainty of ∼2%. An energy calibration with such high accuracy reduces the systematic uncertainty in LAT measurements of, for example, the spectral cutoff in the emission from gamma ray pulsars.     

Energy correction based on fluorescence attenuation of DAMPE

The major scientific goals of the DArk Matter Particle Explorer(DAMPE) are to study cosmicray electrons(including positrons) and gamma rays from 5 GeV to 10 TeV and nuclei from Z = 1 to 26 up to 100 TeV. The deposited energy measured by the Bismuth Germanate Oxide(BGO) calorimeter of DAMPE is affected by fluorescence attenuation in BGO crystals that are 600 mm long. In this work, an in-orbit attenuation calibration method is reported, and energy correction of the sensitive detector unit of the BGO calorimeter is also presented.     

Multi‐Messenger Cosmic Particles

A54 Cosmic Ray Acceleration in Galactic Wind Shocks A71 Detection of Ultra‐High Energy Cosmic Rays and Neutrinos with LOFAR A80 Status of the gravitational‐wave detector GEO600 A87 Recent Results and Future of the MAGIC gamma‐ray telescope A92 Cosmic ray detection with the radio technique A93 Cosmic Ray Physics with IceCube A94 The resonance‐like gamma‐ray absorption processes for use in astrophysics A97 Geometry reconstruction of air shower fluorescence detectors revisited A102 Supermassive Binary Black Holes & Radio Jets A108 Muonic Component of Air Showers Measured by KASCADE‐Grande A110 Towards new frontiers: observation of photons with energies above 10¹⁸ eV A112 The IceCube Neutrino Telescope A114 The ground‐based gamma‐ray observatory CTA A116 IceCube: Recent Results and Prospects A117 Particle Physics with AMANDA and IceCube A118 Altitude dependence of fluorescence light emission by extensive air showers A120 Neutrino‐induced cascades in AMANDA & IceCube A122 Enhancement Telescopes for the Pierre Auger Southern Observatory in Argentina A123 Proton spectra from relativistic shock environments in AGN and GRBs A124 The Baikal Neutrino Telescope – Physics Results A127 Searches for point‐like sources of cosmic neutrinos with IceCube A128 The MAGIC/IceCube Target of Opportunity Programtest run A131 Supernova detection with IceCube: from low to high energy neutrinos A132 Measurement of the UHECR energy spectrum from hybrid data of the Pierre Auger Observatory A133 Extension of IceCube at Lower Energy: the Use of AMANDA as Nested Array and the Future Prospectives A135 Searching for neutrinos with the Pierre Auger Observatory A138 Search for Transient Emission of Neutrinos in IceCube A140 Acoustic Neutrino Detection in Antarctic Ice A159 AMANDA limits on the diffuse muon‐neutrino flux: physics implications A164 Investigation of the Radio Emission of Cosmic Ray Air Showers with LOPES A168 The Northern Site of the Pierre Auger Observatory A170 Shower reconstruction and size spectra with KASCADE‐Grande data A171 Neutrinos from Gamma Ray Bursts: predictions and limits from AMANDA‐II data A172 Simulation study of shower profiles from ultra‐high energy cosmic rays A174 Upper limit to the photon fraction in cosmic rays above 10¹⁹ eV from the Pierre Auger Observatory A176 Astrophysics at MeV energies A180 Study of the Cosmic Ray Composition above 0.4 EeV using the Longitudinal Profiles of Showers observed at the Pierre Auger Observatory A185 Backgrounds for UHE horizontal neutrino showers A186 The Front‐End Cards of the Pierre Auger Surface Detectors: Test Results and Performance in the Field A187 Monte Carlo Studies for MAGIC‐II A194 Measuring the proton‐air cross section from logitudinal air shower profiles A195 The UHECR energy spectrummeasured at the Pierre Auger Observatory A203 Highlights of Observations of Galactic Sources with the MAGIC telescope A207 Adesign study for a 12.5 m ∅︁ Imaging Air Cherenkov Telescope for ground‐based γ ‐ray astronomy A210 The Future of Long‐Wavelengths Radio‐Astronomy in Germany: LOFAR and GLOW A211 Online Monitoring of the Pierre Auger Observatory A216 OPTIMA‐Burst – Catching GRB Afterglows (and other Transients) with High Time Resolution A227 JEM‐EUSO mission A232 Rapid Variations in AGN: Clues on Particle Accelerators A235 Systematic search forVHEgamma‐ray emission from X‐ray bright high‐frequency peaked BL Lac objects A237 Prospects for GeV Astronomy in the Era of GLAST A241 Improvements of the energy reconstruction for the MAGIC telescope by means of analysis and Monte Carlo techniques A265 Discovery of VHE γ ‐rays from BL Lacertae with the MAGIC telescope A266 Results of two observation cycles of LS I+61°303 with the MAGIC telescope A267 Wide Range Multifrequency Observations of Northern TeV Blazars A269 Diffusive and convective cosmic ray transport in elliptical galaxies     

TeV γ-ray fluxes from the long campaigns on Mrk 421 as constraints on the emission of TeV–PeV neutrinos and UHECRs

Long TeV γ-ray campaigns have been carried out to study the spectrum, variability and duty cycle of the BL Lac object Markarian 421. These campaigns have given some evidence of the presence of protons in the jet: (i) Its spectral energy distribution which shows two main peaks; one at low energies (∼1 keV) and the other at high energies (hundreds of GeV), has been described by using synchrotron proton blazar model. (ii) The study of the variability at GeV γ-rays and X-rays has indicated no significant correlation. (iii) TeV γ-ray detections without activity in X-rays, called “orphan flares” have been observed in this object.Recently, The Telescope Array Collaboration reported the arrival of 72 ultra-high-energy cosmic rays with some of them possibly related to the direction of Markarian 421. The IceCube Collaboration reported the detection of 37 extraterrestrial neutrinos in the TeV–PeV energy range collected during three consecutive years. In particular, no neutrino track events were associated with this source. In this paper, we consider the proton–photon interactions to correlate the TeV γ-ray fluxes reported by long campaigns with the neutrino and ultra-high-energy cosmic ray observations around this blazar. Considering the results reported by The IceCube and Telescope Array Collaborations, we found that only from ∼25% to 70% of TeV fluxes described with a power law function with exponential cutoff can come from the proton–photon interactions.     

Improved energy resolution for VHE gamma ray astronomy with systems of Cherenkov telescopes

We present analysis techniques to improve the energy resolution of stereoscopic systems of imaging atmospheric Cherenkov telescopes, using the HEGRA telescope system as an example. The techniques include (i) the determination of the height of the shower maximum, which is then taken into account in the energy determination, and (ii) the determination of the location of the shower core with the additional constraint that the direction of the gamma rays is known a priori. This constraint can be applied for gamma ray point sources, and results in a significant improvement in the localization of the shower core, which translates into better energy resolution. Combining both techniques, the HEGRA telescopes reach an energy resolution between 9% and 12%, over the entire energy range from 1 TeV to almost 100 TeV. Options for further improvements of the energy resolution are discussed.     

The multi-band non-thermal emission from the supernova remnant RX J1713.7−3946

Non-thermal X-rays and very high energy (VHE) γ-rays have been detected from the supernova remnant (SNR) RX J1713.7−3946, and the recent observations with the Suzaku satellite clearly reveal a spectral cut-off in the X-ray spectrum, which directly relates to the cut-off of the energy spectrum of the parent electrons. However, whether the origin of the VHE γ-rays from the SNR is hadronic or leptonic is still in debate. We studied the multi-band non-thermal emission from RX J1713.7−3946 based on a semi-analytical approach towards the non-linear shock acceleration process by including the contribution of the accelerated electrons to the non-thermal radiation. The results show that the multi-band observations on RX J1713.7−3946 can be well explained in the model with appropriate parameters, and the TeV γ-rays have hadronic origin, i.e. they are produced via proton–proton interactions as the relativistic protons accelerated by the shock collide with the ambient matter.     

New constraints on space–time Planck scale fluctuations from established high energy astronomy observations

The space–time metric is widely believed to be subject to stochastic fluctuations induced by quantum gravity at the Planck scale. This work is based on two different phenomenological approaches being currently made to this topic, and theoretical models which describe this phenomenon are not dealt with here. By using the idea developed in one of these two approaches in the framework of the other one, it is shown that the constraints on the nature of Planck scale space–time fluctuations already set by the observation of electrons and gamma-rays with energies above 15 TeV are much stronger than have been shown so far. It is concluded that for the kind of Planck scale fluctuations implied by several models, including the most naive one, to be consistent with the observations, the transformation laws between different reference frames must be modified in order to let the Planck scale be observer independent.     

Polarisation measurements with a CdTe pixel array detector for Laue hard X-ray focusing telescopes

Polarimetry is an area of high energy astrophysics which is still relatively unexplored, even though it is recognized that this type of measurement could drastically increase our knowledge of the physics and geometry of high energy sources. For this reason, in the context of the design of a Gamma-Ray Imager based on new hard-X and soft gamma ray focusing optics for the next ESA Cosmic Vision call for proposals (Cosmic Vision 2015-2025), it is important that this capability should be implemented in the principal on-board instrumentation. For the particular case of wide band-pass Laue optics we propose a focal plane based on a thick pixelated CdTe detector operating with high efficiency between 60–600keV. The high segmentation of this type of detector (1–2mm pixel size) and the good energy resolution (a few keV FWHM at 500keV) will allow high sensitivity polarisation measurements (a few % for a 10mCrab source in 10⁶s) to be performed. We have evaluated the modulation Q factors and minimum detectable polarisation through the use of Monte Carlo simulations (based on the GEANT 4 toolkit) for on and off-axis sources with power law emission spectra using the point spread function of a Laue lens in a feasible configuration.     

TeV γ-rays from old supernova remnants

We study the emission from an old supernova remnant (SNR) with an age of around 10⁵ yr and that from a giant molecular cloud (GMC) encountered by the SNR. When the SNR age is around 10⁵ yr, proton acceleration is efficient enough to emit TeV γ-rays both at the shock of the SNR and that in the GMC. The maximum energy of primarily accelerated electrons is so small that TeV γ-rays and X-rays are dominated by hadronic processes,  π⁰  -decay and synchrotron radiation from secondary electrons, respectively. However, if the SNR is older than several 10⁵ yr, there are few high-energy particles emitting TeV γ-rays because of the energy-loss effect and/or the wave-damping effect occurring at low-velocity isothermal shocks. For old SNRs or SNR–GMC interacting systems capable of generating TeV γ-ray emitting particles, we calculated the ratio of TeV γ-ray (1–10 TeV) to X-ray (2–10 keV) energy flux and found that it can be more than  ∼10²  . Such a source showing large flux ratio may be a possible origin of recently discovered unidentified TeV sources.     

Spiral arms as cosmic ray source distributions

The Milky Way is a spiral galaxy with (or without) a bar-like central structure. There is evidence that the distribution of suspected cosmic ray sources, such as supernova remnants, are associated with the spiral arm structure of galaxies. It is yet not clearly understood what effect such a cosmic ray source distribution has on the particle transport in our Galaxy. We investigate and measure how the propagation of Galactic cosmic rays is affected by a cosmic ray source distribution associated with spiral arm structures.We use the PICARD code to perform high-resolution 3D simulations of electrons and protons in galactic propagation scenarios that include four-arm and two-arm logarithmic spiral cosmic ray source distributions with and without a central bar structure as well as the spiral arm configuration of the NE2001 model for the distribution of free electrons in the Milky Way. Results of these simulation are compared to an axisymmetric radial source distribution. Also, effects on the cosmic ray flux and spectra due to different positions of the Earth relative to the spiral structure are studied.We find that high energy electrons are strongly confined to their sources and the obtained spectra largely depend on the Earth’s position relative to the spiral arms. Similar finding have been obtained for low energy protons and electrons albeit at smaller magnitude. We find that even fractional contributions of a spiral arm component to the total cosmic ray source distribution influences the spectra on the Earth. This is apparent when compared to an axisymmetric radial source distribution as well as with respect to the Earth’s position relative to the spiral arm structure. We demonstrate that the presence of a Galactic bar manifests itself as an overall excess of low energy electrons at the Earth.Using a spiral arm geometry as a cosmic ray source distributions offers a genuine new quality of modeling and is used to explain features in cosmic ray spectra at the Earth that are else-wise attributed to other propagation effects. We show that realistic cosmic ray propagation scenarios have to acknowledge non-axisymmetric source distributions.