In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor (β-)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.

We describe and implement a procedure for determining the couplings of a relativistic mean-field theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers.

Significance: The diffusion of particles and waves through disordered media encompasses a large variety of phenomena, from the motion of insects to the scattering of electrons or light in complex environments. One of the core features of diffusive transport is that the mean length of trajectories traversing a system depends only on the size of the system and of its boundary, which are both independent of the microscopic structure of the underlying medium. Here we show, based on insights from wave-scattering theory, that this fundamental invariance property can be significantly extended beyond the diffusive random walk picture. Our result not only provides an interesting link between all the diverse fields in which wave scattering plays a role but also holds promise for a number of practical applications.

A strong magnetic field enhances the chiral condensate at low temperatures. This so-called magnetic catalysis thus seeks to increase the vacuum mass of nucleons. We employ two relativistic field-theoretical models for nuclear matter, the Walecka model and an extended linear sigma model, to discuss the resulting effect on the transition between vacuum and nuclear matter at zero temperature. In both models we find that the creation of nuclear matter in a sufficiently strong magnetic field becomes energetically more costly due to the heaviness of magnetized nucleons, even though it is also found that nuclear matter is more strongly bound in a magnetic field. Our results are potentially important for dense nuclear matter in compact stars, especially since previous studies in the astrophysical context have always ignored the contribution of the magnetized Dirac sea and thus the effect of magnetic catalysis.

We calculate the nonzero-temperature correction to the beta equilibrium condition in nuclear matter under neutron star merger conditions, in the temperature range 1mEv < T ≲ 5 mEv. We improve on previous work using a consistent description of nuclear matter based on the IUF and SFHo relativistic mean field models. This includes using relativistic dispersion relations for the nucleons, which we show is essential in these models. We find that the nonzero-temperature correction can be of order 10 to 20 MeV, and plays an important role in the correct calculation of Urca rates, which can be wrong by factors of 10 or more if it is neglected.

In this paper, we present a computation of the rates of strangeness-changing processes and the resultant bulk viscosity in matter at the densities and temperatures typical of neutron star mergers. To deal with the high temperature in this environment, we go beyond the Fermi surface approximation in our rate calculations and numerically evaluate the full phase space integral. We include processes where quarks move between baryons via meson exchange: These have generally been omitted in previous analyses but provide the dominant contribution to the rates of strangeness-changing processes and the bulk viscosity. The calculation of these rates is an essential step toward any calculation of dissipation mechanisms in hyperonic matter in mergers. As one application, we calculate the dissipation times for density oscillations at the frequencies seen in merger simulations. We find that hyperon bulk viscosity for temperatures in the MeV regime can probably be neglected in this context but becomes highly relevant for keV-range temperatures.