We investigate the nonlinear growth of large-scale structure in the generalised Proca theory, in which a self-interacting massive vector field plays the role of driving the acceleration of the cosmic expansion. Focusing to the Proca Lagrangian at cubic order -- the cubic vector Galileon model -- we derive the simplified equations for gravity as well as the longitudinal and transverse modes of the vector field under the weak-field and quasi-static approximations, and implement them in a modified version of the ECOSMOG N-body code. Our simulations incorporate the Vainshtein screening effect, which reconciles the fifth force propagated by the longitudinal mode of the cubic vector Galileon model with local tests of gravity. The results confirm that for all scales probed by the simulation, the transverse mode has a negligible impact on structure formation in a realistic cosmological setup. It is well known that in this model the strength of the fifth force is controlled by a free model parameter, which we denote as beta_3. By running a suite of cosmological simulations for different values of beta_3, we show that this parameter also determines the effectiveness of the Vainshtein screening. The model behaves identically to the cubic scalar Galileon for beta_3 going to zero, in which the fifth force is strong in unscreened regions but is efficiently screened in high-density regions. In the opposite limit, beta_3 going to infinity, the model approaches its `quintessence' counterpart, which has a vanishing fifth force but a modified expansion history compared to LambdaCDM. This endows the model with rich phenomenology, which will be investigated in future works.

We carry out ``full-physics'' hydrodynamical simulations of galaxy formation in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld model using a new modified version of the {\sc Arepo} code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) which represent, respectively, weak and moderate departures from GR, in boxes of sizes $62\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$ using $2\times512^3$ dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to $\sim 25$ per cent on large scales. Similar to their impact in the standard cosmological model ($\Lambda$CDM), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to $\sim15$ per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above $\sim0.1$ times the halo radius.

We present an analysis of the matter, halo and galaxy clustering in $f(R)$-gravity employing the SHYBONE full-physics hydrodynamical simulation suite. Our analysis focuses on the interplay between baryonic feedback and $f(R)$-gravity in the matter power spectrum, the matter and halo correlation functions, the halo and galaxy-host-halo mass function, the subhalo and satellite-galaxy count and the correlation function of the stars in our simulations. Our studies of the matter power spectrum in full physics simulations in $f(R)$-gravity show, that it will be very difficult to derive accurate fitting formulae for the power spectrum enhancement in $f(R)$-gravity which include baryonic effects. We find that the enhancement of the halo mass function due to $f(R)$-gravity and its suppression due to feedback effects do not show significant back-reaction effects and can thus be estimated from independent GR-hydro and $f(R)$ dark matter only simulations. Our simulations furthermore show, that the number of subhalos and satellite-galaxies per halo is not significantly affected by $f(R)$-gravity. Low mass halos are nevertheless more likely to be populated by galaxies in $f(R)$-gravity. This suppresses the clustering of stars and the galaxy correlation function in the theory compared to standard cosmology.

Future 21 cm intensity mapping surveys such as SKA can provide precise information on the spatial distribution of the neutral hydrogen (HI) in the post-reionization epoch. This information will allow to test the standard LCDM paradigm and with that the nature of gravity. In this work, we employ the SHYBONE simulations, which model galaxy formation in f(R) modified gravity using the IllustrisTNG model, to study the effects of modified gravity on HI abundance and power spectra. We find that the enhanced growth low-mass dark matter halos experience in f(R) gravity at high redshifts alters the HI power spectrum and can be observable through 21 cm intensity mapping. Our results suggest that the HI power spectrum is suppressed by 13\% on scales $k\lesssim2\,h\,\mathrm{Mpc}^{-1}$ at z=2 for F6, a f(R) model which passes most observational constraints. We show that this suppression can be detectable by SKA1-MID with 1000 hours of exposure time, making HI clustering a novel test of gravity at high redshift.

We have carried out a set of cosmological hydrodynamical simulations that follow galaxy formation in $f(R)$ modified gravity models. Our simulations employ the Illustris-TNG full physics model and a new modified gravity solver in the AREPO code. For the first time we are able to investigate the degeneracy in the matter power spectrum between the effects of $f(R)$-gravity and feedback from active galactic nuclei (AGN), and the imprint of modified gravity on the properties of galaxies and on the distribution of dark matter, gas and stars in the universe. $f(R)$-gravity has an observable effect on the neutral hydrogen power spectrum at high redshift at a level of 20%. For both the F6 and F5 models, this is significantly larger than the predicted errors for the SKA1-MID survey, making this probe a powerful test of gravity on large scales. A similar effect is present in the power spectrum of the stars at high redshift. We also show that rotationally supported disc galaxies can form in $f(R)$-gravity, even in the partially screened regime. Our simulations indicate that there might be more disc galaxies in F6 compared to GR, and fewer in F5. Finally, we show that the back reaction between AGN feedback and modified gravity in the matter power spectrum is not important in the F6 model but has a sizeable effect in F5.

We present a novel fitting formula for the halo concentration enhancement in chameleon $f(R)$ gravity relative to General Relativity (GR). The formula is derived by employing a large set of $N$-body simulations of the Hu-Sawicki $f(R)$ model which cover a wide range of model and cosmological parameters, resolutions and simulation box sizes. The complicated dependence of the concentration on halo mass $M$, redshift $z$, and the $f(R)$ and cosmological parameters can be combined into a simpler form that depends only on a rescaled mass $M/10^{p_2}$, with $p_2\equiv1.5\log_{10}\left[|{\bar{f}_R(z)}|/(1+z)\right]+21.64$ and $\bar{f}_R(z)$ the background scalar field at $z$, irrespective of the $f(R)$ model parameter. Our fitting formula can describe the concentration enhancement well for redshifts $z\leq3$, nearly 7 orders of magnitude in $M/10^{p_2}$ and five decades in halo mass. This is part of a series of works which aims to provide a general framework for self-consistent and unbiased tests of gravity using present and upcoming galaxy cluster surveys. The fitting formula, which is the first quantitative model for the concentration enhancement due to chameleon type modified gravity, is an important part in this framework and will allow continuous exploration of the parameter space. It can also be used to model other statistics such as the matter power spectrum.

Current constraints on gravity are relatively weak on galactic and intergalactic scales. Screened modified gravity models can exhibit complex behaviour there without violating stringent tests of gravity within our Solar System. They might hence provide viable extensions of the theory of gravity. Here, we use galaxy kinematics to constrain screened modified gravity models. We focus on chameleon $f(R)$ gravity and predict its impact on galaxy rotation curves and radial acceleration relations. This is achieved by post-processing state-of-the-art galaxy formation simulations from the \textsc{auriga project}, using the \textsc{mg-gadget} code. For a given galaxy, the surface dividing screened and un-screened regions adopts an oblate shape, reflecting the disc morphology of the galaxy's mass distribution. At the `screening radius'---the radius at which screening is triggered in the disc plane---characteristic `upturns' are present in both rotation curves and radial acceleration relations. The locations of these features depend on various factors, such as the galaxy mass, the concentration of the density profile and the value of the background field amplitude $f_{R0}$. Self-screening of stars and environmental screening also play a role. For Milky Way-size galaxies, we find that a model with $|f_{R0}|=10^{-7}$ results in rotation curves that are indistinguishable from $\Lambda$CDM, while for $|f_{R0}| \geq 2 \times 10^{-6}$ the simulated galaxies are entirely unscreened, violating Solar System constraints. For intermediate values, distinct upturns are present. With a careful statistical analysis of existing samples of observed rotation curves, including lower mass objects, constraints on $f(R)$ gravity with a sensitivity down to $|f_{R0}|\sim10^{-7}$ should be possible.

We introduce a set of four very high resolution cosmological simulations for exploring $f(R)$-gravity, with $2048^3$ particles in $768\,h^{-1}\textrm{Mpc}$ and $1536\,h^{-1}\textrm{Mpc}$ simulation boxes, both for a $|\overline{f_{R0}}| = 10^{-5}$ model and a $\Lambda$CDM comparison universe, making the set the largest simulations of $f(R)$-gravity to date. In order to mimic real observations, the simulations include a continuous 2D and 3D lightcone output which is dedicated to study lensing and clustering statistics in modified gravity. In this work, we present a detailed analysis and resolution study for the matter power spectrum in $f(R)$-gravity over a wide range of scales. We also analyse the angular matter power spectrum and lensing convergence on the lightcone. In addition, we investigate the impact of modified gravity on the halo mass function, matter and halo auto-correlation functions, linear halo bias and the concentration-mass relation. We find that the impact of $f(R)$-gravity is generally larger on smaller scales and increases with decreasing redshift. Comparing our simulations to state-of-the-art hydrodynamical simulations we confirm a degeneracy between $f(R)$-gravity and baryonic feedback in the matter power spectrum on small scales, but also find that scales around $k = 1\, h\, {\rm Mpc}^{-1}$ are promising to distinguish both effects. The lensing convergence power spectrum is increased in $f(R)$-gravity. Interestingly available numerical fits are in good agreement overall with our simulations for both standard and modified gravity, but tend to overestimate their relative difference on non-linear scales by few percent. We also find that the halo bias is lower in $f(R)$-gravity compared to general relativity, whereas halo concentrations are increased for unscreened halos.

We propose a new framework for testing gravity using cluster observations, which aims to provide an unbiased constraint on modified gravity models from Sunyaev Zel'dovich (SZ) and X-ray cluster counts and the cluster gas fraction, among other possible observables. Focusing on a popular $ f(R) $ model of gravity, we propose a novel procedure to recalibrate mass scaling relations from $ \Lambda $CDM to $ f(R) $ gravity for SZ and X-ray cluster observables. We find that the complicated modified gravity effects can be simply modelled as a dependence on a combination of the background scalar field and redshift, $ f_R(z)/(1+z) $, regardless of the $ f(R) $ model parameter. By employing a large suite of N-body simulations, we demonstrate that a theoretically derived tanh fitting formula is in excellent agreement with the dynamical mass enhancement of dark matter haloes for a large range of background field parameters and redshifts. Our framework is sufficiently flexible to allow for tests of other models and inclusion of further observables. The one-parameter description of the dynamical mass enhancement can have important implications on the theoretical modelling of observables and on practical tests of gravity.

The first transiting planetesimal orbiting a white dwarf was recently detected in K2 data of WD1145+017 and has been followed up intensively. The multiple, long, and variable transits suggest the transiting objects are dust clouds, probably produced by a disintegrating asteroid. In addition, the system contains circumstellar gas, evident by broad absorption lines, mostly in the u'-band, and a dust disc, indicated by an infrared excess. Here we present the first detection of a change in colour of WD1145+017 during transits, using simultaneous multi-band fast-photometry ULTRACAM measurements over the u'g'r'i'-bands. The observations reveal what appears to be 'bluing' during transits; transits are deeper in the redder bands, with a u'-r' colour difference of up to ~-0.05 mag. We explore various possible explanations for the bluing. 'Spectral' photometry obtained by integrating over bandpasses in the spectroscopic data in- and out-of-transit, compared to the photometric data, shows that the observed colour difference is most likely the result of reduced circumstellar absorption in the spectrum during transits. This indicates that the transiting objects and the gas share the same line-of-sight, and that the gas covers the white dwarf only partially, as would be expected if the gas, the transiting debris, and the dust emitting the infrared excess, are part of the same general disc structure (although possibly at different radii). In addition, we present the results of a week-long monitoring campaign of the system.

We investigate the impact of f(R) modified gravity on the internal properties of Milky Way sized dark matter halos in a set of cosmological zoom simulations of seven halos from the Aquarius suite, carried out with our code MG-GADGET in the Hu & Sawicki f(R) model. Also, we calculate the fifth forces in ideal NFW-halos as well as in our cosmological simulations and compare them against analytic model predictions for the fifth force inside spherical objects. We find that these theoretical predictions match the forces in the ideal halos very well, whereas their applicability is somewhat limited for realistic cosmological halos. Our simulations show that f(R) gravity significantly affects the dark matter density profile of Milky Way sized objects as well as their circular velocities. In unscreened regions, the velocity dispersions are increased by up to 40% with respect to LCDM for viable f(R) models. This difference is larger than reported in previous works. The Solar circle is fully screened in $f_{R0} = -10^{-6}$ models for Milky Way sized halos, while this location is unscreened for slightly less massive objects. Within the scope of our limited halo sample size, we do not find a clear dependence of the concentration parameter of dark matter halos on $f_{R0}$.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

In this work, we analyze the Lyman-$\alpha$ forest in cosmological hydrodynamical simulations of chameleon-type f(R) gravity with the goal to assess whether the impact of such models is detectable in absorption line statistics. We carry out a set of hydrodynamical simulations with the cosmological simulation code MG-GADGET, including star formation and cooling effects, and create synthetic Lyman-$\alpha$ absorption spectra from the simulation outputs. We statistically compare simulations with f(R) and ordinary general relativity, focusing on flux probability distribution functions (PDFs) and flux power-spectra, an analysis of the column density and line width distributions, as well as the matter power spectrum. We find that the influence of f(R) gravity on the Lyman-$\alpha$ forest is rather small. Even models with strong modifications of gravity, like $|\bar{f}_{R0}| = 10^{-4}$, do not change the statistical Lyman-$\alpha$ properties by more than 10%. The column density and line width distributions are hardly affected at all. It is therefore not possible to get competitive constraints on the background field $f_R$ using current observational data. An improved understanding of systematics in the observations and a more accurate modeling of the baryonic/radiative physics would be required to allow this in the future. The impact of f(R) on the matter power spectrum in our results is consistent with previous works.

We investigate the impact of chameleon-type f(R) gravity models on the properties of galaxy clusters and groups. Our f(R) simulations follow for the first time also the hydrodynamics of the intracluster and intragroup medium. This allows us to assess how f(R) gravity alters the X-ray scaling relations of clusters and how hydrostatic and dynamical mass estimates are biased when modifications of gravity are ignored in their determination. We find that velocity dispersions and intracluster medium temperatures are both increased by up to 1/3 in f(R) gravity in low-mass halos, while the difference disappears in massive objects. The mass scale of the transition depends on the background value f_R0 of the scalar degree of freedom. These changes in temperature and velocity dispersion alter the mass-temperature and X-ray luminosity-temperature scaling relations and bias dynamical and hydrostatic mass estimates that do not explicitly account for modified gravity towards higher values. Recently, a relative enhancement of X-ray compared to weak lensing masses was found by the Planck Collaboration (2013). We demonstrate that an explanation for this offset may be provided by modified gravity and the associated bias effects, which interestingly are of the required size. Finally, we find that the abundance of subhalos at fixed cluster mass is only weakly affected by f(R) gravity.