Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

We study the solar eruptive event on 2017 September 10 that produced long-lasting $>$100 MeV $\gamma$-ray emission and a ground level enhancement (GLE72). The origin of the high-energy ions producing late-phase gamma-ray emission (LPGRE) is still an open question, but a possible explanation is proton acceleration at coronal shocks produced by coronal mass ejections. We examine a common shock acceleration origin for both the LPGRE and GLE72. The $\gamma$-ray emission observed by the Fermi-Large Area Telescope exhibits a weak impulsive phase, consistent with that observed in hard X-and $\gamma$-ray line flare emissions, and what appear to be two distinct stages of LPGRE. From a detailed modeling of the shock wave, we derive the 3D distribution and temporal evolution of the shock parameters, and we examine the shock wave magnetic connection with the visible solar disk. The evolution of shock parameters on field lines returning to the visible disk, mirrors the two stages of LPGRE. We find good agreement between the time history of $>$100 MeV $\gamma$-rays and one produced by a basic shock acceleration model. The time history of shock parameters magnetically mapped to Earth agrees with the rates observed by the Fort Smith neutron monitor during the first hour of the GLE72 if we include a 30% contribution of flare-accelerated protons during the first 10 minutes, having a release time following the time history of nuclear $\gamma$-rays. Our analysis provides compelling evidence for a common shock origin for protons producing the LPGRE and most of the particles observed in GLE72.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

We characterize and provide a catalog of thirty >100 MeV sustained gamma-ray emission (SGRE) events observed by Fermi LAT. These events are temporally and spectrally distinct from the associated solar flares. Their spectra are consistent with decay of pions produced by >300 MeV protons and are not consistent with electron bremsstrahlung. SGRE start times range from CME onset to two hours later. Their durations range from about four minutes to twenty hours and appear to be correlated with durations of >100 MeV SEP proton events. The >300 MeV protons producing SGRE have spectra that can be fit with power laws with a mean index of ~4 and RMS spread of 1.8. Gamma-ray line measurements indicate that SGRE proton spectra are steeper above 300 MeV than they are below 300 MeV. The number of SGRE protons >500 MeV is on average about ten times more than then the number in the associated flare and about fifty to one hundred times less than the number in the accompanying SEP. SGRE can extend tens of degrees from the are site. Sustained bremsstrahlung from MeV electrons was observed in one SGRE event. Flare >100 keV X-ray emission appears to be associated with SGRE and with intense SEPs. From this observation, we provide arguments that lead us to propose that sub-MeV to MeV protons escaping from the flare contribute to the seed population that is accelerated by shocks onto open field lines to produce SEPs and onto field lines returning to the Sun to produce SGRE.

We report abundances of elements from $_{26}$Fe to $_{40}$Zr in the cosmic radiation measured by the SuperTIGER (Trans-Iron Galactic Element Recorder) instrument during 55 days of exposure on a long-duration balloon flight over Antarctica. These observations resolve elemental abundances in this charge range with single-element resolution and good statistics. These results support a model of cosmic-ray origin in which the source material consists of a mixture of 19$^{+11}_{-6}$\% material from massive stars and $\sim$81\% normal interstellar medium (ISM) material with solar system abundances. The results also show a preferential acceleration of refractory elements (found in interstellar dust grains) by a factor of $\sim$4 over volatile elements (found in interstellar gas) ordered by atomic mass (A). Both the refractory and volatile elements show a mass-dependent enhancement with similar slopes.

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.

Secondary nuclear production physics is receiving increased attention given the high-quality measurements of the gamma-ray emissivity of local interstellar gas between ~50 MeV and ~40 GeV, obtained with the Large Area Telescope on board the Fermi space observatory. More than 90% of the gas-related emissivity above 1 GeV is attributed to gamma-rays from the decay of neutral pions formed in collisions between cosmic rays and interstellar matter, with lepton-induced processes becoming increasingly important below 1 GeV. The elementary kinematics of neutral pion production and decay are re-examined in light of two physics questions: does isobaric production follow a scaling behavior? and what is the minimum proton kinetic energy needed to make a gamma-ray of a certain energy formed through intermediate pi0 production? The emissivity spectrum will allow the interstellar cosmic-ray spectrum to be determined reliably, providing a reference for origin and propagation studies as well as input to solar modulation models. A method for such an analysis and illustrative results are presented.

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.

We are in the process of obtaining VLBA polarisation data for the 135 MOJAVE-I Active Galactic Nuclei at four wavelengths in the 18-22cm band. These observations will enable studies of the evolution of the intensity and magnetic-field structures of these AGN jets as they propagate from parsec to kiloparsec scales, as well as studies of the thermal plasma present in the vicinity of the jets on these scales, manifest via Faraday rotation. A wide range of other multi-wavelength studies can also be carried out using these data. Preliminary results for selected sources from the first 3 of 9 observing sessions will be presented. We aim to have 18-cm intensity and polarisation images available via the MOJAVE website within 18 months after the last observing session.

We have studied the radioactive line emission expected from solar active regions after large flares, following the production of long-lived radioisotopes by nuclear interactions of flare-accelerated ions. This delayed X- and gamma-ray line emission can provide unique information on the accelerated particle composition and energy spectrum, as well as on mixing processes in the solar atmosphere. Total cross sections for the formation of the main radioisotopes by proton, 3-He and alpha-particle reactions are evaluated from available data combined with nuclear reaction theory. Thick-target radioisotope yields are provided in tabular form, which can be used to predict fluxes of all of the major delayed lines at any time after a gamma-ray flare. The brightest delayed line for days after the flare is found to be the 511 keV positron-electron annihilation line resulting from the decay of several beta+ radioisotopes. After ~2 days however, the flux of the 511 keV line can become lower than that of the 846.8 keV line from the decay of 56-Co. Our study has revealed other delayed gamma-ray lines that appear to be promising for detection, e.g. at 1434 keV from the radioactivity of both the isomer and the ground state of 52-Mn, 1332 and 1792 keV from 60-Cu, and 931.1 keV from 55-Co. The strongest delayed X-ray line is found to be the Co Kalpha at 6.92 keV, which is produced from both the decay of the isomer of 58-Co by the conversion of a K-shell electron and the decay of 57-Ni by orbital electron capture. Prospects for observation of these lines with RHESSI or future space instruments are discussed.

The Reuven Ramaty High Energy Solar Spectroscopy Imager (RHESSI) has obtained the first high-resolution measurements of nuclear de-excitation lines produced by energetic ions accelerated in a solar flare, a GOES X4.8 event occurring on 23 July, 2002 at a heliocentric angle of 73 degrees. Lines of neon, magnesium, silicon, iron, carbon, and oxygen were resolved for the first time. They exhibit Doppler redshifts of 0.1--0.8% and broadening of 0.1--2.1% (FWHM), generally decreasing with mass. The measured redshifts are larger than expected for a model of an interacting ion distribution isotropic in the downward hemisphere in a radial magnetic field. Possible interpretations of the large redshifts include 1) an inclination of the loop magnetic field to the solar surface so that the ion distribution is oriented more directly away from the observer, and 2) extreme beaming of the ions downward along a magnetic field normal to the solar surface. Bulk downward motion of the plasma in which the accelerated ions interact can be ruled out.

The energies and widths of gamma-ray lines emitted by ambient nuclei excited by flare-accelerated protons and alpha particles provide information on the ions directionality and spectra, and on the characteristics of the interaction region. We have measured the energies and widths of strong lines from de-excitations of 12C, 16O, and 20Ne in solar flares as a function of heliocentric angle. The line energies from all three nuclei exhibit ~1% redshifts for flares at small heliocentric angles, but are not shifted near the limb. The lines have widths of ~3% FWHM. We compare the 12C line measurements for flares at five different heliocentric angles with calculations for different interacting-particle distributions. A downward isotropic distribution (or one with a small upward component) provides a good fit to the line measurements. An angular distribution derived for particles that undergo significant pitch angle scattering by MHD turbulence in coronal magnetic loops provides comparably good fits.

Observations of annihilation radiation in the Galaxy are briefly reviewed. We summarize astrophysical mechanisms leading to positron production, and recent estimates for production rates from nova and supernova nucleosynthesis in the Galaxy. The physical processes involved in the production of annihilation radiation in the interstellar medium are described. These include positron thermalization, charge exchange, radiative recombination, and direct annihilation. Calculations of 2\gamma and 3\gamma spectra and the positronium (Ps) fraction due to the annihilation of positrons in media containing H and He at different temperatures and ionization states are presented. Quenching of Ps by high temperature plasmas or dust could account for differences between 0.511 MeV and 3\gamma Ps continuum maps. These results are presented in the context of the potential of INTEGRAL to map sites of annihilation radiation in the Galaxy. Positron production by compact objects is also considered.

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.