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.     

Introducing the CTA concept

The Cherenkov Telescope Array (CTA) is a new observatory for very high-energy (VHE) gamma rays. CTA has ambitions science goals, for which it is necessary to achieve full-sky coverage, to improve the sensitivity by about an order of magnitude, to span about four decades of energy, from a few tens of GeV to above 100 TeV with enhanced angular and energy resolutions over existing VHE gamma-ray observatories. An international collaboration has formed with more than 1000 members from 27 countries in Europe, Asia, Africa and North and South America. In 2010 the CTA Consortium completed a Design Study and started a three-year Preparatory Phase which leads to production readiness of CTA in 2014. In this paper we introduce the science goals and the concept of CTA, and provide an overview of the project.     

Milagro limits and HAWC sensitivity for the rate-density of evaporating Primordial Black Holes

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ∼5.0 × 10¹⁴ g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV–TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.     

Studies of active galactic nuclei with CTA

In this paper, we review the prospects for studies of active galactic nuclei (AGN) using the envisioned future Cherenkov Telescope Array (CTA). This review focuses on jetted AGN, which constitute the vast majority of AGN detected at gamma-ray energies. Future progress will be driven by the planned lower energy threshold for very high energy (VHE) gamma-ray detections to ∼10 GeV and improved flux sensitivity compared to current-generation Cherenkov Telescope facilities. We argue that CTA will enable substantial progress on gamma-ray population studies by deepening existing surveys both through increased flux sensitivity and by improving the chances of detecting a larger number of low-frequency peaked blazars because of the lower energy threshold. More detailed studies of the VHE gamma-ray spectral shape and variability might furthermore yield insight into unsolved questions concerning jet formation and composition, the acceleration of particles within relativistic jets, and the microphysics of the radiation mechanisms leading to the observable high-energy emission. The broad energy range covered by CTA includes energies where gamma-rays are unaffected from absorption while propagating in the extragalactic background light (EBL), and extends to an energy regime where VHE spectra are strongly distorted. This will help to reduce systematic effects in the spectra from different instruments, leading to a more reliable EBL determination, and hence will make it possible to constrain blazar models up to the highest energies with less ambiguity.     

VAMOS: A pathfinder for the HAWC gamma-ray observatory

VAMOS¹ was a prototype detector built in 2011 at an altitude of 4100 m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.     

Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

Calculation of the Cherenkov light yield from electromagnetic cascades in ice with Geant4

In this work we investigate and parameterize the amount and angular distribution of Cherenkov photons which are generated by electro-magnetic cascades in water or ice. We simulate electromagnetic cascades with Geant4 for primary electrons, positrons and photons with energies ranging from 1 GeV to 10 TeV. We parameterize the total Cherenkov-light yield as a function of energy, the longitudinal evolution of the Cherenkov emission along the cascade-axis and the angular distribution of photons. Furthermore, we investigate the fluctuations of the total light yield, the fluctuations in azimuth and changes of the emission with increasing age of the cascade.     

The major upgrade of the MAGIC telescopes,Part II: A performance study using observations of the Crab Nebula

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ∼ 50  GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220 GeV is (0.66 ± 0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the σ of a 2-dimensional Gaussian distribution, at those energies is ≲ 0.07°, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11%–18% in flux normalization and ± 0.15 for the energy spectrum power-law slope.     

Stellar intensity interferometry: Prospects for sub-milliarcsecond optical imaging

Using kilometric arrays of air Cherenkov telescopes at short wavelengths, intensity interferometry may increase the spatial resolution achieved in optical astronomy by an order of magnitude, enabling images of rapidly rotating hot stars with structures in their circumstellar disks and winds, or mapping out patterns of nonradial pulsations across stellar surfaces. Intensity interferometry (once pioneered by Hanbury Brown and Twiss) connects telescopes only electronically, and is practically insensitive to atmospheric turbulence and optical imperfections, permitting observations over long baselines and through large airmasses, also at short optical wavelengths. The required large telescopes (～10 m) with very fast detectors (～ns) are becoming available as the arrays primarily erected to measure Cherenkov light emitted in air by particle cascades initiated by energetic gamma rays. Planned facilities (e.g., CTA, Cherenkov Telescope Array) envision many tens of telescopes distributed over a few square km. Digital signal handling enables very many baselines (from tens of meters to over a kilometer) to be simultaneously synthesized between many pairs of telescopes, while stars may be tracked across the sky with electronic time delays, in effect synthesizing an optical interferometer in software. Simulated observations indicate limiting magnitudes around m_V = 8, reaching angular resolutions ～30 μarcsec in the violet. The signal-to-noise ratio favors high-temperature sources and emission-line structures, and is independent of the optical passband, be it a single spectral line or the broad spectral continuum. Intensity interferometry directly provides the modulus (but not phase) of any spatial frequency component of the source image; for this reason a full image reconstruction requires phase retrieval techniques. This is feasible if sufficient coverage of the interferometric (u, v)-plane is available, as was verified through numerical simulations. Laboratory and field experiments are in progress; test telescopes have been erected, intensity interferometry has been achieved in the laboratory, and first full-scale tests of connecting large Cherenkov telescopes have been carried out. This paper reviews this interferometric method in view of the new possibilities offered by arrays of air Cherenkov telescopes, and outlines observational programs that should become realistic already in the rather near future.     

Thermospheric/mesospheric temperatures on Venus: Results from ground-based high-resolution spectroscopy of CO2 in 1990/1991 and comparison to results from 2009 and between other techniques

We report temperatures in Venus’ upper mesosphere/lower thermosphere, deduced from reanalyzing very high resolution infrared spectroscopy of CO₂ emission lines acquired in 1990 and 1991. Kinetic temperatures at ～110 km altitude (0.15 Pa) are derived from the Doppler width of fully-resolved single line profiles measured near 10.4 μm wavelength using the NASA GSFC Infrared Heterodyne Spectrometer (IRHS) at the NASA IRTF on Mauna Kea, HI, close to Venus inferior conjunction and two Venus solstices. Measured temperatures range from ～200 to 240 K with uncertainty typically less than 10 K. Temperatures retrieved from similar measurement in 2009 using the Cologne Tuneable Heterodyne Infrared Spectrometer (THIS) at the NOAO McMath Telescope at Kitt Peak, AZ are 10–20 K lower. Temperatures retrieved more recently from the SOIR instrument on Venus EXpress are consistent with these results when the geometry of observation is accounted for. It is difficult to compare ground-based sub-mm retrievals extrapolated to 110 km due to their much larger field of view, which includes the night side regions not accessible to infrared heterodyne observations. Temperature variability appears to be high on day-to-day as well as longer timescales. Observed short term and long term variability may be attributed to atmospheric dynamics, diurnal variability and changes over solar activity and seasons. The Venus International Reference Atmosphere (VIRA) model predicts cooler temperatures at the sampled altitudes in the lower thermosphere/upper mesosphere and is not consistent with these measurements.     

The Cherenkov Telescope Array potential for the study of young supernova remnants

Supernova remnants (SNRs) are among the most important targets for γ-ray observatories. Being prominent non-thermal sources, they are very likely responsible for the acceleration of the bulk of Galactic cosmic rays (CRs). To firmly establish the SNR paradigm for the origin of cosmic rays, it should be confirmed that protons are indeed accelerated in, and released from, SNRs with the appropriate flux and spectrum. This can be done by detailed theoretical models which account for microphysics of acceleration and various radiation processes of hadrons and leptons. The current generation of Cherenkov telescopes has insufficient sensitivity to constrain theoretical models. A new facility, the Cherenkov Telescope Array (CTA), will have superior capabilities and may finally resolve this long standing issue of high-energy astrophysics. We want to assess the capabilities of CTA to reveal the physics of various types of SNRs in the initial 2000 years of their evolution. During this time, the efficiency to accelerate cosmic rays is highest. We perform time-dependent simulations of the hydrodynamics, the magnetic fields, the cosmic-ray acceleration, and the non-thermal emission for type Ia, Ic and IIP SNRs. We calculate the CTA response to the γ-ray emission from these SNRs for various ages and distances, and we perform a realistic analysis of the simulated data. We derive distance limits for the detectability and resolvability of these SNR types at several ages. We test the ability of CTA to reconstruct their morphological and spectral parameters as a function of their distance. Finally, we estimate how well CTA data will constrain the theoretical models.     

Direct wind measurements from November 2007 in Venus’ upper atmosphere using ground-based heterodyne spectroscopy of CO2 at 10 μm wavelength

Between November 23 and 28, 2007, the Cologne Tuneable Heterodyne Infrared Spectrometer THIS was installed at the McMath-Pierce Solar Telescope (Kitt Peak, Arizona, USA) to determine zonal wind velocities and to estimate the subsolar-to-antisolar flow. We investigate dynamics in the upper atmosphere of Venus by measuring the Doppler shift of fully-resolved non-LTE CO₂ emission lines at 959.3917 cm^?1 (10.423 μm), which probe a narrow altitude region in Venus’ atmosphere around 110 ± 10 km (～1 μbar). The results show no significant zonal wind velocity at the equator. An increase with latitude up to 43 ± 13 m/s at a latitude of 33°N was observed. This confirms the deduction of a minor influence of Venus superrotation at an altitude of 110 km from previous measurements in May 2007 (Sornig et al., 2008). The specific observing geometry enables estimating the maximum cross terminator velocity of the subsolar-to-antisolar flow at 72 ± 47 m/s.     

Circulation of the Venus upper mesosphere/lower thermosphere: Doppler wind measurements from 2001–2009 inferior conjunction,sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed within the mesospheric to lower thermospheric altitude (70–120 km) region of the Venus atmosphere, have been mapped across the nightside disk of Venus during 2001–2009 inferior conjunctions, employing the James Clerk Maxwell Telescope (JCMT). Radiative transfer analysis of these thermal line absorptions supports temperature and CO mixing profile retrievals, as described in a companion paper (Clancy et al., 2012). Here, we consider the analysis of the sharp line absorption cores of these CO spectra in terms of accurate Doppler wind profile measurements at 95–115 km altitudes versus local time (～8 pm–4 am) and latitude (～60N–60S). These Doppler wind measurements support determinations of the nightside zonal and subsolar-to-antisolar (SSAS) circulation components over a variety of timescales. The average behavior fitted from 21 retrieved maps of ¹²CO Doppler winds (obtained over hourly, daily, weekly, and interannual intervals) indicates stronger average zonal (85 m/s retrograde) versus SSAS (65 m/s) circulation at the 1 μbar pressure (108–110 km altitude) level. However, the absolute and relative magnitudes of these circulation components exhibit extreme variability over daily to weekly timescales. Furthermore, the individual Doppler wind measurements within each nightside mapping observation generally show significant deviations (20–50 m/s, averaged over 5000 km horizontal scales) from the simple zonal/SSAS solution, with distinct local time and latitudinal characters that are also time variable. These large scale residual circulations contribute 30–70% of the observed nightside Doppler winds at any given time, and may be most responsible for global variations in nightside lower thermospheric trace composition and temperatures, as coincidentally retrieved CO abundance and temperature distributions do not correlate with solution retrograde zonal and SSAS winds (see companion paper, Clancy et al., 2012). Limited comparisons of these nightside submillimeter results with dayside infrared Doppler wind measurements suggest distinct dayside versus nightside circulations, in terms of zonal winds in particular. Combined ¹²CO and ¹³CO Doppler wind mapping observations obtained since 2004 indicate that the average zonal and SSAS wind components increase by 50–100% between altitudes of 100 and 115 km. If gravity waves originating from the cloud levels are responsible for the extension of zonal winds into the thermosphere (Alexander, M.J. [1992]. Geophys. Res. Lett. 19, 2207–2210), such waves deposit substantial momentum (i.e., break) in the lower nightside thermosphere.     

Milagro observations of potential TeV emitters

This paper reports the results from three targeted searches of Milagro TeV sky maps: two extragalactic point source lists and one pulsar source list. The first extragalactic candidate list consists of 709 candidates selected from the Fermi-LAT 2FGL catalog. The second extragalactic candidate list contains 31 candidates selected from the TeVCat source catalog that have been detected by imaging atmospheric Cherenkov telescopes (IACTs). In both extragalactic candidate lists Mkn 421 was the only source detected by Milagro. This paper presents the Milagro TeV flux for Mkn 421 and flux limits for the brighter Fermi-LAT extragalactic sources and for all TeVCat candidates. The pulsar list extends a previously published Milagro targeted search for Galactic sources. With the 32 new gamma-ray pulsars identified in 2FGL, the number of pulsars that are studied by both Fermi-LAT and Milagro is increased to 52. In this sample, we find that the probability of Milagro detecting a TeV emission coincident with a pulsar increases with the GeV flux observed by the Fermi-LAT in the energy range from 0.1 GeV to 100 GeV.     

The distribution of carbon monoxide in the lower atmosphere of Venus

We have obtained spatially resolved near-infrared spectroscopy of the Venus nightside on 15 nights over three observing seasons. We use the depth of the CO absorption band at 2.3 μm to map the two-dimensional distribution of CO across both hemispheres. Radiative transfer models are used to relate the measured CO band depth to the volume mixing ratio of CO. The results confirm previous investigations in showing a general trend of increased CO abundances at around 60° latitude north and south as compared with the equatorial regions. Observations taken over a few nights generally show very similar CO distributions, but significant changes are apparent over longer periods. In past studies it has been assumed that the CO latitudinal variation occurs near 35 km altitude, at which K-band sensitivity to CO is greatest. By modeling the detailed spectrum of the excess CO at high latitudes we show that it occurs at altitudes around 45 km, much higher than has previously been assumed, and that there cannot be significant contribution from levels of 36 km or lower. We suggest that this is most likely due to downwelling of CO-rich gas from the upper atmosphere at these latitudes, with the CO being removed by around 40 km through chemical processes such as the reaction with SO₃.     

The role of sputtering and radiolysis in the generation of Europa exosphere

The exosphere of an atmosphereless icy moon is the result of different surface release processes and subsequent modification of the released particles. At Europa icy moon, water molecules are directly released, but photolysis and radiolysis due to solar UV and Jupiter’s magnetospheric plasma, respectively, can result in OH, H, O and (possibly) H₂ production. These molecules can recombine to reform water and/or new chemical species. As a consequence, Europa’s neutral environment becomes a mixture of different molecules, among which, H₂O dominates in the highest altitudes and O₂, formed mainly by radiolysis of ice and subsequent release of the produced molecules, prevails at lower altitudes. In this work, starting from a previously developed Monte Carlo model for the generation of Europa’s exosphere, where the only considered species was water, we make a first attempt to simulate also the H₂ and O₂ components of the neutral environment around Europa, already observed by the Hubble Space Telescope and the Ultraviolet Imaging Spectrograph on board Cassini, during its flyby of Jupiter. Considering a specific configuration where the leading hemisphere coincides with the sunlit hemisphere, we estimate along the Europa–Sun line an O₂ column density of about 1.5 × 10¹⁹ m^?2 at the dayside and 3 × 10¹⁸ m^?2 at the nightside. In this work we also improve our previous estimation of the sputtered H₂O exosphere of this moon, taking into consideration the trailing–leading asymmetry in the magnetospheric ion bombardment and the energy and temperature dependences of the process yields. We find that a density of 1.5 × 10¹² H₂O/m³ is expected at altitudes ～0.1R_E above the surface of the trailing hemisphere. Additionally, we calculate the escape of H₂O, O₂ and H₂. The total number of neutral atoms in Europa’s neutral torus, is estimated to be in the range 7.8 × 10³²–3.3 × 10³³.