Wavelet imaging cleaning method for atmospheric Cherenkov telescopes

We present a new method of image cleaning for imaging atmospheric Cherenkov telescopes. The method is based on the utilization of wavelets to identify noise pixels in images of gamma-ray and hadronic induced air showers. This method selects more signal pixels with Cherenkov photons than traditional image processing techniques. In addition, the method is equally efficient at rejecting pixels with noise alone. The inclusion of more signal pixels in an image of an air shower allows for a more accurate reconstruction, especially at lower gamma-ray energies that produce low levels of light. We present the results of Monte Carlo simulations of gamma-ray and hadronic air showers which show improved angular resolution using this cleaning procedure. Data from the Whipple Observatory's 10-m telescope are utilized to show the efficacy of the method for extracting a gamma-ray signal from the background of hadronic generated images.     

Using the photons from the Crab Nebula seen by GLAST to calibrate MAGIC and the imaging air Cherenkov telescopes

In this article we discuss the possibility of using the observations by GLAST of standard gamma sources, as the Crab Nebula, to calibrate imaging air Cherenkov detectors, MAGIC in particular, and optimise their energy resolution. We show that at around 100 GeV the absolute energy calibration uncertainty of Cherenkov telescopes can be reduced to 10% by means of such cross-calibration procedure.     

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.     

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.     

Impact of aerosols and adverse atmospheric conditions on the data quality for spectral analysis of the H.E.S.S. telescopes

The Earth’s atmosphere is an integral part of the detector in ground-based imaging atmospheric Cherenkov telescope (IACT) experiments and has to be taken into account in the calibration. Atmospheric and hardware-related deviations from simulated conditions can result in the mis-reconstruction of primary particle energies and therefore of source spectra. During the eight years of observations with the High Energy Stereoscopic System (H.E.S.S.) in Namibia, the overall yield in Cherenkov photons has varied strongly with time due to gradual hardware aging, together with adjustments of the hardware components, and natural, as well as anthropogenic, variations of the atmospheric transparency. Here we present robust data selection criteria that minimize these effects over the full data set of the H.E.S.S. experiment and introduce the Cherenkov transparency coefficient as a new atmospheric monitoring quantity. The influence of atmospheric transparency, as quantified by this coefficient, on energy reconstruction and spectral parameters is examined and its correlation with the aerosol optical depth (AOD) of independent MISR satellite measurements and local measurements of atmospheric clarity is investigated.     

Multi-component study of extended sources with a likelihood method

The spectral and morphological analysis for gamma-ray sources with multiple emission components remains a major challenge for Cherenkov telescopes due to background emission from diffuse gamma rays. Current methods of background suppression, usually based on the bin-by-bin subtraction of OFF-source data do not allow an analysis of the various background components. As an alternative, we present an approach based on an event-by-event likelihood fit of ON-source data, using a combined spectral model for the source emission as well as for the gamma-like background obtained from fits of the OFF-source data. Multiple emission components are separated by successive fits in different energy regimes and spectral variations inside the extended source is derived. The performance of this approach is evaluated with toy Monte-Carlo studies. For the application to real data, two well-studied H.E.S.S. sources are re-examined: the extragalactic point-source PKS 2155-304 and the extended pulsar wind nebula HESS J1825-137. For the latter, radial variation of the emission spectral index was evaluated, confirming earlier findings by the H.E.S.S. collaboration.     

Upgrading and testing the 3D reconstruction of gamma-ray air showers as observed with an array of Imaging Atmospheric Cherenkov Telescopes

Stereoscopic arrays of Imaging Atmospheric Cherenkov Telescopes allow to reconstruct gamma-ray-induced showers in three dimensions, which offers several advantages: direct access to the shower parameters in space and straightforward calorimetric measurement of the incident energy. In addition, correlations between the different images of the same shower are taken into account. An analysis method based on a simple 3D-model of electromagnetic showers was recently implemented in the framework of the H.E.S.S. experiment. In the present article, the method is completed by an additional quality criterion, which reduces the background contamination by a factor of about 2 in the case of extended sources, while keeping gamma-ray efficiency at a high level. On the other hand, the dramatic flares of the blazar PKS 2155-304 in July 2006, which provided H.E.S.S. data with an almost pure gamma-ray sample, offered the unique opportunity of a precision test of the 3D-reconstruction method as well as of the H.E.S.S. simulations used in its calibration. An agreement at a few percent level is found between data and simulations for the distributions of all 3D shower parameters.     

A multivariate analysis approach for the imaging atmospheric Cherenkov telescopes system H.E.S.S.

A multivariate particle classification approach is proposed and applied to the analysis of the data from H.E.S.S. (High Energy Stereoscopic System). The combination of results from the three shower reconstruction methods: Hillas, model and 3D-model, leads to a substantial gain in the discrimination power between photons and hadrons. The construction and use of a combined effective estimator improves by several factors the signal-to-background ratio which is extremely important in case of studies of the faint and extended sources. The results of this approach are presented for a typical set of sources. The consequent gain in the sensitivity is shown through a comparison to the H.E.S.S. published results.     

Instrumentation for very high energy gamma-ray astronomy

Considerable progress has been made in the last half-decade in the field of very high energy (VHE) gamma-ray astronomy (photons with energies between 10¹¹ and 10¹³ eV). The high background level due to the isotropic cosmic ray flux which has bedevilled the field since its inception in the early 1960's can now be reduced to such a degree that significant gamma-ray signals from several sources become visible within a few hours of observation. The instrumentation and methodologies which have made this possible are reviewed. A brief historical introduction is followed by a summary of the salient properties of the atmospheric Cherenkov flash associated with VHE gamma-ray events. The major components of a VHE gamma-ray astronomy telescope are then reviewed. This is followed by a discussion of the different methodologies currently being used to discriminate against the cosmic ray background. Properties of several specific installations are then summarized, and possible future developments in VHE instrumentation are briefly discussed.     

Formation Flying for a Fresnel Lens Observatory Mission

The employment of a large area Phase Fresnel Lens (PFL) in a gamma-ray telescope offers the potential to image astrophysical phenomena with micro-arcsecond (μ^′′) angular resolution [1]. In order to assess the feasibility of this concept, two detailed studies have been conducted of formation flying missions in which a Fresnel lens capable of focussing gamma-rays and the associated detector are carried on two spacecraft separated by up to 10⁶ km. These studies were performed at the NASA Goddard Space Flight Center Integrated Mission Design Center (IMDC) which developed spacecraft, orbital dynamics, and mission profiles. The results of the studies indicated that the missions are challenging but could be accomplished with technologies available currently or in the near term. The findings of the original studies have been updated taking account of recent advances in ion thruster propulsion technology.     

Monte Carlo simulations of radio pulses in atmospheric showers using ZHAireS

We present predictions for the radio pulses emitted by extensive air showers using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full complexity of ultra-high energy cosmic-ray induced shower development in the atmosphere, and allows the calculation of the electric field in both the time and frequency domains. We do not presuppose any emission mechanism and our results are compatible with a superposition of geomagnetic and charge excess radio emission effects. We investigate the polarization of the electric field as well as the effects of the refractive index n and shower geometry on the radio pulses. We show that geometry, coupled to the relativistic effects that appear when using a realistic refractive index n > 1, play a prominent role on the radio emission of air showers.     

Prospects for far infrared arrays

We report initial performance measurements of a 1/8 scale version of a 32×32 pixel array under development for SIRTF. This array demonstrates that we can reach the sensitivity limits set by the natural backgrounds in space while providing good imaging and photometric performance. Based on the achieved performance levels, we project the imaging capabilities of SIRTF in the far infrared to exceed by a factor of more than 10,000 those achieved by any preceding telescope.     

AMICA, an astro-mapper for AMS

The alpha magnetic spectrograph (AMS) is a composite particle detector to be accommodated on the International Space Station (ISS). AMS is mainly devoted to galactic, charged cosmic rays studies, antimatter and dark matter searches. Besides the main, classical physics goals, capabilities in the field of GeV and multi-GeV gamma astrophysics have been established and are under investigation by a number of groups. Due to the unsteadiness of the ISS platform, a star-mapper device is required in order to fully exploit the intrinsic arc-min angular resolution provided by the silicon tracker. A star-mapper is conceptually an imaging, optical instrument able to autonomously recognize a stellar field and to calculate its own orientation with respect to an inertial reference frame. AMICA (Astro Mapper for Instruments Check of Attitude) on AMS is responsible for providing real-time information that is going to be used off-line for compensating the large uncertainties in the ISS flight attitude and the structural degrees of freedom. In this paper, we describe in detail the AMICA sub-system, the accommodation/integration issues and the in-flight alignment procedure adopting identified galactic (Pulsars) and extra-galactic (AGNs) sources.     

A Si/CdTe Compton Camera for gamma-ray lens experiment

Cadmium telluride (CdTe) and cadmium zinc telluride (CdZnTe) have been regarded as promising semiconductor materials for hard X-ray and γ-ray detection. However, a considerable amount of charge loss in these detectors results in a reduced energy resolution. We have achieved a significant improvement in the spectral properties by forming the Schottky junction on the Te side of the CdTe wafer. With the further reduction of leakage current by an adoption of guard ring structure, we have demonstrated a CdTe pixel detector with high energy resolution and full charge collection capabilty. The detector has a pixel size of a few mm and a thickness of 0.5 $-$ 1 mm. We apply this high resolution detector to a new silicon and CdTe Compton Camera which features high angular resolution. We also describe a concept of the stack detector which consists of many thin CdTe layers and provides sufficient efficiency for hard X-rays and gamma-rays up to several hundred keV maintaining good energy resolution. A narrow-FOV Compton telescope can be realized by installing a Si/CdTe Compton Camera inside the deep well of an active shield. This configuration is very suitable as focal plane detector for future focusing gamma-ray missions.     

A new data structure for accelerating kinetic Monte Carlo method

The kinetic Monte Carlo simulation is a rigorous numerical approach to study the chemistry on dust grains in cold dense interstellar clouds. By tracking every single reaction in chemical networks step by step, this approach produces more precise results than other approaches but takes too much computing time. Here we present a method of a new data structure, which is applicable to any physical conditions and chemical networks, to save computing time for the Monte Carlo algorithm. Using the improved structure,the calculating time is reduced by 80 percent compared with the linear structure when applied to the osu-2008 chemical network at 10K. We investigate the effect of the encounter desorption in cold cores using the kinetic Monte Carlo model with an accelerating data structure. We found that the encounter desorption remarkably decreases the abundance of grain-surface H_2 but slightly influences the abundances of other species on the grain.     

A Positioning System based on Communication Satellites and the Chinese Area Positioning System （CAPS）

The Chinese Area Positioning System (CAPS) is a positioning system based on satellite communication that is fundamentally different from the 3"G" (GPS, GLONASS and GALILEO) systems. The latter use special-purpose navigation satellites to broadcast navi-gation information generated on-board to users, while the CAPS transfers ground-generated navigation information to users via the communication satellite. In order to achieve accurate Positioning, Velocity and Time (PVT), the CAPS employs the following strategies to over-come the three main obstacles caused by using the communication satellite: (a) by real-time following-up frequency stabilization to achieve stable frequency; (b) by using a single carrier in the transponder with 36 MHz band-width to gain sufficient power; (c) by incorporating Decommissioned Geostationary Orbit communication satellite (DGEO), barometric pressure and Inclined Geostationary Orbit communication satellite (IGSO) to achieve the 3-D posi-tioning. Furthermore, the abundant transponders available on DGEO can be used to realize the large capacity of communication as well as the integrated navigation and communication. With the communication functions incorporated, five new functions appear in the CAPS: (1) combination of navigation and communication; (2) combination of navigation and high accu-racy orbit measurement; (3) combination of navigation message and wide/local area differen-tial processing; (4) combination of the switching of satellites, frequencies and codes; and (5) combination of the navigation message and the barometric altimetry. The CAPS is thereby labelled a PVT5C system of high accuracy. In order to validate the working principle and the performance of the CAPS, a trial system was established in the course of two years at a cost of about 20 million dollars. The trial constellation consists of two GEO satellites located at E87.5°and E110.5°, two DGEOs located at E130° and E142°, as well as barometric altimetry as a virtual satellite. Static and dynamic performance tests were completed for the Eastern, the Western, the Northern, the Southern and the Middle regions of China. The evaluation results are as follows: (1) land static test, plane accuracy range: C/A code, 15～25 m; P code, 5～10 meters; altitude accuracy range, 1～3 m; (2) land dynamic test, plane accuracy range, C/A code, 15～25 m; P code, 8～10m; (3) velocity accuracy, C/A code, 0.13～0.3 ms-1, P code, 0.15～0.17m s-1; (4) timing accuracy, C/A code, 160ns, P code, 13 ns; (5) timing compared accuracy of Two Way Satellite Time and Frequency Transfer (TWSTFT), average accuracy, 0.068 ns; (6) random error of the satellite ranging, 10.7 mm; (7) orbit determination accuracy, better than 2 m. The above stated random error is 1σ error. At present, this system is used as a preliminary operational system and a complete system with 3 GEO, 3 DGEO and 3 IGSO is being established.     

A Monte Carlo model for the gardening of the lunar regolith

The processes of movement and turnover of the lunar regolith are described by a Monte Carlo model, which includes the effects of collisions by particles from 10^?7 g to 10¹⁰ g. The movement of material by the direct cratering process is the dominant mode, but slumping is also included for angles exceeding the static angle of repose. Using a group of interrelated computer programs a large number of properties are calculated, including topography, formation of layers, depth of the disturbed layer, nuclear track distributions, cosmogenic nuclides and others. In the most complex program, the history of a 36 point square array is followed for times up to 4 × 10⁸ yr. As expected the histories generated are complex and exhibit great variety. Because a crater covers much less area than its ejecta blanket, there is a tendency for the height change at a test point to exhibit the ‘gambler's ruin’ phenomenon: periods of slow accumulation followed by sudden excavation. In general the agreement with experiment and observation seems good. Two areas of disagreement stand out. First, the calculated surface is rougher than that observed. This problem is understood, and will not occur in a newer version of the model. Second, the observed bombardment ages, of the order of 4 × 10⁸ yr, are shorter than expected (by perhaps a factor of 5). We cannot accept Fireman's (1974) explanation; this remains an important puzzle.