Instrumentation and software for the VERITAS Project and the detection of TeV Gamma-Ray Emission from the BL Lacertae Object 1ES1959+650

Transmission of Cherenkov light through the atmosphere is strongly influenced by the optical clarity of the atmosphere and the prevailing weather conditions. The performance of telescopes measuring this light is therefore dependent on atmospheric effects. This thesis presents software and hardware developed to implement a prototype sky monitoring system for use on the proposed next-generation gamma-ray telescope array, VERITAS. The system, consisting of a CCD camera and a far-infrared pyrometer, was successfully installed and tested on the ten metre atmospheric Cherenkov imaging telescope operated by the VERITAS Collaboration at the F.L. Whipple Observatory in Arizona. The thesis also presents the results of observations of the BL Lacertae object, 1ES1959+650, made with the Whipple ten metre telescope. The observations provide evidence for TeV gamma-ray emission from the BL Lacertae object, 1ES1959+650, at a level of more than 15 standard deviations above background. This represents the first unequivocal detection of this object at TeV energies, making it only the third extragalactic source seen at such levels of significance in this energy range. The flux variability of the source on a number of timescales is also investigated.

A numerical model for the gamma-ray emission of the microquasar LS 5039

The possible association between the microquasar LS 5039 and the EGRET source 3EG J1824-1514 suggests that microquasars could also be sources of high energy gamma-rays. In this paper, we explore, with a detailed numerical model, if this system can produce the emission detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. The computed spectrum reproduces the observed spectral characteristics at very high energy.

High-energy gamma-ray emission from microquasars: LS 5039 and LS I+61303

The possible associations between the microquasars LS 5039 and LS I +61 303 and the EGRET sources 3EG J1824-1514 and 3EG J0241+6103 suggest that microquasars could also be sources of high-energy gamma-rays. In this work, we present a detailed numerical inverse Compton (IC) model, based on a microquasar scenario, that reproduces the high-energy gamma-ray spectra and variability observed by EGRET for the mentioned sources. Our model considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet that interact through IC scattering with both the radiation and the magnetic fields.

Search for very high energy gamma-ray emission from pulsar-pulsar wind nebula systems with the MAGIC telescope

The MAGIC collaboration has searched for high-energy gamma-ray emission of some of the most promising pulsar candidates above an energy threshold of 50 GeV, an energy not reachable up to now by other ground-based instruments. Neither pulsed nor steady gamma-ray emission has been observed at energies of 100 GeV from the classical radio pulsars PSR J0205+6449 and PSR J2229+6114 (and their nebulae 3C58 and Boomerang, respectively) and the millisecond pulsar PSR J0218+4232. Here, we present the flux upper limits for these sources and discuss their implications in the context of current model predictions.

Variable Very-High-Energy Gamma-Ray Emission from the Microquasar LS I +61 303

Microquasars are binary star systems with relativistic radio-emitting jets. They are potential sources of cosmic rays and can be used to elucidate the physics of relativistic jets. We report the detection of variable gamma-ray emission above 100 gigaelectron volts from the microquasar LS I 61 + 303. Six orbital cycles were recorded. Several detections occur at a similar orbital phase, which suggests that the emission is periodic. The strongest gamma-ray emission is not observed when the two stars are closest to one another, implying a strong orbital modulation of the emission or absorption processes.

Gamma-ray emission from microquasars: a numerical model for LSI+61º303

We explore the possible association between the microquasar LSI +61°303 and the EGRET source 2CG 135+01/3EG J0241+6103 by studying, with a detailed numerical model, whether this system can produce the emission and the variability detected by EGRET (>100 MeV) through inverse Compton (IC) scattering. Our numerical approach considers a population of relativistic electrons entrained in a cylindrical inhomogeneous jet, interacting with both the radiation and the magnetic fields, taking into account the Thomson and Klein-Nishina regimes of interaction. Our results reproduce the observed spectral characteristics and variability at γ-rays, thus strengthening the identification of LSI +61°303 as a high-energy γ-ray source.

Gamma-ray emission from massive young stellar objects

Context.Massive stars form in dense and massive molecular cores. The exact formation mechanism is unclear, but it is possible that some massive stars are formed by processes similar to those that produce the low-mass stars, with accretion/ejection phenomena occurring at some point of the evolution of the protostar. This picture seems to be supported by the detection of a collimated stellar wind emanating from the massive protostar IRAS 16547-4247. A triple radio source is associated with the protostar: a compact core and two radio lobes. The emission of the southern lobe is clearly non-thermal. Such emission is interpreted as synchrotron radiation produced by relativistic electrons locally accelerated at the termination point of a thermal jet. Since the ambient medium is determined by the properties of the molecular cloud in which the whole system is embedded, we can expect high densities of particles and infrared photons. Because of the confirmed presence of relativistic electrons, inverse Compton and relativistic Bremsstrahlung interactions are unavoidable. Aims.We aim to make quantitative predictions of the spectral energy distribution of the non-thermal spots generated by massive young stellar objects, with emphasis on the particular case of IRAS 16547-4247. Methods.We study the high-energy emission generated by the relativistic electrons which produce the non-thermal radio source in IRAS 16547-4247. We also study the result of proton acceleration at the terminal shock of the thermal jet and make estimates of the secondary gamma rays and electron-positron pairs produced by pion decay. Results.We present spectral energy distributions for the southern lobe of IRAS 16547-4247, for a variety of conditions. We show that high-energy emission might be detectable from this object in the gamma-ray domain. The source may also be detectable in X-rays through long exposures with current X-ray instruments. Conclusions.Gamma-ray telescopes such as GLAST, and even ground-based Cherenkov arrays of new generation can be used to study non-thermal processes occurring during the formation of massive stars.

Measuring energy spectra of TeV gamma-ray emission from the Cygnus region of our galaxy with Milagro

High energy gamma rays can provide fundamental clues to the origins of cosmic rays. In this thesis, TeV gamma-ray emission from the Cygnus region is studied. Previously the Milagro experiment detected five TeV gamma-ray sources in this region and a significant excess of TeV gamma rays whose origin is still unclear. To better understand the diffuse excess the separation of sources and diffuse emission is studied using the latest and most sensitive data set of the Milagro experiment. In addition, a newly developed technique is applied that allows the energy spectrum of the TeV gamma rays to be reconstructed using Milagro data. No conclusive statement can be made about the spectrum of the diffuse emission from the Cygnus region because of its low significance of 2.2 σ above the background in the studied data sample. The entire Cygnus region emission is best fit with a power law with a spectral index of α=2.40 (68% confidence interval: 1.35-2.92) and a exponential cutoff energy of 31.6 TeV (10.0-251.2 TeV). In the case of a simple power law assumption without a cutoff energy the best fit yields a spectral index of α=2.97 (68% confidence interval: 2.83-3.10). Neither of these best fits are in good agreement with the data. The best spectral fit to the TeV emission from MGRO J2019+37, the brightest source in the Cygnus region, yields a spectral index of α=2.30 (68% confidence interval: 1.40-2.70) with a cutoff energy of 50.1 TeV (68% confidence interval: 17.8-251.2 TeV) and a spectral index of α=2.75 (68% confidence interval: 2.65-2.85) when no exponential cutoff energy is assumed. According to the present analysis, MGRO J2019+37 contributes 25% to the differential flux from the entire Cygnus at 15 TeV.

Deep observation of the NGC1275 region with MAGIC: search of diffuse gamma-ray emission from cosmic rays in the Perseus cluster

Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse γ-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E ≥ 100 GeV) for a total of 253 hr from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6º from the centre, have been detected as point-like VHE γ-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index Γ = 3.6 ± 0.2_(stat) ± 0.2_(syst) between 90 GeV and 1200 GeV. We do not detect any diffuse γ-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes α. For α = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be ≤ 1 − 2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to ≤ 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B_(0), that depend on the CR distribution. For α = 2.2, B_(0) ≥ 5 − 8 µG, which is below the ∼25 µG inferred from Faraday rotation measurements, whereas for α ≤ 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our γ-ray flux upper limits independently of the magnetic field strength.

Search for VHE gamma-ray emission from Geminga pulsar and nebula with the MAGIC telescopes

The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula have been claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 hours of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 × 10^(−13) TeV cm^(−2)s^(−1) for the Geminga pulsar and 3.5 × 10^(−12) TeV cm^(−2)s^(−1) for the surrounding nebula at 50 GeV are the most constraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 years of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.

The star forming region Monoceros R2 as a gamma-ray source

After the release of the gamma-ray source catalog produced by the Fermi satellite during its first two years of operation, a significant fraction of sources still remain unassociated at lower energies. In addition to well-known high-energy emitters (pulsars, blazars, supernova remnants, etc.), theoretical expectations predict new classes of gamma-ray sources. In particular, gamma-ray emission could be associated with some of the early phases of stellar evolution, but this interesting possibility is still poorly understood. Aims: The aim of this paper is to assess the possibility of the Fermi gamma-ray source 2FGL J0607.5-0618c being associated with the massive star forming region Monoceros R2. Methods: A multi-wavelength analysis of the Monoceros R2 region is carried out using archival data at radio, infrared, X-ray, and gamma-ray wavelengths. The resulting observational properties are used to estimate the physical parameters needed to test the different physical scenarios. Results: We confirm the 2FGL J0607.5-0618c detection with improved confidence over the Fermi two-year catalog. We find that a combined effect of the multiple young stellar objects in Monoceros R2 is a viable picture for the nature of the source.

High mass microquasars and low latitude gamma-ray sources

Population studies of unidentified EGRET sources suggest that there exist at least three different populations of galactic gamma-ray sources. One of these populations is formed by young objects distributed along the galactic plane with a strong concentration toward the inner spiral arms of the Galaxy. Variability, spectral and correlation analysis indicate that this population is not homogeneous. In particular, there is a subgroup of sources that display clear variability in their gamma-ray fluxes on timescales from days to months. Following the proposal by Kaufman Bernad\'o et al. (2002), we suggest that this group of sources might be high-mass microquasars, i.e. accreting black holes or neutron stars with relativistic jets and early-type stellar companions. We present detailed inhomogeneous models for the gamma-ray emission of these systems that include both external and synchrotron self-Compton interactions. We have included effects of interactions between the jet and all external photon fields to which it is exposed: companion star, accretion disk, and hot corona. We make broadband calculations to predict the spectral energy distribution of these objects from radio up to GeV energies. The results and predictions can be tested by present and future gamma-ray instruments like INTEGRAL, AGILE, and GLAST.

Massive protostars as gamma-ray sources.

Massive protostars have associated bipolar outflows with velocities of hundreds of km s-1. Such outflows can produce strong shocks when they interact with the ambient medium leading to regions of nonthermal radio emission. Aims. We aim at exploring under which conditions relativistic particles are accelerated at the terminal shocks of the protostellar jets and whether they can produce significant gamma-ray emission. Methods. We estimate the conditions necessary for particle acceleration up to very high energies and gamma-ray production in the nonthermal hot spots of jets associated with massive protostars embedded in dense molecular clouds. Results. We show that relativistic bremsstrahlung and proton-proton collisions can make molecular clouds with massive young stellar objects detectable by the Fermi satellite at MeV-GeV energies and by Cherenkov telescope arrays in the GeV-TeV range. Conclusions. Gamma-ray astronomy can be used to probe the physical conditions in star-forming regions and particle acceleration processes in the complex environment of massive molecular clouds.

A Search for Muon Neutrinos in Coincidence with Gamma-Ray Bursts in the Southern Hemisphere Sky Using the IceCube Neutrino Observatory

The origin of observed ultra-high energy cosmic rays (UHECRs, energies in excess of $10^{18.5}$ eV) remains unknown, as extragalactic magnetic fields deflect these charged particles from their true origin. Interactions of these UHECRs at their source would invariably produce high energy neutrinos. As these neutrinos are chargeless and nearly massless, their propagation through the universe is unimpeded and their detection can be correlated with the origin of UHECRs. Gamma-ray bursts (GRBs) are one of the few possible origins for UHECRs, observed as short, immensely bright outbursts of gamma-rays at cosmological distances. The energy density of GRBs in the universe is capable of explaining the measured UHECR flux, making them promising UHECR sources. Interactions between UHECRs and the prompt gamma-ray emission of a GRB would produce neutrinos that would be detected in coincidence with the GRB’s gamma-ray emission. The IceCube Neutrino Observatory can be used to search for these neutrinos in coincidence with GRBs, detecting neutrinos through the Cherenkov radiation emitted by secondary charged particles produced in neutrino interactions in the South Pole glacial ice. Restricting these searches to be in coincidence with GRB gamma-ray emis- sion, analyses can be performed with very little atmospheric background. Previous searches have focused on detecting muon tracks from muon neutrino interactions fromthe Northern Hemisphere, where the Earth shields IceCube’s primary background of atmospheric muons, or spherical cascade events from neutrinos of all flavors from the entire sky, with no compelling neutrino signal found. Neutrino searches from GRBs with IceCube have been extended to a search for muon tracks in the Southern Hemisphere in coincidence with 664 GRBs over five years of IceCube data in this dissertation. Though this region of the sky contains IceCube’s primary background of atmospheric muons, it is also where IceCube is most sensitive to neutrinos at the very highest energies as Earth absorption in the Northern Hemisphere becomes relevant. As previous neutrino searches have strongly constrained neutrino production in GRBs, a new per-GRB analysis is introduced for the first time to discover neutrinos in coincidence with possibly rare neutrino-bright GRBs. A stacked analysis is also performed to discover a weak neutrino signal distributed over many GRBs. Results of this search are found to be consistent with atmospheric muon backgrounds. Combining this result with previously published searches for muon neutrino tracks in the Northern Hemisphere, cascade event searches over the entire sky, and an extension of the Northern Hemisphere track search in three additional years of IceCube data that is consistent with atmospheric backgrounds, the most stringent limits yet can be placed on prompt neutrino production in GRBs, which increasingly disfavor GRBs as primary sources of UHECRs in current GRB models.

The relationship between gamma Cassiopeiae's X-ray emission and its circumstellar environment

gamma Cas is the prototypical classical Be star and is recently best known for its variable hard X-ray emission. To elucidate the reasons for this emission, we mounted a multiwavelength campaign in 2010 centered around four XMM-Newton observations. The observational techniques included long baseline optical interferometry (LBOI) from two instruments at CHARA, photometry carried out by an automated photometric telescope and H alpha observations. Because gamma Cas is also known to be in a binary, we measured radial velocities from the H alpha line and redetermined its period as 203.55 +/- 0.20 days and its eccentricity as near zero. The LBOI observations suggest that the star's decretion disk was axisymmetric in 2010, has an system inclination angle near 45 degrees, and a larger radius than previously reported. In addition, the Be star began an "outburst" at the beginning of our campaign, made visible by a brightening and reddening of the disk during our campaign and beyond. Our analyses of the new high resolution spectra disclosed many attributes also found from spectra obtained in 2001 (Chandra) and 2004 (XMM-Newton). As well as a dominant hot (approximate to 14 keV) thermal component, the familiar attributes included: (i) a fluorescent feature of Fe K even stronger than observed at previous times; (ii) strong lines of N VII and Ne XI lines indicative of overabundances; and (iii) a subsolar Fe abundance from K-shell lines but a solar abundance from L-shell ions. We also found that two absorption columns are required to fit the continuum. While the first one maintained its historical average of 1 x 10(21) cm(-2), the second was very large and doubled to 7.4 x 10(23) cm(-2) during our X-ray observations. Although we found no clear relation between this column density and orbital phase, it correlates well with the disk brightening and reddening both in the 2010 and earlier observations. Thus, the inference from this study is that much (perhaps all?) of the X-ray emission from this source originates behind matter ejected by gamma Cas into our line of sight.