We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened f(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z=0.5 and z=1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.

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.

We compare analytic predictions for real and Fourier space two-point statistics for biased tracers from a variety of Lagrangian Perturbation Theory approaches against those from state of the art N-body simulations in $f(R)$ Hu-Sawicki and the nDGP braneworld modified gravity theories. We show that the novel physics of gravitational collapse in scalar tensor theories with the chameleon or the Vainshtein screening mechanism can be effectively factored in with bias parameters analytically predicted using the Peak-Background Split formalism when updated to include the environmental sensitivity of modified gravity theories as well as changes to the halo mass function. We demonstrate that Convolution Lagrangian Perturbation Theory (CLPT) and Standard Perturbation Theory (SPT) approaches provide accurate analytic methods to predict the correlation function and power spectra, respectively, for biased tracers in modified gravity models and are able to characterize both the BAO, power-law and small scale regimes needed for upcoming galaxy surveys such as DESI, Euclid, LSST and WFIRST.

Models that seek to explain cosmic acceleration through modifications to General Relativity (GR) evade stringent Solar System constraints through a restoring, screening mechanism. Down-weighting the high density, screened regions in favor of the low density, unscreened ones offers the potential to enhance the amount of information carried in such modified gravity models. In this work, we assess the performance of a new "marked" transformation and perform a systematic comparison with the clipping and logarithmic transformations, in the context of $\Lambda$CDM and the symmetron and $f(R)$ modified gravity models. Performance is measured in terms of the fractional boost in the Fisher information and the signal-to-noise ratio (SNR) for these models relative to the statistics derived from the standard density distribution. We find that all three statistics provide improved Fisher boosts over the basic density statistics. The model parameters for the "marked" and clipped transformation that best enhance signals and the Fisher boosts are determined. We also show that the mark is useful both as a Fourier and real space transformation; a marked correlation function also enhances the SNR relative to the standard correlation function, and can on mildly non-linear scales show a significant difference between the $\Lambda$CDM and the modified gravity models. Our results demonstrate how a series of simple analytical transformations could dramatically increase the predicted information extracted on deviations from GR, from large-scale surveys, and give the prospect for a potential detection much more feasible.

We implement an adaptation of the COLA approach, a hybrid scheme that combines Lagrangian perturbation theory with an N-body approach, to model non-linear collapse in chameleon and symmetron modified gravity models. Gravitational screening is modeled effectively through the attachment of a suppression factor to the linearized Klein-Gordon equations. The adapted COLA approach is benchmarked, with respect to an N-body code both for the $\Lambda$CDM scenario and for the modified gravity theories. It is found to perform well in the estimation of the dark matter power spectra, with consistency of 1 % to $k\sim2.5$ h/Mpc. Redshift space distortions are shown to be effectively modeled through a Lorentzian parameterization with a velocity dispersion fit to the data. We find that COLA performs less well in predicting the halo mass functions, but has consistency, within $1\sigma$ uncertainties of our simulations, in the relative changes to the mass function induced by the modified gravity models relative to $\Lambda$CDM. The results demonstrate that COLA, proposed to enable accurate and efficient, non-linear predictions for $\Lambda$CDM, can be effectively applied to a wider set of cosmological scenarios, with intriguing properties, for which clustering behavior needs to be understood for upcoming surveys such as LSST, DESI, Euclid and WFIRST.