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.     

Io: Volcanic thermal sources and global heat flow

We have examined thermal emission from 240 active or recently-active volcanic features on Io and quantified the magnitude and distribution of their volcanic heat flow during the Galileo epoch. We use spacecraft data and a geological map of Io to derive an estimate of the maximum possible contribution from small dark areas not detected as thermally active but which nevertheless appear to be sites of recent volcanic activity. We utilize a trend analysis to extrapolate from the smallest detectable volcanic heat sources to these smallest mapped dark areas. Including the additional heat from estimates for “outburst” eruptions and for a multitude of very small (“myriad”) hot spots, we account for ～62 × 10¹² W (～59 ± 7% of Io’s total thermal emission). Loki Patera contributes, on average, 9.6 × 10¹² W (～9.1 ± 1%). All dark paterae contribute 45.3 × 10¹² W (～43 ± 5%). Although dark flow fields cover a much larger area than dark paterae, they contribute only 5.6 × 10¹² W (～5.3 ± 0.6%). Bright paterae contribute ～2.6 × 10¹² W (～2.5 ± 0.3%). Outburst eruption phases and very small hot spots contribute no more than ～4% of Io’s total thermal emission: this is probably a maximum value. About 50% of Io’s volcanic heat flow emanates from only 1.2% of Io’s surface. Of Io’s heat flow, 41 ± 7.0% remains unaccounted for in terms of identified sources. Globally, volcanic heat flow is not uniformly distributed. Power output per unit surface area is slightly biased towards mid-latitudes, although there is a stronger bias toward the northern hemisphere when Loki Patera is included. There is a slight favoring of the northern hemisphere for outbursts where locations were well constrained. Globally, we find peaks in thermal emission at ～315°W and ～105°W (using 30° bins). There is a minimum in thermal emission at around 200°W (almost at the anti-jovian longitude) which is a significant regional difference. These peaks and troughs suggest a shift to the east from predicted global heat flow patterns resulting from tidal heating in an asthenosphere. Global volcanic heat flow is dominated by thermal emission from paterae, especially from Loki Patera (312°W, 12°N). Thermal emission from dark flows maximises between 165°W and 225°W. Finally, it is possible that a multitude of very small hot spots, smaller than the present angular resolution detection limits, and/or cooler, secondary volcanic processes involving sulphurous compounds, may be responsible for at least part of the heat flow that is not associated with known sources. Such activity should be sought out during the next mission to Io.     

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.     

Two-photon positron–electron annihilation in a strong magnetic field

We consider the two-photon positron and electron annihilation in flight, it means the annihilating particles exhibit the middly relativistic momenta in a super strong magnetic field. Such particles are present in the corona of pulsars and magnetars. The paper presents how the total emission rate for the two-photon process is affected not only by magnetic field but also by the relativistic momentum of the annihilating particles. We found that the momenta influence significantly the total emission rate and the directions of the emitted photons. Additionally, the total emission for the two-photon process is comparable to that for the one-photon process at the momentum of annihilating particles of about m₀, where m₀ is the electron mass, and the magnetic field close to the critical Schwinger value of 4.41 × 10¹³ G. The latter is reported as a main annihilation channel in a super strong magnetic field. We calculated also the energetic spectra of annihilating photons emitted, which are also affected by the magnetic field and the momenta of the annihilating particles.     

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.     

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.     

Characteristics of planetary-scale waves simulated by a new venusian mesosphere and thermosphere general circulation model

We have developed a new general circulation model (GCM) for the venusian mesosphere and thermosphere (80-about 180 km). Our GCM simulations show that winds in the subsolar-to-antisolar direction (SS–AS) are predominant above about 90 km. A weak return flow of the SS–AS is seen below about 90 km. We performed GCM simulations imposing the planetary-scale waves (thermal tides, Rossby wave, and Kelvin wave) at the lower boundary. Although the diurnal and semidiurnal tides are damped below 95 km, the Rossby wave propagates up to around 130 km. However, the amplitude of the Rossby wave is too small (<1 m/s) to affect the general circulation. On the other hand, the Kelvin wave propagates up to about 130 km with a maximum zonal wind fluctuation of approximately 5.9 m/s on average. The amplitude of the Kelvin wave sometimes exceeds 10 m/s around the terminator. The Kelvin wave causes a temporal variation in the wind velocity at the altitude of the O₂-1.27 μm nightglow emission (about 95 km). Using a newly developed 1-D nightglow model and the composition distribution calculated from our GCM, we investigated the impact of the Kelvin wave on the nightglow distribution. Our results suggest that the Kelvin wave would cause temporal variations in the nightglow emission in the 23:50–00:20 LT region with an intensity of 1.1–1.3 MR and a period of approximately 4 days.     

Rheologies and ages of lava flows on Elysium Mons,Mars

We present results of our study of the rheologies and ages of lava flows in the Elysium Mons region of Mars. Previous studies have shown that the geometric dimensions of lava flows reflect rheological properties such as yield strength, effusion rate and viscosity. In this study the rheological properties of lava flows in the Elysium Mons region were determined and compared to the rheologies of the Ascraeus Mons lava flows. We also derived new crater size-frequency distribution measurements (CSFDs) for the Elysium lava flows to identify possible changes in the rheological properties with time. In addition, possible changes in the rheological properties with the distance from the caldera of Elysium Mons were analyzed.In total, 35 lava flows on and around Elysium Mons were mapped, and divided into three groups, lava flows on the flanks of Elysium Mons, in the plains between the three volcanoes Elysium Mons, Hecates and Albor Tholus and lava flows south of Albor Tholus. The rheological properties of 32 of these flows could be determined. Based on our morphometric measurements of each individual lava flow, estimates for the yield strengths, effusion rates, viscosities, and eruption duration of the studied lava flows were made. The yield strengths of the investigated lava flows range from ～3.8 × 10² Pa to ～1.5 × 10⁴ Pa, with an average of ～3.0 × 10³ Pa. These yield strengths are in good agreement with estimates for terrestrial basaltic lava flows. The effusion rates are on average ～747 m³ s^?1, ranging from ～99 to 4450 m³ s^?1. The viscosities are on average ～4.1 × 10⁶ Pa s, with a range of 1.2 × 10⁵ Pa s to 3.1 × 10⁷ Pa s. The eruption durations of the flows were calculated to be between 6 and 183 days, with an average of ～51 days. The determined rheological properties are generally very similar to those of other volcanic regions on Mars, such as on Ascraeus Mons in the Tharsis region. Calculated yield strengths and viscosities point to a basaltic/andesitic composition of the lava flows, similar to basaltic or andesitic a’a lava flows on Earth.Absolute model ages of all 35 lava flows on Elysium Mons were derived from crater size-frequency distribution measurements (CSFD). The derived model ages show a wide variation from about 632 Ma to 3460 Ma. Crater size-frequency distribution measurements of the Elysium Mons caldera show an age of ～1640 Ma, which is consistent with the resurfacing age of Werner (2009). Significant changes of the rheologies with time could not be observed. Similarly, we did not observe systematic changes in ages with increasing distances of lava flows from the Elysium Mons caldera.     

Preliminary results of the spectral analysis of Suzaku data of SW Ursae Majoris and BZ Ursae Majoris

Herein, we present the preliminary results of the spectral analysis of the Suzaku data of SW Ursae Majoris (SW UMa) and BZ Ursae Majoris (BZ UMa) which were obtained in their quiescent states. The aim of this study was to analyse 0.5–6.0 keV energy range X-ray photons of the Suzaku satellite and to determine their emission mechanisms and the areas where they were radiated from these two dwarf novae. No study on this energy range for BZ UMa has been reported in the literature. We analysed the Suzaku data in the (0.5–6.0) keV and in (0.5–2.5) keV energy bands of the systems separatelly to search for possible differences between the emission mechanisms in the energy ranges; however, none were found. We obtained the best-fitted spectral models, flux, luminosity, temperature and mass accretion rate values for the systems. The white dwarf radius of BZ UMa, which has not been reported in the literature, was found to be 6.88 × 10⁸ cm. From this study, we concluded that the mass accretion rate values for both systems were lower than the critical mass accretion rate value in the boundary layers of cataclysmic variables, which indicates that the boundary layers are composed of optically thin thermal X-ray emitting gas in their quiescent state; these results support those of previous studies regarding SW UMa and some studies regarding BZ UMa. Another result was that there can be a coronal structure above the boundary layer which act as soft X-ray emitter for SW UMa.     

High energy neutrinos from the TeV Blazar 1ES 1959 + 650

The AMANDA neutrino telescope has recently reported the detection of high-energy neutrinos spatially and temporally coincident with the flaring of the TeV blazar 1ES 1959 + 650. If high-energy neutrinos are in fact generated by this blazar, it would be the first strong evidence for the hadronic acceleration of cosmic rays. At present, the statistical significance of this observation cannot be reliably assessed, however. In this letter, we investigate whether circumstances exist where the source can produce the flux implied by the coincident events. We show that if the TeV gamma-ray emission observed from 1ES 1959 + 650 or other nearby TeV blazars is the result of accelerated protons interacting with nucleons, it is reasonable that AMANDA could detect several events during a flaring period. Such rates require that the spectral index of the source be rather high (for instance ∼2.8 for 1ES 1959 + 650) and that the Lorentz factor of the jet be fairly small (Γ ∼ 1).     

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.     

The abundance and vertical distribution of the unknown ultraviolet absorber in the venusian atmosphere from analysis of Venus Monitoring Camera images

Observations of Venus using the ultraviolet filter of the Venus Monitoring Camera (VMC) on ESA’s Venus Express Spacecraft (VEX) provide the best opportunity for study of the spatial and temporal distribution of the venusian unknown ultraviolet absorber since the Pioneer Venus (PV) mission. We compare the results of two sets of 125 radiative transfer models of the upper atmosphere of Venus to each pixel in a subset of VMC UV channel images. We use a quantitative best fit criterion based upon the notion that the distribution of the unknown absorber should be independent of the illumination and observing geometry. We use the product of the cosines of the incidence and emission angles and search for absorber distributions that are uncorrelated with this geometric parameter, finding that two models can describe the vertical distribution of the unknown absorber. One model is a well-mixed vertical profile above a pressure level of roughly 120 mb (～63 km). This is consistent with the altitude of photochemical formation of sulfuric acid. The second model describes it as a thin layer of pure UV absorber at a pressure level roughly around 24 mb (～71 km) and this altitude is consistent with the top of upper cloud deck. We find that the average abundance of unknown absorber in the equatorial region is 0.21 ± 0.04 optical depth and it decreases in the polar region to 0.08 ± 0.05 optical depth at 365 nm.     

Atomic oxygen on the Venus nightside: Global distribution deduced from airglow mapping

The Visible and Infra-Red Thermal Imaging Spectrometer (VIRTIS) instrument on board the Venus Express spacecraft has measured the O₂(a¹Δ) nightglow distribution at 1.27 μm in the Venus mesosphere for more than two years. Nadir observations have been used to create a statistical map of the emission on Venus nightside. It appears that the statistical 1.6 MR maximum of the emission is located around the antisolar point. Limb observations provide information on the altitude and on the shape of the emission layer. We combine nadir observations essentially covering the southern hemisphere, corrected for the thermal emission of the lower atmosphere, with limb profiles of the northern hemisphere to generate a global map of the Venus nightside emission at 1.27 μm. Given all the O₂(a¹Δ) intensity profiles, O₂(a¹Δ) and O density profiles have been calculated and three-dimensional maps of metastable molecular and atomic oxygen densities have been generated. This global O density nightside distribution improves that available from the VTS3 model, which was based on measurements made above 145 km. The O₂(a¹Δ) hemispheric average density is 2.1 × 10⁹ cm^?3, with a maximum value of 6.5 × 10⁹ cm^?3 at 99.2 km. The O density profiles have been derived from the nightglow data using CO₂ profiles from the empirical VTS3 model or from SPICAV stellar occultations. The O hemispheric average density is 1.9 × 10¹¹ cm^?3 in both cases, with a mean altitude of the peak located at 106.1 km and 103.4 km, respectively. These results tend to confirm the modeled values of 2.8 × 10¹¹ cm^?3 at 104 km and 2.0 × 10¹¹ cm^?3 at 110 km obtained by Brecht et al. [Brecht, A., Bougher, S.W., Gérard, J.-C., Parkinson, C.D., Rafkin, S., Foster, B., 2011a. J. Geophys. Res., in press] and Krasnopolsky [Krasnopolsky, V.A., 2010. Icarus 207, 17–27], respectively. Comparing the oxygen density map derived from the O₂(a¹Δ) nightglow observations, it appears that the morphology is very different and that the densities obtained in this study are about three times higher than those predicted by the VTS3 model.     

Robust limits on Lorentz violation from gamma-ray bursts

We constrain the possibility of a non-trivial refractive index in free space corresponding to an energy-dependent velocity of light: c(E) ≃ c₀(1 − E/M), where M is a mass scale that might represent effect of quantum-gravitational space-time foam, using the arrival times of sharp features observed in the intensities of radiation with different energies from a large sample of gamma-ray bursters (GRBs) with known redshifts. We use wavelet techniques to identify genuine features, which we confirm in simulations with artificial added noise. Using the weighted averages of the time-lags calculated using correlated features in all the GRB light curves, we find a systematic tendency for more energetic photons to arrive earlier. However, there is a very strong correlation between the parameters characterizing an intrinsic time-lag at the source and a distance-dependent propagation effect. Moreover, the significance of the earlier arrival times is less evident for a subsample of more robust spectral structures. Allowing for intrinsic stochastic time-lags in these features, we establish a statistically robust lower limit: M > 0.9 × 10¹⁶ GeV on the scale of violation of Lorentz invariance.     

Sodium atoms in the lunar exotail: Observed velocity and spatial distributions

The lunar sodium tail extends long distances due to radiation pressure on sodium atoms in the lunar exosphere. Our earlier observations measured the average radial velocity of sodium atoms moving down the lunar tail beyond Earth (i.e., near the anti-lunar point) to be ～12.5 km/s. Here we use the Wisconsin H-alpha Mapper to obtain the first kinematically resolved maps of the intensity and velocity distribution of this emission over a 15° × 15 ° region on the sky near the anti-lunar point. We present both spatially and spectrally resolved observations obtained over four nights bracketing new Moon in October 2007. The spatial distribution of the sodium atoms is elongated along the ecliptic with the location of the peak intensity drifting 3° east along the ecliptic per night. Preliminary modeling results suggest the spatial and velocity distributions in the sodium exotail are sensitive to the near surface lunar sodium velocity distribution. Future observations of this sort along with detailed modeling offer new opportunities to describe the time history of lunar surface sputtering over several days.     

The XMM-Newton view of the yellow hypergiant IRC + 10420 and its surroundings

Among evolved massive stars likely in transition to the Wolf–Rayet phase, IRC + 10420 is probably one of the most enigmatic. It belongs to the category of yellow hypergiants and it is characterized by quite high mass loss episodes. Even though IRC + 10420 benefited of many observations in several wavelength domains, it has never been a target for an X-ray observatory. We report here on the very first dedicated observation of IRC + 10420 in X-rays, using the XMM-Newton satellite. Even though the target is not detected, we derive X-ray flux upper limits of the order of 1–3 × 10⁻¹⁴ erg cm⁻² s⁻¹ (between 0.3 and 10.0 keV), and we discuss the case of IRC + 10420 in the framework of emission models likely to be adequate for such an object. Using the Optical/UV Monitor on board XMM-Newton, we present the very first upper limits of the flux density of IRC + 10420 in the UV domain (between 1800 and 2250 Å and between 2050 and 2450 Å). Finally, we also report on the detection in this field of 10 X-ray and 7 UV point sources, and we briefly discuss their properties and potential counterparts at longer wavelengths.     

A photometric study of the W UMa-type binary DF Canum Venaticorum

We present new photometric observations for the eclipsing binary DF CVn, and determined five light minimum times. By using the Wilson–Devinney code, two sets of photometric solutions were deduced from our observations in 2009. The asymmetric light curves obtained on 2009 March 5 were modeled by a dark spot on the more massive component. The results indicate that DF CVn is a W-type weak-contact binary, with a mass ratio of q ～ 0.28 and an overcontact degree of f ～ 20%. From the O ? C curve of minimum times, it is found that there exists a cyclic variation, whose period and semi-amplitude are P₃ = 17.2(±0.9) year and A = 0.^d0070(±0.^d0008), respectively. This kind of cyclic oscillation may possibly result from the light-time effect due to the presence of an unseen third body. This kind of additional body may extract angular momentum from the central system. The low-amplitude changes of the light curves on a short-time scale (e.g., half a month) may be attributed to the dark spot activity, which may result in angular momentum loss via magnetic breaking. With angular momentum loss, the weak-contact binary DF CVn will evolve into a deep-contact configuration.