We consider how upcoming photometric large scale structure surveys can be optimized to measure the properties of dark energy and possible cosmic scale modifications to General Relativity in light of realistic astrophysical and instrumental systematic uncertainities. In particular we include flexible descriptions of intrinsic alignments, galaxy bias and photometric redshift uncertainties in a Fisher Matrix analysis of shear, position and position-shear correlations, including complementary cosmological constraints from the CMB. We study the impact of survey tradeoffs in depth versus breadth, and redshift quality. We parameterise the results in terms of the Dark Energy Task Force figure of merit, and deviations from General Relativity through an analagous Modified Gravity figure of merit. We find that intrinsic alignments weaken the dependence of figure of merit on area and that, for a fixed observing time, a fiducial Stage IV survey plateaus above roughly 10,000deg2 for DE and peaks at about 5,000deg2 as the relative importance of IAs at low redshift penalises wide, shallow surveys. While reducing photometric redshift scatter improves constraining power, the dependence is shallow. The variation in constraining power is stronger once IAs are included and is slightly more pronounced for MG constraints than for DE. The inclusion of intrinsic alignments and galaxy position information reduces the required prior on photometric redshift accuracy by an order of magnitude for both the fiducial Stage III and IV surveys, equivalent to a factor of 100 reduction in the number of spectroscopic galaxies required to calibrate the photometric sample.

We consider the impact of key astrophysical and measurement systematics on constraints on dark energy and modifications to gravity on cosmic scales. We focus on upcoming photometric "Stage III" and "Stage IV" large scale structure surveys such as DES, SuMIRe, Euclid, LSST and WFIRST. We illustrate the different redshift dependencies of gravity modifications compared to intrinsic alignments, the main astrophysical systematic. The way in which systematic uncertainties, such as galaxy bias and intrinsic alignments, are modelled can change dark energy equation of state and modified gravity figures of merit by a factor of four. The inclusion of cross-correlations of cosmic shear and galaxy position measurements helps reduce the loss of constraining power from the lensing shear surveys. When forecasts for Planck CMB and Stage IV surveys are combined, constraints on the dark energy equation of state and modified gravity model are recovered, relative to those from shear data with no systematic uncertainties, if fewer than 36 free parameters in total are used to describe the galaxy bias and intrinsic alignment models as a function of scale and redshift. To facilitate future investigations, we also provide a fitting function for the matter power spectrum arising from the phenomenological modified gravity model we consider.

We consider the cosmological constraints on theories in which there exists a nontrivial coupling between the dark matter sector and the sector responsible for the acceleration of the universe, in light of the most recent supernovae, large scale structure and cosmic microwave background data. For a variety of models, we show that the strength of the coupling of dark matter to a quintessence field is constrained to be less than 7% of the coupling to gravity. We also show that long range interactions between fermionic dark matter particles mediated by a light scalar with a Yukawa coupling are constrained to be less than 5% of the strength of gravity at a distance scale of 10 Mpc. We show that all of the models we consider are quantum mechanically weakly coupled, and argue that some other models in the literature are ruled out by quantum mechanical strong coupling.

Theoretical differences in the growth of structure offer the possibility that we might distinguish between modified gravity theories of dark energy and \LambdaCDM. A significant impediment to applying current and prospective large scale galaxy and weak lensing surveys to this problem is that, while the mildly nonlinear regime is important, there is a lack of numerical simulations of nonlinear growth in modified gravity theories. A major question exists as to whether existing analytical fits, created using simulations of standard gravity, can be confidently applied. In this paper we address this, presenting results of N-body simulations of a variety of models where gravity is altered including the Dvali, Gabadadze and Porrati model. We consider modifications that alter the Poisson equation and also consider the presence of anisotropic shear stress that alters how particles respond to the gravitational potential gradient. We establish how well analytical fits of the matter power spectrum by Peacock and Dodds and Smith et al. are able to predict the nonlinear growth found in the simulations from z=50 up to today, and also consider implications for the weak lensing convergence power spectrum. We find that the analytical fits provide good agreement with the simulations, being within 1\sigma of the simulation results for cases with and without anisotropic stress and for scale-dependent and independent modifications of the Poisson equation. No strong preference for either analytical fit is found.