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.     

Compton scattered transition radiation from very high energy particles

X-ray transition radiation can be used to measure the Lorentz factor of relativistic particles. At energies approaching γ=E/mc²=10⁵, transition radiation detectors can be optimized by using thick (5–10 mil) foils with large (5–10 mm) spacings. This implies X-ray energies 100 keV and the use of scintillators as the X-ray detectors. Compton scattering of the X-rays out of the particle beam then becomes an important effect. We discuss the design of very high energy detectors, the use of metal radiator foils rather than the standard plastic foils, inorganic scintillators for detecting Compton scattered transition radiation, and the application to the ACCESS cosmic ray experiment.     

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.     

The configuration of the loop-like structure of the solar atmosphere

An integral, governing steady flows in an isolated thin magnetic flux tube in the hydrostatic plane-stratified atmosphere, has been obtained. The integral, that we named as the shape integral, is expressed as (1 − M_A²)B cos θ = const. Here M_A² is the Alfven Mach number, B is the magnetic field strength and θ is the flux tube inclination to the horizontal. The shape integral should hold for most loop models because it represents just the momentum balance laws and has no relation to any energy balance mode. Its application to the isothermal and static cases is discussed and illustrated.     

A theoretical study of the direction and polarization of solar radio millisecond spike radiation

We investigated the angular direction and polarization of the solar radio millisecond spike emission in the model in which the spike emission is due to the second harmonic instability modes driven by electron cyclotron maser of loss cone distributed electrons during the propagation of a nonlinear plasma density wave near the magnetic mirror. We found that, when the angle θ between the wave vector and the magnetic field is > 60 °, the emission is in 100% X-mode polarization; when 40 ° < θ<60 °, the emission is in 100% O-mode polarization provided the amplitude of the density wave is below a certain limit; above that limit, the polarization will fall from 100% O-mode to even the X-mode. We also found that only 0.1% of the free energy of energy carrying electrons in the source region is converted into radiation wave energy.     

Mangeto-optical effects in sunspots and their vector magnetic field

The system of transfer equations of the four Stokes parameters I, Q, U, V under the action of the magneto-optical effect (i.e. the Unno-Beckers equations) are numerically solved in this paper for the magneto-sensitive lines FeI λλ 6302.499 and 5324.191 using an appropriate sunspot model. The errors in the expressions for the coefficients _r and _W in Beckers' paper [2] have been corrected for. From the results of calculations, features of the profiles of the Stokes parameters dependent on the magnetic vector have been isolated. Our computations also show that the magneto-optical effect should be taken into consideration in the measurement of the vector magnetic fields.

In the fourth section of this paper we have established a simple and convenient method for obtaining-information on the magnetic vector (including the field strength B, its inclination to the line of sight γ and its azimuth χ) from the profiles of the Stokes parameters. It consists of three steps: (1) The value of B is determined from the distance of the highest point in the V-profile from the central line. (2) γ is then found from V_max, i.e maximum value of V. (3) Lastly, the angle χ is found from Q₀, i.e. the value of Q at line centre.     

Implications of modified AGN spectra due to absorption by infrared photons

Recent observations of high energy photons from active galactic nuclei (AGN) have lead to a renewed interest in the effect of gamma-ray cascading in background radiation fields. In recent years, numerous authors have explored possible modifications to the observed high energy photon spectra from AGN due to propagation through such fields. This paper will re-examine a number of these issues. The major conclusions of this paper are: (1) The reaction γ + γ → e⁺ + e⁻ provides a highly sensitive probe of background infrared (IR) fields. However little can be inferred about the specific nature of the IR background or the effect on spectra from more distant AGN from a study of nearby sources alone. Currently, only upper limits to the background IR density can be established. (2) The contribution of secondary photons from pair-cascading off of background microwave and infrared radiation is most likely unobservable in the regime of 1 TeV unless the strength of extragalactic magnetic fields are much less than 10⁻¹³ G and the inherent source spectra continue to much higher energies. The possible contribution from this process may be conclusively ruled out through further spectral measurements of AGN in the TeV regime during high and low states. (3) There is little hope of unambiguously extracting the values of either the Hubble constant or Ω from γ-γ attenuation measurements. (4) The sensitivity of the attenuation to the density of the background IR produced prior to the epoch of the observed source suggests a future possibility to probe directly the evolution of the background IR radiation through multiple source observations.     

Prospects of solar axion searches with crystal detectors

A calculation of the expected signal due to Primakoff coherent conversion of solar axions into photons via Bragg scattering in several crystal detectors is presented. The results are confronted with the experimental sensitivities of present and future experiments concluding that the sensitivity of crystal detectors does not challenge the globular cluster limit on the axion-photon coupling g_aγγ. In particular, in the axion mass window m_a 0.03 eV explored with this technique (not accessible at present by other methods), g_aγγ might be constrained down to 10⁻⁹ GeV⁻¹ (the recent helioseismological bound) provided that significant improvements in the parameters and performances of these detectors be achieved and large statistics accumulated. This bound should be considered as a minimal goal for the sensitivity of future crystal experiments. Consequently, finding a positive signal at this level of sensitivity would necessarily imply revisiting other more stringent astrophysical limits derived for the same range of m_a values.     

General relativistic effects in the γ-ray radiation from pulsars

General relativistic and quantum-mechanical effects in the γ-ray radiation from neutron stars are studied in this paper. For three values of the magnetic dipole moment, we calculated the upper limit of the escaping energy of the γ-photon as a function of the angle between the line of sight and the magnetic axis, with and without these effects. We conclude that for low magnetic moments, the quantum-mechanical effects are unimportant but the general relativistic effects are to reduce significantly the upper limit. For high magnetic moments, the two effects nearly cancel out and the net effect is small. We discuss the implications of these results on the γ-ray spectrum.     

Probing pulsar winds using inverse Compton scattering

We investigate the effects of inverse Compton scattering by electrons and positrons in the unshocked winds of rotationally-powered binary pulsars. This process can scatter low energy target photons to produce gamma rays with energies from MeV to TeV. The binary radio pulsars PSR B1259−63 and PSR J0045−73 are both in close eccentric orbits around bright main sequence stars which provide a huge density of low energy target photons. The inverse Compton scattering process transfers momentum from the pulsar wind to the scattered photons, and therefore provides a drag which tends to decelerate the pulsar wind. We present detailed calculations of the dynamics of a pulsar wind which is undergoing inverse Compton scattering, showing that the deceleration of the wind of PSR B1259−63 due to ‘inverse Compton drag' is small, but that this process may confine the wind of PSR J0045−73 before it attains pressure balance with the outflow of its companion star. We calculate the spectra and light curves of the resulting inverse Compton emission from PSR B1259−63 and show that if the size of the pulsar wind nebula is comparable to the binary separation, then the γ-ray emission from the unshocked wind may be detectable by atmospheric Cherenkov detectors or by the new generation of satellite-borne γ-ray detectors such as INTEGRAL and GLAST. This mechanism may therefore provide a direct probe of the freely-expanding regions of pulsar winds, previously thought to be invisible.     

Gamma ray spectrum of the crab nebula in the multi TeV region

The THEMISTOCLE array of 18 Cherenkov detectors which has a 3 TeV gamma energy threshold, has detected a signal from the Crab nebula at a 5.8 standard deviation level. Information on the energy spectrum is obtained in the range 3–15 TeV. The integrated flux can be fitted with the form, Φ (> E) = (3.7 ± 0.5) × 10^-12 (E/5)-^{−1.5 ± 0.20} cm⁻² s⁻¹ (E in TeV) compatible with the extrapolation of results at lower energies. The Crab signal is used to measure the angular resolution of the multi-telescope technique. The value obtained is 2.3 mr (0.15°) in agreement with the results of simulations, and confirms the interest of this new method for multi-TeV gamma-ray detection.     

Neutrinos from SN1987A: flavor conversion and interpretation of results

After recent results from solar neutrino experiments and KamLAND we can definitely say that neutrinos from SN1987A underwent flavor conversion, and the conversion effects must be taken into account in the analysis of the data. Assuming the normal mass hierarchy of neutrinos we calculate the permutation factors p for the Kamiokande-2, IMB and Baksan detectors. The conversion inside the star leads to p=0.28–0.32; the oscillations in the matter of the Earth give partial (and different for different detectors) regeneration of the original signal, reducing this factor down to 0.15–0.20 (at E=40 MeV). We study in details the influence of conversion on the observed signal depending on the parameters of the original neutrino spectra. For a given set of these parameters, the conversion could lead to an increase of the average energy of the observed events up to 50% and of the number of events by a factor of 2 at Kamiokande-2 and by a factor of 3–5 at IMB. Inversely, we find that neglecting the conversion effects can lead up to 50% error in the determination of the average energy of the original spectrum and about 50% error in the original luminosity. Comparing our calculations with experimental data we conclude that the Kamiokande-2 data alone do not favor strong conversion effect, which testifies for small difference of the original and spectra. In contrast, the combined analysis of the Kamiokande and IMB results slightly favors strong conversion effects (that is, large difference of the original spectra). In comparison with the no-oscillation case, the latter requires lower average energy and higher luminosity of the original flux.     

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).     

Are diffuse high energy neutrinos and γ-rays from starburst galaxies observable?

Loeb and Waxman have argued that high energy neutrinos from the decay of pions produced in interactions of cosmic rays with interstellar gas in starburst galaxies would be produced with a large enough flux to be observable. Their model is reexamined here and we obtain an upper limit to the diffuse neutrino flux from starburst galaxies. The upper limit obtained here is a factor of 5 lower than the flux which they predict. Our predicted neutrino flux would be below the atmospheric neutrino foreground flux at energies below 300 TeV and therefore would be unobservable. Compared with predicted fluxes from other extragalactic high energy neutrino sources, starburst neutrinos with PeV energies would have a flux considerably below that predicted for AGN models.

We also estimate an upper limit for the diffuse GeV γ-ray flux from starbust galaxies to be of the observed γ-ray background, much less than the component from unresolved blazars and more than an order of magnitude below the estimate of Thompson et al.     

The self-gravity of Uranian rings

We derive an approximate, analytic expression for the self-gravity of Uranian elliptical rings at any point in the N-streamline formalism (N → ∞), and numerically evaluate by a Monte Carlo calculation, the correction factor γ that corrects for the approximations used in the derivation. We confirm that the 2-streamline (N = 2) model overestimates the mass of the ring: for the ε, , ßrings, the mass given by the N-streamline model (N → ∞) is only 2/5 that by the 2-streamline model.     

The pionic contribution to diffuse γ-rays: upper limits

Diffuse γ-rays probe the highest-energy processes at the largest scales. Here we derive model-independent constraints on the hadronic contribution to the Galactic and extragalactic γ-ray spectra at in the energy range 50 MeVE_γ10 GeV. The hadronic component is dominated by emission from neutral pions, with a characteristic spectrum symmetric about m_π⁰/2. We exploit the well-defined properties of the pion decay spectrum to quantify the maximum pionic fraction of the observed γ-ray intensity. We find that the Galactic spectrum above 30 MeV can be atmost about 50% pionic. The maximum pionic contribution to the extragalactic spectrum is energy dependent; it also depends on the redshift range over which the sources are distributed, ranging from as low as about 20% for pions generated very recently, to as much as 90% if the pions are generated around redshift 10. The implications of these constraints for models of γ-ray and neutrino emission are briefly discussed.     

Density profiles of dark matter halos with anisotropic velocity tensors

We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ³ ∝ r⁻, and (ii) the velocity anisotropy parameter is given by the relation β_a(r)=β₁+2β₂ (r/r_*)/[1+(r/r_*)²] where β₁, β₂ are parameters and r_* equals twice the virial radius, r_vir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10^−2.5r_vir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β₁ and for increasing concentration of the system. On the other hand, slopes at r=r_vir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10⁻³r_vir and r_vir and connects kinematic and structural characteristics of spherical systems is described.     

Separation of gamma and hadron initiated air showers with energies between 20 and 500 TeV

The discrimination between air showers initiated by γ rays and by hadrons is one of the fundamental problems in experimental cosmic-ray physics. The physics of this ‘γ/hadron separation’ is discussed in this paper. We restrict ourselves to the energy range from about 20 to 500 TeV, and take only the information contained in the lateral Čerenkov light distribution and the number of electrons at the detector level into consideration. An understanding of the differences between air showers generated by γ rays and those due to hadrons leads us to formulate suitable observables for the separation process. Angle integrating Čerenkov arrays (AICA) offer a promising new approach to ground-based γ-ray astronomy in the energy region from about 20 to 500 TeV. In order to establish this technique, an efficient suppression of the overwhelming hadronic background radiation is required. As an example for our general discussion, we present one method for γ/hadron separation in AICAs called ‘LES’. It is based on the simultaneous determination of the shower size and some characteristic parameters of the lateral distribution of the Čerenkov light. The potential inherent within this technique is demonstrated in quantitative detail for the existing ‘AIROBICC’ AICA. We also propose an objective measure of the intrinsic sensitivity of a detection scheme in ground-based γ-ray astronomy, the ‘reduced quality factor’. It is shown that AICAs may reach a sensitivity to γ-ray point sources in the high VHE range similar to that of the Čerenkov-telescope imaging technique in the low VHE region.     

Diagnostics and error analysis of the magnetic field strength in solar microwave burst regions

This paper points out that the previous method of finding the magnetic field strength using the approximate formula of the maximum frequency of microwave burst of Dulk and Marsh gives non-unique and self-inconsistent results, and with relative errors as high as −3.0±5.3, whereas my recently proposed three-parameter (ν_p, β, T) method will give unique and self-consistent solutions of the field strength, with greatly reduced errors (±0.48).     

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae will be detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, Λ, on the spectrum, since some recent cosmological observations strongly suggest the existence of Λ. We find following results. (1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimates (6–10 MeV). (2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae are dominant, is about 25h₅₀²(R_SN/0.1 yr⁻¹)(n_Gh₅₀⁻³/0.02 Mpc⁻³) events per year, where R_SN is the supernova rate in a galaxy, n_G is the number density of galaxies, and h₅₀ = H₀/(50 km/s Mpc), where H₀ is the Hubble constant. (3) The event rate is almost insensitive to Λ. The flux increases in the low energy side (< 10 MeV) with increasing Λ, but decreases in the high energy side (> 10 MeV) in models in which the integrated number of supernovae in one galaxy is fixed.