We present an effective model for the one-dimensional Lyman-$alpha$ flux power spectrum far above the baryonic Jeans scale. The main new ingredient is constituted by a set of two parameters that encode the impact of small, highly non-linear scales on
the one-dimensional power spectrum on large scales, where it is measured by BOSS. We show that, by marginalizing over the model parameters that capture the impact of the intergalactic medium, the flux power spectrum from both simulations and observations can be described with high precision. The model displays a degeneracy between the neutrino masses and the (unknown, in our formalism) normalization of the flux power spectrum. This degeneracy can be lifted by calibrating one of the model parameters with simulation data, and using input from Planck CMB data. We demonstrate that this approach can be used to extract bounds on the sum of neutrino masses with comparably low numerical effort, while allowing for a conservative treatment of uncertainties from the dynamics of the intergalactic medium. An explorative analysis yields an upper bound of $0.16,$eV at $95%$ C.L. when applied to BOSS data at $3leq zleq 4.2$. We also forecast that if the systematic and statistical errors will be reduced by a factor two the upper bound will become $0.1,$eV at $95%$ C.L., and $0.056,$eV when assuming a $1%$ error.
The main features of the gravitational dynamics of binary neutron star systems are now well established. While the inspiral can be precisely described in the post-Newtonian approximation, fully relativistic magneto-hydrodynamical simulations are requ
ired to model the evolution of the merger and post-merger phase. However, the interpretation of the numerical results can often be non-trivial, so that toy models become a very powerful tool. Not only do they simplify the interpretation of the post-merger dynamics, but also allow to gain insights into the physics behind it. In this work, we construct a simple toy model that is capable of reproducing the whole angular momentum evolution of the post-merger remnant, from the merger to the collapse. We validate the model against several fully general-relativistic numerical simulations employing a genetic algorithm, and against additional constraints derived from the spectral properties of the gravitational radiation. As a result, from the remarkably close overlap between the model predictions and the reference simulations within the first milliseconds after the merger, we are able to systematically shed light on the currently open debate regarding the source of the low-frequency peaks of the gravitational wave power spectral density. Additionally, we also present two original relations connecting the angular momentum of the post-merger remnant at merger and collapse to initial properties of the system.
We probe the self-interactions of dark matter using observational data of relaxed galaxy groups and clusters. Our analysis uses the Jeans formalism and considers a wider range of systematic effects than in previous work, including adiabatic contracti
on and stellar anisotropy, to robustly constrain the self-interaction cross section. For both groups and clusters, our results show a mild preference for a nonzero cross section compared with cold collisionless dark matter. Our groups result, $sigma/m=0.5pm0.2~mathrm{cm}^2/mathrm{g}$, places the first constraint on self-interacting dark matter (SIDM) at an intermediate scale between galaxies and massive clusters. Our clusters result is $sigma/m=0.19pm0.09~mathrm{cm}^2/mathrm{g}$, with an upper limit of $sigma / m < 0.35~mathrm{cm}^2/mathrm{g}$ (95% CL). Thus, our results disfavor a velocity-independent cross section of order $1~mathrm{cm}^2/mathrm{g}$ or larger needed to address small scale structure problems in galaxies, but are consistent with a velocity-dependent cross section that decreases with increasing scattering velocity. Comparing the cross sections with and without the effect of adiabatic contraction, we find that adiabatic contraction produces slightly larger values for our data sample, but they are consistent at the $1sigma$ level. Finally, to validate our approach, we apply our Jeans analysis to a sample of mock data generated from SIDM-plus-baryons simulations with $sigma/m = 1~mathrm{cm}^2/mathrm{g}$. This is the first test of the Jeans model at the level of stellar and lensing observables directly measured from simulations. We find our analysis gives a robust determination of the cross section, as well as consistently inferring the true baryon and dark matter density profiles.
Although the main features of the evolution of binary neutron star systems are now well established, many details are still subject to debate, especially regarding the post-merger phase. In particular, the lifetime of the hyper-massive neutron stars
formed after the merger is very hard to predict. In this work, we provide a simple analytic relation for the lifetime of the merger remnant as function of the initial mass of the neutron stars. This relation results from a joint fit of data from observational evidence and from various numerical simulations. In this way, a large range of collapse times, physical effects and equation of states is covered. Finally, we apply the relation to the gravitational wave event GW170817 to constrain the equation of state of dense matter.
The observation of gravitational waves from a binary neutron star merger by LIGO/VIRGO and the associated electromagnetic counterpart provides a high precision test of orbital dynamics, and therefore a new and sensitive probe of extra forces and new
radiative degrees of freedom. Axions are one particularly well-motivated class of extensions to the Standard Model leading to new forces and sources of radiation, which we focus on in this paper. Using an effective field theory (EFT) approach, we calculate the first post-Newtonian corrections to the orbital dynamics, radiated power, and gravitational waveform for binary neutron star mergers in the presence of an axion. This result is applicable to many theories which add an extra massive scalar degree of freedom to General Relativity. We then perform a detailed forecast of the potential for Advanced LIGO to constrain the free parameters of the EFT, and map these to the mass $m_a$ and decay constant $f_a$ of the axion. At design sensitivity, we find that Advanced LIGO can potentially exclude axions with $m_a lesssim 10^{-11} {rm eV}$ and $f_a sim (10^{14} - 10^{17}) {rm GeV}$. There are a variety of complementary observational probes over this region of parameter space, including the orbital decay of binary pulsars, black hole superradiance, and laboratory searches. We comment on the synergies between these various observables.
The Lyman-$alpha$ forest is a valuable probe of dark matter models featuring a scale-dependent suppression of the power spectrum as compared to $Lambda$CDM. In this work, we present a new estimator of the Lyman-$alpha$ flux power spectrum that does n
ot rely on hydrodynamical simulations. Our framework is characterized by nuisance parameters that encapsulate the complex physics of the intergalactic medium and sensitivity to highly non-linear small-scale modes. After validating the approach based on high-resolution hydrodynamical simulations for $Lambda$CDM, we derive conservative constraints on interacting dark matter models from BOSS Lyman-$alpha$ data on large scales, k<0.02(km/s)^(-1), with the relevant nuisance parameters left free in the model fit. The estimator yields lower bounds on the mass of cannibal dark matter, where freeze-out occurs through 3-to-2 annihilation, in the MeV range. Furthermore, we find that models of dark matter interacting with dark radiation, which have been argued to address the $H_0$ and $sigma_8$ tensions, are compatible with BOSS Lyman-$alpha$ data.
Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation
from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degrees of freedom are a characteristic feature of many extensions of General Relativity. For concreteness, we work in the context of metric $f(R)$ gravity, which is equivalent to General Relativity and a universally coupled scalar field with a non-linear potential whose form is fixed by the choice of $f(R)$. In theories where neutron stars (or other compact objects) obtain a significant scalar charge, the resulting attractive finite-range scalar force has implications for both the inspiral and merger phases of binary systems. We first present an analysis of the inspiral dynamics in Newtonian limit, and forecast the constraints on the mass of the scalar and charge of the compact objects for the Advanced LIGO gravitational wave observatory. We then perform a comparative study of binary neutron star mergers in General Relativity with those of a one-parameter model of $f(R)$ gravity using fully relativistic hydrodynamical simulations. These simulations elucidate the effects of the scalar on the merger and post-merger dynamics. We comment on the utility of the full waveform (inspiral, merger, post-merger) to probe different regions of parameter space for both the particular model of $f(R)$ gravity studied here and for finite-range scalar forces more generally.
