Search

Stage IV large scale structure surveys are promising probes of gravity on cosmological scales. Due to the vast model-space in the modified gravity literature, model-independent parameterisations represent useful and scalable ways to test extensions of $\Lambda$CDM. In this work we use a recently validated approach of computing the non-linear $3\times 2$pt observables in modified gravity models with a time-varying effective gravitational constant $\mu$ and a gravitational slip $\eta$ that is binned in redshift to produce Fisher forecasts for an LSST Y10-like survey. We also include in our modelling an effective nulling scheme for weak-lensing by applying the BNT transformation that localises the weak-lensing kernel enabling well-informed scale cuts. We show that the combination of improved non-linear modelling and better control of the scales that are modelled/cut yields high precision constraints on the cosmological and modified gravity parameters. We find that 4 redshift bins for $\mu$ of width corresponding to equal incremental $\Lambda$CDM growth is optimal given the state-of-the-art modelling and show how the BNT transformation can be used to mitigate the impact of small-scale systematic effects, such as baryonic feedback., Comment: 27 pages, 12 figures A few typos corrected and a couple of small changes made to the text to improve presentation of results, added missing reference. Comments welcome!

The Parameterised Post-Newtonian (PPN) approach is the default framework for performing precision tests of gravity in nearby astrophysical systems. In recent works we have extended this approach for cosmological applications, and in this paper we use observations of the anisotropies in the Cosmic Microwave Background to constrain the time variation of the PPN parameters $\alpha$ and $\gamma$ between last scattering and the present day. We find their time-averages over cosmological history should be within $\sim 20\%$ of their values in GR, with $\bar{\alpha}=0.89^{+0.08}_{-0.09}$ and $\bar{\gamma}=0.90^{+0.07}_{-0.08}$ at the $68\%$ confidence level. We also constrain the time derivatives of these parameters, and find that their present-day values should be within a factor of two of the best Solar System constraints. Many of these results have no counter-part from Solar System observations, and are entirely new constraints on the gravitational interaction. In all cases, we find that the data strongly prefer $\bar{\alpha}\simeq \bar{\gamma}$, meaning that observers would typically find local gravitational physics to be compatible with GR, despite considerable variation of $\alpha$ and $\gamma$ being allowed over cosmic history. This study lays the groundwork for future precision tests of gravity that combine observations made over all cosmological and astrophysical scales of length and time., Comment: Comments welcome

Several prominent forthcoming Cosmic Microwave Background polarisation experiments will employ a Continuously Rotating Half-Wave Plate (CRHWP), the primary purpose of which is to mitigate instrumental systematic effects on relatively large angular scales, where the $B$-mode polarisation signal generated by primordial gravitational waves is expected to peak. The use of a CRHWP necessitates demodulating the time-ordered data during the early stages of data processing. The standard approach is to ``lock in'' on the polarisation signal using the known polarisation modulation frequency and then use Fourier techniques to filter out the remaining unwanted components in the data. However, an alternative less well-studied option is to incorporate the demodulation directly into the map-making step. Using simulations, we compare the performance of these two approaches to determine which is most effective for $B$-mode signal recovery. Testing the two techniques in multiple experimental scenarios, we find that the lock-in technique performs best over the full multipole range explored. However, for the recovery of the largest angular scales ($\ell < 100$) we find essentially no difference in the recovery of the signal between the lock-in and map-making approaches, suggesting that a parallel analysis based on the latter approach could represent a powerful consistency check for primordial $B$-mode experiments employing a CRHWP. We also investigate the impact of a detector-differencing step, implemented prior to demodulation, finding that in most scenarios it makes no difference whether differencing is used or not. However, analysing detectors individually allows the point at which information from multiple detectors is combined to be moved to later stages in the analysis pipeline. This presents alternative options for dealing with instrumental systematic effects that are not mitigated by the CRHWP.

Model-independent tests of gravity with cosmology are important when testing extensions to the standard cosmological model. To maximise the impact of these tests one requires predictions for the matter power spectrum on non-linear scales. In this work we validate the \texttt{ReACT} approach to the non-linear matter power spectrum against a suite of phenomenological modified gravity N-body simulations with a time-varying gravitational constant, covering a wider range of parameter space than previously examined. This vanilla application of \texttt{ReACT} has limited range and precision due to the different concentration-mass relation $c(M)$ that occurs when gravity is modified. We extend this approach with a fitting function for a modified concentration-mass relation, allowing for accurate (1$\%$) computation of the matter power spectrum up $k=2\,h\,{\rm Mpc}^{-1}$ across a substantial range of parameter space. This fitting function allows precision model-independent tests of modified gravity to be carried out using the data from upcoming large scale structure surveys., Comment: 28 pages, 13 figures Comments Welcome!

General relativity manifests very similar equations in different regimes, notably in large scale cosmological perturbation theory, non-linear cosmological structure formation, and in weak field galactic dynamics. The same is not necessarily true in alternative gravity theories, in particular those that possess MONDian behaviour ("relativistic extensions" of MOND). In these theories different regimes are typically studied quite separately, sometimes even with the freedom in the theories chosen differently in different regimes. If we wish to properly and fully test complete cosmologies containing MOND against the $\Lambda$CDM paradigm then we need to understand cosmological structure formation on all scales, and do so in a coherent and consistent manner. We propose a method for doing so and apply it to generalised Einstein-Aether theories as a case study. We derive the equations that govern cosmological structure formation on all scales in these theories and show that the same free function (which may contain both Newtonian and MONDian branches) appears in the cosmological background, linear perturbations, and non-linear cosmological structure formation. We show that MONDian behaviour on galactic scales does not necessarily result in MONDian behaviour on cosmological scales, and for MONDian behaviour to arise cosmologically, there will be no modification to the Friedmann equations governing the evolution of the homogeneous cosmological background. We comment on how existing N-body simulations relate to complete and consistent generalised Einstein-Aether cosmologies. The equations derived in this work allow consistent cosmological N-body simulations to be run in these theories whether or not MONDian behaviour manifests on cosmological scales., Comment: Typo corrected on page 14; submitted to JCAP

Parameterised Post-Newtonian Cosmology (PPNC) is a theory-agnostic framework for testing gravity in cosmology, which connects gravitational physics on small and large scales in the Universe. It is a direct extension of the Parameterised Post-Newtonian (PPN) approach to testing gravity in isolated astrophysical systems, and therefore allows constraints on gravity from vastly different physical regimes to be compared and combined. We investigate the application of this framework to a class of example scalar-tensor theories of gravity in order to verify theoretical predictions, and to investigate for the first time the scale-dependence of the gravitational couplings that appear within its perturbation equations. In doing so, we evaluate the performance of some simple interpolating functions in the transition region between small and large cosmological scales, as well as the uncertainties that using such functions would introduce into the calculation of observables. We find that all theoretical predictions of the PPNC framework are verified to high accuracy in the relevant regimes, and that simple interpolating functions perform well (but not perfectly) between these regimes. This study is an important step towards being able to use the PPNC framework to analyse cosmological datasets, and to thereby test if/how the gravitational interaction has changed as the Universe has evolved., Comment: 28 pages, 16 figures

We present approaches to quickly simulate systematics affecting CMB observations, including the effects of the scanning strategy. Using summary properties of the scan we capture features of full time ordered data (TOD) simulations, allowing maps and power spectra to be generated at much improved speed for a number of systematics - the cases we present experienced speed ups of 3-4 orders of magnitude when implementing the map-based approaches. We demonstrate the effectiveness of the approaches at capturing the salient features of the scan by directly comparing to full TOD simulations - seeing agreement at sub-percent levels of accuracy. We simulate the effects of differential gain, pointing, and ellipticity to show the effectiveness of the approaches, but note that one could extend these techniques to other systematics. We finally show how to apply these fast map-based simulations of systematic effects to a full focal plane showing their ability to incorporate thousands of detectors as seen in modern CMB experiments., Comment: 17 pages, 13 figures, 1 table, prepared for submission to MNRAS

Upcoming cosmological intensity mapping surveys will open new windows on the Universe, but they must first overcome a number of significant systematic effects, including polarization leakage. We present a formalism that uses scan strategy information to model the effect of different instrumental systematics on the recovered cosmological intensity signal for `single-dish' (autocorrelation) surveys. This modelling classifies different systematics according to their spin symmetry, making it particularly relevant for dealing with polarization leakage. We show how to use this formalism to calculate the expected contamination from different systematics as a function of the scanning strategy. Most importantly, we show how systematics can be disentangled from the intensity signal based on their spin properties via map-making. We illustrate this, using a set of toy models, for some simple instrumental systematics, demonstrating the ability to significantly reduce the contamination to the observed intensity signal. Crucially, unlike existing foreground removal techniques, this approach works for signals that are non-smooth in frequency, e.g. polarized foregrounds. These map-making approaches are simple to apply and represent an orthogonal and complementary approach to existing techniques for removing systematics from upcoming 21cm intensity mapping surveys., Comment: 19 pages, 14 Figures, 2 Tables, published in MNRAS

Model-independent constraints on modified gravity models hitherto exist mainly on linear scales. A recently developed formalism presented a consistent parameterisation that is valid on all scales. Using this approach, we perform model-independent modified gravity $N$-body simulations on all cosmological scales with a time-dependent $\mu$. We present convergence tests of our simulations, and we examine how well existing fitting functions reproduce the non-linear matter power spectrum of the simulations. We find that although there is a significant variation in the accuracy of all of the fitting functions over the parameter space of our simulations, the ReACT framework delivers the most consistent performance for the matter power spectrum. We comment on how this might be improved to the level required for future surveys such as Euclid and the Vera Rubin Telescope (LSST). We also show how to compute weak-lensing observables consistently from the simulated matter power spectra in our approach, and show that ReACT also performs best when fitting the weak-lensing observables. This paves the way for a full model-independent test of modified gravity using all of the data from such upcoming surveys., Comment: 32 pages, 14 Figures; Comments welcome

Instrumental systematics need to be controlled to high precision for upcoming Cosmic Microwave Background (CMB) experiments. The level of contamination caused by these systematics is often linked to the scan strategy, and scan strategies for satellite experiments can significantly mitigate these systematics. However, no detailed study has been performed for ground-based experiments. Here we show that under the assumption of constant elevation scans (CESs), the ability of the scan strategy to mitigate these systematics is strongly limited, irrespective of the detailed structure of the scan strategy. We calculate typical values and maps of the quantities coupling the scan to the systematics, and show how these quantities vary with the choice of observing elevations. These values and maps can be used to calculate and forecast the magnitude of different instrumental systematics without requiring detailed scan strategy simulations. As a reference point, we show that inclusion of even a single boresight rotation angle significantly improves over sky rotation alone for mitigating these systematics. A standard metric for evaluating cross-linking is related to one of the parameters studied in this work, so a corollary of our work is that the cross-linking will suffer from the same CES limitations and therefore upcoming CMB surveys will unavoidably have poorly cross-linked regions if they use CESs, regardless of detailed scheduling choices. Our results are also relevant for non-CMB surveys that perform constant elevation scans and may have scan-coupled systematics, such as intensity mapping surveys., Comment: Narrative adjusted to make scope of paper and results clearer; version published in the Open Journal of Astrophysics

The CMB $B$-mode polarisation signal -- both the primordial gravitational wave signature and the signal sourced by lensing -- is subject to many contaminants from systematic effects. Of particular concern are systematics that result in mixing of signals of different ``spin'', particularly leakage from the much larger spin-0 intensity signal to the spin-2 polarisation signal. We present a general formalism, which can be applied to arbitrary focal plane setups, that characterises signals in terms of their spin. We provide general expressions to describe how spin-coupled signals observed by the detectors manifest at map-level, in the harmonic domain, and in the power spectra, focusing on the polarisation spectra -- the signals of interest for upcoming CMB surveys. We demonstrate the presence of a previously unidentified cross-term between the systematic and the intrinsic sky signal in the power spectrum, which in some cases can be the dominant source of contamination. The formalism is not restricted to intensity to polarisation leakage but provides a complete elucidation of all leakage including polarisation mixing, and applies to both full and partial (masked) sky surveys, thus covering space-based, balloon-borne, and ground-based experiments. Using a pair-differenced setup, we demonstrate the formalism by using it to completely characterise the effects of differential gain and pointing systematics, incorporating both intensity leakage and polarisation mixing. We validate our results with full time ordered data simulations. Finally, we show in an Appendix that an extension of simple binning map-making to include additional spin information is capable of removing spin-coupled systematics during the map-making process., Comment: Minimal revisions - some additional references added, typos fixed, etc. 31 pages, 8 figures, 1 table. Accepted by MNRAS

The cosmological phenomenology of gravity is typically studied in two limits: relativistic perturbation theory (on large scales) and Newtonian gravity (required for smaller, non-linear, scales). Traditional approaches to model-independent modified gravity are based on perturbation theory, so do not apply on non-linear scales. Future surveys such as Euclid will produce significant data on both linear and non-linear scales, so a new approach is required to constrain model-independent modified gravity by simultaneously using all of the data from these surveys. We use the higher order equations from the post-Friedmann approach to derive a single set of "simple 1PF" (first post-Friedmann) equations that apply in both the small scale and large scale limits, and we examine the required conditions for there to be no intermediate regime, meaning that these simple equations are valid on all scales. We demonstrate how the simple 1PF equations derived here can be used as a model-independent framework for modified gravity that applies on all cosmological scales, and we present an algorithm for determining which modified gravity theories are subsumed under this approach. This modified gravity framework provides a rigorous approach to phenomenological N-body simulations, and paves the way to consistently using all of the data from upcoming surveys to constrain modified gravity in a model-independent fashion., Comment: Minor updates to match version published in PRD, primarily the addition of section V.C.6. Chosen as editor's suggestion

We perform the first test of dark matter (DM) stress-energy evolution through cosmic history, using cosmic microwave background measurements supplemented with baryon acoustic oscillation data and the Hubble Space Telescope key project data. We constrain the DM equation of state (EoS) in 8 redshift bins, and its sound speed and (shear) viscosity in 9 redshift bins, finding no convincing evidence for non-$\Lambda$CDM values in any of the redshift bins. Despite this enlarged parameter space, the sound speed and viscosity are constrained relatively well at late times (due to the inclusion of CMB lensing), whereas the EoS is most strongly constrained around recombination. These results constrain for the first time the level of "coldness" required of DM across various cosmological epochs at both the background and perturbative levels. We show that simultaneously allowing time dependence for both the EoS and sound speed parameters shifts the posterior of the DM abundance before recombination to a higher value, while keeping the present day DM abundance similar to the $\Lambda$CDM value. This shifts the posterior for the present day Hubble constant compared to $\Lambda$CDM, suggesting that DM with time-dependent parameters is well-suited to explore possible solutions to persistent tensions within the $\Lambda$CDM model. We perform a detailed comparison with our previous study involving a vanishing sound speed and viscosity using the same datasets in order to explain the physical mechanism behind these shifts., Comment: 20 pages, 10 figures, 5 tables; updated to match version published in PRD

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers., Comment: 20 pages, 14 figures, 4 tables. Accepted in MNRAS. Links to Data Release catalogues will be updated when available upon publication. Before this, catalogues are available on request from the authors

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured., Comment: 25 pages, 18 figures, 3 tables. Accepted in MNRAS

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4., Comment: Astro2020 Decadal Project Whitepaper. arXiv admin note: text overlap with arXiv:1808.07445

Future CMB experiments will require exquisite control of systematics in order to constrain the $B$-mode polarisation power spectrum. One class of systematics that requires careful study is instrumental systematics. The potential impact of such systematics is most readily understood by considering analysis pipelines based on pair differencing. In this case, any differential gain, pointing or beam ellipticity between the two detectors in a pair can result in intensity leakage into the $B$-mode spectrum, which needs to be controlled to a high precision due to the much greater magnitude of the total intensity signal as compared to the $B$-mode signal. One well known way to suppress such systematics is through careful design of the scan-strategy, in particular making use of any capability to rotate the instrument about its pointing (boresight) direction. Here, we show that the combination of specific choices of such partial boresight rotation angles with redundancies present in the scan strategy is a powerful approach for suppressing systematic effects. This mitigation can be performed in analysis in advance of map-making and, in contrast to other approaches (e.g. deprojection or filtering), results in no signal loss. We demonstrate our approach explicitly with time ordered data simulations relevant to next-generation ground-based CMB experiments, using deep and wide scan strategies appropriate for experiments based in Chile. These simulations show a reduction of multiple orders of magnitude in the spurious $B$-mode signal arising from differential gain and differential pointing systematics., Comment: 6 figures. Updated to match published version with a few minor clarifications; the main change is a brief discussion of bandpass mismatch

Constraints on the properties of the cosmological dark matter have previously been obtained in a model-independent fashion using the Generalised Dark Matter (GDM) framework. Here we extend that work in several directions: We consider the inclusion of WiggleZ matter power spectrum data, and show that this improves the constraints on the two perturbative GDM parameters, $c^2_s$ and $c^2_\text{vis}$, by a factor of 3, for a conservative choice of wavenumber range. A less conservative choice can yield an improvement of up to an order of magnitude compared to previous constraints. In order to examine the robustness of this result we develop a GDM halo model to explore how non-linear structure formation could proceed in this framework, since currently GDM has only been defined perturbatively and only linear theory has been used when generating constraints. We then examine how the halo model affects the constraints obtained from the matter power spectrum data. The less-conservative wavenumber range shows a significant difference between linear and non-linear modelling, with the latter favouring GDM parameters inconsistent with $\Lambda$CDM, underlining the importance of careful non-linear modelling when using this data. We also use this halo model to establish the robustness of previously obtained constraints, particularly those that involve weak gravitational lensing of the cosmic microwave background. Additionally, we show how the inclusion of neutrino mass as a free parameter affects previous constraints on the GDM parameters., Comment: 27 pages, 10 figures. Updated to match published version: a few minor clarifications on the choices in the halo model construction and relation to WDM and an additional plot showing the mass function

Recently, an extension to the Parameterized Post-Newtonian (PPN) formalism has been proposed. This formalism, the Parameterized Post-Newtonian-Vainshteinian (PPNV) formalism, is well suited to theories which exhibit Vainshtein screening of scalar fields. In this paper we apply the PPNV formalism to the Quartic and Quintic Galileon theories for the first time. As simple generalizations of standard scalar-tensor field theories they are important guides for the generalization of parameterized approaches to the effects of gravity beyond General Relativity. In the Quartic case, we find new PPNV potentials for both screened and un-screened regions of spacetime, showing that in principle these theories can be tested. In the Quintic case we show that Vainshtein screening does not occur to Newtonian order, meaning that the theory behaves as Brans-Dicke to this order, and we discuss possible higher order effects., Comment: 18 pages, 2 figure

In this proceeding, we present studies of instrumental systematic effects for the Simons Obsevatory (SO) that are associated with the detector system and its interaction with the full SO experimental systems. SO will measure the Cosmic Microwave Background (CMB) temperature and polarization anisotropies over a wide range of angular scales in six bands with bandcenters spanning from 27 GHz to 270 GHz. We explore effects including intensity-to-polarization leakage due to coupling optics, bolometer nonlinearity, uncalibrated gain variations of bolometers, and readout crosstalk. We model the level of signal contamination, discuss proposed mitigation schemes, and present instrument requirements to inform the design of SO and future CMB projects., Comment: Proceeding from SPIE Astronomical Telescopes+Instrumentation 2018 (27 pages, 13 figures) v2: Added HEALPix reference

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources., Comment: This paper presents an overview of the Simons Observatory science goals, details about the instrument will be presented in a companion paper. The author contribution to this paper is available at https://simonsobservatory.org/publications.php (Abstract abridged) -- matching version published in JCAP

New telescopes are being built to measure the Cosmic Microwave Background (CMB) with unprecedented sensitivity, including Simons Observatory (SO), CCAT-prime, the BICEP Array, SPT-3G, and CMB Stage-4. We present observing strategies for telescopes located in Chile that are informed by the tools used to develop recent Atacama Cosmology Telescope (ACT) and Polarbear surveys. As with ACT and Polarbear, these strategies are composed of scans that sweep in azimuth at constant elevation. We explore observing strategies for both small (0.42 m) aperture telescopes (SAT) and a large (6 m) aperture telescope (LAT). We study strategies focused on small sky areas to search for inflationary gravitational waves as well as strategies spanning roughly half the low-foreground sky to constrain the effective number of relativistic species and measure the sum of neutrino masses via the gravitational lensing signal due to large scale structure. We present these strategies specifically considering the telescope hardware and science goals of the SO, located at 23 degrees South latitude, 67.8 degrees West longitude. Observations close to the Sun and the Moon can introduce additional systematics by applying additional power to the instrument through telescope sidelobes. Significant side lobe contamination in the data can occur even at tens of degrees or more from bright sources. Therefore, we present several strategies that implement Sun and Moon avoidance constraints into the telescope scheduling. Strategies for resolving conflicts between simultaneously visible fields are discussed. We focus on maximizing telescope time spent on science observations. It will also be necessary to schedule calibration measurements, however that is beyond the scope of this work. The outputs of this study are algorithms that can generate specific schedule commands for the Simons Observatory instruments., Comment: 15 pages, 9 figures, 3 tables, SPIE Astronomical Telescopes + Instrumentation 2018

Cold Dark Matter (CDM) is a crucial constituent of the current concordance cosmological model. Having a vanishing equation of state (EoS), its energy density scales with the inverse cosmic volume and is thus uniquely described by a single number, its present abundance. We test the inverse cosmic volume law for Dark Matter (DM) by allowing its EoS to vary independently in eight redshift bins in the range $z=10^5$ and $z=0$. We use the latest measurements of the Cosmic Microwave Background radiation from the Planck satellite and supplement them with Baryon Acoustic Oscillation (BAO) data from the 6dF and SDSS-III BOSS surveys, and with the Hubble Space Telescope (HST) key project data. We find no evidence for nonzero EoS in any of the eight redshift bins. With Planck data alone, the DM abundance is most strongly constrained around matter-radiation equality $\omega^{\rm eq}_g = 0.1193^{+0.0036}_{-0.0035}$ (95% c.l.), whereas its present day value is more weakly constrained $\omega^{(0)}_g = 0.16^{+0.12}_{-0.10}$ (95% c.l.). Adding BAO or HST data does not significantly change the $\omega^{\rm eq}_g$ constraint, while $\omega^{(0)}_g$ tightens to $0.160^{+0.069}_{-0.065} $ (95% c.l.) and $0.124^{+0.081}_{-0.067}$ (95% c.l.) respectively. Our results constrain for the first time the level of "coldness" required of the DM across various cosmological epochs and show that the DM abundance is strictly positive at all times., Comment: 6 pages, 4 figures, changed color scheme for figures

Weak lensing has become an increasingly important tool in cosmology and the use of galaxy shapes to measure cosmic shear has become routine. The weak-lensing distortion tensor contains two other effects in addition to the two components of shear: the convergence and rotation. The rotation mode is not measurable using the standard cosmic shear estimators based on galaxy shapes, as there is no information on the original shapes of the images before they were lensed. Due to this, no estimator has been proposed for the rotation mode in cosmological weak-lensing surveys, and the rotation mode has never been constrained. Here, we derive an estimator for this quantity, which is based on the use of radio polarisation measurements of the intrinsic position angles of galaxies. The rotation mode can be sourced by physics beyond $\Lambda$CDM, and also offers the chance to perform consistency checks of $\Lambda$CDM and of weak-lensing surveys themselves. We present simulations of this estimator and show that, for the pedagogical example of cosmic string spectra, this estimator could detect a signal that is consistent with the constraints from Planck. We examine the connection between the rotation mode and the shear $B$-modes and thus how this estimator could help control systematics in future radio weak-lensing surveys.

We examine how the properties of dark matter, parameterised by an equation of state parameter $w$ and two perturbative Generalised Dark Matter (GDM) parameters $c^2_s$ (the sound speed) and $c^2_\text{vis}$ (the viscosity), are constrained by existing cosmological data, particularly the Planck 2015 data release. We find that the GDM parameters are consistent with zero, and are strongly constrained, showing no evidence for extending the dark matter model beyond the Cold Dark Matter (CDM) paradigm. The dark matter equation of state is constrained to be within $-0.000896

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field., Comment: Updated to match version published in JCAP

Newtonian simulations are routinely used to examine the matter dynamics on non-linear scales. However, even on these scales, Newtonian gravity is not a complete description of gravitational effects. A post-Friedmann approach shows that the leading order correction to Newtonian theory is a vector potential in the metric. This vector potential can be calculated from N-body simulations, requiring a method for extracting the velocity field. Here, we present the full details of our calculation of the post-Friedmann vector potential, using the Delauney Tesselation Field Estimator (DTFE) code. We include a detailed examination of the robustness of our numerical result, including the effects of box size and mass resolution on the extracted fields. We present the power spectrum of the vector potential and find that the power spectrum of the vector potential is $\sim 10^5$ times smaller than the power spectrum of the fully non-linear scalar gravitational potential at redshift zero. Comparing our numerical results to perturbative estimates, we find that the fully non-linear result can be more than an order of magnitude larger than the perturbative estimate on small scales. We extend the analysis of the vector potential to multiple redshifts, showing that this ratio persists over a range of scales and redshifts. We also comment on the implications of our results for the validity and interpretation of Newtonian simulations., Comment: Matches version published in MNRAS; substantially reduced number of figures compared to v1. Supplementary figures available at MNRAS website or in v1. Detailed follow up to http://arxiv.org/abs/1306.1562

In this paper we examine cosmological weak lensing on non-linear scales and show that there are Newtonian and relativistic contributions and that the latter can also be extracted from standard Newtonian simulations. We use the post-Friedmann formalism, a post-Newtonian type framework for cosmology, to derive the full weak-lensing deflection angle valid on non-linear scales for any metric theory of gravity. We show that the only contributing term that is quadratic in the first order deflection is the expected Born correction and lens-lens coupling term. We use this deflection angle to analyse the vector and tensor contributions to the E- and B- mode cosmic shear power spectra. In our approach, once the gravitational theory has been specified, the metric components are related to the matter content in a well-defined manner. Specifying General Relativity, we write down a complete set of equations for a GR$+\Lambda$CDM universe for computing all of the possible lensing terms from Newtonian N-body simulations. We illustrate this with the vector potential and show that, in a GR$+\Lambda$CDM universe, its contribution to the E-mode is negligible with respect to that of the conventional Newtonian scalar potential, even on non-linear scales. Thus, under the standard assumption that Newtonian N-body simulations give a good approximation of the matter dynamics, we show that the standard ray tracing approach gives a good description for a $\Lambda$CDM cosmology., Comment: Updated to match version published in JCAP

We present the first calculation of an intrinsically relativistic quantity in fully non-linear cosmolog- ical large-scale structure studies. Traditionally, non-linear structure formation in standard {\Lambda}CDM cosmology is studied using N-body simulations, based on Newtonian gravitational dynamics on an expanding background. When one derives the Newtonian regime in a way that is a consistent ap- proximation to the Einstein equations, a gravito-magnetic vector potential - giving rise to frame dragging - is present in the metric in addition to the usual Newtonian scalar potential. At leading order, this vector potential does not affect the matter dynamics, thus it can be computed from Newtonian N-body simulations. We explain how we compute the vector potential from simulations in {\Lambda}CDM and examine its magnitude relative to the scalar potential. We also discuss some possible observable effects., Comment: 5 pages, 3 figurs

Model-independent parametrisations for examining departures from General Relativity have been increasingly studied over the past few years. Various observables have been used to constrain the parameters and forecasts for future surveys have been carried out. In one such forecast, galaxy cluster counts were used to constrain the parameters. Here, we carry out a limited set of $N$-body simulations, with a modified Poisson equation, to examine the accuracy of existing mass functions for modified gravity cosmologies. As well as altering the gravitational calculation, we include the effect of a screening scale to ensure consistency of the theory with solar system tests. Our results suggest that if a screening scale exists its effect can be taken into account in the cluster count calculation through its effect on the linear matter power spectrum. If this is done, the accuracy of the standard mass function formalism in modified gravity theories with reasonably small departures from General Relativity, as tested in this work, is comparable to the standard case., Comment: 7 pages, 3 figures

Model-independent parametrisations of modified gravity have attracted a lot of attention over the past few years and numerous combinations of experiments and observables have been suggested to constrain the parameters used in these models. Galaxy clusters have been mentioned, but not looked at as extensively in the literature as some other probes. Here we look at adding galaxy clusters into the mix of observables and examine how they could improve the constraints on the modified gravity parameters. In particular, we forecast the constraints from combining Planck satellite Cosmic Microwave Background (CMB) measurements and Sunyaev-Zeldovich (SZ) cluster catalogue with a DES-like weak lensing survey. We find that cluster counts significantly improve the constraints over those derived using CMB and WL. We then look at surveys further into the future, to see how much better it may be feasible to make the constraints., Comment: 14 pages, 8 figures, Updated to match version published in JCAP

Vector perturbations sourced by topological defects can generate rotations in the lensing of background galaxies. This is a potential smoking gun for the existence of defects since rotation generates a curl-like component in the weak lensing signal which is not generated by standard density perturbations at linear order. This rotation signal is calculated as generated by cosmic strings. Future large scale weak lensing surveys should be able to detect this signal even for string tensions an order of magnitude lower than current constraints.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form (‘SO-Nominal’) consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating (“Stage 3”) experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4. Construction of SO-Nominal is fully funded, and operations and data analysis are funded for part of the planned five-year observations. We will seek federal funding to complete the observations and analysis of SO-Nominal, at the $25M level. The SO has a low risk and cost efficient upgrade path – the 6 m LAT can accommodate almost twice the baseline number of detectors and the SATs can be duplicated at low cost. We will seek funding at the $75M level for an expansion of the SO (‘SO-Enhanced’) that fills the remaining focal plane in the LAT, adds three SATs, and extends operations by five years, substantially improving our science return. By this time SO may be operating as part of the larger CMB-S4 project. This white paper summarizes and extends material presented in, which describes the science goals of SO-Nominal, and which describe the instrument design.