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 contraction 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.5\pm0.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.19\pm0.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 $1\sigma$ 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 not 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.

We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an `anisotropic separate universe' framework. Our results suggest that the importance of these coefficient(s) is a $\sim 10 \%$ effect, and plausibly smaller.

We study soft limits of correlation functions for the density and velocity fields in the theory of structure formation. First, we re-derive the (resummed) consistency conditions at unequal times using the eikonal approximation. These are solely based on symmetry arguments and are therefore universal. Then, we explore the existence of equal-time relations in the soft limit which, on the other hand, depend on the interplay between soft and hard modes. We scrutinize two approaches in the literature: The time-flow formalism, and a background method where the soft mode is absorbed into a locally curved cosmology. The latter has been recently used to set up (angular averaged) `equal-time consistency relations'. We explicitly demonstrate that the time-flow relations and `equal-time consistency conditions' are only fulfilled at the linear level, and fail at next-to-leading order for an Einstein de-Sitter universe. While applied to the velocities both proposals break down beyond leading order, we find that the `equal-time consistency conditions' quantitatively approximates the perturbative results for the density contrast. Thus, we generalize the background method to properly incorporate the effect of curvature in the density and velocity fluctuations on short scales, and discuss the reasons behind this discrepancy. We conclude with a few comments on practical implementations and future directions.