Tests of Statistical Significance and Background Estimation in Gamma Ray Air Shower Experiments

In this paper we discuss several methods of significance calculation and point out the limits of their applicability. We then introduce a self consistent scheme for source detection and discuss some of its properties. The method allows incorporating background anisotropies by vetoing existing small scale regions on the sky and compensating for known large scale anisotropies. By giving an example using Milagro gamma ray observatory we demonstrate how the method can be employed to relax the detector stability assumption. Two practical implementations of the method are discussed. The method is universal and can be used with any large field-of-view detector, where the object of investigation, steady or transient, point or extended, traverses its field of view.Comment: 12 pages, 7 figures. Submitted to ApJ, submission number: 5833

TeV Gamma-Ray Sources from a Survey of the Galactic Plane with Milagro

A survey of Galactic gamma-ray sources at a median energy of ~20 TeV has been performed using the Milagro Gamma Ray Observatory. Eight candidate sources of TeV emission are detected with pre-trials significance $>4.5\sigma$ in the region of Galactic longitude $l\in[30^\circ,220^\circ]$ and latitude $b\in[-10^\circ,10^\circ]$. Four of these sources, including the Crab nebula and the recently published MGRO J2019+37, are observed with significances $>4\sigma$ after accounting for the trials involved in searching the 3800 square degree region. All four of these sources are also coincident with EGRET sources. Two of the lower significance sources are coincident with EGRET sources and one of these sources is Geminga. The other two candidates are in the Cygnus region of the Galaxy. Several of the sources appear to be spatially extended. The fluxes of the sources at 20 TeV range from ~25% of the Crab flux to nearly as bright as the Crab.Comment: Submitted to Ap

Discovery of Localized Regions of Excess 10-TeV Cosmic Rays

An analysis of 7 years of Milagro data performed on a 10-degree angular scale has found two localized regions of excess of unknown origin with greater than 12 sigma significance. Both regions are inconsistent with gamma-ray emission with high confidence. One of the regions has a different energy spectrum than the isotropic cosmic-ray flux at a level of 4.6 sigma, and it is consistent with hard spectrum protons with an exponential cutoff, with the most significant excess at ~10 TeV. Potential causes of these excesses are explored, but no compelling explanations are found.Comment: Submitted to PhysRevLet

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. 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.Comment: Accepted for publication in the Astrophysical Journa

A Measurement of the Spatial Distribution of Diffuse TeV Gamma Ray Emission from the Galactic Plane with Milagro

Diffuse $\gamma$-ray emission produced by the interaction of cosmic-ray particles with matter and radiation in the Galaxy can be used to probe the distribution of cosmic rays and their sources in different regions of the Galaxy. With its large field of view and long observation time, the Milagro Gamma Ray Observatory is an ideal instrument for surveying large regions of the Northern Hemisphere sky and for detecting diffuse $\gamma$-ray emission at very high energies. Here, the spatial distribution and the flux of the diffuse $\gamma$-ray emission in the TeV energy range with a median energy of 15 TeV for Galactic longitudes between 30$^\circ$ and 110$^\circ$ and between 136$^\circ$ and 216$^\circ$ and for Galactic latitudes between -10$^\circ$ and 10$^\circ$ are determined. The measured fluxes are consistent with predictions of the GALPROP model everywhere except for the Cygnus region ($l\in[65^\circ,85^\circ]$). For the Cygnus region, the flux is twice the predicted value. This excess can be explained by the presence of active cosmic ray sources accelerating hadrons which interact with the local dense interstellar medium and produce gamma rays through pion decay.Comment: 15 pages, 3 figures, accepted by Ap

The Study of TeV Variability and Duty Cycle of Mrk 421 from 3 Years of Observations with the Milagro Observatory

TeV flaring activity with time scales as short as tens of minutes and an orphan TeV flare have been observed from the blazar Markarian 421 (Mrk 421). The TeV emission from Mrk 421 is believed to be produced by leptonic synchrotron self-Compton (SSC) emission. In this scenario, correlations between the X-ray and the TeV fluxes are expected, TeV orphan flares are hardly explained and the activity (measured as duty cycle) of the source at TeV energies is expected to be equal or less than that observed in X-rays if only SSC is considered. To estimate the TeV duty cycle of Mrk 421 and to establish limits on its variability at different time scales, we continuously observed Mrk 421 with the Milagro observatory. Mrk 421 was detected by Milagro with a statistical significance of 7.1 standard deviations between 2005 September 21 and 2008 March 15. The observed spectrum is consistent with previous observations by VERITAS. We estimate the duty cycle of Mrk 421 for energies above 1 TeV for different hypothesis of the baseline flux and for different flare selections and we compare our results with the X-ray duty cycle estimated by Resconi et al. 2009. The robustness of the results is discussed.Comment: 27 pages, 6 figures, ApJ accepte

The Intrinsic Origin of Spin Echoes in Dipolar Solids Generated by Strong Pi Pulses

In spectroscopy, it is conventional to treat pulses much stronger than the linewidth as delta-functions. In NMR, this assumption leads to the prediction that pi pulses do not refocus the dipolar coupling. However, NMR spin echo measurements in dipolar solids defy these conventional expectations when more than one pi pulse is used. Observed effects include a long tail in the CPMG echo train for short delays between pi pulses, an even-odd asymmetry in the echo amplitudes for long delays, an unusual fingerprint pattern for intermediate delays, and a strong sensitivity to pi-pulse phase. Experiments that set limits on possible extrinsic causes for the phenomena are reported. We find that the action of the system's internal Hamiltonian during any real pulse is sufficient to cause the effects. Exact numerical calculations, combined with average Hamiltonian theory, identify novel terms that are sensitive to parameters such as pulse phase, dipolar coupling, and system size. Visualization of the entire density matrix shows a unique flow of quantum coherence from non-observable to observable channels when applying repeated pi pulses.Comment: 24 pages, 27 figures. Revised from helpful referee comments. Added new Table IV, new paragraphs on pages 3 and 1

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.Comment: 11 pages, 10 figure