Tim Dietrich, Katy Clough

The QCD axion is a hypothetical particle motivated by the Strong CP problem of particle physics. One of the primary ways in which its existence can be inferred is via its function as an additional cooling channel in stars, with some of the strongest constraints coming from the supernova observation SN1987A. Multimessenger observations of binary neutron star mergers (such as those of GW170817, AT2017gfo, and GRB170817A) may provide another scenario in which such constraints could be obtained. In particular, the axion could potentially alter the lifetime, the ejection of material, and the emitted gravitational wave signal of the postmerger remnant. In this article, we perform numerical relativity simulations of a binary neutron star merger, including a phenomenological description of the nucleon-nucleon-axion bremsstrahlung to quantify the effects of such a cooling channel on the dynamical evolution. While our simulations show a difference in the temperature profile of the merger remnant, the imprint of the axion via nucleon-nucleon-axion bremsstrahlung on the emitted gravitational wave signal and the ejecta mass is too small to improve constraints on the axion mass with current or future planned detectors. Whilst we consider a limited number of cases, and a simplified cooling model, these broadly represent the "best case" scenario, thus, a more thorough investigation is unlikely to change the conclusions, at least for this particular interaction channel.

Johannes U. Lange, Frank C. van den Bosch, Andrew R. Zentner, Kuan Wang, Andrew P. Hearin, Hong Guo

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

Claudia de Rham, Andrew J. Tolley

Within the standard effective field theory of General Relativity, we show that the speed of gravitational waves deviates, ever so slightly, from luminality on cosmological and other spontaneously Lorentz-breaking backgrounds. This effect results from loop contributions from massive fields of any spin, including Standard Model fields, or from tree level effects from massive higher spins $s \ge 2$. We show that for the choice of interaction signs implied by S-matrix and spectral density positivity bounds suggested by analyticity and causality, the speed of gravitational waves is in general superluminal at low-energies on NEC preserving backgrounds, meaning gravitational waves travel faster than allowed by the metric to which photons and Standard Model fields are minimally coupled. We show that departure of the speed from unity increases in the IR and argue that the speed inevitably returns to luminal at high energies as required by Lorentz invariance. Performing a special tuning of the EFT so that renormalization sensitive curvature-squared terms are set to zero, we find that finite loop corrections from Standard Model fields still lead to an epoch dependent modification of the speed of gravitational waves which is determined by the precise field content of the lightest particles with masses larger than the Hubble parameter today. Depending on interpretation, such considerations could potentially have far-reaching implications on light scalar models, such as axionic or fuzzy cold dark matter.