Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

After an introduction to the QCD phase diagram, the nuclear equations of state, and neutron star mergers, I discuss three projects related to transport and nuclear matter in neutron star mergers. The first is the nature of beta equilibrium in the portion of a merger that is transparent to neutrinos. We calculate the weak interaction (Urca) rates and find that the beta equilibrium condition needs to be modified by adding an additional chemical potential, which changes slightly the particle content in neutrino-transparent beta equilibrium. Secondly, we calculate the bulk viscosity in neutrino-transparent nuclear matter in conditions encountered in neutron star mergers. Bulk viscosity arises from a phase lag between the pressure and density in the nuclear matter, which is due to the finite rate of beta equilibration. When bulk viscosity is sufficiently strong, which happens when the equilibration rate nearly matches the frequency of the density oscillation, it can noticeably dampen the oscillation. We find that in certain thermodynamic conditions likely encountered in mergers, oscillations in nuclear matter can be damped on timescales on the order of 10 milliseconds, so we conclude that bulk viscosity should be included in merger simulations. Finally, we study thermal transport due to axions in neutron star mergers. We conclude that axions are never trapped in mergers, but instead escape, carrying energy away from the merger. We calculate the cooling time due to the energy carried away by axions and find that within current constraints on the axion-nucleon coupling, axions could cool fluid elements in mergers on timescales which could affect the dynamics of the merger.

We calculate the bulk-viscous dissipation time for adiabatic density oscillations in nuclear matter at densities of 1-7 times nuclear saturation density and at temperatures ranging from 1 MeV, where corrections to previous low-temperature calculations become important, up to 10 MeV, where the assumption of neutrino transparency is no longer valid. Under these conditions, which are expected to occur in neutron star mergers, damping of density oscillations arises from beta equilibration via weak interactions. We find that for 1 kHz oscillations the shortest dissipation times are in the 5 to 20 ms range, depending on the equation of state, which means that bulk viscous damping could affect the dynamics of a neutron star merger. For higher frequencies the dissipation time can be even shorter.

We show that the commonly used criterion for beta equilibrium in neutrino-transparent dense nuclear matter becomes invalid as temperatures rise above 1 MeV. Such temperatures are attained in neutron star mergers. By numerically computing the relevant weak interaction rates we find that the correct criterion for beta equilibrium requires an isospin chemical potential that can be as large as 10-20 MeV, depending on the temperature at which neutrinos become trapped.

The Dark Energy Spectroscopic Instrument (DESI) is under construction to measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 square degrees will be measured during the life of the experiment. A new prime focus corrector for the KPNO Mayall telescope will deliver light to 5000 fiber optic positioners. The fibers in turn feed ten broad-band spectrographs. We present an overview of the instrumentation, the main technical requirements and challenges, and the current status of the project.

The Weak Equivalence Principle states that the geodesics of a test particle in a gravitational field are independent of the particle's constitution. To constrain violations of the Weak Equivalence Principle, we use the one-meter telescope at Table Mountain Observatory near Los Angeles to monitor the relative arrival times of pulses from the Crab Pulsar in the optical ($\lambda \approx 585$ nm) and near-infrared ($\lambda \approx 814$ nm) using an instrument which detects single photons with nanosecond-timing resolution in those two bands. The infrared pulse arrives slightly before the visible pulse. Our three analysis methods give delays with statistical errors of $\Delta t_{obs} = 7.41 \pm 0.58$, $0.4 \pm 3.6$, and $7.35 \pm 4.48$ microseconds (at most 1/4000 of the pulsar period). We attribute this discrepancy to systematic error from the fact that the visible and infrared pulses have slightly different shapes. Whether this delay emerges from the pulsar, is caused by passing through wavelength-dependent media, or is caused by a violation of the equivalence principle, unless there is a fine-tuned cancellation among these, we set the first upper limit on the differential post-Newtonian parameter at these wavelengths of $\Delta \gamma < 1.07 \times 10^{-10}~(3\sigma)$. This result falls in an unexplored region of parameter space and complements existing limits on equivalence-principle violation from fast radio bursts, gamma ray bursts, as well as previous limits from the Crab.

We investigate the stability of stars with a density discontinuity between a high-density core and a very low density mantle. Previous work on "strange dwarfs" suggested that such a discontinuity could stabilize stars that would have been classified as unstable by the conventional criteria based on extrema in the mass-radius relation. We investigate the stability of such stars by numerically solving the Sturm-Liouville equations for the lowest-energy modes of the star. We find that the conventional criteria are correct, and strange dwarfs are not stable.