Stereoscopic observations of gamma rays at the Whipple observatory

Observations of photons at E ≥ 550 GeV from the Crab Nebula are presented and used to assess the potential of multi-telescope systems for γ-ray astronomy.

The Whipple observatory 10 m and 8 m imaging atmospheric erenkov telescopes have been used to provide a stereoscopic view of air showers to make a more complete measurement of air shower parameters. Here we present a measurement of the spread in the arrival direction of primary γ-rays originating from a point source. The data show that the shower arrival direction can be reconstructed with an accuracy of σ = 0.°14.     

The potential of the ground based arrays of imaging atmospheric Cherenkov telescopes. II. Gamma ray flux sensitivities

We discuss the capability of ‘100 GeV’ class imaging atmospheric Cherenkov telescope (IACT) arrays as future powerful instruments of ground-based gamma-ray astronomy. It is assumed that the array is gathered from individually triggered quadrangular 4-IACT ‘cells’ with a linear size of about 100 m. The multi-cell concept allows coverage of large detection areas economically, and at the same time the effective exploitation of the stereoscopic approach of determination of the shower parameters using information obtained by several IACTs simultaneously. Determination of arrival directions of γ-ray primaries on an event-by-event basis with accuracy δθ ≤ 0.1° combined with high suppression efficiency (at both the hardware and software levels) of the background hadronic showers by a factor of ≈ 10³, and large, up to 1 km² collection areas, can provide minimum detectable energy fluxes of ≥ 100 GeV γ-rays from point sources down to 10⁻¹³ erg/cm² s which is about 3 orders of magnitude lower than the current sensitivities achieved by the satellite-borne detectors at MeV and GeV energies. High sensitivities of multi-IACT arrays would partially compensate the limited efficiency of the technique for all-sky surveys, as well as allow study of moderately extended (≤ 1°) γ-ray sources. IACT arrays with minimum detectable fluence of ≥ 100 GeV γ-rays S_γ < 10⁻⁸ erg/cm² are well suited for effective exploration of highly sporadic nonthermal phenomena from different classes of astrophysical objects on time-scales from ≤ 1 s to several minutes.     

Limits on the TeV flux of diffuse gamma rays as measured with the HEGRA air shower array

Using data from the HEGRA air shower array, taken in the period from April 1998 to March 2000, upper limits on the ratio I_γ/I_CR of the diffuse photon flux I_γ to the hadronic cosmic ray flux I_CR are determined for the energy region 20–100 TeV. The analysis uses a gamma–hadron discrimination which is based on differences in the development of photon- and hadron-induced air showers after the shower maximum. A method which is sensitive only to the non-isotropic component of the diffuse photon flux yields an upper limit of I_γ/I_CR (at 54 TeV) <2.0×10⁻³ (at the 90% confidence level) for a sky region near the inner galaxy (20°< galactic longitude <60° and |galactic latitude |<5°). A method which is sensitive to both the isotropic and the non-isotropic component yields global upper limits of I_γ/I_CR (at 31 TeV) <1.2×10⁻² and I_γ/I_CR (at 53 TeV) <1.4×10⁻² (at the 90% confidence level).     

Methods to determine the angular resolution of the HEGRA extended air shower scintillator array

A precise knowledge of the angular resolution of scintillator arrays used to observe extended air showers (EAS) is of key importance in the search for VHE/UHE γ point sources. Four independent methods have been used to determine the mean resolution for which a value of ΔΘ₆₃ of 0.8°(1.0°) at a proton threshold of 50 (40) TeV has been found for the HEGRA EAS-array.     

Improvement of gamma-hadron discrimination at TeV energies using a new parameter, image Surface Brightness

The recent advances in TeV γ-ray astronomy using imaging Cerenkov telescopes is mainly due to the good separation between the properties of γ-ray and hadronic showers. A new method for this discrimination is described. The method is mainly based on the change of the image's Size and Surface Brightness (Size/Area) with energy and the differences in these parameters between those initiated by γ-rays and hadrons. The application of the method to the Whipple's 10 m telescope database from the Crab Nebula and Markarian 421 demonstrates its sensitivity to gamma-hadron separation. This study shows the importance of the image Surface Brightness, in addition to the previously used Shape and Orientation parameters.     

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.     

Constraining the power spectrum using clusters

We analyze an extended redshift sample of Abell/ACO clusters and compare the results with those coming from numerical simulations of the cluster distribution, based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, so that we generate a set of 48 independent mock Abell/ACO cluster samples per model, on which we estimate cosmic variance effects. Other than the standard CDM model, we consider (a) Ω₀ = 1 CDM models based on lowering the Hubble parameter and/or on tilting the primordial spectrum; (b) Ω₀ = 1 Cold + Hot DM models with 0.1 ≤Ω_ν ≤0.5; (c) low-density flat ΛCDM models with 0.3 ≤Ω₀ ≤0.5. We compare real and simulated cluster distributions by analysing correlation statistics, the probability density function, and supercluster properties from percolation analysis. We introduce a generalized definition of the spectrum shape parameter Γ in terms of σ₂₅/σ₈, where σ_ris the rms fluctuation amplitude within a sphere of radius r. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 Γ 0.25 and an effective spectral index at the 20 h⁻¹ Mpc scale −1.1 n_eff −0.9 are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press-Schechter approach, whose reliability is explicitly tested against N-body simulations. By combining results from the analysis of the distribution and the abundance of clusters we conclude that, of the cosmological models considered here, the only viable models are either Cold + Hot DM ones with 0.2 Ω_ν 0.3, better if shared between two massive ν species, and ΛCDM ones with 0.3 Ω₀0.5.     

Gamma-ray bursters as sources of cosmic rays

From the little we know of the physical conditions in γ-ray bursters, it seems that they are potentially effective in the acceleration of high-energy cosmic rays (CRs), especially if the bursters are at cosmological distances. We find that, with the observed statistics and fluxes of γ-ray bursts, cosmological bursters may be an important source of cosmic rays in two regions of the observed spectrum: (1) At the very-high-energy end (E > 10¹⁹ eV), where CRs must be of extragalactic origin. (2) Around and above the spectral feature that has been described as a bump and/or a knee, which occurs around 10¹⁵ eV. The occasional bursters that occur inside the Galaxy — about once in a few hundred thousand years if burst emission is isotropic; more often, if it is beamed — could maintain the density of galactic cosmic rays at the observed level in this range. These two energy ranges might correspond to two typical CR energy scales characteristic of bursters: one pertinent to CR acceleration due to interaction of a magnetized-fireball front with an ambient medium; the other to acceleration in the fireball itself (e.g. shock acceleration).     

Discussion of trignometric parallaxes determined along right ascension and declination

We discuss the data of trignometric parallaxes determined along and δ published by the US Naval Observatory and Van Vleck Observatory. We found that the absolute difference δπ = π_x−π_y shows a very similar periodic variation with for both observations. This variation is shown to be due to errors in π_y and to be related to the current practice of maximising the parallax factor for parallax determination in right ascension. Therefore, if possible, we should increase observations where the parallax factor in declination is maximised.     

Calculation of the TeV prompt muon component in very high energy cosmic ray showers

HEMAS-DPM is a Monte Carlo for the simulation of very high energy cosmic ray showers, which includes the DPMJET-II code based on the two component Dual Parton Model. DPMJET-II provides also charm production in agreement with data and, for p exceeding 5 GeV/c, with perturbative QCD results in proton-proton interactions. In this respect, a new scheme has been considered for the inclusive production of D mesons at large p in hadronic collisions in the framework of perturbative fragmentation functions, allowing an analysis at the Next to Leading Order (NLO) level which goes beyond the fixed I(_s³) perturbative theory of open charm production. We have applied HEMAS-DPM to the calculation of the prompt muon component for E_μ ≥ 1 TeV in air showers considering the two extreme cases of primary protons and Fe nuclei.     

AGN at very high energies: Cosmic accelerators and probes of space-time

During the last decade, very high energy astrophysics emerged as a new branch of astronomy with major discoveries achieved by the present ground-based gamma-ray Cherenkov telescopes. The sample of cosmic sources firmly detected at very high energy (VHE) now exceeds two hundred objects, including active galactic nuclei (AGN), pulsar wind nebulae, and several other types of sources of which a significant number are unidentified ones. The scientific return from recent VHE data is particularly interesting for AGN science, shedding new light on particle acceleration and emission processes around supermassive black holes, and probing the intergalactic space by the analysis of VHE photons propagating from bright remote sources to the Earth. The perspectives of this research field are promising with new generation VHE instruments such as CTA, a project of open observatory at extreme energies at the horizon 2023, allowing a deep analysis of the sky in the highest part of the electromagnetic spectrum, from 20 GeV to 300 TeV.     

Triple splitting with the F2-region of the ionosphere at high and mid-latitudes

The occurrence of the third (z-ray) component of the F2-trace on ionograms is investigated at high- and mid-latitudes. Diurnal variations show a systematic shift, with magnetic inclination, of the time of maximum occurrence. Seasonal variations show a winter maximum, and an inverse sunspot-cycle relationship exists. Maximum occurrence appears between a magnetic inclination of 70° and 80° with a fall-off either side.

Evidence is presented to suggest a z-ray association with “Spread-F” fronts, and a possible mechanism for the recording of the z-ray trace at the transmitter site is described. This involves longitudinal propagation of the o-mode at its normal reflection level, coupling at this point, and ultimate reflection for the z-ray mode as a result of sloping ionization contours belonging to “Spread-F” fronts extending in directions perpendicular to the magnetic meridian.

An association with V.L.F. emissions (“dawn-chorus”) is discussed.     

Measurable difference in Cherenkov light between gamma and hadron induced EAS

We describe the possible measurable difference in the Cherenkov light component of EAS induced by an electromagnetic particle (i.e., e⁺, e⁻ or γ) and induced by a hadron (i.e., proton or heavier nuclei) in the TeV range. The method can be applied in experiments which use wave front sampling method of EAS Cherenkov light detection (e.g., THEMISTOCLE, ASGAT).     

Gamma-Hadron Separation using Čerenkov Photon Density Fluctuations

In the atmospheric Čerenkov technique γ-rays are detected against the abundant background produced by hadronic showers. In order to improve the signal to noise ratio of theexperiment, it is necessary to reject a significant fraction of hadronic showers. Traditional background rejection methods based on image shape parameters have been extensively used for the data from imaging telescopes. However, non-imaging Čerenkov telescopes have to develop very different means of statistically identifying and removing cosmic ray events. Some of the parameters, which could be potentially important for non-imaging arrays, are the temporal and spectral differences, the lateral distributions and density fluctuations of Čerenkov photons generated by γ-ray and hadron primaries. Here we study the differences in fluctuations of Čerenkov photon density in the light pool at the observation level from showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the PACT array has been used to evaluate the efficiency of the new technique. Various types of density fluctuations like the short range and medium range fluctuations as well as flatness parameter are studied. The estimated quality factors reflect the efficiencies with which the hadrons can be rejected from the data. Since some of these parameters are independent, the cuts may be applied in tandem and we demonstrate that the proton rejection efficiency of ∼90% can be achieved. Use of density fluctuations is particularly suited for wavefront sampling observations and it seems to be a good technique to improve the signal to noise ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

Changes of the cosmic-ray mass composition in the 10–10 eV energy range

Data taken with ten Cosmic Ray Tracking (CRT) detectors and the HEGRA air-shower array on La Palma, Canary Islands, have been analysed to investigate changes of the cosmic ay mass composition at the ‘knee’ of the cosmic-ray flux spectrum near 10¹⁵ eV energy. The analysis is based on the angular distributions of particles in air showers. HEGRA data provided the shower size, direction, and core position and CRT data the particle track information. It is shown that the angular distribution of muons in air showers is sensitive to the composition over a wide range of shower sizes and, thus, primary cosmic-ray energies with little systematic uncertainties. Results can be easily expressed in terms of ln A of primary cosmic rays. In the lower part of the energy range covered, we have considerable overlap with direct composition measurements by the JACEE collaboration and find compatible results in the observed rise of ln A. Above about 10¹⁵ eV energy we find no or at most a slow further rise of ln A. Simple cosmic-ray composition models are presented which are fully consistent with our results as well as the JACEE flux and composition measurements and the flux measurements of the Tibet ASγ collaboration. Minimal three-parameter composition models defined by the same power-law slope of all elements below the knee and a common change in slope at a fixed rigidity are inconsistent with these data.     

Application of the complete linearization method to the analysis of solar prominence emission spectrum

We analysed the emission spectra of solar prominences using the complete linearization method [5] and found simultaneously the optical depth at the line centre τ₀, the doppler width of the line Δλ_D and the damping width a. The results show 1) that the complete linearization method has a larger radius of convergence, 2) that we must consider the variation of the source function with depth, when determining τ₀, and 3) that the calculated values of the damping constant for the H, Hβ of hydrogen and H and K lines of Calcium are all much greater than the theoretical values from doppler broadening and radiation damping, showing that other mechanisms besides these two also contribute to the broadening of prominence lines.     

The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Electron, muon, and hadron lateral distributions measured in air showers by the KASCADE experiment

Measurements of electron, muon and hadron lateral distributions of extensive air showers as recorded in the Karlsruhe shower core and array detector experiment are presented. The data cover the energy range from 5×10¹⁴ eV up to almost 10¹⁷ eV and extend from the inner core region to distances of 200 m. The electron and muon distributions are corrected for mutual contaminations by taking into account the detector properties in the experiment. All distributions are well described by NKG functions. The scale radii describing the electron and hadron data best are 30 and 10 m, respectively. We discuss the correlation between the scale radii and the ‘age' parameter as well as their dependence on shower size, zenith angle, and particle energy threshold.     

Optical characteristics of the whipple observatory TeV gamma-ray imaging telescope

In TeV -ray astronomy, large mirrors are used to collect erenkov light from electromagnetic cascades in the atmosphere in order to obtain low energy thresholds. The flux sensitivity of TeV -ray detectors is limited by background due to erenkov light bursts from isotropic, cosmic-ray showers which are much more numerous than -ray showers. It has recently been established that most of this background can be eliminated on the basis of the shapes of erenkov light images on the focal plane of a telescope. In order for this technique to work, the light collector must have adequate resolution over a relatively wide field of view. In this paper, the optical characteristics of the 10 m reflector used in the imaging detection of the Crab Nebula are examined and contrasted with those of a standard parabolic design. This 10 m reflector has a unique (Davies-Cotton) design with small spherical facet mirrors placed on spherical support structure with radius equal to exactly 1/2 the curvature radius of the facet mirrors. The off-axis focusing properties of this type of telescope have not been examined previously.     

On the reliability of the so-far performed tests for measuring the Lense–Thirring effect with the LAGEOS satellites

In this paper, we critically discuss the so-far performed attempts aimed at the detection of the general relativistic gravitomagnetic Lense–Thirring effect in the gravitational field of the Earth with the existing LAGEOS satellites. In the latest reported measurement of the gravitomagnetic shift with the nodes of the LAGEOS satellites and the second generation GRACE-only EIGEN-GRACE02S Earth gravity model over an observational time span of 11 years a 5–10% total accuracy is claimed at 1–3σ, respectively. We will show that, instead, it might be 15–45% (1–3σ) if the impact of the secular variations of the even zonal harmonics is considered. Possible strategies in order both to make more robust and reliable the tests with the node-only LAGEOS–LAGEOS II combination used and to overcome the problems affecting it with other alternative combinations are presented.