Milagro Data Research Articles

Reconcilement of VHE γ -ray/X-ray correlation studies in Mrk 421 and break-down at high fluxes

Multi-wavelength campaigns have been carried out to study the correlation between the very high energy (VHE) $\gamma$-ray and the X-ray emissions in blazars but, no conclusive results have been achieved yet. In this paper, we add Milagro data to the existing VHE $\gamma$-ray data from HEGRA-CT1 and Whipple and test the consistency and robustness of the reported correlation between VHE $\gamma$-ray and X-ray fluxes in Mrk 421. We found that at monthly time scale the correlation is robust, consistent between instruments and described as a linear function. Furthermore, most of the fluxes on shorter time scales are consistent with the correlation within 3 $\sigma_A$ even, where $\sigma_A$ is an estimated intrinsic scatter. However, a break-down of the correlation becomes clearly evident at high states of activity with fluxes $\rm \gtrsim 2.5\times 10^{-10}\, cm^{-2}s^{-1}$ at energies above 400 GeV independently of the time scale, observational period or instrument, even for single flares, the X-ray and VHE $\gamma$-ray emissions lie on the correlation until the VHE $\gamma$-ray flux reaches values higher than the one mentioned above. We have interpreted our results within the one-zone synchrotron self-Compton model. We found that describing a single and unique $\gamma$-ray/X-ray correlation strongly narrows the range of possible values of the magnetic field $B$ when a constant value of the spectral index along the correlation is assumed.

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×1014g 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 100TeV 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.

SPECTRUM AND MORPHOLOGY OF THE TWO BRIGHTEST MILAGRO SOURCES IN THE CYGNUS REGION: MGRO J2019+37 AND MGRO J2031+41

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. MGRO J2019+37 is detected with a significance of 12.3 standard deviations ($\sigma$), and is better fit by a power law with an exponential cutoff than by a simple power law, with a probability $>98$% (F-test). The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of $7^{+5}_{-2}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$, a spectral index of $2.0^{+0.5}_{-1.0}$ and a cutoff energy of $29^{+50}_{-16}$ TeV. MGRO J2031+41 is detected with a significance of 7.3$\sigma$, with no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of $2.4^{+0.6}_{-0.5}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$ and a spectral index of $3.08^{+0.19}_{-0.17}$. The overall flux is subject to an $\sim$30% systematic uncertainty. The systematic uncertainty on the power law indices is $\sim$0.1. A comparison with previous results from TeV J2032+4130, MGRO J2031+41 and MGRO J2019+37 is also presented.

CONSTRAINTS ON THE EMISSION MODEL OF THE “NAKED-EYE BURST” GRB 080319B

On 2008 March 19, one of the brightest gamma-ray bursts (GRBs) ever recorded was detected by several groundand space-based instruments spanning the electromagnetic spectrum from radio to gamma rays. With a peak visual magnitude of 5.3, GRB 080319B was dubbed the “naked-eye” GRB, as an observer under dark skies could have seen the burst without the aid of an instrument. Presented here are results from observations of the prompt phase of GRB 080319B taken with the Milagro TeV observatory. The burst was observed at an elevation angle of 47 ◦ . Analysis of the data is performed using both the standard air shower method and the scaler or single-particle technique, which results in a sensitive energy range that extends from ∼ 5G eV to>20 TeV. These observations provide the only direct constraints on the properties of the high-energy gamma-ray emission from GRB 080319B at these energies. No evidence for emission is found in the Milagro data, and upper limits on the gamma-ray flux above 10 GeV are derived. The limits on emission between ∼25 and 200 GeV are incompatible with the synchrotron self-Compton model of gamma-ray production and disfavor a corresponding range (2 eV–16 eV) of assumed synchrotron peak energies. This indicates that the optical photons and soft (∼650 keV) gamma rays may not be produced by the same electron population.

OBSERVATION AND SPECTRAL MEASUREMENTS OF THE CRAB NEBULA WITH MILAGRO

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV - 100 TeV gamma ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between ~1 and ~100 TeV. The TeV emission is believed to be caused by inverse-Compton up-scattering scattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fit. The fit is consistent with values measured by IACTs between 1 and 20 TeV. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4 to 2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20 to 40 TeV.

OBSERVATION OF TeV GAMMA RAYS FROM THE CYGNUS REGION WITH THE ARGO-YBJ EXPERIMENT

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

MILAGRO OBSERVATIONS OF MULTI-TeV EMISSION FROM GALACTIC SOURCES IN THE FERMI BRIGHT SOURCE LIST

We present the result of a search of the Milagro sky map for spatial correlations with sources from a subset of the recent Fermi Bright Source List (BSL). The BSL consists of the 205 most significant sources detected above 100 MeV by the Fermi Large Area Telescope. We select sources based on their categorization in the BSL, taking all confirmed or possible Galactic sources in the field of view of Milagro. Of the 34 Fermi sources selected, 14 are observed by Milagro at a significance of 3 standard deviations or more. We conduct this search with a new analysis which employs newly optimized gamma-hadron separation and utilizes the full eight-year Milagro data set. Milagro is sensitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10 to 50 TeV depending on the source spectrum and declination. These results extend the observation of these sources far above the Fermi energy band. With the new analysis and additional data, multi-TeV emission is definitively observed associated with the Fermi pulsar, J2229.0+6114, in the Boomerang pulsar wind nebula (PWN). Furthermore, an extended region of multi-TeV emission is associated with the Fermi pulsar, J0634.0+1745, the Geminga pulsar.

Evaluation of WRF mesoscale simulations and particle trajectory analysis for the MILAGRO field campaign

Abstract. Accurate numerical simulations of the complex wind flows in the Mexico City Metropolitan Area (MCMA) can be an invaluable tool for interpreting the MILAGRO field campaign results. This paper uses three methods to evaluate numerical simulations of basin meteorology using the MM5 and WRF models: statistical comparisons with observations, "Concentration Field Analysis" (CFA) using measured air pollutant concentrations, and comparison of flow features using cluster analysis. CFA is shown to be a better indication of simulation quality than statistical metrics, and WRF simulations are shown to be an improvement on the MM5 ones. Comparisons with clusters identifies an under-representation of the drainage flows into the basin and an over-representation of wind shear in the boundary layer. Particle trajectories simulated with WRF-FLEXPART are then used to analyse the transport of the urban plume and show rapid venting and limited recirculation during MILAGRO. Lagrangian impacts were identified at the campaign supersites, and age spectra of the pollutants evaluated at those same sites. The evaluation presented in the paper show that mesoscale meteorological simulations are of sufficient accuracy to be useful for MILAGRO data analysis.

Milagro Constraints on Very High Energy Emission from Short‐Duration Gamma‐Ray Bursts

Recent rapid localizations of short, hard gamma-ray bursts (GRBs) by the Swift and HETE satellites have led to the observation of the first afterglows and the measurement of the first redshifts from this type of burst (Fox et al. 2005; Gehrels et al. 2005; Villasenor et al. 2005; Berger et al. 2005; Barthelmy et al. 2005). Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Seventeen short-duration (<5 s) GRBs detected by satellites occurred within the field of view of the Milagro gamma-ray observatory between 2000 January and 2006 December. We have searched the Milagro data for >100 GeV counterparts to these GRBs and find no significant emission correlated with these bursts. Due to the absorption of high-energy gamma rays by the extragalactic background light (EBL), detections are only expected for redshifts less than ~0.5. While most long-duration GRBs occur at redshifts higher than 0.5, the opposite is thought to be true of short GRBs. Lack of a detected VHE signal thus allows setting meaningful fluence limits. One GRB in the sample (050509b) has a likely association with a galaxy at a redshift of 0.225, while another (051103) has been tentatively linked to the nearby galaxy M81. Fluence limits are corrected for EBL absorption, either using the known measured redshift, or computing the corresponding absorption for a redshift of 0.1 and 0.5, as well as for the case of z = 0.

Search for very high energy gamma rays from WIMP annihilations near the Sun with the Milagro detector

The neutralino, the lightest stable supersymmetric particle, is a strong theoretical candidate for the missing astronomical ``dark matter''. A profusion of such neutralinos can accumulate near the Sun when they lose energy upon scattering and are gravitationally captured. Pair-annihilations of those neutralinos may produce very high energy (VHE, above $100 GeV$) gamma-rays. Milagro is an air shower array which uses the water Cherenkov technique to detect extensive air showers and is capable of observing VHE gamma-rays from the direction of the Sun with an angular resolution of $0.75^{\circ}$. Analysis of Milagro data with an exposure to the Sun of 1165 hours presents the first attempt to detect TeV gamma-rays produced by annihilating neutralinos captured by the Solar system and shows no statistically significant signal. Resulting limits that can be set on gamma-ray flux due to near-Solar neutralino annihilations and on neutralino cross-section are presented.