Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Several neutrino observatories have searched for coincident neutrino signals associated with gravitational waves induced by the merging of two black holes.No statistically significant neutrino signal in excess of background level was observed.These experiments use different neutrino detection technologies and are sensitive to various neutrino types.A combined analysis was performed on the KamLAND, Super-Kamiokande and Borexino experimental data with a frequentist statistical approach to achieve a global picture of the associated neutrino fluence.Both monochromatic and Fermi-Dirac neutrino spectra were assumed in the calculation.The final results are consistent with null neutrino signals associated with a binary of black holes merging process.The derived 90\% confidence level upper limits on the fluence and luminosity of various neutrino types are presented for neutrino energy less than 110 MeV.

A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $\alpha_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.).

A combined BATSE-GBM hard X-ray catalog is presented based on Earth Occultation Imaging applied to a reanalysis of BATSE data. An imaging approach has been developed for the reanalysis of Earth Occultation analysis of BATSE data. The standard occultation analysis depends on a predetermined catalog of potential sources, so that a real source not present in the catalog may induce systematic errors when source counts associated with an uncatalogued source are incorrectly attributed to catalog sources. The goal of the imaging analysis is to find a complete set of hard X-ray sources, including sources not in the original BATSE occultation catalog. Using the imaging technique, we have identified 15 known sources and 17 unidentified sources and added them to the BATSE occultation catalog. The resulting expanded BATSE catalog of sources observed during 1991-2000 is compared to the ongoing GBM survey.

We describe an improved method of mapping the gamma-ray sky by applying the Linear Radon Transform to data from BATSE on NASA's CGRO. Based on a method similar to that used in medical imaging, we use the relatively sharp (~0.25 deg) limb of the Earth to collimate BATSE's eight Large Area Detectors (LADs). Coupling this to the ~51-day precession cycle of the CGRO orbit, we can complete a full survey of the sky, localizing point sources to < 1 deg accuracy. This technique also uses a physical model for removing many sources of gamma-ray background, which allows us to image strong gamma-ray sources such as the Crab up to ~2 MeV with only a single precession cycle. We present the concept of the Radon Transform technique as applied to the BATSE data for imaging the gamma-ray sky and show sample images in three broad energy bands (23-98 keV, 98-230 keV, and 230-595 keV) centered on the positions of selected sources from the catalog of 130 known sources used in our Enhanced BATSE Occultation Package (EBOP) analysis system. Any new sources discovered during the sky survey will be added to the input catalog for EBOP allowing daily light curves and spectra to be generated. We also discuss the adaptation of tomographic imaging to the Fermi GBM occultation project.

We report BATSE earth-occultation observations of the long-term soft gamma-ray (30 keV - 1.7 MeV) flux and spectral variability of the transient source GRO J1719-24 covering a 1000-day period between 13 January 1993 and 10 October 1995. During this period, the source underwent a major outburst in the fall of 1993 when the 35-100 keV flux rose from a quiescent state of less than 16 mCrab before 17 September 1993 to a level of 1.5 Crab on 3 October and remained approximately at this level over the next ~70 days. It then decreased sharply to the pre-transition quiescent level of ~44 mCrab on 21 December where it remained until 5 September 1994. During a 400-day period between 5 September 1994 and 10 October 1995, the source again underwent a series of five transitions when the 35-100 keV flux increased to low-intensity levels of \~200-400 mCrab. The low-state spectra were generally characterized by a power law with spectral index of ~2. The high-state spectra on the other hand have two components: a thermal Comptonized shape below ~200 keV with electron temperature kTe of ~37 keV and optical depth of ~2.8, and a soft power-law tail with photon index of ~3.4 above 200 keV that extends to ~500 keV. The softer high-intensity spectrum and the harder low-intensity spectrum intercept at ~400 keV. The non-thermal power-law gamma-ray component seen in both the high and low-intensity spectra suggests that the persistent non-thermal emission source is coupled to the hot and variable thermal emission source in the system. This is similar to that seen in two other gamma-ray emitting black-hole candidates GRO J0422+32 and Cygnus X-1, suggesting that perhaps similar system configurations and processes are occurring in these systems.

Cygnus X-1 and the gamma-ray transients GROJ0422+32 and GROJ1719-24 displayed similar spectral properties when they underwent transitions between the high and low gamma-ray (30 keV to few MeV) intensity states. When these sources were in the high gamma-ray intensity state (gamma-2, for Cygnus X-1), their spectra featured two components: a Comptonized shape below 200-300 keV with a soft power-law tail (photon index >3) that extended to ~1 MeV or beyond. When the sources were in the low-intensity state (gamma-0, for Cygnus X-1), the Comptonized spectral shape below 200 keV typically vanished and the entire spectrum from 30 keV to ~1 MeV can be characterized by a single power law with a relatively harder photon index ~2-2.7. Consequently the high- and low-intensity gamma-ray spectra intersect, generally in the ~400 keV - 1 MeV range, in contrast to the spectral pivoting seen previously at lower (~10 keV) energies. The presence of the power-law component in both the high- and low-intensity gamma-ray spectra strongly suggests that the non-thermal process is likely to be at work in both the high and the low-intensity situations. We have suggested a possible scenario (Ling & Wheaton, 2003), by combining the ADAF model of Esin et al. (1998) with a separate jet region that produces the non-thermal gamma-ray emission, and which explains the state transitions. Such a scenario will be discussed in the context of the observational evidence, summarized above, from the database produced by EBOP, JPL's BATSE earth occultation analysis system.

The BATSE earth-occultation database provides nine years of coverage for 75 $\gamma$-ray sources in the energy range 35-1700 keV. For transient sources, this long time-base dataset makes it possible to study the repeated outbursts from individual objects. We have used the JPL data analysis package EBOP (Enhanced BATSE Occultation Package) to derive the light curves and the time evolution of the spectra for the black hole candidate and microquasar sources GRO J1655-40 and GRS 1915+105. We find that GRO J1655-40, during high-intensity flaring periods, is characterized by a single power-law spectrum up to 500 keV with a spectral index consistent with that observed by OSSE. During one flare observed contemporaneously with OSSE and HEXTE, the GRO J1655-40 spectrum was observed to steepen as the $\gamma$-ray intensity increased. For GRS 1915+105, the spectrum during high intensity flaring periods can be characterized by a broken power law with a time-varying high-energy component. The spectra of these microquasars differ from the black hole candidates Cygnus X-1, GRO J0422+32, and GRO J1719-24, which have thermal contributions to their spectra when in high $\gamma$-ray states. This suggests that there may be two different classes of Galactic black hole candidates.

We report ~ 600 days of BATSE earth-occultation observations of the total gamma-ray (30 keV to 1.7 MeV) emission from the Crab nebula, between 1991 May 24 (TJD 8400) and 1994 October 2 (TJD 9627). Lightcurves from 35-100, 100-200, 200-300, 300-400, 400-700, and 700-1000 keV, show that positive fluxes were detected by BATSE in each of these six energy bands at significances of approximately 31, 20, 9.2, 4.5, 2.6, and 1.3 sigma respectively per day. We also observed significant flux and spectral variations in the 35-300 keV energy region, with time scales of days to weeks. The spectra below 300 keV, averaged over typical CGRO viewing periods of 6-13 days, can be well described by a broken power law with average indices of ~ 2.1 and ~ 2.4 varying around a spectral break at ~ 100 keV. Above 300 keV, the long-term averaged spectra, averaged over three 400 d periods (TJD 8400-8800, 8800-9200, and 9200-9628, respectively) are well represented by the same power law with index of ~ 2.34 up to ~ 670 keV, plus a hard spectral component extending from ~ 670 keV to ~ 1.7 MeV, with a spectral index of ~ 1.75. The latter component could be related to a complex structure observed by COMPTEL in the 0.7-3 MeV range. Above 3 MeV, the extrapolation of the power-law continuum determined by the low-energy BATSE spectrum is consistent with fluxes measured by COMPTEL in the 3-25 MeV range, and by EGRET from 30-50 MeV. We interpret these results as synchrotron emission produced by the interaction of particles ejected from the pulsar with the field in different dynamical regions of the nebula system, as observed recently by HST, XMM-Newton, and Chandra.

We report results of a comprehensive study of the soft gamma-ray (30 keV to 1.7 MeV) emission of GROJ0422+32 during its first known outburst in 1992. These results were derived from the BATSE earth-occultation database with the JPL data analysis package, EBOP (Enhanced BATSE Occultation Package). Results presented here focus primarily on the long-term temporal and spectral variability of the source emission associated with the outburst. The light curves with 1-day resolution in six broad energy-bands show the high-energy flux (>200 keV) led the low-energy flux (<200 keV) by ~5 days in reaching the primary peak, but lagged the latter by ~7 days in starting the declining phase. We confirm the "secondary maximum" of the low-energy (<200 keV) flux at TJD 8970-8981, ~120 days after the first maximum. Our data show that the "secondary maximum" was also prominent in the 200-300 keV band, but became less pronounced at higher energies. During this 200-day period, the spectrum evolved from a power-law with photon index of 1.75 on TJD 8839, to a shape that can be described by a Comptonized model or an exponential power law below 300 keV, with a variable power-law tail above 300 keV. The spectrum remained roughly in this two-component shape until ~9 November (TJD 8935) and then returned to the initial power-law shape with an index of ~2 until the end of the period. The correlation of the two spectral shapes with the high and low luminosities of the soft gamma-ray emission is strongly reminiscent of that seen in Cygnus X-1. We interpret these results in terms of the Advection Dominated Accretion Flow (ADAF) model with possibly a "jet-like" region that persistently produced the non-thermal power-law gamma rays observed throughout the event.

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.