We use HII starburst galaxy apparent magnitude measurements to constrain cosmological parameters in six cosmological models. A joint analysis of HII galaxy, quasar angular size, baryon acoustic oscillations peak length scale, and Hubble parameter measurements result in relatively model-independent and restrictive estimates of the current values of the non-relativistic matter density parameter $\Omega_{\rm m_0}$ and the Hubble constant $H_0$. These estimates favor a 2.0$\sigma$ to 3.4$\sigma$ (depending on cosmological model) lower $H_0$ than what is measured from the local expansion rate. The combined data are consistent with dark energy being a cosmological constant and with flat spatial hypersurfaces, but do not strongly rule out mild dark energy dynamics or slightly non-flat spatial geometries.

The radio source 3C 264, hosted by the giant elliptical galaxy NGC 3862, was observed with VERITAS between February 2017 and May 2019. These deep observations resulted in the discovery of very-high-energy (VHE; E $>100$ GeV) $\gamma$-ray emission from this active galaxy. An analysis of $\sim$57 hours of quality-selected live time yields a detection at the position of the source, corresponding to a statistical significance of 7.8 standard deviations above background. The observed VHE flux is variable on monthly time scales, with an elevated flux seen in 2018 observations. The VHE emission during this elevated state is well-characterized by a power-law spectrum with a photon index $\Gamma = 2.20 \pm 0.27$ and flux F($>315$ GeV) = ($7.6\pm 1.2_{\mathrm stat} \pm 2.3_{\mathrm syst})\times 10^{-13}$ cm$^{-2}$ s$^{-1}$, or approximately 0.7% of the Crab Nebula flux above the same threshold. 3C 264 ($z = 0.0217$) is the most distant radio galaxy detected at VHE, and the elevated state is thought to be similar to that of the famously outbursting jet in M 87. Consequently, extensive contemporaneous multi-wavelength data were acquired in 2018 at the time of the VHE high state. An analysis of these data, including VLBA, VLA, HST, Chandra and Swift observations in addition to the VERITAS data, is presented, along with a discussion of the resulting spectral energy distribution.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

In this paper I examine how well Ryskin's model of emergent cosmic acceleration fits several sets of cosmological observations. I find that while Ryskin's model is somewhat compatible with the standard model of cosmological acceleration ($\Lambda$CDM) for low redshift (z < 1) measurements, its predictions diverge considerably from those of the standard model for measurements made at high redshift (for which z > 1), and it is therefore not a compelling substitute for the standard model.

Characteristic variability timescales in blazar gamma-ray light curves can provide insight into the physical processes responsible for the gamma-ray variability. The power spectral density (PSD) is capable of revealing such timescales, which may appear as breaks or periodicities. Continuous-time autoregressive moving-average (CARMA) models can be used to accurately estimate a light curve's PSD. Through a lightcurve simulation study, we develop a methodology to identify PSD breaks using CARMA models. Using this methodology, we study the gamma-ray light curves of 13 bright blazars observed with the Fermi Large Area Telescope in the 0.1-300 GeV band over 9.5 years. We present the blazar gamma-ray PSDs, which provide evidence for low-frequency breaks on timescales ~1 year in four sources, and an additional high-frequency break on a timescale ~9 days in one source.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

The High Altitude Water Cherenkov (HAWC) array is sensitive to high energy Cosmic Rays (CR) in the $\sim 10$ to $\sim 200$ TeV energy range, making it possible to construct maps of the so called "Sun Shadow" ($SS$), i. e. of the deficit of CR coming from the direction of the Sun. In this work, we present the variation of the Relative Intensity of the deficit ($SS_{RI}$) for three years of HAWC observations form 2016 to 2018 in which we found a clear decreasing trend of the ($SS_{RI}$) over the studied period, corresponding to the declining phase of the solar cycle 24. By comparing the $SS_{RI}$ with the photospheric magnetic field evolution, we show that there is a linear relationship between the $SS_{RI}$ and the median photospheric magnetic field of the Active Region belt (-40$^\circ \le$ lat $\le$ 40$^\circ$) and a inverse linear relationship with the polar photospheric magnetic field (lat $\ge \pm$ 60$^\circ$). The former relationship is due to the magnetic field causing a deviation of the CR, whereas the latter reflects the change of the heliospheric field topology from multipolar to dipolar configurations. These relationships are valid only when the median magnetic field is lower than 8 G, during the declining and minimum phases of the solar cycle 24. Finally, we show that relativistic charged particles, in the 10 to 200 TeV energy range, are deflected a few degrees.

Strong electric fields in thunderclouds have long been known to accelerate secondary charged particles. We investigate this effect using three years (2015-2017) of data from the scalar system of the High Altitude Water Cherenkov (HAWC) observatory, which is an air shower array deployed 4100 m a.s.l. in central Mexico. The experimental site is frequently affected by strong thunderstorms, and the detector's high altitude, large area, and high sensitivity to cosmic-ray air showers make it ideal for investigating particle acceleration due to the electric fields present inside the thunder storm clouds. In particular, the scaler system of HAWC records the output of each one of the 1200 PMTs as well as the 2, 3, and 4-fold multiplicities (logic AND in a time window of 30 ns) of each water Cherenkov detectors (WCD) with a sampling rate of 40 Hz. Using data from this scaler system, we identify approximately 100 increases in the scaler rate which is in time coincidence with thunderstorms. These events show high cross correlation between the scaler rate and the electric field, hence can be produced by the acceleration of secondary particle by the thunderstorm electric fields. In this work we present the method of identification of these events and their general characteristics.

The variation in atmospheric pressure is due to changes in mass of the air column above, which in turn resembles the density variation of atmosphere and will affect the decay of secondary particles of cosmic rays. The ground based cosmic ray detectors observe pressure dependent variation in their flux. The High Altitude Water Cherenkov (HAWC) gamma ray observatory is a great detector of secondary particles because of its high altitude, high uptime, and large area (including total photo-cathode area), which makes the HAWC scaler system an ideal instrument for solar modulation studies. Although, in order to perform these studies it is necessary to isolate and remove the atmospheric modulations. The observed rate in each PMT has signatures of both the solar and atmospheric modulations, which makes it difficult to measure the pressure coefficient ($\beta_P$ ). The pressure at the HAWC site shows a periodic behavior ($\sim$ 12 hours), which also reflects in the scalar rates. This periodic property was used to isolate the pressure modulation and $\beta_P$ were estimated with accuracy. Since the pressure dependence is a physical phenomenon, the estimated coefficients for PMTs should be identical, any deviation from this can be due to malfunction of the PMT. This make this method a useful tool to identify the malfunctioning PMTs and help us to isolate them from the analysis. In this analysis we are presenting the method of estimation of the pressure coefficients and its usage to correct the HAWC scalar data to make it suitable for the solar modulations studies.

The General Antiparticle Spectrometer (GAPS) will carry out a sensitive dark matter search by measuring low-energy ($\mathrm{E} < 0.25 \mathrm{GeV/nucleon}$) cosmic ray antinuclei. The primary targets are low-energy antideuterons produced in the annihilation or decay of dark matter. At these energies antideuterons from secondary/tertiary interactions are expected to have very low fluxes, significantly below those predicted by well-motivated, beyond the standard model theories. GAPS will also conduct low-energy antiproton and antihelium searches. Combined, these observations will provide a powerful search for dark matter and provide the best observations to date on primordial black hole evaporation on Galactic length scales. The GAPS instrument detects antinuclei using the novel exotic atom technique. It consists of a central tracker with a surrounding time-of-flight (TOF) system. The tracker is a one cubic meter volume containing 10 cm-diameter lithium-drifted silicon (Si(Li)) detectors. The TOF is a plastic scintillator system that will both trigger the Si(Li) tracker and enable better reconstruction of particle tracks. After coming to rest in the tracker, antinuclei will form an excited exotic atom. This will then de-excite via characteristic X-ray transitions before producing a pion/proton star when the antiparticle annihilates with the nucleus. This unique event topology will give GAPS the nearly background-free detection capability required for a rare-event search. Here we present the scientific motivation for the GAPS experiment, its design and its current status as it prepares for flight in the austral summer of 2021-22.

In this paper we use all available baryon acoustic oscillation, Hubble parameter, and quasar angular size data to constrain six dark energy cosmological models, both spatially flat and non-flat. Depending on the model and data combination considered, these data mildly favor closed spatial hypersurfaces (by as much as $1.7\sigma$) and dark energy dynamics (up to a little over $2\sigma$) over a cosmological constant $\Lambda$. The data also favor, at $1.8\sigma$ to $3.4\sigma$, depending on the model and data combination, a lower Hubble constant than what is measured from the local expansion rate.

Little is known about the origin of the high-energy and sustained emission from solar Long-Duration Gamma-Ray Flares (LDGRFs), identified with the Compton Gamma Ray Observatory (CGRO), the Solar Maximum Mission (SMM), and now Fermi. Though Fermi/Large Area Space Telescope (LAT) has identified dozens of flares with LDGRF signature, the nature of this phenomenon has been a challenge to explain both due to the extreme energies and long durations. The highest-energy emission has generally been attributed to pion production from the interaction of >300 MeV protons with the ambient matter. The extended duration suggests that particle acceleration occurs over large volumes extending high in the corona, either from stochastic acceleration within large coronal loops or from back precipitation from coronal mass ejection driven shocks. It is possible to test these models by making direct comparison between the properties of the accelerated ion population producing the gamma-ray emission derived from the Fermi/LAT observations, and the characteristics of solar energetic particles (SEPs) measured by the Payload for Matter-Antimatter Exploration and Light Nuclei Astrophysics (PAMELA) spacecraft in the energy range corresponding to the pion-related emission detected with Fermi. For fourteen of these events we compare the two populations -- SEPs in space and the interacting particles at the Sun -- and discuss the implications in terms of potential sources. Our analysis shows that the two proton numbers are poorly correlated, with their ratio spanning more than five orders of magnitude, suggesting that the back precipitation of shock-acceleration particles is unlikely the source of the LDGRF emission.

The General AntiParticle Spectrometer (GAPS) is a balloon-borne instrument designed to detect cosmic-ray antimatter using the novel exotic atom technique, obviating the strong magnetic fields required by experiments like AMS, PAMELA, or BESS. It will be sensitive to primary antideuterons with kinetic energies of $\approx0.05-0.2$ GeV/nucleon, providing some overlap with the previously mentioned experiments at the highest energies. For $3\times35$ day balloon flights, and standard classes of primary antideuteron propagation models, GAPS will be sensitive to $m_{\mathrm{DM}}\approx10-100$ GeV c$^{-2}$ WIMPs with a dark-matter flux to astrophysical flux ratio approaching 100. This clean primary channel is a key feature of GAPS and is crucial for a rare event search. Additionally, the antiproton spectrum will be extended with high statistics measurements to cover the $0.07 \leq E \leq 0.25 $ GeV domain. For $E>0.2$ GeV GAPS data will be complementary to existing experiments, while $E<0.2$ GeV explores a new regime. The first flight is scheduled for late 2020 in Antarctica. These proceedings will describe the astrophysical processes and backgrounds relevant to the dark matter search, a brief discussion of detector operation, and construction progress made to date.

Despite the significant progress achieved in recent years, the physical mechanisms underlying the origin of solar energetic particles (SEPs) are still a matter of debate. The complex nature of both particle acceleration and transport poses challenges to developing a universal picture of SEP events that encompasses both the low-energy (from tens of keV to a few hundreds of MeV) observations made by space-based instruments and the GeV particles detected by the worldwide network of neutron monitors in ground-level enhancements (GLEs). The high-precision data collected by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment offer a unique opportunity to study the SEP fluxes between $\sim$80 MeV and a few GeV, significantly improving the characterization of the most energetic events. In particular, PAMELA can measure for the first time with good accuracy the spectral features at moderate and high energies, providing important constraints for current SEP models. In addition, the PAMELA observations allow the relationship between low and high-energy particles to be investigated, enabling a clearer view of the SEP origin. No qualitative distinction between the spectral shapes of GLE, sub-GLE and non-GLE events is observed, suggesting that GLEs are not a separate class, but are the subset of a continuous distribution of SEP events that are more intense at high energies. While the spectral forms found are to be consistent with diffusive shock acceleration theory, which predicts spectral rollovers at high energies that are attributed to particles escaping the shock region during acceleration, further work is required to explore the relative influences of acceleration and transport processes on SEP spectra.

We use all available baryon acoustic oscillation distance measurements and Hubble parameter data to constrain the cosmological constant $\Lambda$, dynamical dark energy, and spatial curvature in simple cosmological models. We find that the consensus spatially flat $\Lambda$CDM model provides a reasonable fit to the data, but depending on the Hubble constant prior and cosmological model, it can be a little more than 1$\sigma$ away from the best-fit model, which can favor mild dark energy dynamics or non-flat spatial hypersurfaces.

New results on the short-term galactic cosmic ray (GCR) intensity variation (Forbush decrease) in December 2006 measured by the PAMELA instrument are presented. Forbush decreases are sudden suppressions of the GCR intensities which are associated with the passage of interplanetary transients such as shocks and interplanetary coronal mass ejections (ICMEs). Most of the past measurements of this phenomenon were carried out with ground-based detectors such as neutron monitors or muon telescopes. These techniques allow only the indirect detection of the overall GCR intensity over an integrated energy range. For the first time, thanks to the unique features of the PAMELA magnetic spectrometer, the Forbush decrease commencing on 2006 December 14, following a CME at the Sun on 2006 December 13 was studied in a wide rigidity range (0.4 - 20 GV) and for different species of GCRs detected directly in space. The daily averaged GCR proton intensity was used to investigate the rigidity dependence of the amplitude and the recovery time of the Forbush decrease. Additionally, for the first time, the temporal variations in the helium and electron intensities during a Forbush decrease were studied. Interestingly, the temporal evolutions of the helium and proton intensities during the Forbush decrease were found in good agreement, while the low rigidity electrons (< 2 GV) displayed a faster recovery. This difference in the electron recovery is interpreted as a charge-sign dependence introduced by drift motions experienced by the GCRs during their propagation through the heliosphere.

We report on the measurement of the $^{7}$Be($n, p$)$^{7}$Li cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux experimental area (EAR2) of the n_TOF facility at CERN. This reaction plays a key role in the lithium yield of the Big Bang Nucleosynthesis (BBN) for standard cosmology. The only two previous time-of-flight measurements performed on this reaction did not cover the energy window of interest for BBN, and showed a large discrepancy between each other. The measurement was performed with a Si-telescope, and a high-purity sample produced by implantation of a $^{7}$Be ion beam at the ISOLDE facility at CERN. While a significantly higher cross section is found at low-energy, relative to current evaluations, in the region of BBN interest the present results are consistent with the values inferred from the time-reversal $^{7}$Li($p, n$)$^{7}$Be reaction, thus yielding only a relatively minor improvement on the so-called Cosmological Lithium Problem (CLiP). The relevance of these results on the near-threshold neutron production in the p+$^{7}$Li reaction is also discussed.

The High Altitude Water Cherenkov observatory (HAWC) is an air shower array devised for TeV gamma-ray astronomy. HAWC is located at an altitude of 4100 m a.s.l. in Sierra Negra, Mexico. HAWC consists of 300 Water Cherenkov Detectors, each instrumented with 4 photomultiplier tubes (PMTs). HAWC re-uses the Front-End Boards from the Milagro experiment to receive the PMT signals. These boards are used in combination with Time to Digital Converters (TDCs) to record the time and the amount of light in each PMT hit (light flash). A set of VME TDC modules (128 channels each) is operated in a continuous (dead time free) mode. The TDCs are read out via the VME bus by Single-Board Computers (SBCs), which in turn are connected to a gigabit Ethernet network. The complete system produces ~ 500 MB/s of raw data. A high-throughput data processing system has been designed and built to enable real-time data analysis. The system relies on off-the-shelf hardware components, an open-source software technology for data transfers (ZeroMQ) and a custom software framework for data analysis (AERIE). Multiple trigger and reconstruction algorithms can be combined and run on blocks of data in a parallel fashion, producing a set of output data streams which can be analyzed in real time with minimal latency (< 5 s). This paper provides an overview of the hardware set-up and an in-depth description of the software design, covering both the TDC data acquisition system and the real-time data processing system. The performance of these systems is also discussed.

List of proceedings from the HAWC Collaboration presented at the 35th International Cosmic Ray Conference, 12 July - 20 July 2017, Bexco, Busan, Korea.

We present the results of the first high-altitude balloon flight test of a concept for an advanced Compton telescope making use of modern scintillator materials with silicon photomultiplier (SiPM) readouts. There is a need in the fields of high-energy astronomy and solar physics for new medium-energy gamma-ray (~0.4 - 10 MeV) detectors capable of making sensitive observations. A fast scintillator- based Compton telescope with SiPM readouts is a promising solution to this instrumentation challenge, since the fast response of the scintillators permits the rejection of background via time-of-flight (ToF) discrimination. The Solar Compton Telescope (SolCompT) prototype was designed to demonstrate stable performance of this technology under balloon-flight conditions. The SolCompT instrument was a simple two-element Compton telescope, consisting of an approximately one-inch cylindrical stilbene crystal for a scattering detector and a one-inch cubic LaBr3:Ce crystal for a calorimeter detector. Both scintillator detectors were read out by 2 x 2 arrays of Hamamatsu S11828-3344 MPPC devices. Custom front-end electronics provided optimum signal rise time and linearity, and custom power supplies automatically adjusted the SiPM bias voltage to compensate for temperature-induced gain variations. A tagged calibration source, consisting of ~240 nCi of Co-60 embedded in plastic scintillator, was placed in the field of view and provided a known source of gamma rays to measure in flight. The SolCompT balloon payload was launched on 24 August 2014 from Fort Sumner, NM, and spent ~3.75 hours at a float altitude of ~123,000 feet. The instrument performed well throughout the flight. After correcting for small (~10%) residual gain variations, we measured an in-flight ToF resolution of ~760 ps (FWHM). Advanced scintillators with SiPM readouts continue to show great promise for future gamma-ray instruments.

Data from the PAMELA satellite experiment were used to perform a detailed measurement of under-cutoff protons at low Earth orbits. On the basis of a trajectory tracing approach using a realistic description of the magnetosphere, protons were classified into geomagnetically trapped and re-entrant albedo. The former include stably-trapped protons in the South Atlantic Anomaly, which were analyzed in the framework of the adiabatic theory, investigating energy spectra, spatial and angular distributions; results were compared with the predictions of the AP8 and the PSB97 empirical trapped models. The albedo protons were classified into quasi-trapped, concentrating in the magnetic equatorial region, and un-trapped, spreading over all latitudes and including both short-lived (precipitating) and long-lived (pseudo-trapped) components. Features of the penumbra region around the geomagnetic cutoff were investigated as well. PAMELA observations significantly improve the characterization of the high energy proton populations in near Earth orbits.

The PAMELA satellite-borne experiment is providing first direct measurements of Solar Energetic Particles (SEPs) with energies from $\sim$80 MeV to several GeV in near-Earth space. Its unique observational capabilities include the possibility of measuring the flux angular distribution and thus investigating possible anisotropies related to SEP events. This paper focuses on the analysis methods developed to estimate SEP energy spectra as a function of the particle pitch angle with respect to the Interplanetary Magnetic Field (IMF). The crucial ingredient is provided by an accurate simulation of the asymptotic exposition of the PAMELA apparatus, based on a realistic reconstruction of particle trajectories in the Earth's magnetosphere. As case study, the results of the calculation for the May 17, 2012 event are reported.

The PAMELA satellite experiment is providing first direct measurements of Solar Energetic Particles (SEPs) with energies from about 80 MeV to several GeV in near-Earth space, bridging the low energy data by other space-based instruments and the Ground Level Enhancement (GLE) data by the worldwide network of neutron monitors. Its unique observational capabilities include the possibility of measuring the flux angular distribution and thus investigating possible anisotropies. This work reports the analysis methods developed to estimate the SEP energy spectra as a function of the particle pitch-angle with respect to the Interplanetary Magnetic Field (IMF) direction. The crucial ingredient is provided by an accurate simulation of the asymptotic exposition of the PAMELA apparatus, based on a realistic reconstruction of particle trajectories in the Earth's magnetosphere. As case study, the results for the May 17, 2012 event are presented.

Data from the PAMELA satellite experiment were used to measure the geomagnetic cutoff for high-energy ($\gtrsim$ 80 MeV) protons during the solar particle events on 2006 December 13 and 14. The variations of the cutoff latitude as a function of rigidity were studied on relatively short timescales, corresponding to single spacecraft orbits (about 94 minutes). Estimated cutoff values were cross-checked with those obtained by means of a trajectory tracing approach based on dynamical empirical modeling of the Earth's magnetosphere. We find significant variations in the cutoff latitude, with a maximum suppression of about 6 deg for $\sim$80 MeV protons during the main phase of the storm. The observed reduction in the geomagnetic shielding and its temporal evolution were compared with the changes in the magnetosphere configuration, investigating the role of IMF, solar wind and geomagnetic (Kp, Dst and Sym-H indexes) variables and their correlation with PAMELA cutoff results.

List of proceedings from the HAWC Collaboration presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

The nature of particle acceleration at the Sun, whether through flare reconnection processes or through shocks driven by coronal mass ejections (CMEs), is still under scrutiny despite decades of research. The measured properties of solar energetic particles (SEPs) have long been modeled in different particle-acceleration scenarios. The challenge has been to disentangle to the effects of transport from those of acceleration. The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instrument, enables unique observations of SEPs including composition and the angular distribution of the particles about the magnetic field, i.e. pitch angle distribution, over a broad energy range (>80 MeV) -- bridging a critical gap between space-based measurements and ground-based. We present high-energy SEP data from PAMELA acquired during the 2012 May 17 SEP event. These data exhibit differential anisotropies and thus transport features over the instrument rigidity range. SEP protons exhibit two distinct pitch angle distributions; a low-energy population that extends to 90{\deg} and a population that is beamed at high energies (>1 GeV), consistent with neutron monitor measurements. To explain a low-energy SEP population that exhibits significant scattering or redistribution accompanied by a high-energy population that reaches the Earth relatively unaffected by dispersive transport effects, we postulate that the scattering or redistribution takes place locally. We believe these are the first comprehensive measurements of the effects of solar energetic particle transport in the Earth's magnetosheath.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

The High-Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays and charged cosmic rays at TeV energies. The detector is still under construction, but data acquisition with the partially deployed detector started in 2013. An analysis of the cosmic-ray arrival direction distribution based on $4.9\times 10^{10}$ events recorded between June 2013 and February 2014 shows anisotropy at the $10^{-4}$ level on angular scales of about $10^\circ$. The HAWC cosmic-ray sky map exhibits three regions of significantly enhanced cosmic-ray flux; two of these regions were first reported by the Milagro experiment. A third region coincides with an excess recently reported by the ARGO-YBJ experiment. An angular power spectrum analysis of the sky shows that all terms up to $\ell=15$ contribute significantly to the excesses.

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 x 10^14 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.

We argue that the hour-long neutron transient detected by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) Neutron Spectrometer beginning at 15:45 UT on 2011 June 4 is due to secondary neutrons from energetic protons interacting in the spacecraft. The protons were probably accelerated by a shock that passed the spacecraft about thirty minutes earlier. We reach this conclusion after a study of data from the MESSENGER neutron spectrometer, gamma-ray spectrometer, X-ray Spectrometer, and Energetic Particle Spectrometer, and from the particle spectrometers on STEREO A. Our conclusion differs markedly from that given by Lawrence et al. [2014] who claimed that there is "strong evidence" that the neutrons were produced by the interaction of ions in the solar atmosphere. We identify significant faults with the authors' arguments that led them to that conclusion.

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.

We describe the science motivation and development of a pair production telescope for medium-energy gamma-ray polarimetry. Our instrument concept, the Advanced Energetic Pair Telescope (AdEPT), takes advantage of the Three-Dimensional Track Imager, a low-density gaseous time projection chamber, to achieve angular resolution within a factor of two of the pair production kinematics limit (~0.6 deg at 70 MeV), continuum sensitivity comparable with the Fermi-LAT front detector (<3x10-6 MeV cm-2 s-1 at 70 MeV), and minimum detectable polarization less than 10% for a 10 millicrab source in 106 seconds.

This paper reports the results from three targeted searches of Milagro TeV sky maps: two extragalactic point source lists and one pulsar source list. The first extragalactic candidate list consists of 709 candidates selected from the Fermi-LAT 2FGL catalog. The second extragalactic candidate list contains 31 candidates selected from the TeVCat source catalog that have been detected by imaging atmospheric Cherenkov telescopes (IACTs). In both extragalactic candidate lists Mkn 421 was the only source detected by Milagro. This paper presents the Milagro TeV flux for Mkn 421 and flux limits for the brighter Fermi-LAT extragalactic sources and for all TeVCat candidates. The pulsar list extends a previously published Milagro targeted search for Galactic sources. With the 32 new gamma-ray pulsars identified in 2FGL, the number of pulsars that are studied by both Fermi-LAT and Milagro is increased to 52. In this sample, we find that the probability of Milagro detecting a TeV emission coincident with a pulsar increases with the GeV flux observed by the Fermi-LAT in the energy range from 0.1 GeV to 100 GeV.

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.

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.

We present the sensitivity of HAWC to Gamma Ray Bursts (GRBs). HAWC is a very high-energy gamma-ray observatory currently under construction in Mexico at an altitude of 4100 m. It will observe atmospheric air showers via the water Cherenkov method. HAWC will consist of 300 large water tanks instrumented with 4 photomultipliers each. HAWC has two data acquisition (DAQ) systems. The main DAQ system reads out coincident signals in the tanks and reconstructs the direction and energy of individual atmospheric showers. The scaler DAQ counts the hits in each photomultiplier tube (PMT) in the detector and searches for a statistical excess over the noise of all PMTs. We show that HAWC has a realistic opportunity to observe the high-energy power law components of GRBs that extend at least up to 30 GeV, as it has been observed by Fermi LAT. The two DAQ systems have an energy threshold that is low enough to observe events similar to GRB 090510 and GRB 090902b with the characteristics observed by Fermi LAT. HAWC will provide information about the high-energy spectra of GRBs which in turn could help to understanding about e-pair attenuation in GRB jets, extragalactic background light absorption, as well as establishing the highest energy to which GRBs accelerate particles.

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.

Spacecraft observations in the inner heliosphere offer the first opportunity to measure 1-10 MeV solar neutrons. We discuss the physics of low-energy neutron production in solar flares and show that, even at interacting-particle energies of 2 MeV/nucleon, neutrons with energies >10 MeV are produced. On the other hand, a significant fraction of 1-10 MeV neutrons result from interactions of >10 MeV/nucleon ions in typical flare spectra. We calculate the escaping neutron spectra for mono-energetic and power-law particle spectra at the Sun for the location and observation angle of MESSENGER at the time of its reported detection of low-energy neutrons associated with the 2007 December 31 solar flare. We detail concerns about this questionable observation of solar neutrons: 1. the inferred number of accelerated protons at the Sun for this modest M2-class flare was 10X larger than any flare observed to date, 2. the onset and duration of the 'solar' neutron count rate was similar to that of the solar energetic particles (SEPs), and 3. the authors' argument that the SEPs were dominated by electrons and so could not have produced the neutron counts locally in the spacecraft. In contrast we argue that solar energetic protons and alpha particles, through local neutron production and accidental coincidences, were the source of most of the reported 'solar-neutron' counts.

The coming decade will see the routine use of solar data of unprecedented spatial and spectral resolution, time cadence, and completeness. To capitalize on the new (or soon to be available) facilities such as SDO, ATST and FASR, and the challenges they present in the visualization and synthesis of multi-wavelength datasets, we propose that realistic, sophisticated, 3D active region and flare modeling is timely and critical, and will be a forefront of coronal studies over the coming decade. To make such modeling a reality, a broad, concerted effort is needed to capture the wealth of information resulting from the data, develop a synergistic modeling effort, and generate the necessary visualization, interpretation and model-data comparison tools to accurately extract the key physics.

Gamma-ray astronomy presents an extraordinary scientific potential for the study of the most powerful sources and the most violent events in the Universe. In order to take full advantage of this potential, the next generation of instrumentation for this domain will have to achieve an improvement in sensitivity over present technologies of at least an order of magnitude. The DUAL mission concept takes up this challenge in two complementary ways: a very long observation of the entire sky, combined with a large collection area for simultaneous observations of Type Ia SNe. While the Wide-Field Compton Telescope (WCT) accumulates data from the full gamma-ray sky (0.1-10 MeV) over the entire mission lifetime, the Laue-Lens Telescope (LLT) focuses on 56Co emission from SNe Ia (0.8-0.9 MeV), collecting gamma-rays from its large area crystal lens onto the WCT. Two separated spacecraft flying in formation will maintain the DUAL payloads at the lens' focal distance.

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 8-year Milagro dataset. Milagro is sensitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10-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.

Results are presented of a harmonic analysis of the large scale cosmic-ray anisotropy as observed by the Milagro observatory. We show a two-dimensional display of the sidereal anisotropy pro jections in right ascension generated by the fitting of three harmonics to 18 separate declination bands. The Milagro observatory is a water Cherenkov detector located in the Jemez mountains near Los Alamos, New Mexico. With a high duty cycle and large field-of-view, Milagro is an excellent instrument for measuring this anisotropy with high sensitivity at TeV energies. The analysis is conducted using a seven year data sample consisting of more than 95 billion events, the largest such data set in existence. We observe an anisotropy with a magnitude around 0.1% for cosmic rays with a median energy of 6 TeV. The dominant feature is a deficit region of depth (2.49 +/- 0.02 stat. +/- 0.09 sys.)x10^(-3) in the direction of the Galactic North Pole centered at 189 degrees right ascension. We observe a steady increase in the magnitude of the signal over seven years.

We propose to advance investigations of electromagnetic radiation originating in atomic nuclei beyond its current infancy to a true astronomy. This nuclear emission is independent from conditions of gas, thus complements more traditional stronomical methods used to probe the nearby universe. Radioactive gamma-rays arise from isotopes which are made in specific locations inside massive stars, their decay in interstellar space traces an otherwise not directly observable hot and tenuous phase of the ISM, which is crucial for feedback from massive stars. Its intrinsic clocks can measure characteristic times of processes within the ISM. Frontier questions that can be addressed with studies in this field are the complex interiors of massive stars and supernovae which are key agents in galactic dynamics and chemical evolution, the history of star-forming and supernova activity affecting our solar-system environment, and explorations of occulted and inaccessible regions of young stellar nurseries in our Galaxy.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

The Gamma-RAy Polarimeter Experiment (GRAPE) is a concept for an astronomical hard X-ray Compton polarimeter operating in the 50 - 500 keV energy band. The instrument has been optimized for wide-field polarization measurements of transient outbursts from energetic astrophysical objects such as gamma-ray bursts and solar flares. The GRAPE instrument is composed of identical modules, each of which consists of an array of scintillator elements read out by a multi-anode photomultiplier tube (MAPMT). Incident photons Compton scatter in plastic scintillator elements and are subsequently absorbed in inorganic scintillator elements; a net polarization signal is revealed by a characteristic asymmetry in the azimuthal scattering angles. We have constructed a prototype GRAPE module containing a single CsI(Na) calorimeter element, at the center of the MAPMT, surrounded by 60 plastic elements. The prototype has been combined with custom readout electronics and software to create a complete "engineering model" of the GRAPE instrument. This engineering model has been calibrated using a nearly 100% polarized hard X-ray beam at the Advanced Photon Source at Argonne National Laboratory. We find modulation factors of 0.46 +/- 0.06 and 0.48 +/- 0.03 at 69.5 keV and 129.5 keV, respectively, in good agreement with Monte Carlo simulations. In this paper we present details of the beam test, data analysis, and simulations, and discuss the implications of our results for the further development of the GRAPE concept.

POET (Polarimeters for Energetic Transients) is a Small Explorer mission concept proposed to NASA in January 2008. The principal scientific goal of POET is to measure GRB polarization between 2 and 500 keV. The payload consists of two wide FoV instruments: a Low Energy Polarimeter (LEP) capable of polarization measurements in the energy range from 2-15 keV and a high energy polarimeter (Gamma-Ray Polarimeter Experiment -- GRAPE) that will measure polarization in the 60-500 keV energy range. Spectra will be measured from 2 keV up to 1 MeV. The POET spacecraft provides a zenith-pointed platform for maximizing the exposure to deep space. Spacecraft rotation will provide a means of effectively dealing with systematics in the polarization response. POET will provide sufficient sensitivity and sky coverage to measure statistically significant polarization for up to 100 GRBs in a two-year mission. Polarization data will also be obtained for solar flares, pulsars and other sources of astronomical interest.

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.

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.

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.

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.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.

The Advanced Compton Telescope (ACT), the next major step in gamma-ray astronomy, will probe the fires where chemical elements are formed by enabling high-resolution spectroscopy of nuclear emission from supernova explosions. During the past two years, our collaboration has been undertaking a NASA mission concept study for ACT. This study was designed to (1) transform the key scientific objectives into specific instrument requirements, (2) to identify the most promising technologies to meet those requirements, and (3) to design a viable mission concept for this instrument. We present the results of this study, including scientific goals and expected performance, mission design, and technology recommendations.

This paper reviews the development status of GRAPE (the Gamma-Ray Polarimeter Experiment), a hard X-ray Compton Polarimeter. The purpose of GRAPE is to measure the polarization of hard X-rays in the 50-300 keV energy range. We are particularly interested in X-rays that are emitted from solar flares and gamma-ray bursts (GRBs), although GRAPE could also be employed in the study of other astrophysical sources. Accurately measuring the polarization of the emitted radiation will lead to a better understating of both emission mechanisms and source geometries. The GRAPE design consists of an array of plastic scintillators surrounding a central high-Z crystal scintillator. The azimuthal distribution of photon scatters from the plastic array into the central calorimeter provides a measure of the polarization fraction and polarization angle of the incident radiation. The design of the detector provides sensitivity over a large field-of-view (>pi steradian). The design facilitates the fabrication of large area arrays with minimal deadspace. This paper presents the latest design concept and the most recent results from laboratory tests of a GRAPE science model.

This paper discusses the latest progress in the development of GRAPE (Gamma-Ray Polarimeter Experiment), a hard X-ray Compton Polarimeter. The purpose of GRAPE is to measure the polarization of hard X-rays in the 50-300 keV energy range. We are particularly interested in X-rays that are emitted from solar flares and gamma-ray bursts (GRBs). Accurately measuring the polarization of the emitted radiation from these sources will lead, to a better understating of both the emission mechanisms and source geometries. The GRAPE design consists of an array of plastic scintillators surrounding a central high-Z crystal scintillator. We can monitor individual Compton scatters that occur in the plastics and determine whether the photon is photo absorbed by the high-Z crystal or not. A Compton scattered photon that is immediately photo absorbed by the high-Z crystal constitutes a valid event. These valid events provide us with the interaction locations of each incident photon and ultimately produces a modulation pattern for the Compton scattering of the polarized radiation. Comparing with Monte Carlo simulations of a 100% polarized beam, the level of polarization of the measured beam can then be determined. The complete array is mounted on a flat-panel multi-anode photomultiplier tube (MAPMT) that can measure the deposited energies resulting from the photon interactions. The design of the detector allows for a large field-of-view (>pi steradian), at the same time offering the ability to be close-packed with multiple modules in order to reduce deadspace. We plan to present in this paper the latest laboratory results obtained from GRAPE using partially polarized radiation sources.

The primary scientific mission of the Black Hole Finder Probe (BHFP), part of the NASA Beyond Einstein program, is to survey the local Universe for black holes over a wide range of mass and accretion rate. One approach to such a survey is a hard X-ray coded-aperture imaging mission operating in the 10--600 keV energy band, a spectral range that is considered to be especially useful in the detection of black hole sources. The development of new inorganic scintillator materials provides improved performance (for example, with regards to energy resolution and timing) that is well suited to the BHFP science requirements. Detection planes formed with these materials coupled with a new generation of readout devices represent a major advancement in the performance capabilities of scintillator-based gamma cameras. Here, we discuss the Coded Aperture Survey Telescope for Energetic Radiation (CASTER), a concept that represents a BHFP based on the use of the latest scintillator technology.

The Milagro gamma-ray observatory employs a water Cherenkov detector to observe extensive air showers produced by high energy particles interacting in the Earth's atmosphere. Milagro has a wide field of view and high duty cycle, monitoring the northern sky almost continuously in the 100 GeV to 100 TeV energy range. Milagro is, thus, uniquely capable of searching for very high-energy emission from gamma-ray bursts (GRBs) during the prompt emission phase. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Twenty-five satellite-triggered GRBs occurred within the field of view of Milagro between January 2000 and December 2001. We have searched for counterparts to these GRBs and found no significant emission from any of the burst positions. Due to the absorption of high-energy gamma rays by the extragalactic background light, detections are only expected to be possible for redshifts less than ~0.5. Three of the GRBs studied have measured redshifts. GRB 010921 has a redshift low enough (0.45) to allow an upper limit on the fluence to place an observational constraint on potential GRB models.

The Milagro telescope monitors the northern sky for 100 GeV to 100 TeV transient emission through continuous very high energy wide-field observations. The large effective area and ~100 GeV energy threshold of Milagro allow it to detect very high energy (VHE) gamma-ray burst emission with much higher sensitivity than previous instruments and a fluence sensitivity at VHE energies comparable to that of dedicated gamma-ray burst satellites at keV to MeV energies. Even in the absence of a positive detection, VHE observations can place important constraints on gamma-ray burst (GRB) progenitor and emission models. We present limits on the VHE flux of 40 s -- 3 h duration transients nearby to earth, as well as sensitivity distributions which have been corrected for gamma-ray absorption by extragalactic background light and cosmological effects. The sensitivity distributions suggest that the typical intrinsic VHE fluence of GRBs is similar or weaker than the keV -- MeV emission, and we demonstrate how these sensitivity distributions may be used to place observational constraints on the absolute VHE luminosity of gamma-ray bursts for any GRB emission and progenitor model.

We present a concept for an Advanced Compton Telescope (ACT) based on the use of pixelized gas micro-well detectors to form a three-dimensional electron track imager. A micro-well detector consists of an array of individual micro-patterned proportional counters opposite a planar drift electrode. When combined with thin film transistor array readouts, large gas volumes may be imaged with very good spatial and energy resolution at reasonable cost. The third dimension is determined from the drift time of the ionization electrons. The primary advantage of this approach is the excellent tracking of the Compton recoil electron that is possible in a gas volume. Such good electron tracking allows us to reduce the point spread function of a single incident photon dramatically, greatly improving the imaging capability and sensitivity. The polarization sensitivity, which relies on events with large Compton scattering angles, is particularly enhanced. We describe a possible ACT implementation of this technique, in which the gas tracking volume is surrounded by a CsI calorimeter, and present our plans to build and test a small prototype over the next three years.

Polarization measurements of hard X-rays from solar flares can provide a direct handle on the extent to which the accelerated electrons are beamed, which, in turn, has important implications for particle acceleration models. Past efforts to measured X-ray polarization have met with only limited success. Although designed primarily as a hard X-ray imager and spectrometer, the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is also capable of measuring the polarization of hard X-rays (20-100 keV) from solar flares. This capability arises from the inclusion of a small unobstructed Be scattering element that is strategically located within the cryostat that houses the array of nine germanium detectors. The Ge detectors are segmented, with both a front and rear active volume. Low energy photons (below about 100 keV) can reach a rear segment of a Ge detector only indirectly, by scattering. Low energy photons from the Sun have a direct path to the Be and have a high probability of Compton scattering into a rear segment of a Ge detector. The azimuthal distribution of these scattered photons carries with it a signature of the linear polarization of the incident flux. Sensitivity estimates, based on Monte Carlo simulations and in-flight background measurements, indicate that a 20-100 keV polarization sensitivity of less than a few percent can be achieved for X-class flares.

Supernova 1998bu in the galaxy M96 was observed by COMPTEL for a total of 88 days starting 17 days after the explosion. We searched for a signal in the 847 keV and 1238 keV lines of radioactive 56Co from this type Ia supernova. Using several different analysis methods, we did not detect SN1998bu. Our measurements should have been sensitive enough to detect 60Co gamma-rays as predicted from supernova models. Our 2-sigma flux limit is 2.3 10^{-5} photons cm^{-2} s^{-1}; this would correspond to 0.35 solar mass of ejected 56Ni, if SN1998bu were at a distance of 11.3 Mpc and transparent to MeV gamma rays for the period of our measurements. We discuss our measurements in the context of common supernova models, and conclude disfavoring a supernova event with large mixing and major parts of the freshly-generated radioactivity in outer layers.

Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.

We have used observations of Cygnus X-1 from the Compton Gamma-Ray Observatory (CGRO) and BeppoSAX to study the variation in the MeV gamma-ray emission between the hard and soft spectral states, using spectra that cover the energy range 20 keV up to 10 MeV. These data provide evidence for significant spectral variability at energies above 1 MeV. In particular, whereas the hard X-ray flux decreases during the soft state, the flux at energies above 1 MeV increases, resulting in a significantly harder gamma-ray spectrum at energies above 1 MeV. This behavior is consistent with the general picture of galactic black hole candidates having two distinct spectral forms at soft gamma-ray energies. These data extend this picture, for the first time, to energies above 1 MeV. We have used two different hybrid thermal/non-thermal Comptonization models to fit broad band spectral data obtained in both the hard and soft spectral states. These fits provide a quantitative estimate of the electron distribution and allow us to probe the physical changes that take place during transitions between the low and high X-ray states. We find that there is a significant increase (by a factor of ~4) in the bolometric luminosity as the source moves from the hard state to the soft state. Furthermore, the presence of a non-thermal tail in the Comptonizing electron distribution provides significant constraints on the magnetic field in the source region.

The development of CdZnTe detectors with an orthogpnal coplanar anode structure has important implications for future astrophysical instrumentation. As electron-only devices, like pixel detectors, coplanar anode strip detectors can be fabricated in the thickness required to be effective for photons with energies in excess of 500 keV. Unlike conventional double-sided strip detectors, the coplanar anode strip detectors require segmented contacts and signal processing electronics on only one surface. This facilitates the fabrication of closely-packed large area arrays that will be required for the next generation of coded aperture imagers. These detectors provide both very good energy resolution and sub-millimeter spatial resolution (in all three spatial dimensions) with far fewer electronic channels than are required for pixel detectors. Here we summarize results obtained from prototype detectors having a thickness of 5 mm and outline a concept for large area applications.

We have used observations from CGRO to study the variation in the MeV emission of Cygnus X-1 between its low and high X-ray states. These data provide a measurement of the spectral variability above 1 MeV. The high state MeV spectrum is found to be much harder than that of the low state MeV spectrum. In particular, the power-law emission seen at hard X-ray energies in the high state spectrum (with a photon spectral index of 2.6) is found to extend out to at least 5 MeV, with no evidence for any cutoff. Here we present the data and describe our efforts to model both the low state and high state spectra using a hybrid thermal/nonthermal model in which the emission results from the Comptonization of an electron population that consists of both a thermal and nonthermal component.

PKS 1622-297 was detected as a source of gamma-rays by the EGRET experiment at energies above 100 MeV during a gamma-ray outburst June and July 1995. We analyzed the COMPTEL data (0.75 - 30 MeV) of this time period to investigate the behaviour of PKS 1622-297 at lower gamma-ray energies. The blazar is significantly (5.7 sigma) detected by COMPTEL at energies above 10 MeV. Below 10 MeV the source is only marginally (3-10 MeV band) or not (below 3 MeV) detected. The summed MeV spectrum shows a 'hard' (alpha < 2, ~ E^(-alpha)) shape and, if combined with the simultaneous EGRET spectrum, a spectral break at MeV energies is indicated. We present the COMPTEL results (light curves, spectra) and compare them to results derived in neighbouring energy bands, in particular to the EGRET one.

We present the latest update of the 1.809 MeV sky survey obtained with COMPTEL. Based on all observations taken since the launch of CGRO in spring 1991 to early summer this year we obtain 1.809 MeV all sky maps using different imaging methods. The background is modelled on the basis of an adjacent energy approach. We confirm the previously reported characteristics of the galactic 1.809 MeV emission, specifically excesses in regions away from the inner Galaxy. The observed 1.8 MeV gamma-ray line is ascribed to the radioactive decay of 26Al in the interstellar medium. 26Al has been found to be predominantly synthesised in massive stars and their subsequent core-collapse supernovae, which is confirmed in tracer comparisons. Due to this, one anticipates flux enhancements aligned with regions of recent star formation, such as apparently observed in the Cygnus and Vela regions.

The type Ia supernova SN1998bu in M96 was observed by COMPTEL for a total of 88 days starting 17 days after the first detection of the supernova. The accumulated effective observation time was 4.14 Msec. The COMPTEL observations were performed in a special instrument mode improving the low-energy sensitivity. We generated images in the 847 keV and 1238 keV lines of 56Co, using improved point spread functions for the low-energy mode. We do not detect SN1998bu. A spectral analysis of our data also confirms the non-detection of the supernova. We discuss the event for which our upper limits constrain the standard supernova models.

The instrumental line background of the Compton telescope COMPTEL onboard the Compton Gamma-Ray Observatory is due to the activation and/or decay of many isotopes. The major components of this background can be attributed to eight individual isotopes, namely 2D, 22Na, 24Na, 28Al, 40K, 52Mn, 57Ni, and 208Tl. The identification of instrumental lines with specific isotopes is based on the line energies as well as on the variation of the event rate with time, cosmic-ray intensity, and deposited radiation dose during passages through the South-Atlantic Anomaly. The characteristic variation of the event rate due to a specific isotope depends on its life-time, orbital parameters such as the altitude of the satellite above Earth, and the solar cycle. A detailed understanding of the background contributions from instrumental lines is crucial at MeV energies for measuring the cosmic diffuse gamma-ray background and for observing gamma-ray line emission in the interstellar medium or from supernovae and their remnants. Procedures to determine the event rate from each background isotope are described, and their average activity in spacecraft materials over the first seven years of the mission is estimated.

The Compton Gamma-Ray Observatory (CGRO) has observed the Cygnus region on several occasions since its launch in 1991. The data collected by the COMPTEL experiment on CGRO represent the most sensitive observations to date of Cygnus X-1 in the 0.75-30 MeV range. A spectrum accumulated by COMPTEL over 10 weeks of observation time shows significant evidence for emission extending out to several MeV. We have combined these data with contemporaneous data from both BATSE and OSSE to produce a broad-band gamma-ray spectrum, corresponding to the low X-ray state of Cygnus X-1, extending from 50 keV up to approximately 5 MeV. Although there is no evidence for any broad line-like emissions in the MeV region, these data further confirm the presence of a hard tail at energies above several hundred keV. In particular, the spectrum at MeV energies can be described as a power-law with a photon spectral index of $\alpha$ = -3.2, with no evidence for a cutoff at high energies. For the 200 keV to 5 MeV spectrum, we provide a quantitative description of the underlying electron spectrum, in the context of a hybrid thermal/non-thermal model for the emission. The electron spectrum can be described by a thermal Maxwellian with a temperature of $kT_e$ = 86 keV and a non-thermal power-law component with a spectral index of $p_e$ = 4.5. The spectral data presented here should provide a useful basis for further theoretical modeling.

We report first observational results of multifrequency campaigns on the prominent Virgo blazars 3C 273 and 3C 279 which were carried out in January and February 1999. Both blazars are detected from radio to gamma-ray energies. We present the measured X- to gamma-ray spectra of both sources, and for 3C 279 we compare the 1999 broad-band (radio to gamma-ray) spectrum to measured previous ones.

In previous work, we have used data from the first three years of the CGRO mission to assemble a broad-band $\gamma$-ray spectrum of the galactic black hole candidate Cygnus X-1. Contemporaneous data from the COMPTEL, OSSE and BATSE experiments on CGRO were selected on the basis of the hard X-ray flux (45--140 keV) as measured by BATSE. This provided a spectrum of Cygnus X-1 in its canonical low X-ray state (as measured at energies below 10 keV), covering the energy range from 50 keV to 5 MeV. Here we report on a comparison of this spectrum to a COMPTEL-OSSE spectrum collected during a high X-ray state of Cygnus X-1 (May, 1996). These data provide evidence for significant spectral variability at energies above 1 MeV. In particular, whereas the hard X-ray flux {\it decreases} during the high X-ray state, the flux at energies above 1 MeV {\it increases}, resulting in a significantly harder high energy spectrum. This behavior is consistent with the general picture of galactic black hole candidates having two distinct spectral forms at soft $\gamma$-ray energies. These data extend this picture, for the first time, to energies above 1 MeV.

The Compton Gamma-Ray Observatory (CGRO) experiments EGRET and COMPTEL observed the Virgo sky region continuously for 7 weeks between December 10, 1996 and January 28, 1997. The prominent quasar 3C~273 was found to be the brightest source in gamma-rays and was significantly detected by EGRET and COMPTEL. The EGRET experiment observed a time-variable flux at energies above 100 MeV, which reached in a 2-week flaring period (December 30, 1996 to January 14, 1997) its highest flux level observed during the CGRO-era. COMPTEL, however, does not observe obvious time variability at energies below ~30 MeV contemporaneous to EGRET. In particular, no flare was observed, indicating that this outburst is solely a high-energy (>100 MeV) phenomenon. The energy spectrum between 3 MeV and 10 GeV is well represented by a simple power-law model. Below 3 MeV a spectral turnover is indicated. Performing spectral analysis for different time periods, we found evidence for a spectral hardening during the flaring period, which is consistent with the flare occurring mainly at the higher energies and with its absence at COMPTEL energies of a few MeV. This may be interpreted as an indication that the emission in the EGRET energy range is dominated by a different radiation mechanism than the MeV emission. We argue that the most likely mechanism for the high-energy flare is inverse-Compton scattering of reprocessed accretion-disk radiation.

Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.

We use all COMPTEL data from the beginning of the CGRO mission (April '91) up to the end of CGRO Cycle 6 (November '97) to carry out all-sky point source analyses in the four standard COMPTEL energy bands for different time periods. We apply our standard maximum-likelihood method to generate all-sky significance and flux maps for point sources by subtracting off the diffuse emission components via model fitting. In addition, fluxes of known sources have been determined for individual CGRO Phases/Cycles to generate lightcurves with a time resolution of the order of one year. The goal of the analysis is to derive quantitative results -- significances, fluxes, light curves -- of our brightest and most significant sources such as 3C 273, and to search for additional new COMPTEL sources, showing up in time-averaged maps only.

The COMPTEL experiment aboard the COMPTON Gamma-Ray Observatory (CGRO) has observed the prominent TeV-blazars Mkn 421 and Mkn 501 many times between the start of its mission in April '91 and December '98. This paper reports first COMPTEL results from mainly time-averaged (CGRO Cycles) data. No evidence for both sources is found up to the end Cycle VI. However, the sum of all 10-30 MeV Cycle VII data shows a weak (3.2 sigma detection) MeV-source being positionally consistent with Mkn 421. During Cycle VII Mkn 421 was rather active at TeV-energies. Due to the lack of other known Gamma-ray sources in this sky region, we consider Mkn 421 as the most likely counterpart for this Gamma-ray emission. However, its connection cannot be proven by COMPTEL.

Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.

GRB 990123 was the first burst from which simultaneous optical, X-ray and gamma-ray emission was detected; its afterglow has been followed by an extensive set of radio, optical and X-ray observations. We have studied the gamma-ray burst itself as observed by the CGRO detectors. We find that gamma-ray fluxes are not correlated with the simultaneous optical observations, and the gamma-ray spectra cannot be extrapolated simply to the optical fluxes. The burst is well fit by the standard four-parameter GRB function, with the exception that excess emission compared to this function is observed below ~15 keV during some time intervals. The burst is characterized by the typical hard-to-soft and hardness-intensity correlation spectral evolution patterns. The energy of the peak of the nu f_nu spectrum, E_p, reaches an unusually high value during the first intensity spike, 1470 +/- 110 keV, and then falls to \~300 keV during the tail of the burst. The high-energy spectrum above ~MeV is consistent with a power law with a photon index of about -3. By fluence, GRB 990123 is brighter than all but 0.4% of the GRBs observed with BATSE, clearly placing it on the -3/2 power-law portion of the intensity distribution. However, the redshift measured for the afterglow is inconsistent with the Euclidean interpretation of the -3/2 power-law. Using the redshift value of >= 1.61 and assuming isotropic emission, the gamma-ray fluence exceeds 10E54 ergs.

We present a new algorithm, called Multiresolution Regularized Expectation Maximization (MREM), for the reconstruction of gamma-ray intensity maps from COMPTEL data. The algorithm is based on the iterative Richardson-Lucy scheme to which we added a wavelet thresholding step in order to eliminate image-noise in the reconstruction. The wavelet thresholding explicitly accounts for spatial correlations in the data, and adapts the angular resolution locally, depending on the significance of the signal in the data. We compare the performance of MREM to that of the maximum entropy and the Richardson-Lucy algorithms by means of Monte-Carlo simulations of COMPTEL 1.809 MeV gamma-ray line observations. The simulations demonstrate that the maximum entropy and Richardson-Lucy algorithms provide virtually identical reconstructions which are heavily disturbed by image noise. MREM largely suppresses this noise in the reconstructions, showing only the significant structures that are present in the data. Application of MREM to COMPTEL 1.8 MeV gamma-ray line data results in a 1.809 MeV sky map that is much smoother than the maximum entropy or Richardson-Lucy reconstructions presented previously. The essential features of this map are (1) an asymmetric galactic ridge emission reaching from l=45 deg to l=240 deg, (2) a bright localised emission feature in the Cygnus region around (l,b)=(80 deg,0 deg), (3) two emission spots at l=317 deg and l=332 deg situated in the galactic plane, and (4) an extended emission region around (l,b)=(160 deg, 0 deg). Comparison of the MREM map to the simulated reconstructions demonstrates that the 1.809 MeV emission is confined to the galactic plane.

Deconvolving COMPTEL gamma-ray data into images presents a major methodological challenge. We developed a new algorithm called Multiresolution Regularized Expectation Maximization (MREM), which explicitly accounts for spatial correlations in the image by using wavelets (Dixon et al. 1998). We demonstrate that MREM largely suppresses image noise in the reconstruction, showing only significant structures that are present in the data. Application to 1.8 MeV data results in a sky map that is much smoother than the maximum entropy reconstructions presented previously, but it shows the same characteristic emission features which have been established earlier.

We examined the question of a possible correlation between 1.8 MeV emission and intensity distributions observed at other wavelengths by comparing the 1.8 MeV data to an extended set of all-sky maps, covering the entire explored wavelength range from the radio band up to high-energy gamma-rays. This analysis revealed that tracers of the old stellar population or the local interstellar medium provide only a poor description of the data. Tracers of the young stellar population considerably improve the fit. Residuals are minimized for the free-free emission map obtained by the DMR at microwave wavelengths and the 158 um fine-structure line map of C+ observed by FIRAS. Within the statistics of the present data, both maps provide an overall satisfactory fit. We therefore conclude that adequate tracers of 26Al sources have now been identified. Implications for 26Al source types are discussed (see also Knodlseder, these proceedings).

The COMPTEL experiment aboard CGRO, exploring the previously unknown sky at MeV-energies, has so far detected 10 Active Galactic Nuclei (AGN): 9 blazars and the radio galaxy Centaurus A. No Seyfert galaxy has been found yet. With these results COMPTEL has opened the field of extragalactic Gamma-ray astronomy in the MeV-band.

The COMPTEL experiment aboard the Compton Gamma-Ray Observatory (CGRO) has observed the gamma-ray blazar 3C 279 several times between April 1991 and September 1995. This paper reports on a consistent analysis of these observations using the most recent COMPTEL data analysis tools. Detections and non-detections of 3C 279 along the CGRO-mission indicate a time-variable MeV-flux. Spectral variability is indicated as well, however can not be significantly proven by spectral fitting. The average MeV-spectrum of 3C 279, as measured by COMPTEL over the four-year period, is consistent with a photon power-law slope of about -1.9. This spectrum smoothly connects to the simultaneous 30 MeV to 10 GeV spectrum obtained from an analysis of the combined EGRET observations. No spectral break is required to fit the time-averaged MeV- to GeV-spectrum of 3C 279.

The COMPTEL observations of the blazar-type quasar PKS 0528+134 in the energy range 0.75 MeV to 30 MeV carried out between April 1991 and September 1994 have been analyzed. During the first two years PKS 0528+134 was most significantly detected at energies above 3 MeV. During the last year there is only evidence for the quasar at energies below 3 MeV indicating a spectral change. The time-averaged COMPTEL energy spectrum between 0.75 MeV and 30 MeV is well represented by a power-law shape. Spectra collected from different observational periods reveal different power-law shapes: a hard state during flaring observations reported by EGRET, and a soft state otherwise. The combined simultaneous EGRET and COMPTEL spectra indicate these two spectral states as well. During low intensisty gamma-ray phases no spectral break is obvious from the combined COMPTEL and EGRET measurements. For the gamma-ray flaring phases however, the combined COMPTEL and EGRET data require a spectral bending at MeV-energies. By fitting broken power-law functions the best-fit values for the break in photon index range between 0.6 and 1.7, and for the break energy between ~5 MeV and ~20 MeV. Because the flux values measured by COMPTEL below 3 MeV in both states are roughly equal, the observations would be consistent with an additional spectral component showing up during gamma-ray flaring phases of PKS 0528+134. Such a component could be introduced by e.g. a high-energy electron-positron population with a low-energy cutoff in their bulk Lorentz factor distribution. The multiwavelength spectrum of PKS 0528+134 for gamma-ray flaring phases shows that the major energy release across the entire electro-magnetic spectrum is measured at MeV-energies.

We report on a study of the gamma-ray continuum emission from the Crab supernova nebula and on a search for nuclear de-excitation gamma-ray lines. Crab is the brightest continuum source in the 1-10 MeV gamma-ray sky, and its continuum radiation is most likely of synchrotron origin. It is a likely source of cosmic rays through shock acceleration and thus a potential candidate for gamma-ray line emission from nuclear interactions. Five years of COMPTEL observations enable a fine spectral binning to investigate the behaviour of the 0.75-30 MeV emission in detail and to search for nuclear de-excitation lines on top of the continuum. The nebular spectrum shows a break at the edge of the COMPTEL energy range and connects well to the EGRET spectrum, probably reflecting electron energy losses in the synchrotron emission scenario. Such a smooth continuum model alone may not be sufficient to explain the observations. A weak bump in the spectrum at 1-2 MeV may be present. No significant evidence for distinct line emission is seen, but the presence of a blend of line features or another synchrotron component cannot be excluded.

The COMPTEL instrument on the Compton Gamma Ray Observatory is used to measure the locations of gamma-ray bursts through direct imaging of MeV photons. In a comprehensive search, we have detected and localized 29 bursts observed between 1991 April 19 and 1995 May 31. The average location accuracy of these events is 1.25$\arcdeg$ (1$\sigma$), including a systematic error of $\sim0.5\arcdeg$, which is verified through comparison with Interplanetary Network (IPN) timing annuli. The combination of COMPTEL and IPN measurements results in locations for 26 of the bursts with an average ``error box'' area of only $\sim$0.3 deg$^2$ (1$\sigma$). We find that the angular distribution of COMPTEL burst locations is consistent with large-scale isotropy and that there is no statistically significant evidence of small-angle auto-correlations. We conclude that there is no compelling evidence for burst repetition since no more than two of the events (or $\sim$7% of the 29 bursts) could possibly have come from the same source. We also find that there is no significant correlation between the burst locations and either Abell clusters of galaxies or radio-quiet quasars. Agreement between individual COMPTEL locations and IPN annuli places a lower limit of $\sim$100~AU (95% confidence) on the distance to the stronger bursts.

We estimate the 1.8 MeV luminosity of the Sco-Cen association due to radioactive decay of 26Al to (4-15) 10e-5 ph cm**-2 s**-1. We propose a low surface brightness, limb brightened bubble for the 1.8 MeV intensity distribution. The detectibility of this distribution with existing gamma-ray telescopes is discussed.

COMPTEL imaging analysis revealed a patchy, asymmetric distribution of diffuse 1.8 MeV emission along the Galactic plane, which is attributed to the decay of radioactive 26Al in the ISM. If massive stars were the major source of Galactic 26Al, the 1.8 MeV emission should be asymmetric and trace the spiral arms of the Galaxy, presumed site of massive star formation. Using model fits, we indeed find weak evidence in the COMPTEL data that the observed 1.8 MeV emission is at least partly confined to spiral arms. We derive a total Galactic 26Al mass of 2.5 Msol from which at least 0.7 Msol can be attributed to massive stars.

Significant 1.8 MeV emission from the Carina region has been detected by COMPTEL. The emission is concentrated within 6 degrees or less near the Carina nebula NGC 3372, one of the brightest HII regions known in our Galaxy. This region contains a wealth of extreme young open clusters whose massive stars possibly contributed to an enrichment of 26Al in the ISM within the last few million years. The relation of these clusters and the peculiar object Eta Carinae with the observed emission is discussed. The 26Al yield of the clusters is estimated using current theoretical nucleosynthesis models.

We have analyzed the Ulysses, BATSE, and COMPTEL spectral data from the $\gamma$-ray burst of June 22, 1992 (GB 920622). COMPTEL data reveal a hard to soft evolution within the first pulse of the burst, while the mean hardness ratios of the three pulses are the same. Unlike the single instrument spectra, the composite spectrum of GB 920622 averaged over the total burst duration ranging from 20 keV up to 10 MeV cannot be fit by a single power law. Instead, the spectrum shows continuous curvature across the full energy range. COMPTEL imaging and BATSE/Ulysses triangulation constrain the source location of GB 920622 to a ring sector 1.1 arcmin wide and 2 degrees long. This area has been searched for quiescent X-ray sources using \ros survey data collected about two years before the burst. After the optical identification of the X-ray sources in and near the GRB location we conclude that no quiescent X-ray counterpart candidate for GB 920622 has been found.