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.

All three components of the current density are required to compute the heating rate due to free magnetic energy dissipation. Here we present a first test of a new model developed to determine if the times of increases in the resistive heating rate in active region (AR) photospheres are correlated with the subsequent occurrence of M and X flares in the corona. A data driven, 3 D, non-force-free magnetohydrodynamic model restricted to the near-photospheric region is used to compute time series of the complete current density and the resistive heating rate per unit volume $(Q(t))$ in each pixel in neutral line regions (NLRs) of 14 ARs. The model is driven by time series of the magnetic field ${\bf B}$ measured by the Helioseismic & Magnetic Imager on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series for ${\bf B}$ in every AR pixel. For each AR, the cumulative distribution function (CDF) of the values of the NLR area integral $Q_i(t)$ of $Q(t)$ is found to be a scale invariant power law distribution essentially identical to the observed CDF for the total energy released in coronal flares. This suggests that coronal flares and the photospheric $Q_i$ are correlated, and powered by the same process. The model predicts spikes in $Q_i$ with values orders of magnitude above background values. These spikes are driven by spikes in the non-force free component of the current density. The times of these spikes are plausibly correlated with times of subsequent M or X flares a few hours to a few days later. The spikes occur on granulation scales, and may be signatures of heating in horizontal current sheets. It is also found that the times of relatively large values of the rate of change of the NLR unsigned magnetic flux are also plausibly correlated with the times of subsequent M and X flares, and spikes in $Q_i$.

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.

The origin of ultra-compact dwarf galaxies (UCDs) is not yet clear. One possible formation path of UCDs is the threshing of a nucleated elliptical dwarf galaxy (dE, N), however, it remains unclear how such massive nuclear stellar systems were formed in dwarf galaxies. To better establish the early history of UCDs, we investigate the formation of UCD progenitor clusters from super giant molecular clouds (SGMCs), using hydrodynamical simulations. In this study we focus on SGMCs with masses $10^{7} - 10^{8} \rm M_{\odot}$ that can form massive star clusters that display physical properties similar to UCDs. We find that the clusters have extended star formation histories with two phases, producing multiple distinct stellar populations, and that the star formation rate is dependent on the feedback effects of SNe and AGB stars. The later generations of stars formed in these clusters are more compact, leading to a clearly nested structure, and these stars will be more He-rich than those of the first generation, leading to a slight colour gradient. The simulated clusters demonstrate scaling relations between $R_{\rm eff}$ and $M$ and $\sigma_{v}$ and $M$ consistent with those observed in UCDs and strongly consistent with those of the original SGMC. We discuss whether SGMCs such as these can be formed through merging of self-gravitating molecular clouds in galaxies at high-z.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

A magnetohydrodynamic model that includes a complete electrical conductivity tensor is used to estimate conditions for photospherically driven, linear, non-plane Alfvenic oscillations extending from the photosphere to the lower corona to drive a chromospheric heating rate due to Pedersen current dissipation that is comparable to the net chromospheric net radiative loss of $\sim 10^7$ ergs-cm$^{-2}$-sec$^{-1}$. The heating rates due to electron current dissipation in the photosphere and corona are also computed. The wave amplitudes are computed self-consistently as functions of an inhomogeneous background (BG) atmosphere. The effects of the conductivity tensor are resolved numerically using a resolution of 3.33 m. The oscillations drive a chromospheric heating flux $F_{Ch} \sim 10^7 - 10^8$ ergs-cm$^{-2}$-sec$^{-1}$ at frequencies $\nu \sim 10^2 - 10^3$ mHz for BG magnetic field strengths $B \gtrsim 700$ G and magnetic field perturbation amplitudes $\sim 0.01 - 0.1$ $B$. The total resistive heating flux increases with $\nu$. Most heating occurs in the photosphere. Thermalization of Poynting flux in the photosphere due to electron current dissipation regulates the Poynting flux into the chromosphere, limiting $F_{Ch}$. $F_{Ch}$ initially increases with $\nu$, reaches a maximum, and then decreases with increasing $\nu$ due to increasing electron current dissipation in the photosphere. The resolution needed to resolve the oscillations increases from $\sim 10$ m in the photosphere to $\sim 10$ km in the upper chromosphere, and is proportional to $\nu^{-1/2}$. Estimates suggest that these oscillations are normal modes of photospheric flux tubes with diameters $\sim 10-20$ km, excited by magnetic reconnection in current sheets with thicknesses $\sim 0.1$ km.

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.).

An MHD model of a Hydrogen plasma with flow, an energy equation, NLTE ionization and radiative cooling, and an Ohm's law with anisotropic electrical conduction and thermoelectric effects is used to self-consistently generate atmospheric layers over a $50$ km height range. A subset of these solutions contain current sheets, and have properties similar to those of the lower and middle chromosphere. The magnetic field profiles are found to be close to Harris sheet profiles, with maximum field strengths $\sim 25-150$ G. The radiative flux $F_R$ emitted by individual sheets is $\sim 4.9 \times 10^5 - 4.5 \times 10^6$ ergs-cm$^{-2}$-s$^{-1}$, to be compared with the observed chromospheric emission rate of $\sim 10^7$ ergs-cm$^{-2}$-s$^{-1}$. Essentially all emission is from regions with thicknesses $\sim 0.5 - 13$ km containing the neutral sheet. About half of $F_R$ comes from sub-regions with thicknesses 10 times smaller. A resolution $\lesssim 5-130$ m is needed to resolve the properties of the sheets. The sheets have total H densities $\sim 10^{13}-10^{15}$ cm$^{-3}$. The ionization fraction in the sheets is $\sim 2-20$ times larger, and the temperature is $\sim 2000-3000$ K higher than in the surrounding plasma. The Joule heating flux $F_J$ exceeds $F_R$ by $\sim 4-34 \%$, the difference being balanced in the energy equation mainly by a negative compressive heating flux. Proton Pedersen current dissipation generates $\sim 62-77\%$ of the positive contribution to $F_J$. The remainder of this contribution is due to electron current dissipation near the neutral sheet where the plasma is weakly magnetized.

A magnetohydrodynamic model is used to determine conditions under which the Lorentz force accelerates plasma to type 2 spicule speeds in the chromosphere. The model generalizes a previous model to include a more realistic pre-spicule state, and the vertical viscous force. Two cases of acceleration under upper chromospheric conditions are considered. The magnetic field strength for these cases is <= 12.5 and 25 G. Plasma is accelerated to terminal vertical speeds of 66 and 78 km/s in 100 s, compared with 124 and 397 km/s for the case of zero viscosity. The flows are localized within horizontal diameters ~ 80 and 50 km. The total thermal energy generated by viscous dissipation is ~ 10 times larger than that due to Joule dissipation, but the magnitude of the total cooling due to rarefaction is >~ this energy. Compressive heating dominates during the early phase of acceleration. The maximum energy injected into the corona by type 2 spicules, defined as the energy flux in the upper chromosphere, may largely balance total coronal energy losses in quiet regions, possibly also in coronal holes, but not in active regions. It is proposed that magnetic flux emergence in inter-granular regions drives type 2 spicules.

Observed since the 1950's, the seasonal effect on underground muons is a well studied phenomenon. The interaction height of incident cosmic rays changes as the temperature of the atmosphere changes, which affects the production height of mesons (mostly pions and kaons). The decay of these mesons produces muons that can be detected underground. The production of muons is dominated by pion decay, and previous work did not include the effect of kaons. In this work, the methods of Barrett and MACRO are extended to include the effect of kaons. These efforts give rise to a new method to measure the atmospheric K/$\pi$ ratio at energies beyond the reach of current fixed target experiments. These methods were applied to data from the MINOS far detector. A method is developed for making these measurements at other underground detectors, including OPERA, Super-K, IceCube, Baksan and the MINOS near detector.

MINOS is the first large magnetic detector deep underground and is the first to measure the muon charge ratio with high statistics in the region near 1 TeV.\cite{bib:adamson} An approximate formula for the muon charge ratio can be expressed in terms of $\epsilon_\pi$ = 115 GeV, $\epsilon_K$ = 850 GeV and $\ec$. The implications for K production in the atmosphere will be discussed.

Although formal proposals have not yet been made, the UNO and Hyper-Kamiokande projects are being developed to follow-up the tremendously successful program at Super-Kamiokande using a detector that is 20-50 times larger. The potential of such a detector to continue the study of astrophysical neutrinos is considered and contrasted with the program for cubic kilometer neutrino observatories.