GeV emission from gamma-ray burst afterglows

We calculate the GeV afterglow emission expected from a few mechanisms related to gamma-ray bursts (GRBs) and their afterglows. Given the brightness of the early X-ray afterglow emission measured by Swift /X-Ray Telescope, Gamma-ray Large Area Space Telescope (GLAST)/Large Area Telescope (LAT) should detect the self-Compton emission from the forward shock driven by the GRB ejecta into the circumburst medium. Novel features discovered by Swift in X-ray afterglows (plateaus and chromatic light-curve breaks) indicate the existence of a pair-enriched, relativistic outflow located behind the forward shock. Bulk and inverse-Compton upscattering of the prompt GRB emission by such outflows provide another source of GeV afterglow emission detectable by LAT. The large-angle burst emission and synchrotron forward-shock emission are, most likely, too dim at high photon energy to be observed by LAT. The spectral slope of the high-energy afterglow emission and its decay rate (if it can be measured) allow the identification of the mechanism producing the GeV transient emission following GRBs.     

Comparison of Fermi-LAT and CTA in the region between 10–100 GeV

The past decade has seen a dramatic improvement in the quality of data available at both high (HE: 100 MeV to 100 GeV) and very high (VHE: 100 GeV to 100 TeV) gamma-ray energies. With three years of data from the Fermi Large Area Telescope (LAT) and deep pointed observations with arrays of Cherenkov telescope, continuous spectral coverage from 100 MeV to ∼10 TeV exists for the first time for the brightest gamma-ray sources. The Fermi-LAT is likely to continue for several years, resulting in significant improvements in high energy sensitivity. On the same timescale, the Cherenkov Telescope Array (CTA) will be constructed providing unprecedented VHE capabilities. The optimisation of CTA must take into account competition and complementarity with Fermi, in particularly in the overlapping energy range 10–100 GeV. Here we compare the performance of Fermi-LAT and the current baseline CTA design for steady and transient, point-like and extended sources.     

The comparative spectra of cosmic-ray protons and helium nuclei

We have re-examined and extended the measurements of the primary cosmic ray proton and helium nuclei intensities in the range from a few MeV nuc^–1 to 100 GeV nuc^–1 using a considerable body of recently published data. The differential spectra obtained from this data are determined as a function of both energy and rigidity. The exponents of the energy spectra of both protons and helium nuclei are found to be different at the same energy/nucleon and to increase with increasing energy between 1 and 100 GeV nuc^–1 reaching a value=–2.70 at higher energies and in addition, theP/He ratio changes from a value 5 at 1 GeV nuc^–1 and below to a value 30 at 100 GeV nuc^–1. On a rigidity representation the spectral exponent for each species is nearly identical and remains virtually constant above several GV at a value of –2.70, and in addition, theP/He ratio is also a constant 7 above 3 GeV. The changingP/He ratio and spectral exponent on an energy representation occur at energies well above those at which interplanetary modulation effects or interstellar ionization energy loss effects can significantly affect the spectra. In effect by comparing energy spectra and rigidity spectra in the intermediate energy range above the point where solar modulation effects and interstellar energy loss effects are important, but in the range where there are significant differences between energy and rigidity spectra, we deduce that the cosmic ray source spectra are effectively rigidity spectra. This fact has important implications regarding the mechanism of acceleration of this radiation and also with regard to the form of the assumed galactic spectrum at low energies. The relationship between the proton and helium spectra derived here and the heavier nuclei spectral differences recently reported in the literature is also examined.If rigidity spectra are adopted for protons and helium nuclei, then the source abundance ratio of these two components is determined to be 7:1. Some cosmological implications of this ratio are discussed.     

IACT observations of gamma-ray bursts: prospects for the Cherenkov Telescope Array

Gamma rays at rest frame energies as high as 90 GeV have been reported from gamma-ray bursts (GRBs) by the Fermi Large Area Telescope (LAT). There is considerable hope that a confirmed GRB detection will be possible with the upcoming Cherenkov Telescope Array (CTA), which will have a larger effective area and better low-energy sensitivity than current-generation imaging atmospheric Cherenkov telescopes (IACTs). To estimate the likelihood of such a detection, we have developed a phenomenological model for GRB emission between 1 GeV and 1 TeV that is motivated by the high-energy GRB detections of Fermi-LAT, and allows us to extrapolate the statistics of GRBs seen by lower energy instruments such as the Swift-BAT and BATSE on the Compton Gamma-ray Observatory. We show a number of statistics for detected GRBs, and describe how the detectability of GRBs with CTA could vary based on a number of parameters, such as the typical observation delay between the burst onset and the start of ground observations. We also consider the possibility of using GBM on Fermi as a finder of GRBs for rapid ground follow-up. While the uncertainty of GBM localization is problematic, the small field-of-view for IACTs can potentially be overcome by scanning over the GBM error region. Overall, our results indicate that CTA should be able to detect one GRB every 20–30 months with our baseline instrument model, assuming consistently rapid pursuit of GRB alerts, and provided that spectral breaks below ~100 GeV are not a common feature of the bright GRB population. With a more optimistic instrument model, the detection rate can be as high as 1 to 2 GRBs per year.     

Pulsar simulations for the Fermi Large Area Telescope

Pulsars are among the prime targets for the Large Area Telescope (LAT) aboard the recently launched Fermi observatory. The LAT will study the gamma-ray Universe between 20 MeV and 300 GeV with unprecedented detail. Increasing numbers of gamma-ray pulsars are being firmly identified, yet their emission mechanisms are far from being understood. To better investigate and exploit the LAT capabilities for pulsar science, a set of new detailed pulsar simulation tools have been developed within the LAT collaboration. The structure of the pulsar simulator package (PulsarSpectrum) is presented here. Starting from photon distributions in energy and phase obtained from theoretical calculations or phenomenological considerations, gamma-rays are generated and their arrival times at the spacecraft are determined by taking into account effects such as barycentric effects and timing noise. Pulsars in binary systems also can be simulated given orbital parameters. We present how simulations can be used for generating a realistic set of gamma-rays as observed by the LAT, focusing on some case studies that show the performance of the LAT for pulsar observations.     

Influence of the geomagnetic field on the IACT detection technique for possible sites of CTA observatories

We investigate the influence of the geomagnetic field (GF) on the Imaging Air Cherenkov Telescope technique for two northern (Tenerife and San Pedro Martir) and three southern (Salta, Leoncito and Namibia (the H.E.S.S.-site)) site candidates for Cherenkov Telescope Array (CTA) observatories. We use the CORSIKA and sim_telarray programs for Monte Carlo simulations of gamma ray showers, hadronic background and the telescope response. We focus here on gamma ray measurements in the low energy, sub-100 GeV, range. Therefore, we only consider the performance of arrays of several large telescopes. Neglecting the GF effect, we find (in agreement with previous studies) that such arrays have lower energy thresholds, and larger collection areas below 30 GeV, when located at higher altitudes. We point out, however, that in the considered ranges of altitudes and magnetic field intensities, 1800–3600 m a.s.l. and 0–40 μT, respectively, the GF effect has a similar magnitude to this altitude effect. We provide the trigger-level performance parameters of the observatory affected by the GF effect, in particular the collection areas, detection rates and the energy thresholds for all five locations, which information may be useful in the selection of sites for CTA. We also find simple scaling of these parameters with the magnetic field strength, which can be used to assess the magnitude of the GF effect for other sites; in this work we use them to estimate the performance parameters for five sites: South Africa-Beaufort West, USA-Yavapai Ranch, Namibia-Calapanzi, Chile-La Silla and India-Hanle. We roughly investigate the impact of the geophysical conditions on gamma/hadron separation procedures involving image shape and direction cuts. We note that the change of altitude has an opposite effect at the trigger and analysis levels, i.e. gains in triggering efficiency at higher altitudes are partially balanced by losses in the separation efficiency. In turn, a stronger GF spoils both the shape and the direction discrimination of gamma rays, thus its effects at the trigger and analysis levels add up resulting in a significant reduction of the observatory performance. Overall, our results indicate that the local GF strength at a site can be equally important as its altitude for the low-energy performance of CTA.     

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.     

High energy γ-radiation from the core of radio galaxy Centaurus A

The results of analysis of approximately 3 years of gamma-ray observations (August 2008-July 2011) of the radio galaxy Centaurus A with the Fermi Large Area Telescope (Fermi LAT) are presented. By modeling the surrounding (background) sources including the giant lobes of Centaurus A, and using the standard binned likelihood analysis method, the energy spectrum of the core is derived. In the energy range below several GeV it is described by a single power-law with photon index Γ = 2.73 ± 0.06 in agreement with the report of the Fermi LAT collaboration based on the first 10 months observations of the source. However, at higher energies the new data show significant excess above the extrapolation of the energy spectrum from low energies. The total flux between 200 MeV to 50 GeV is estimated to be (1.63 ± 0.14) × 10⁻⁷ ph cm⁻² s⁻¹. The comparison of the corresponding Spectral Energy Distribution (SED) at GeV energies with the SED in the TeV energy band reported by the H.E.S.S. collaboration shows that we deal with two or perhaps even three components of gamma-radiation originating from different regions located within the central 10 kpc of Centaurus A.     

Performance of the MAGIC stereo system obtained with Crab Nebula data

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50 − 60 GeV. Advanced stereo analysis techniques allow MAGIC to achieve a sensitivity as good as (0.76 ± 0.03)% of the Crab Nebula flux in 50 h of observations above 290 GeV. The angular resolution at those energies is better than ∼0.07°. We also perform a detailed study of possible systematic effects which may influence the analysis of the data taken with the MAGIC telescopes.     

Implications of the detection of pulsed TeV gamma-rays from Geminga

The implications of the recently reported 237 ms pulsed TeV gamma-ray signal from Geminga with respect to system energetics and gamma-ray emission mechanisms are discussed. Using the reported flux values in the MeV-GeV range from the EGRET and COMPTEL experiments and the published results on TeV emission, the derived emission spectrum of Geminga in the MeV-TeV energy range is shown to suggest a spectral break at 2 GeV, with the differential exponent changing from a value of –1.50 upto 2 GeV to a value of –2.7 in the TeV range.     

Fine structure in the cosmic ray spectrum: Further analysis and the next step

An analysis is made of the fine structure in the cosmic ray energy spectrum: new facets of present observations and their interpretation and the next step. It is argued that less than about 10% of the intensity of the helium ‘peak’ at the knee at ≈5 PeV is due to just a few sources (SNR) other than the single source. The apparent concavity in the rigidity spectra of protons and helium nuclei which have maximum curvature at about 200 GV is confirmed by a joint analysis of the PAMELA, CREAM and ATIC experiments. The spectra of heavier nuclei also show remarkable structure in the form of ‘ankles’ at several hundred GeV/nucleon. Possible mechanisms are discussed. The search for ‘pulsar peaks’ has not yet proved successful.     

甚大面积伽马射线空间望远镜计划

高能伽马射线探测是研究极端天体物理的主要途径之一.空间高能伽马射线探测具有覆盖波段宽、时间连续性好、能量分辨率高等突出优势.在成功研发并运行我国首颗天文卫星—“悟空”号(DArk Matter Particle Explorer, DAMPE)的基础上,紫金山天文台联合国内的多家单位提议研制甚大面积伽马射线空间望远镜(Very Large Area gamma-ray Space Telescope, VLAST),该望远镜在GeV–TeV能段接受度高达10 m2·sr,并具有强的MeV–GeV波段探测能力,其综合性能预期比费米卫星的大面积伽马望远镜(Fermi-LAT (Large Area Telescope))提升10倍之上.重点介绍了VLAST的主要科学目标,探测器的初步配置及预期性能指标.     

A new component of cosmic rays?

Recent direct measurements of the energy spectra of the major mass components of cosmic rays have indicated the presence of a ‘kink’ in the region of 200 GeV per nucleon. The kink, which varies in magnitude from one element to another, is much sharper than predicted by our cosmic ray origin model in which supernova remnants are responsible for cosmic ray acceleration and it appears as though a new, steeper component is responsible.The component amounts to about 20 percent of the total at 30 GeV/nucleon for protons and helium nuclei and its magnitude varies with nuclear charge; the unweighted fraction for all cosmic rays being 36%.The origin of the new component is subject to doubt but the contenders include O, B, A, supergiant and Wolf-Rayet stars, by way of their intense stellar winds. Another explanation is also in terms of these particles as the sources but then being trapped, and even further accelerated, in the Local Bubble.     

Multiband fitting to three long GRBs with Fermi/LAT data:structured ejecta sweeping up a density-jump medium

We present broadband (radio,optical,X-ray and GeV) fits to the afterglow light curves and spectra of three long-duration gamma-ray bursts (GRBs 080916C,090902B,and 090926A) detected by the Gamma-Ray Burst Monitor and Large Area Telescope (LAT) instruments on the Fermi satellite.Using the observed broadband data,we study the origin of the high energy emission,and suggest that the early-time GeV emission and the late-time radio,optical,and X-ray afterglows can be understood as being due to synchrotron emissio...     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

Study of X-ray and γ-ray flare from Mkn-421 during February 2010 using Swift and Fermi observations

In this work, we report on the intense flaring activity from Mkn-421 in X-ray and γ-ray regimes simultaneously observed by Swift-XRT/BAT and Fermi-LAT satellite telescopes in February 2010. With the aim of understanding the underlying physics of the flaring state in Mkn-421, we have performed a detailed spectral analysis of Swift/XRT and Fermi/LAT observations of Mkn-421 during February 12–25, 2010 (MJD 55239–55252). Over this period, we study the daily light curves and spectral variability of the source in 1–10 keV, 0.1–1 GeV and 1–100 GeV energy bands. We have performed the spectral analysis of Swift-XRT and Fermi/LAT observations to study the spectral evolution in the X-ray and gamma-ray energy domains respectively. We also compute the fractional variability amplitude in both the energy bands during the above period. We study trends between spectral parameters and physical insights provided by the parameter responsible for X-ray and γ-ray emission from the source. We search for energetic features phenomenologically linked to the single zone SSC model for blazar emission. We also produce the broad band SED with a leptonic single zone SSC model for the source.     

A tool to estimate the Fermi Large Area Telescope background for short-duration observations

The proper estimation of the background is a crucial component of data analyses in astrophysics, such as source detection, temporal studies, spectroscopy, and localization. For the case of the Large Area Telescope (LAT) on board the Fermi spacecraft, approaches to estimate the background for short (≲1000 s duration) observations fail if they ignore the strong dependence of the LAT background on the continuously changing observational conditions. We present a (to be) publicly available background-estimation tool created and used by the LAT Collaboration in several analyses of Gamma Ray Bursts. This tool can accurately estimate the expected LAT background for any observational conditions, including, for example, observations with rapid variations of the Fermi spacecraft’s orientation occurring during automatic repointings.     

Keck and VLT AO observations and models of the uranian rings during the 2007 ring plane crossings

We present observations of the uranian ring system at a wavelength of 2.2 μm, taken between 2003 and 2008 with NIRC2 on the W.M. Keck telescope in Hawaii, and on 15–17 August 2007 with NaCo on the Very Large Telescope (VLT) in Chile. Of particular interest are the data taken around the time of the uranian ring plane crossing with Earth on 16 August 2007, and with the Sun (equinox) on 7 December 2007. We model the data at the different viewing aspects with a Monte Carlo model to determine: (1) the normal optical depth τ₀, the location, and the radial extent of the main rings, and (2) the parameter Aτ₀ (A is the particle geometric albedo), the location, and the radial plus vertical extent of the dusty rings. Our main conclusions are: (i) The brightness of the ? ring is significantly enhanced at small phase and ring inclination angles; we suggest this extreme opposition effect to probably be dominated by a reduction in interparticle shadowing. (ii) A broad sheet of dust particles extends inwards from the λ ring almost to the planet itself. This dust sheet has a vertical extent of ∼140 km, and Aτ₀ = 2.2 × 10⁻⁶. (iii) The dusty rings between ring 4 and the α ring and between the α and β rings are vertically extended with a thickness of ∼300 km. (iv) The ζ ring extends from ∼41,350 km almost all the way inwards to the planet. The main ζ ring, centered at ∼39,500 km from the planet, is characterized by Aτ₀ = 3.7 × 10⁻⁶; this parameter decreases closer to the planet. The ζ ring has a full vertical extent of order 800–900 km, with a pronounced density enhancement in the mid-plane. (v) The η_c ring is optically thin and less than several tens of km in the vertical direction. This ring may be composed of macroscopic material, surrounded by clumps of dust.     

Determination of physical properties of the Asteroid (41) Daphne from interferometric observations in the thermal infrared

We describe interferometric observations of the Asteroid (41) Daphne in the thermal infrared obtained with the Mid-Infrared Interferometric Instrument (MIDI) and the Auxiliary Telescopes (ATs) of the European Southern Observatory (ESO) Very Large Telescope Interferometer (VLTI). We derived the size and the surface thermal properties of (41) Daphne by means of a thermophysical model (TPM), which is used for the interpretation of interferometric data for the first time. From our TPM analysis, we derived a volume equivalent diameter for (41) Daphne of 189 km, using a non-convex 3-D shape model derived from optical lightcurves and adaptive optics images (B. Carry, private communication). On the other hand, when using the convex shape of Kaasalainen et al. (Kaasalainen, M., Mottola, S., Fulchignoni, M. [2002]. Icarus 159, 369-395) in our TPM analysis, the resulting volume equivalent diameter of (41) Daphne is between 194 and 209 km, depending on the surface roughness. The shape of the asteroid is used as an a priori information in our TPM analysis. No attempt is made to adjust the shape to the data. Only the size of the asteroid and its thermal parameters such as, albedo, thermal inertia and roughness are adjusted to the data. We estimated our model systematic uncertainty to be of 4% and of 7% on the determination of the asteroid volume equivalent diameter depending on whether the non-convex or the convex shape is used, respectively. In terms of thermal properties, we derived a value of the surface thermal inertia smaller than 50 J m⁻² s^−0.5 K⁻¹ and preferably in the range between 0 and ∼30 J m⁻² s^−0.5 K⁻¹. Our TPM analysis also shows that Daphne has a moderate macroscopic surface roughness.     

Ground-based calibration and characterization of the Fermi gamma-ray burst monitor detectors

One of the scientific objectives of NASA’s Fermi Gamma-ray Space Telescope is the study of Gamma-Ray Bursts (GRBs). The Fermi Gamma-Ray Burst Monitor (GBM) was designed to detect and localize bursts for the Fermi mission. By means of an array of 12 NaI(Tl) (8 keV to 1 MeV) and two BGO (0.2 to 40 MeV) scintillation detectors, GBM extends the energy range (20 MeV to > 300 GeV) of Fermi’s main instrument, the Large Area Telescope, into the traditional range of current GRB databases. The physical detector response of the GBM instrument to GRBs is determined with the help of Monte Carlo simulations, which are supported and verified by on-ground individual detector calibration measurements. We present the principal instrument properties, which have been determined as a function of energy and angle, including the channel-energy relation, the energy resolution, the effective area and the spatial homogeneity.