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.

The events of multiple neutron production under 2000g/cm$^2$ thick rock absorber were studied at the Tien~Shan mountain cosmic ray station, at the altitude of 3340m above the sea level. From comparison of the experimental and Geant4 simulated neutron multiplicity spectra it follows that the great bulk of these events can be explained by interaction of cosmic ray muons with internal material of the neutron detector. In synchronous operation of the underground neutron monitor with the Tien~Shan shower detector system it was found that the characteristics of the muonic component of extensive air showers which is seemingly responsible for generation of the neutron events underground do change noticeably within the energy range of the knee of primary cosmic ray spectrum. Some peculiar shower events were detected when the neutron signal reveals itself only $\sim$(100--1000)\,$\mu$s after the passage of the shower particles front which probably means an existence of corresponding delay of the muon flux in such events.

We present a description of the new complex installation for the study of extensive air showers which was created at the Tien Shan mountain cosmic ray station, as well as the results of the test measurements made there in 2014-2016. At present, the system for registration of electromagnetic shower component consists of $\sim$100 detector points built on the basis of plastic scintillator plates with the sensitive area of 0.25m$^2$ and 1m$^2$, spread equidistantly over $\sim$10$^4$m$^2$ space. The dynamic range of scintillation amplitude measurements is currently about $(3-7)\cdot 10^4$, and there is a prospect of it being extended up to $\sim$10$^6$. The direction of shower arrival is defined by signal delays from a number of the scintillators placed cross-wise at the periphery of the detector system. For the investigation of nuclear active shower components there was created a multi-tier 55m$^2$ ionization-neutron calorimeter with a sum absorber thickness of $\sim$1000g/cm$^2$, typical spatial resolution of the order of 10cm, and dynamic range of ionization measurement channel about $\sim$10$^5$. Also, the use of saturation-free neutron detectors is anticipated for registration of the high- and low-energy hadron components in the region of shower core. A complex of underground detectors is designed for the study of muonic and penetrative nuclear-active components of the shower. The full stack of data acquisition, detector calibration, and shower parameters restoration procedures are now completed, and the newly obtained shower size spectrum and lateral distribution of shower particles occur in agreement with conventional data. Future studies in the field of $10^{14}-10^{17}$eV cosmic ray physics to be held at the new shower installation are discussed.

Purposeful investigation of radiation fluxes strongly delayed in relation to the main particles front of extensive air shower (EAS) was undertaken at the Tien Shan Mountain Cosmic Ray Station. It was found that the passage of the EAS can be accompanied by the delayed thermal neutrons and by the soft $(30-50)$ keV gamma rays, mostly concentrated within a region of about $(5-10)$ m around shower axis, where the integral radiation fluence can vary in the limits of $(10^{-4}-1)$ cm$^{-2}$ for neutrons, and of $(0.1-1000)$ cm$^{-2}$ for gamma rays. The dependence of signal multiplicity on the shower size $N_e$ has a power shape both for the neutron and gamma ray components, with a sharp increase of its power index around the value of $N_e\approx 10^6$, which corresponds to the position of the $3\cdot10^{15}$ eV knee in the primary cosmic ray spectrum. Total duration of detectable radiation signal after the EAS passage can be of some tens of milliseconds in the case of neutron component, and up to a few whole seconds for gamma rays. The delayed accompaniment of low-energy radiation particles can be an effective probe to study the interaction of the hadronic component of EAS.

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.

By analysing a database of 26 years of observations of Jupiter with the Nancay Decameter Array, we unambiguously identify the radio emissions caused by the Ganymede-Jupiter interaction. We study the energetics of these emissions via the distributions of their intensities, duration, and power, and compare them to the energetics of the Io-Jupiter radio emissions. This allows us to demonstrate that the average emitted radio power is proportional to the Poynting flux from the rotating Jupiter's magnetosphere intercepted by the obstacle. We then generalize this result to the radio-magnetic scaling law that appears to apply to all plasma interactions between a magnetized flow and an obstacle, magnetized or not. Extrapolating this scaling law to the parameter range corresponding to hot Jupiters, we predict large radio powers emitted by these objects, that should result in detectable radio flux with new-generation radiotelescopes. Comparing the distributions of the durations of Ganymede-Jupiter and Io-Jupiter emission events also suggests that while the latter results from quasi-permanent Alfven wave excitation by Io, the former likely results from sporadic reconnection between magnetic fields Ganymede and Jupiter, controlled by Jupiter's magnetic field geometry and modulated by its rotation.

Analyzing a database of 26 years of observations of Jupiter from the Nan\c{c}ay Decameter Array, we study the occurrence of Io-independent emissions as a function of the orbital phase of the other Galilean satellites and Amalthea. We identify unambiguously the emissions induced by Ganymede and characterize their intervals of occurrence in CML and Ganymede phase and longitude. We also find hints of emissions induced by Europa and, surprisingly, by Amalthea. The signature of Callisto-induced emissions is more tenuous.

This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. This paper describes digital radio astronomical receivers used for decameter and meter wavelength observations. Since 1998, digital receivers performing on-the-fly dynamic spectrum calculations or waveform data recording without data loss have been used at the UTR-2 radio telescope, the URAN VLBI system, and the GURT new generation radio telescope. Here we detail these receivers developed for operation in the strong interference environment that prevails in the decameter wavelength range. Data collected with these receivers allowed us to discover numerous radio astronomical objects and phenomena at low frequencies, a summary of which is also presented.