Search

Current and future large scale structure surveys aim to constrain the neutrino mass and the equation of state of dark energy. We aim to construct accurate and interpretable symbolic approximations to the linear and nonlinear matter power spectra as a function of cosmological parameters in extended $\Lambda$CDM models which contain massive neutrinos and non-constant equations of state for dark energy. This constitutes an extension of the syren-halofit emulators to incorporate these two effects, which we call syren-new (SYmbolic-Regression-ENhanced power spectrum emulator with NEutrinos and $W_0-w_a$). We also obtain a simple approximation to the derived parameter $\sigma_8$ as a function of the cosmological parameters for these models. Our results for the linear power spectrum are designed to emulate CLASS, whereas for the nonlinear case we aim to match the results of EuclidEmulator2. We compare our results to existing emulators and $N$-body simulations. Our analytic emulators for $\sigma_8$, the linear and nonlinear power spectra achieve root mean squared errors of 0.1%, 0.3% and 1.3%, respectively, across a wide range of cosmological parameters, redshifts and wavenumbers. We verify that emulator-related discrepancies are subdominant compared to observational errors and other modelling uncertainties when computing shear power spectra for LSST-like surveys. Our expressions have similar accuracy to existing (numerical) emulators, but are at least an order of magnitude faster, both on a CPU and GPU. Our work greatly improves the accuracy, speed and range of applicability of current symbolic approximations to the linear and nonlinear matter power spectra. We provide publicly available code for all symbolic approximations found., Comment: 18 pages, 15 figures

We show that a universe with a nonminimally coupled scalar field can fit current measurements of the expansion rate of the Universe better than the standard $\Lambda$-Cold Dark Matter model or other minimally coupled dark energy models. While we find a clear improvement in the goodness of fit for this dark energy model with respect to others that have been considered in the recent literature, using various information theoretic criteria, we show that the evidence for it is still inconclusive., Comment: Updated version. Minor changes to the discussion and references added. 5 pages, 2 figures

Gravitational waves can provide crucial insights about the environments in which black holes live. In this work, we use numerical relativity simulations to study the behaviour of self-interacting scalar (wave-like) dark matter clouds accreting onto isolated and binary black holes. We find that repulsive self-interactions smoothen the ``spike'' of an isolated black hole and saturate the density. Attractive self-interactions enhance the growth and result in more cuspy profiles, but can become unstable and undergo explosions akin to the superradiant bosenova that reduce the local cloud density. We quantify the impact of self-interactions on an equal-mass black hole merger by computing the dephasing of the gravitational-wave signal for a range of couplings. We find that repulsive self-interactions saturate the density of the cloud, thereby reducing the dephasing. For attractive self-interactions, the dephasing may be larger, but if these interactions dominate prior to the merger, the dark matter can undergo bosenova during the inspiral phase, disrupting the cloud and subsequently reducing the dephasing., Comment: 8 pages, 7 figures, 1 movie: https://youtu.be/ZeoMlUHknUM. Matches version published in PRD

New constraints on the expansion rate of the Universe seem to favor evolving dark energy in the form of thawing quintessence models, i.e., models for which a canonical, minimally coupled scalar field has, at late times, begun to evolve away from potential energy domination. We scrutinize the evidence for thawing quintessence by exploring what it predicts for the equation of state. We show that, in terms of the usual Chevalier-Polarski-Linder parameters, ($w_0$, $w_a$), thawing quintessence is, in fact, only marginally consistent with a compilation of the current data. Despite this, we embrace the possibility that thawing quintessence is dark energy and find constraints on the microphysics of this scenario. We do so in terms of the effective mass $m^2$ and energy scale $V_0$ of the scalar field potential. We are particularly careful to enforce un-informative, flat priors on these parameters so as to minimize their effect on the final posteriors. While the current data favors a large and negative value of $m^2$, when we compare these models to the standard $\Lambda$CDM model we find that there is scant evidence for thawing quintessence., Comment: Accepted and forthcoming in PRD. Minor changes to discussions. 18 pages, 14 figures

Physics, as a fundamental science, aims to understand the laws of Nature and describe them in mathematical equations. While the physical reality manifests itself in a wide range of phenomena with varying levels of complexity, the equations that describe them display certain statistical regularities and patterns, which we begin to explore here. By drawing inspiration from linguistics, where Zipf's law states that the frequency of any word in a large corpus of text is roughly inversely proportional to its rank in the frequency table, we investigate whether similar patterns for the distribution of operators emerge in the equations of physics. We analyse three corpora of formulae and find, using sophisticated implicit-likelihood methods, that the frequency of operators as a function of their rank in the frequency table is best described by an exponential law with a stable exponent, in contrast with Zipf's inverse power-law. Understanding the underlying reasons behind this statistical pattern may shed light on Nature's modus operandi or reveal recurrent patterns in physicists' attempts to formalise the laws of Nature. It may also provide crucial input for symbolic regression, potentially augmenting language models to generate symbolic models for physical phenomena. By pioneering the study of statistical regularities in the equations of physics, our results open the door for a meta-law of Nature, a (probabilistic) law that all physical laws obey., Comment: 9 pages, 5 figures

Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data allows us to place constraints on the (bias-weighted) production rate density of gravitational waves by astrophysical sources as a function of cosmic time. Our constraints range from $\langle b\dot{\Omega}_{\rm GW}\rangle<3.0\times10^{-9}\,{\rm Gyr}^{-1}$ at $z\sim0.06$ to $\langle b\dot{\Omega}_{\rm GW}\rangle<2.7\times10^{-7}\,{\rm Gyr}^{-1}$ at $z\sim1.5$ (95\% C.L.), assuming a frequency spectrum of the form $f^{2/3}$ (corresponding to an astrophysical background of binary mergers), and a reference frequency $f_{\rm ref}=25\,{\rm Hz}$. Although these constraints are $\sim2$ orders of magnitude higher than the expected signal, we show that a detection may be possible with future experiments., Comment: 14 pages, 11 figures

Vacuum decay and symmetry breaking play an important role in the fundamental structure of the matter and the evolution of the universe. In this work we study how the purely classical effect of accretion of fundamental fields onto black holes can lead to shells of symmetry restoration in the midst of a symmetry broken phase. We also show how it can catalyze vacuum decay, forming a bubble that expands asymptotically at the speed of light. These effects offer an alternative, purely classical mechanism to quantum tunnelling for seeding phase transitions in the universe., Comment: 4 pages, 4 figures, 1 movie: https://youtu.be/c2vcNoX8fVs Comments welcome!

Rapid and accurate evaluation of the nonlinear matter power spectrum, $P(k)$, as a function of cosmological parameters and redshift is of fundamental importance in cosmology. Analytic approximations provide an interpretable solution, yet current approximations are neither fast nor accurate relative to numerical emulators. We use symbolic regression to obtain simple analytic approximations to the nonlinear scale, $k_\sigma$, the effective spectral index, $n_{\rm eff}$, and the curvature, $C$, which are required for the halofit model. We then re-optimise the coefficients of halofit to fit a wide range of cosmologies and redshifts. We explore the space of analytic expressions to fit the residuals between $P(k)$ and the optimised predictions of halofit. Our results are designed to match the predictions of EuclidEmulator2, but are validated against $N$-body simulations. Our symbolic expressions for $k_\sigma$, $n_{\rm eff}$ and $C$ have root mean squared fractional errors of 0.8%, 0.2% and 0.3%, respectively, for redshifts below 3 and a wide range of cosmologies. The re-optimised halofit parameters reduce the root mean squared fractional error (compared to EuclidEmulator2) from 3% to below 2% for wavenumbers $k=9\times10^{-3}-9 \, h{\rm Mpc^{-1}}$. We introduce syren-halofit (symbolic-regression-enhanced halofit), an extension to halofit containing a short symbolic correction which improves this error to 1%. Our method is 2350 and 3170 times faster than current halofit and hmcode implementations, respectively, and 2680 and 64 times faster than EuclidEmulator2 (which requires running class) and the BACCO emulator. We obtain comparable accuracy to EuclidEmulator2 and BACCO when tested on $N$-body simulations. Our work greatly increases the speed and accuracy of symbolic approximations to $P(k)$, making them significantly faster than their numerical counterparts without loss of accuracy., Comment: 11 pages, 8 figures. Accepted for publication in A&A

For dark matter to be detectable with gravitational waves from binary black holes, it must reach higher than average densities in their vicinity. In the case of light (wave-like) dark matter, the density of dark matter between the binary can be significantly enhanced by accretion from the surrounding environment. Here we show that the resulting dephasing effect on the last ten orbits of an equal mass binary is maximized when the Compton wavelength of the scalar particle is comparable to the orbital separation, $2\pi/\mu\sim d$. The phenomenology of the effect is different from the channels that are usually discussed, where dynamical friction (along the orbital path) and radiation of energy and angular momentum drive the dephasing, and is rather dominated by the radial force (the spacetime curvature in the radial direction) towards the overdensity between the black holes. Whilst our numerical studies limit us to scales of the same order, this effect may persist at larger separations and/or particle masses, playing a significant role in the merger history of binaries., Comment: 4 pages, 4 figures, 1 appendix, 1 movie: https://youtu.be/2VJIfqCp7D8. Matches version published in PRL

Computing the matter power spectrum, $P(k)$, as a function of cosmological parameters can be prohibitively slow in cosmological analyses, hence emulating this calculation is desirable. Previous analytic approximations are insufficiently accurate for modern applications, so black-box, uninterpretable emulators are often used. We utilise an efficient genetic programming based symbolic regression framework to explore the space of potential mathematical expressions which can approximate the power spectrum and $\sigma_8$. We learn the ratio between an existing low-accuracy fitting function for $P(k)$ and that obtained by solving the Boltzmann equations and thus still incorporate the physics which motivated this earlier approximation. We obtain an analytic approximation to the linear power spectrum with a root mean squared fractional error of 0.2% between $k = 9\times10^{-3} - 9 \, h{\rm \, Mpc^{-1}}$ and across a wide range of cosmological parameters, and we provide physical interpretations for various terms in the expression. Our analytic approximation is 950 times faster to evaluate than camb and 36 times faster than the neural network based matter power spectrum emulator BACCO. We also provide a simple analytic approximation for $\sigma_8$ with a similar accuracy, with a root mean squared fractional error of just 0.1% when evaluated across the same range of cosmologies. This function is easily invertible to obtain $A_{\rm s}$ as a function of $\sigma_8$ and the other cosmological parameters, if preferred. It is possible to obtain symbolic approximations to a seemingly complex function at a precision required for current and future cosmological analyses without resorting to deep-learning techniques, thus avoiding their black-box nature and large number of parameters. Our emulator will be usable long after the codes on which numerical approximations are built become outdated., Comment: 9 pages, 5 figures. Accepted for publication in A&A

Inflation is a highly favoured theory for the early Universe. It is compatible with current observations of the cosmic microwave background and large scale structure and is a driver in the quest to detect primordial gravitational waves. It is also, given the current quality of the data, highly under-determined with a large number of candidate implementations. We use a new method in symbolic regression to generate all possible simple scalar field potentials for one of two possible basis sets of operators. Treating these as single-field, slow-roll inflationary models we then score them with an information-theoretic metric ("minimum description length") that quantifies their efficiency in compressing the information in current data. We explore two possible priors on the parameter space of potentials, one related to the functions' structural complexity and one that uses a Katz back-off language model to prefer functions that may be theoretically motivated. This enables us to identify the inflaton potentials that optimally balance simplicity with accuracy at explaining current data, which may subsequently find theoretical motivation. Our exploratory study opens the door to extraction of fundamental physics directly from data, and may be augmented with more refined theoretical priors in the quest for a complete understanding of the early Universe., Comment: 13+4 pages, 4 figures; Accepted for publication in Physical Review D

There is compelling evidence that the Universe is undergoing a late phase of accelerated expansion. One of the simplest explanations for this behaviour is the presence of dark energy. A plethora of microphysical models for dark energy have been proposed. The hope is that, with the ever increasing precision of cosmological surveys, it will be possible to precisely pin down the model. We show that this is unlikely and that, at best, we will have a phenomenological description for the microphysics of dark energy. Furthermore, we argue that the current phenomenological prescriptions are ill-equipped for shedding light on the fundamental theory of dark energy., Comment: Updated the figures to have higher resolution. 14 pages, 8 figures

We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle masses (relative to the BH size). In the case of a small dark matter cloud, the effect of accretion is a non-negligible contribution to the total force on the black hole, even for small scalar masses. We confirm that this momentum accretion transitions between two regimes (wave- and particle-like) and we identify the mass of the scalar at which the transition between regimes occurs., Comment: 11 pages, 5 figures. Minor corrections and references added to match published version

When choosing between competing symbolic models for a data set, a human will naturally prefer the "simpler" expression or the one which more closely resembles equations previously seen in a similar context. This suggests a non-uniform prior on functions, which is, however, rarely considered within a symbolic regression (SR) framework. In this paper we develop methods to incorporate detailed prior information on both functions and their parameters into SR. Our prior on the structure of a function is based on a $n$-gram language model, which is sensitive to the arrangement of operators relative to one another in addition to the frequency of occurrence of each operator. We also develop a formalism based on the Fractional Bayes Factor to treat numerical parameter priors in such a way that models may be fairly compared though the Bayesian evidence, and explicitly compare Bayesian, Minimum Description Length and heuristic methods for model selection. We demonstrate the performance of our priors relative to literature standards on benchmarks and a real-world dataset from the field of cosmology., Comment: 8+2 pages, 2 figures. Accepted for The Genetic and Evolutionary Computation Conference (GECCO) 2023 Workshop on Symbolic Regression

As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte-Carlo sampling methods. A particularly relevant example is the redshift distribution, $p(z)$, of the source samples, which may require tens of parameters to describe fully. However, relatively tight priors can be usually placed on these parameters through calibration of the associated systematics. In this paper we show, quantitatively, that a linearisation of the theoretical prediction with respect to these calibratable systematic parameters allows us to analytically marginalise over these extra parameters, leading to a factor $\sim30$ reduction in the time needed for parameter inference, while accurately recovering the same posterior distributions for the cosmological parameters that would be obtained through a full numerical marginalisation over 160 $p(z)$ parameters. We demonstrate that this is feasible not only with current data and current achievable calibration priors but also for future Stage-IV datasets., Comment: 11 pages, 8 figures, prepared for submission to MNRAS, comments welcome

We apply a new method for learning equations from data -- Exhaustive Symbolic Regression (ESR) -- to late-type galaxy dynamics as encapsulated in the radial acceleration relation (RAR). Relating the centripetal acceleration due to baryons, $g_\text{bar}$, to the total dynamical acceleration, $g_\text{obs}$, the RAR has been claimed to manifest a new law of nature due to its regularity and tightness, in agreement with Modified Newtonian Dynamics (MOND). Fits to this relation have been restricted by prior expectations to particular functional forms, while ESR affords an exhaustive and nearly prior-free search through functional parameter space to identify the equations optimally trading accuracy with simplicity. Working with the SPARC data, we find the best functions typically satisfy $g_\text{obs} \propto g_\text{bar}$ at high $g_\text{bar}$, although the coefficient of proportionality is not clearly unity and the deep-MOND limit $g_\text{obs} \propto \sqrt{g_\text{bar}}$ as $g_\text{bar} \to 0$ is little evident at all. By generating mock data according to MOND with or without the external field effect, we find that symbolic regression would not be expected to identify the generating function or reconstruct successfully the asymptotic slopes. We conclude that the limited dynamical range and significant uncertainties of the SPARC RAR preclude a definitive statement of its functional form, and hence that this data alone can neither demonstrate nor rule out law-like gravitational behaviour., Comment: 12+4 pages, 4 figures, 3 tables; minor revision to match MNRAS published version

Symbolic Regression (SR) algorithms attempt to learn analytic expressions which fit data accurately and in a highly interpretable manner. Conventional SR suffers from two fundamental issues which we address here. First, these methods search the space stochastically (typically using genetic programming) and hence do not necessarily find the best function. Second, the criteria used to select the equation optimally balancing accuracy with simplicity have been variable and subjective. To address these issues we introduce Exhaustive Symbolic Regression (ESR), which systematically and efficiently considers all possible equations -- made with a given basis set of operators and up to a specified maximum complexity -- and is therefore guaranteed to find the true optimum (if parameters are perfectly optimised) and a complete function ranking subject to these constraints. We implement the minimum description length principle as a rigorous method for combining these preferences into a single objective. To illustrate the power of ESR we apply it to a catalogue of cosmic chronometers and the Pantheon+ sample of supernovae to learn the Hubble rate as a function of redshift, finding $\sim$40 functions (out of 5.2 million trial functions) that fit the data more economically than the Friedmann equation. These low-redshift data therefore do not uniquely prefer the expansion history of the standard model of cosmology. We make our code and full equation sets publicly available., Comment: 15 pages, 7 figures, 2 tables. Accepted for publication in the IEEE Transactions on Evolutionary Computation

An accurate calibration of the source redshift distribution $p(z)$ is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the $p(z)$ statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the ``clumpiness'' parameter to $S_8 = 0.768^{+0.021}_{-0.017}$. This corresponds to a $\sim\sqrt{2}$ improvement in uncertainties with respect to either DES or HSC alone., Comment: 32 pages, 12 figures, comments welcome; v2: minor changes, accepted version in JCAP

The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behaviour of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasi-stationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. Whilst different initial configurations converge to the preferred quasi-stationary one after several orbits, unwanted transient oscillations are generated in the process, which may impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way., Comment: 9 pages, 5 figures, 1 appendix, 2 movies: https://youtube.com/playlist?list=PLSkfizpQDrcb8Zd7C9yEfnMpnKQk7I9ts Comments welcome!

In the context of modified gravity, at linear level, the growth of structure in the Universe will be affected by modifications to the Poisson equation and by the background expansion rate of the Universe. It has been shown that these two effects lead to a degeneracy which must be properly accounted for if one is to place reliable constraints on new forces on large scales or, equivalently, modifications to General Relativity. In this paper we show that current constraints are such that assumptions about the background expansion have little impact on constraints on modifications to gravity. We do so by considering the background of a $\Lambda$ Cold Dark Matter ($\Lambda$CDM) universe, a universe with a more general equation of state for the dark energy, and finally, a general, model-independent, expansion rate. We use Gaussian Processes to model modifications to Poisson's equation and, in the case of a general expansion rate, to model the redshift dependent Hubble rate. We identify a degeneracy between modifications to Poisson's equation and the background matter density, $\Omega_M$, which can only be broken by assuming a model-dependent expansion rate. We show that, with current data, the constraints on modifications to the Poisson equation via measurements of the growth rate range between $10-20\%$ depending on the strength of our assumptions on the Universe's expansion rate., Comment: Prepared for submission to PRD. Comments welcome

We study the generation and propagation of gravitational waves in scalar-tensor gravity using numerical relativity simulations of scalar field collapses beyond spherical symmetry. This allows us to compare the tensor and additional massive scalar waves that are excited. As shown in previous work in spherical symmetry, massive propagating scalar waves decay faster than 1/r and disperse, resulting in an inverse chirp. These effects obscure the ringdown in any extracted signal by mixing it with the transient responses of the collapse during propagation. In this paper we present a simple method to rewind the extracted signals to horizon formation, which allows us to clearly identify the ringdown phase and extract the amplitudes of the scalar quasinormal modes, quantifying their excitation in strong gravity events and verifying the frequencies to perturbative calculations. The effects studied are relevant to any theories in which the propagating waves have a dispersion relation, including the tensor case., Comment: 10 pages, 8 figures. Movies: https://youtu.be/YYBwyAbH5Fk. Matches version published in PRD

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas., Comment: Accepted in: Living Reviews in Relativity

We constrain the expansion history of the Universe and the cosmological matter density fraction in a model-independent way by exclusively making use of the relationship between background and perturbations under a minimal set of assumptions. We do so by employing a Gaussian process to model the expansion history of the Universe from present time to the recombination era. The expansion history and the cosmological matter density are then constrained using recent measurements from cosmic chronometers, Type-Ia supernovae, baryon acoustic oscillations, and redshift-space distortion data. Our results show that the evolution in the reconstructed expansion history is compatible with the \textit{Planck} 2018 prediction at all redshifts. The current data considered in this study can constrain a Gaussian process on $H(z)$ to an average $9.4 \%$ precision across redshift. We find $\Omega_m = 0.224 \pm 0.066$, lower but statistically compatible with the \textit{Planck} 2018 cosmology. Finally, the combination of future DESI measurements with the CMB measurement considered in this work holds the promise of $8 \%$ average constraints on a model-independent expansion history as well as a five-fold tighter $\Omega_m$ constraint using the methodology developed in this work., Comment: Accepted to MNRAS

GRChombo is an open-source code for performing Numerical Relativity time evolutions, built on top of the publicly available Chombo software for the solution of PDEs. Whilst GRChombo uses standard techniques in NR, it focusses on applications in theoretical physics where adaptability, both in terms of grid structure, and in terms of code modification, are key drivers., Comment: JOSS submission for first official release and review of the GRChombo GitHub repository: https://github.com/GRChombo/GRChombo, https://doi.org/10.21105/joss.03703

Lorentz invariance violation in quantum gravity (QG) models or a nonzero photon mass, $m_\gamma$, would lead to an energy-dependent propagation speed for photons, such that photons of different energies from a distant source would arrive at different times, even if they were emitted simultaneously. By developing source-by-source, Monte Carlo-based forward models for such time delays from gamma ray bursts, and marginalising over empirical noise models describing other contributions to the time delay, we derive constraints on $m_\gamma$ and the QG length scale, $\ell_{\rm QG}$, using spectral lag data from the BATSE satellite. We find $m_\gamma < 4.0 \times 10^{-5} \, h \, {\rm eV}/c^2$ and $\ell_{\rm QG} < 5.3 \times 10^{-18} \, h \, {\rm \, GeV^{-1}}$ at 95% confidence, and demonstrate that these constraints are robust to the choice of noise model. The QG constraint is among the tightest from studies which consider multiple gamma ray bursts and the constraint on $m_\gamma$, although weaker than from using radio data, provides an independent constraint which is less sensitive to the effects of dispersion by electrons., Comment: 6 pages, 1 figure, 1 table, accepted for publication in Physical Review D

We analyse Starobinsky inflation in the presence of the Brout Englert Higgs (BEH) boson with a non-minimal coupling to the Ricci scalar, $R$. The latter induces a coupling of the massive scaleron associated with the $R^2$ term to the BEH boson and this leads to a radiative correction to the BEH mass that must be fine tuned to keep the scalar light. For the case of $R^2$ driven inflation this requires a high level of fine tuning of order 1 part in $10^{8}$; for the case of Higgs inflation it is very much greater. We consider a scale invariant extension of the $R^2$/Higgs model and find that for $R^2$ driven inflation but not for Higgs inflation the required fine tuning is significantly reduced to one part in $10^{3-4}$. We consider the vacuum stability of the fine tuned model and its reheating and dilaton abundance after inflation. We also discuss possible gravitational wave signals associated with the model and the constraint on the mass of scalar or fermion dark matter candidates if they are produced by the gravitational couplings of the scalaron.

Theories of gravity that obey the Weak Equivalence Principle have the same Parametrised Post-Newtonian parameter $\gamma$ for all particles at all energies. The large Shapiro time delays of extragalactic sources allow us to put tight constraints on differences in $\gamma$ between photons of different frequencies from spectral lag data, since a non-zero $\Delta \gamma$ would result in a frequency-dependent arrival time. The majority of previous constraints have assumed that the Shapiro time delay is dominated by a few local massive objects, although this is a poor approximation for distant sources. In this work we consider the cosmological context of these sources by developing a source-by-source, Monte Carlo-based forward model for the Shapiro time delays by combining constrained realisations of the local density field using the Bayesian origin reconstruction from galaxies algorithm with unconstrained large-scale modes. Propagating uncertainties in the density field reconstruction and marginalising over an empirical model describing other contributions to the time delay, we use spectral lag data of Gamma Ray Bursts from the BATSE satellite to constrain $\Delta \gamma < 2.1 \times 10^{-15}$ at $1 \sigma$ confidence between photon energies of $25 {\rm \, keV}$ and $325 {\rm \, keV}$., Comment: 14 pages, 7 figures, 3 tables, accepted for publication in Physical Review D

Light bosonic scalars (e.g. axions) may form clouds around black holes via superradiant instabilities, or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we study numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. We show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar the rapid force oscillations could in principle leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral., Comment: 17 pages, 13 figures, 1 table; Version as accepted in PRD - minor changes

In order to investigate the origin of the ongoing tension between the amplitude of matter fluctuations measured by weak lensing experiments at low redshifts and the value inferred from the cosmic microwave background anisotropies, we reconstruct the evolution of this amplitude from $z\sim2$ using existing large-scale structure data. To do so, we decouple the linear growth of density inhomogeneities from the background expansion, and constrain its redshift dependence making use of a combination of 6 different data sets, including cosmic shear, galaxy clustering and CMB lensing. We analyze these data under a consistent harmonic-space angular power spectrum-based pipeline. We show that current data constrain the amplitude of fluctuations mostly in the range $0.2

Attempts at constraining theories of late time accelerated expansion often assume broad priors for the parameters in their phenomenological description. Focusing on shift-symmetric scalar-tensor theories with standard gravitational wave speed, we show how a more careful analysis of their dynamical evolution leads to much narrower priors. In doing so, we propose a simple and accurate parametrisation of these theories, capturing the redshift dependence of the equation of state, $w(z)$, and the kinetic braiding parameter, $\alpha_{\rm B}(z)$, with only two parameters each, and derive their statistical distribution (a.k.a. theoretical priors) that fit the cosmology of the underlying model. We have considered two versions of the shift-symmetric model, one where the energy density of dark energy is given solely by the scalar field, and another where it also has a contribution from the cosmological constant. By including current data, we show how theoretical priors can be used to improve constraints by up to an order of magnitude. Moreover, we show that shift-symmetric theories without a cosmological constant are observationally viable. We work up to quartic order in first derivatives of the scalar in the action and our results suggest this truncation is a good approximation to more general shift-symmetric theories. This work establishes an actionable link between phenomenological parameterisations and Lagrangian-based theories, the two main approaches to test cosmological gravity and cosmic acceleration., Comment: 18 pages, 13 figures; Version as accepted in PRD - minor changes

Galactic scale tests have proven to be powerful tools in constraining fundamental physics in previously under-explored regions of parameter space. The astrophysical regime which they probe is inherently complicated, and the inference methods used to make these constraints should be robust to baryonic effects. Previous analyses have assumed simple empirical models for astrophysical noise without detailed calibration or justification. We outline a framework for assessing the reliability of such methods by constructing and testing more advanced baryonic models using cosmological hydrodynamical simulations. As a case study, we use the Horizon-AGN simulation to investigate warping of stellar disks and offsets between gas and stars within galaxies, which are powerful probes of screened fifth forces. We show that the degree of `U'-shaped warping of galaxies is well modelled by Gaussian random noise, but that the magnitude of the gas-star offset is correlated with the virial radius of the host halo. By incorporating this correlation we confirm recent results ruling out astrophysically relevant Hu-Sawicki $f(R)$ gravity, and identify a $\sim 30\%$ systematic uncertainty due to baryonic physics. Such an analysis must be performed case-by-case for future galactic tests of fundamental physics., Comment: 17 pages, 11 figures, 2 tables, accepted for publication in Physical Review D

The ringdown of a perturbed black hole contains fundamental information about space-time in the form of Quasi Normal Modes (QNM). Modifications to general relativity, or extended profiles of other fields surrounding the black hole, so called "black hole hair", can perturb the QNM frequencies. Previous works have examined the QNM frequencies of spherically symmetric "dirty" black holes - that is black holes surrounded by arbitrary matter fields. Such analyses were restricted to static systems, making the assumption that the metric perturbation was independent of time. However, in most physical cases such black holes will actually be growing dynamically due to accretion of the surrounding matter. Here we develop a perturbative analytic method that allows us to compute for the first time the time dependent QNM deviations of such growing dirty black holes. Whilst both are small, we show that the change in QNM frequency due to the accretion can be of the same order or larger than the change due to the static matter distribution itself, and therefore should not be neglected in such calculations. We present the case of spherically symmetric accretion of a complex scalar field as an illustrative example, but the method has the potential to be extended to more complicated cases., Comment: 14 pages, 1 figure

Scalar fields around compact objects are of interest for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions. We study the behaviour of a scalar field around a Kerr black hole with non trivial asymptotic boundary conditions - both non zero density and non zero angular momentum. Starting from an initial radially homogeneous configuration, a scalar cloud is accreted, which asymptotes to known stationary configurations over time. We study the cloud growth for different parameters including black hole spin, scalar field mass, and the scalar field density and angular momentum far from the black hole. We characterise the transient growth of the mass and angular momentum in the cloud, and the spatial profile of the scalar around the black hole, and relate the results of fully non-linear simulations to an analytic perturbative expansion. We also highlight the potential for these accreted clouds to create monochromatic gravitational wave signals - similar to the signals from superradiant clouds, although significantly weaker in amplitude., Comment: 21 pages, 24 figures

We provide an end-to-end exploration of a distinct modified gravitational theory in Jordan-Brans-Dicke (JBD) gravity, from an analytical and numerical description of the background expansion and linear perturbations, to the nonlinear regime captured with a hybrid suite of $N$-body simulations, to the parameter constraints from existing cosmological probes. The nonlinear corrections to the matter power spectrum due to baryons, massive neutrinos, and modified gravity are simultaneously modeled and propagated in the cosmological analysis for the first time. In the combined analysis of the Planck CMB temperature, polarization, and lensing reconstruction, Pantheon supernova distances, BOSS measurements of BAO distances, the Alcock-Paczynski effect, and the growth rate, along with the joint ($3\times2$pt) dataset of cosmic shear, galaxy-galaxy lensing, and overlapping redshift-space galaxy clustering from KiDS and 2dFLenS, we constrain the JBD coupling constant, $\omega_{\rm BD}>1540$ (95% CL), the effective gravitational constant, $G_{\rm matter}/G=0.997\pm0.029$, the sum of neutrino masses, $\sum m_{\nu}<0.12$ eV (95% CL), and the baryonic feedback amplitude, $B<2.8$ (95% CL), all in agreement with the standard model expectation. We show that the uncertainty in the gravitational theory alleviates the tension between KiDS$\times$2dFLenS and Planck to below $1\sigma$ and the tension in the Hubble constant between Planck and the direct measurement of Riess et al. (2019) down to ~$3\sigma$; however, we find no substantial model selection preference for JBD gravity relative to $\Lambda$CDM. We further show that the neutrino mass bound degrades by up to a factor of three as the $\omega_{\rm BD}$ parameterization becomes more restrictive and that a positive shift in $G_{\rm matter}/G$ suppresses the CMB damping tail in a way that might complicate future inferences of small-scale physics. (Abridged), Comment: 53 pages, 26 figures, results unchanged, version accepted for publication by PRD. Abstract abridged. Our MCMC chains, likelihood code, and data products are public at https://github.com/sjoudaki/CosmoJBD

Cosmic shear is one of the most powerful probes of Dark Energy, targeted by several current and future galaxy surveys. Lensing shear, however, is only sampled at the positions of galaxies with measured shapes in the catalog, making its associated sky window function one of the most complicated amongst all projected cosmological probes of inhomogeneities, as well as giving rise to inhomogeneous noise. Partly for this reason, cosmic shear analyses have been mostly carried out in real-space, making use of correlation functions, as opposed to Fourier-space power spectra. Since the use of power spectra can yield complementary information and has numerical advantages over real-space pipelines, it is important to develop a complete formalism describing the standard unbiased power spectrum estimators as well as their associated uncertainties. Building on previous work, this paper contains a study of the main complications associated with estimating and interpreting shear power spectra, and presents fast and accurate methods to estimate two key quantities needed for their practical usage: the noise bias and the Gaussian covariance matrix, fully accounting for survey geometry, with some of these results also applicable to other cosmological probes. We demonstrate the performance of these methods by applying them to the latest public data releases of the Hyper Suprime-Cam and the Dark Energy Survey collaborations, quantifying the presence of systematics in our measurements and the validity of the covariance matrix estimate. We make the resulting power spectra, covariance matrices, null tests and all associated data necessary for a full cosmological analysis publicly available., Comment: 42 pages, 16 figures, 4 tables, to be submitted to JCAP

Galileons are scalar field theories which obey the Galileon symmetry $\varphi \to \varphi + b + c_\mu x^\mu$ and are capable of self-acceleration if they have an inverted sign for the kinetic term. These theories violate the Strong Equivalence Principle, such that black holes (BHs) do not couple to the Galileon field, whereas non-relativistic objects experience a fifth force with strength $\Delta G / G_{\rm N}$ relative to gravity. For galaxies falling down a gradient in the Galileon field, this results in an offset between the centre of the galaxy and its host supermassive BH. We reconstruct the local gravitational and Galileon fields through a suite of constrained N-body simulations (which we dub CSiBORG) and develop a Monte Carlo-based forward model for these offsets on a galaxy-by-galaxy basis. Using the measured offset between the optical centre and active galactic nucleus of 1916 galaxies from the literature, propagating uncertainties in the input quantities and marginalising over an empirical noise model describing astrophysical and observational noise, we constrain the Galileon coupling to be $\Delta G / G_{\rm N} < 0.16$ at $1\sigma$ confidence for Galileons with crossover scale $r_{\rm C} \gtrsim H_0^{-1}$., Comment: 18 pages, 6 figures, 2 tables, accepted for publication in Physical Review D

$f(R)$ is a paradigmatic modified gravity theory that typifies extensions to General Relativity with new light degrees of freedom and hence screened fifth forces between masses. These forces produce observable signatures in galaxy morphology, caused by a violation of the weak equivalence principle due to a differential impact of screening among galaxies' mass components. We compile statistical datasets of two morphological indicators -- offsets between stars and gas in galaxies and warping of stellar disks -- and use them to constrain the strength and range of a thin-shell-screened fifth force. This is achieved by applying a comprehensive set of upgrades to past work (Desmond et al 2018a,b): we construct a robust galaxy-by-galaxy Bayesian forward model for the morphological signals, including full propagation of uncertainties in the input quantities and marginalisation over an empirical model describing astrophysical noise. Employing more stringent data quality cuts than previously we find no evidence for a screened fifth force of any strength $\Delta G/G_\text{N}$ in the Compton wavelength range $0.3-8$ Mpc, setting a $1\sigma$ bound of $\Delta G/G_\text{N}<0.8$ at $\lambda_C=0.3$ Mpc that strengthens to $\Delta G/G_\text{N}<3\times10^{-5}$ at $\lambda_C=8$ Mpc. These are the tightest bounds to date beyond the Solar System by over an order of magnitude. For the Hu-Sawicki model of $f(R)$ with $n=1$ we require a background scalar field value $f_{R0} < 1.4 \times 10^{-8}$, forcing practically all astrophysical objects to be screened. We conclude that this model can have no relevance to astrophysics or cosmology., Comment: 15 pages, 6 figures; minor revision, matches PRD accepted version

We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

We construct a model for the angular power spectrum of the instrumental noise in interferometer networks mapping gravitational wave backgrounds (GWBs) as a function of detector noise properties, network configuration and scan strategy. We use the model to calculate the noise power spectrum for current and future ground-based experiments, as well as for planned space missions. We present our results in a language similar to that used in cosmic microwave background and intensity mapping experiments, and connect the formalism with the sensitivity curves that are common lore in GWB analyses. Our formalism is implemented in a lightweight python module that we make publicly available at https://github.com/damonge/schNell., Comment: Comments welcome. Code available at https://github.com/damonge/schNell

The spatial and temporal discreteness of gravitational wave sources leads to shot noise that may, in some regimes, swamp any attempts at measuring the anisotropy of the gravitational wave background. Cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can alleviate this issue, and potentially recover some of the underlying properties of the gravitational wave background. We quantify the shot noise level and we explicitly show that cross-correlating the gravitational wave background and a galaxy catalog improves the chances of a first detection of the background anisotropy with a gravitational wave observatory operating in the frequency range (10Hz,100Hz), given sufficient sensitivity., Comment: 8 pages, 1 figure

The decay timescales of the quasinormal modes of a massive scalar field have an intriguing behavior: they either grow or decay with increasing angular harmonic numbers $\ell$, depending on whether the mass of the scalar field is small or large. We identify the properties of the effective potential of the scalar field that leads to this behavior and characterize it in detail. If the scalar field is non-minimally coupled, considered here, the scalar quasinormal modes will leak into the gravitational wave signal and will have decaying times that are comparable or smaller than those typical in General Relativity. Hence, these modes could be detectable in the future. Finally, we find that the anomalous behavior in the decay timescales of quasinormal modes is present in a much larger class of models beyond a simple massive scalar field., Comment: Updated to published version

Self-interacting dark matter (SIDM) is an able alternative to collisionless dark matter. If dark matter does have self-interactions, we would expect this to cause a separation between the collisionless stars and the dark matter halo of a galaxy as it falls through a dark matter medium. For stars arranged in a disk, this would generate a U-shaped warp. The magnitude of this warping depends on the SIDM cross section, type of self-interaction, relative velocity of galaxy and background, halo structure, and density of the dark matter medium. In this paper, we set constraints on long-range (light mediator) dark matter self-interaction by means of this signal. We begin by measuring U-shaped warps in $3,213$ edge-on disk galaxies within the Sloan Digital Sky Survey. We then forward-model the expected warp from SIDM on a galaxy-by-galaxy basis by combining models of halo structure, density and velocity field reconstructions, and models for the dark matter interactions. We find no evidence for a contribution to the warps from SIDM. Our constraints are highly dependent on the uncertain velocities of our galaxies: for a normalized Rutherford-like cross section we find $\tilde{\sigma}/m_{\rm{DM}} \lesssim 3\times 10^{-13}~\rm{cm}^2/\rm{g}$ at fixed velocity $v = 300~\rm{km/s}$ -- a bound that scales roughly linearly with increasing $v$. In the appendix we translate these bounds into limits on the momentum transfer cross section, finding $\sigma_T(300~\rm{km/s})/m_{\rm{DM}} \lesssim 0.1~\rm{cm}^2/\rm{g}$. [abridged], Comment: v2: minor modifications + appendix with results in terms of the momentum transfer cross section

The late time acceleration of the Universe can be characterized in terms of an extra, time dependent, component of the universe -- dark energy. The simplest proposal for dark energy is a scalar-tensor theory -- quintessence -- which consists of a scalar field, $\phi$, whose dynamics is solely dictated by its potential, $V(\phi)$. Such a theory can be uniquely characterized by the equation of state of the scalar field energy momentum-tensor. We find the time dependence of the equation of state for a broad family of potentials and, using this information, we propose an analytic prior distribution for the most commonly used parametrization. We show that this analytic prior can be used to accurately predict the distribution of observables for the next generation of cosmological surveys. Including the theoretical priors in the comparison with observations considerably improves the constraints on the equation of state., Comment: V1: 13 pages, 11 figures, comments welcome; V2: 14 pages, 11 figures, minor clarifications added, accepted for publication in PRD

Scale invariant theories of gravity give a compelling explanation to the early and late time acceleration of the Universe. Unlike most scalar-tensor theories, fifth forces are absent and it would therefore seem impossible to distinguish scale invariant gravity from general relativity. We show that the ringdown of a Schwarschild-de Sitter black hole may have a set of massive modes which are characteristic of scale invariant gravity. In principle these new modes can be used to distinguish scale invariant gravity from general relativity. In practice, we discuss the obstacles to generating these new massive modes and their detectability with future gravitational wave experiments but also speculate on their role in Kerr black holes., Comment: 10 pages, submitted to PRD

Weyl (scale) invariant theories of scalars and gravity can generate all mass scales spontaneously. In this paper we study a particularly simple version -- scale invariant $R^2$ gravity -- and show that, during an inflationary period, it leads to fluctuations which, for a particular parameter choice, are almost indistinguishable from normal $R^2$ inflation. Current observations place tight constraints on the primary coupling constant of this theory and predict a tensor to scalar ratio, $0.0033 > r > 0.0026$, which is testable with the next generation of cosmic microwave background experiments.

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value., Comment: 61 pages. Some acknowledgments and references added. This is version-1.1 of an internal collaboration document of LSST-DESC that is being made public and is not planned for submission to a journal

Through numerical simulations of a minimally coupled massive Klein-Gordon scalar field, we show that it is possible to grow hair on a Schwarzschild black hole if one assumes an initial periodically time-varying but spatially homogeneous scalar background. By "hair", we mean a non-trivial profile in the scalar field. We find that this profile emerges on a timescale related to the mass of the black hole, with features related to the mass of the scalar particle. We undertake simulations with and without backreaction on the metric and see that the essential, qualitative features remain consistent. We also contrast the results from higher mass scalars to the case of a low mass with a large Compton wavelength. The results are particularly relevant for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions., Comment: 9 pages, 8 figures, updated to incorporate low mass cases and referee suggestions