The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning (ML) formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two mono-energetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results intimate the promise of this \lq\lq variational annealing\rq\rq\ methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

We calculate the rate of production of dark photons from electron-positron pair annihilation in hot and dense matter characteristic of supernova progenitors. Given the non-linear dependence of the emission rate on the dark photon mass and current astrophysical constraints on the dark photon parameter space, we focus on the mass range of 1--10 MeV. For the conditions under consideration both mixing with the in-medium photon and plasma effects on the electron dispersion relation are non-negligible and are explored in detail. We perform our calculations to the leading order in the fine-structure constant. Transverse and longitudinal photon modes are treated separately given their different dispersion relations. We consider the implications for the evolution of massive stars when dark photons decay either into particles of the standard model or of the dark sector.

If the thermal evolution of the hot young neutron star in the supernova remnant HESS J1731-347 is driven by neutrino emission, it provides a stringent constraint on the coupling of light (mass $\ll 10$ keV) axion-like particles to neutrons. Using Markov-Chain Monte Carlo we find that for the values of axion-neutron coupling $g_{ann}^2 > 7.7 \times 10^{-20}$ (90% c.l.) the axion cooling from the bremsstrahlung reaction $n+n\rightarrow n+n +a$ is too rapid to account for the high observed surface temperature. This implies that the Pecci-Quinn scale or axion decay constant $f_a > 6.7 \times 10^7$ GeV for KSVZ axions and $f_a > 1.7 \times 10^9$ GeV for DFSZ axions. The high temperature of this neutron star also allows us to tighten constraints on the size of the nucleon pairing gaps.

We present a comprehensive and pedagogic discussion of the properties of photons in cold and dense nuclear matter based on the resummed one-loop photon self energy. Correlations between electrons, muons, protons and neutrons in beta equilibrium that arise due to electromagnetic and strong interactions are consistently taken into account within the random phase approximation. Screening effects and damping as well as collective excitations are systematically studied in a fully relativistic setup. Our study is relevant to linear response theory of dense nuclear matter, calculations of transport properties of cold dense matter and to investigations of the production and propagation of hypothetical vector bosons such as the dark photons.

We calculate the rate of production of hypothetical light vector bosons (LVBs) from nucleon-nucleon bremsstrahlung reactions in the soft radiation limit directly in terms of the measured nucleon-nucleon elastic cross sections. We use these results and the observation of neutrinos from supernova SN1987a to deduce constraints on the couplings of vector bosons with masses $\lesssim 200$ MeV to either electric charge (dark photons) or to baryon number. We establish for the first time strong constraints on LVB that couple only to baryon number, and revise earlier constraints on the dark photon. For the latter, we find that the excluded region of parameter space is diminished by about a factor of 10.

We calculate neutrino absorption rates due to charged-current reactions $\nu_e+n \rightarrow e^- + p $ and $\bar{\nu}_e+p \rightarrow e^+ + n $ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the $\nu_e$ cross-section and suppress the $\bar{\nu}_e$ cross-section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudo-potential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and the charged-current rates are studied in detail. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density. We find that for typical ambient conditions in the neutrino-sphere, $T=5-10$ MeV and $\rho =10^{11}-10^{13}$ g/cm$^3$, the difference between the $\nu_{e}$ and $\bar{\nu}_e$ absorption rates is much larger than predicted earlier. Our results have important implications for heavy-element nucleosynthesis in supernovae and for supernova neutrino detection.

We show that electrons can couple to the neutron excitations in neutron stars and find that this can limit their contribution to the transport properties of dense matter, especially the shear viscosity. The coupling between electrons and neutrons is induced by protons in the core, and by ions in the crust. We calculate the effective electron-neutron interaction for the kinematics of relevance to the scattering of degenerate electrons at high density. We use this interaction to calculate the electron thermal conductivity, electrical conductivity, and shear viscosity in the neutron star inner crust, and in the core where we consider both normal and superfluid phases of neutron-rich matter. In some cases, particularly when protons are superconducting and neutrons are in their normal phase, we find that electron-neutron scattering can be more important than the other scattering mechanisms considered previously.