Measurements of gravitational lensing of the cosmic microwave background (CMB) hold the promise of yielding unique insights into cosmology at high redshift. Uncertainties due to baryonic effects associated with galaxy formation and evolution, including gas cooling, star formation, and feedback from active galactic nuclei (AGN) and supernovae, have typically been neglected when forecasting the sensitivity of future CMB surveys. In this paper, we determine the impact of these effects using four suites of hydrodynamical simulations which incorporate various prescriptions for baryonic processes, namely OWLS, BAHAMAS, Horizon, and IllustrisTNG. Our analysis shows characteristic power suppressions of several percent in CMB lensing due to baryonic effects, compared to dark-matter only simulations, at experimentally observable angular scales. We investigate the associated bias in the inferred neutrino mass for experiments like the upcoming Simons Observatory and CMB-S4. Depending on the experimental precision and the strength of the baryonic feedback within the simulations, biases in the neutrino mass sum show significant dispersion, ranging from very small to an over-estimation by 1.1$\sigma$. We conclude that baryonic effects will likely be non-negligible for a detection of neutrino mass using CMB lensing.

The cosmic microwave background (CMB) is sensitive to the recent phase of accelerated cosmic expansion through the late-time integrated Sachs-Wolfe (ISW) effect, which manifests as secondary temperature fluctuations on large angular scales. However, the large cosmic variance from primary CMB fluctuations limits the usefulness of this effect in constraining dark energy or modified gravity. In this paper, we propose a novel method to separate the ISW signal from the primary signal using gravitational lensing, based on the fact that the ISW signal is, to a good approximation, not gravitationally lensed. We forecast how well we can isolate the ISW signal for different experimental configurations, and discuss various applications, including modified gravity, large-scale CMB anomalies, and measurements of local-type primordial non-Gaussianity. Although not within reach of current experiments, the proposed method is a unique way to remove the cosmic variance of the primary signal, allowing for better CMB-based constraints on late-time phenomena than previously thought possible.

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project.

Wave effects are often neglected in microlensing studies; however, for coherent point-like sources, such as pulsars and fast radio bursts (FRBs), wave effects will become important in their gravitational lensing. In this paper, we describe the wave optics formalism, its various limits, and the conditions for which these limits hold. Using the simple point lens as an example, we will show that the frequency dependence of wave effects breaks degeneracies that are present in the usual geometric optics limit, and constructive interference results in larger magnifications further from the lens. This latter fact leads to a generic increase in cross section for microlensing events in the wave-optics regime compared to the geometric optics regime. For realistic percent-level spectral sensitivities, this leads to a relative boost in lensing cross section of more than an order of magnitude. We apply the point-lens model to the lensing of FRBs and pulsars and find that these radio sources will be lensed in the full wave-optics regime by isolated masses in the range of $0.1-100\,M_\oplus$, which includes free-floating planets (FFPs), whose Einstein radius is smaller than the Fresnel scale. More generally, the interference pattern allows an instantaneous determination of lens masses, unlike traditional microlensing techniques which only yield a mass inference from the event timescale.

The large-scale clustering of matter is impacted by baryonic physics, particularly AGN feedback. Modelling or mitigating this impact will be essential for making full use of upcoming measurements of cosmic shear and other large-scale structure probes. We study baryonic effects on the matter bispectrum, using measurements from a selection of state-of-the-art hydrodynamical simulations: IllustrisTNG, Illustris, EAGLE, and BAHAMAS. We identify a low-redshift enhancement of the bispectrum, peaking at $k\sim 3h\,{\rm Mpc}^{-1}$, that is present in several simulations, and discuss how it can be associated to the evolving nature of AGN feedback at late times. This enhancement does not appear in the matter power spectrum, and therefore represents a new source of degeneracy breaking between two- and three-point statistics. In addition, we provide physical interpretations for other aspects of these measurements, and make initial comparisons to predictions from perturbation theory, empirical fitting formulas, and the response function formalism. We publicly release our measurements (including estimates of their uncertainty due to sample variance) and bispectrum measurement code as resources for the community.

Line-intensity mapping of the 21cm line is a powerful probe of large scale structure at z<6, tracing large-scale structure via neutral hydrogen content that is found within galaxies. In principle, it enables cost-efficient surveys of the matter distribution up to z~6, unlocking orders of magnitude more modes for observational cosmology. Canada has been a traditional leader in this field, having led the first detections of the cosmological 21cm signal via cross-correlations with optical galaxy surveys and having constructed the Canadian Hydrogen Intensity Mapping Experiment (CHIME). The field is now entering a new era where data is abundant, allowing studies in how to overcome systematics to be tackled in an empirical, head-on fashion. In the next few years, this will produce the first detection of the 21cm auto power spectrum, which will pave the way towards a large suite of scientific possibilities. These potentially include precision measurements on the dark energy equation of state and other LCDM parameters, constraints on how HI mass traces dark matter, a detection of neutrino effects on large-scale structure, and the use of 21cm lensing to further constrain cosmology. To turn these promising directions into reality, we recommend a sustained program of investment in 21cm cosmology, starting with funding for the Canadian Hydrogen Observatory and Radio transient Detector (CHORD), followed by small-scale development efforts targeting next-generation hardware and sustained support for theory and technical staff support. Additionally, Canada should invest in complementary line-intensity mapping efforts (such as with CO or [CII] lines) and maintain participation in next-generation international efforts such as the Packed Ultra-wideband Mapping Array (PUMA) and the Square Kilometre Array (SKA).

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

We investigate the detection prospects for gravitational lensing of three-dimensional maps from upcoming line intensity surveys, focusing in particular on the impact of gravitational nonlinearities on standard quadratic lensing estimators. Using perturbation theory, we show that these nonlinearities can provide a significant contaminant to lensing reconstruction, even for observations at reionization-era redshifts. However, we show how this contamination can be mitigated with the use of a "bias-hardened" estimator. Along the way, we present an estimator for reconstructing long-wavelength density modes, in the spirit of the "tidal reconstruction" technique that has been proposed elsewhere, and discuss the dominant biases on this estimator. After applying bias-hardening, we find that a detection of the lensing potential power spectrum will still be challenging for the first phase of SKA-Low, CHIME, and HIRAX, with gravitational nonlinearities decreasing the signal to noise by a factor of a few compared to forecasts that ignore these effects. On the other hand, cross-correlations between lensing and galaxy clustering or cosmic shear from a large photometric survey look promising, provided that systematics can be sufficiently controlled. We reach similar conclusions for a single-dish survey inspired by CII measurements planned for CCAT-prime, suggesting that lensing is an interesting science target not just for 21cm surveys, but also for intensity maps of other lines.

We evaluate the potential for current and future cosmic shear measurements from large galaxy surveys to constrain the impact of baryonic physics on the matter power spectrum. We do so using a model-independent parameterization that describes deviations of the matter power spectrum from the dark-matter-only case as a set of principal components that are localized in wavenumber and redshift. We perform forecasts for a variety of current and future datasets, and find that at least ~90% of the constraining power of these datasets is contained in no more than nine principal components. The constraining power of different surveys can be quantified using a figure of merit defined relative to currently available surveys. With this metric, we find that the final Dark Energy Survey dataset (DES Y5) and the Hyper Suprime Cam Survey will be roughly an order of magnitude more powerful than existing data in constraining baryonic effects. Upcoming Stage IV surveys (LSST, Euclid, and WFIRST) will improve upon this by a further factor of a few. We show that this conclusion is robust to marginalization over several key systematics. The ultimate power of cosmic shear to constrain galaxy formation is dependent on understanding systematics in the shear measurements at small (sub-arcminute) scales. If these systematics can be sufficiently controlled, cosmic shear measurements from DES Y5 and other future surveys have the potential to provide a very clean probe of galaxy formation and to strongly constrain a wide range of predictions from modern hydrodynamical simulations.

We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

The most effective use of data from current and upcoming large scale structure~(LSS) and CMB observations requires the ability to predict the clustering of LSS with very high precision. The Effective Field Theory of Large Scale Structure (EFTofLSS) provides an instrument for performing analytical computations of LSS observables with the required precision in the mildly nonlinear regime. In this paper, we develop efficient implementations of these computations that allow for an exploration of their dependence on cosmological parameters. They are based on two ideas. First, once an observable has been computed with high precision for a reference cosmology, for a new cosmology the same can be easily obtained with comparable precision just by adding the difference in that observable, evaluated with much less precision. Second, most cosmologies of interest are sufficiently close to the Planck best-fit cosmology that observables can be obtained from a Taylor expansion around the reference cosmology. These ideas are implemented for the matter power spectrum at two loops and are released as public codes. When applied to cosmologies that are within 3$\sigma$ of the Planck best-fit model, the first method evaluates the power spectrum in a few minutes on a laptop, with results that have 1\% or better precision, while with the Taylor expansion the same quantity is instantly generated with similar precision. The ideas and codes we present may easily be extended for other applications or higher-precision results.

We study the prediction of the dark matter power spectrum at two-loop order in the Effective Field Theory of Large Scale Structures (EFTofLSS) using high precision numerical simulations. In our universe, short distance non-linear perturbations, not under perturbative control, affect long distance fluctuations through an effective stress tensor that needs to be parametrized in terms of counterterms that are functions of the long distance fluctuating fields. We find that at two-loop order it is necessary to include three counterterms: a linear term in the over density, $\delta$, a quadratic term, $\delta^2$, and a higher derivative term, $\partial^2\delta$. After the inclusion of these three terms, the EFTofLSS at two-loop order matches simulation data up to $k\simeq 0.34 \,h\, {\rm Mpc}^{-1}$ at redshift $z=0$, up to $k\simeq 0.55\,h\, {\rm Mpc}^{-1}$ at $z=1$, and up to $k\simeq 1.1\,h\, {\rm Mpc}^{-1}$ at $z=2$. At these wavenumbers, the cosmic variance of the simulation is at least as small as $10^{-3}$, providing a high precision comparison between theory and data. The actual reach of the theory is affected by theoretical uncertainties associated to not having included higher order terms in perturbation theory, for which we provide an estimate, and by potentially overfitting the data, which we also try to address. Since in the EFTofLSS the coupling constants associated with the counterterms are unknown functions of time, we show how a simple parametrization gives a sensible description of their time-dependence. Overall, the $k$-reach of the EFTofLSS is much larger than previous analytical techniques, showing that the amount of cosmological information amenable to high-precision analytical control might be much larger than previously believed.

Given the importance of future large scale structure surveys for delivering new cosmological information, it is crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbative scheme to compute the clustering of dark matter in the weakly nonlinear regime in an expansion in $k/k_{\rm NL}$, where $k$ is the wavenumber of interest and $k_{\rm NL}$ is the wavenumber associated to the nonlinear scale. It has been recently shown that the EFTofLSS matches to $1\%$ level the dark matter power spectrum at redshift zero up to $k\simeq 0.3 h\,$Mpc$^{-1}$ and $k\simeq 0.6 h\,$Mpc$^{-1}$ at one and two loops respectively, using only one counterterm that is fit to data. Similar results have been obtained for the momentum power spectrum at one loop. This is a remarkable improvement with respect to former analytical techniques. Here we study the prediction for the equal-time dark matter bispectrum at one loop. We find that at this order it is sufficient to consider the same counterterm that was measured in the power spectrum. Without any remaining free parameter, and in a cosmology for which $k_{\rm NL}$ is smaller than in the previously considered cases ($\sigma_8=0.9$), we find that the prediction from the EFTofLSS agrees very well with $N$-body simulations up to $k\simeq 0.25 h\,$Mpc$^{-1}$, given the accuracy of the measurements, which is of order a few percent at the highest $k$'s of interest. While the fit is very good on average up to $k\simeq 0.25 h\,$Mpc$^{-1}$, the fit performs slightly worse on equilateral configurations, in agreement with expectations that for a given maximum $k$, equilateral triangles are the most nonlinear.

We study the prediction of the Effective Field Theory of Large Scale Structures (EFTofLSS) for the matter power spectrum at different redshifts. In previous work, we found that the two-loop prediction can match the nonlinear power spectrum measured from $N$-body simulations at redshift zero within approximately 2% up to $k\sim 0.6\,h\, {\rm Mpc}^{-1}$ after fixing a single free parameter, the so-called "speed of sound". We determine the time evolution of this parameter by matching the EFTofLSS prediction to simulation output at different redshifts, and find that it is well-described by a fitting function that only includes one additional parameter. After the two free parameters are fixed, the prediction agrees with nonlinear data within approximately 2% up to at least $k\sim 1\,h\, {\rm Mpc}^{-1}$ at $z\geq 1$, and also within approximately 5% up to $k\sim 1.2\,h\, {\rm Mpc}^{-1}$ at $z=1$ and $k\sim 2.3\,h\, {\rm Mpc}^{-1}$ at $z=3$, a major improvement with respect to other perturbative techniques. We also develop an accurate way to estimate where the EFTofLSS predictions at different loop orders should fail, based on the sizes of the next-order terms that are neglected, and find agreement with the actual comparisons to data. Finally, we use our matter power spectrum results to perform analytical calculations of lensing potential power spectra corresponding to both CMB and galaxy lensing. This opens the door to future direct applications of the EFTofLSS to observations of gravitational clustering on cosmic scales.

Large scale structure surveys promise to be the next leading probe of cosmological information. It is therefore crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbation theory for the weakly non-linear regime of dark matter, where correlation functions are computed in an expansion of the wavenumber k of a mode over the wavenumber associated with the non-linear scale k_nl. Since most of the information is contained at high wavenumbers, it is necessary to compute higher order corrections to correlation functions. After the one-loop correction to the matter power spectrum, we estimate that the next leading one is the two-loop contribution, which we compute here. At this order in k/k_nl, there is only one counterterm in the EFTofLSS that must be included, though this term contributes both at tree-level and in several one-loop diagrams. We also discuss correlation functions involving the velocity and momentum fields. We find that the EFTofLSS prediction at two loops matches to percent accuracy the non-linear matter power spectrum at redshift zero up to k~0.6 h/Mpc, requiring just one unknown coefficient that needs to be fit to observations. Given that Standard Perturbation Theory stops converging at redshift zero at k~0.1 h/Mpc, our results demonstrate the possibility of accessing a factor of order 200 more dark matter quasi-linear modes than naively expected. If the remaining observational challenges to accessing these modes can be addressed with similar success, our results show that there is tremendous potential for large scale structure surveys to explore the primordial universe.

Large scale structure surveys are likely the next leading probe of cosmological information. It is therefore crucial to reliably predict their observables. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a manifestly convergent perturbation theory for the weakly non-linear regime, where dark matter correlation functions are computed in an expansion of the wavenumber k over the wavenumber associated to the non-linear scale knl. To push the predictions to higher wavenumbers, it is necessary to compute the 2-loop matter power spectrum. For equal-time correlators, exactly as with standard perturturbation theory, there are IR divergences present in each diagram that cancel completely in the final result. We develop a method by which all 2-loop diagrams are computed as one integral, with an integrand that is manifestly free of any IR divergences. This allows us to compute the 2-loop power spectra in a reliable way that is much less numerically challenging than standard techniques. We apply our method to scaling universes where the linear power spectrum is a single power law of k, and where IR divergences can particularly easily interfere with accurate evaluation of loop corrections if not handled carefully. We show that our results are independent of IR cutoff and, after renormalization, of the UV cutoff, and comment how the method presented here naturally generalizes to higher loops.

It has recently been claimed that measurements of the baryonic Tully-Fisher relation (BTFR), a power-law relationship between the observed baryonic masses and outer rotation velocities of galaxies, support the predictions of modified Newtonian dynamics for the slope and scatter in the relation, while challenging the cold dark matter (CDM) paradigm. We investigate these claims, and find that: 1) the scatter in the data used to determine the BTFR is in conflict with observational uncertainties on the data; 2) these data do not make strong distinctions regarding the best-fit BTFR parameters; 3) the literature contains a wide variety of measurements of the BTFR, many of which are discrepant with the recent results; and 4) the claimed CDM "prediction" for the BTFR is a gross oversimplification of the complex galaxy-scale physics involved. We conclude that the BTFR is currently untrustworthy as a test of CDM.

We perform a Fisher matrix calculation of the predicted uncertainties on estimates of the cosmic string tension Gmu from upcoming observational data (namely, cosmic microwave background power spectra from the Planck satellite and an idealized future polarization experiment). We employ simulations that are more general than others commonly used in the literature, leaving the mean velocity of strings, correlation length of the string network, and "wiggliness" (which parametrizes smaller-scale structure along the strings) as free parameters that can be observationally measured. In a new code, StringFast, we implement a method for efficient computation of the C_l spectra induced by a network of strings, which is fast enough to be used in Markov Chain Monte Carlo analyses of future data. Performing a calculation with the string parameters left free results in projected constraints on Gmu that are larger than those obtained by fixing their values a priori, typically by a factor of ~2-7. We also find that if Gmu is equal to the current observational maximum, Planck will be able to make a confident detection of strings. However, if Gmu is two orders of magnitude smaller, even a perfect, lensing-free measurement of polarization power spectra will not be able to detect a nonzero string tension at better than 2 sigma confidence.

There has been considerable interest in recent years in cosmological models in which we inhabit a very large, underdense void as an alternative to dark energy. A longstanding objection to this proposal is that observations limit our position to be very close to the void centre. By selecting from a family of void profiles that fit supernova luminosity data, we carefully determine how far from the centre we could be. To do so, we use the observed dipole component of the cosmic microwave background, as well as an additional stochastic peculiar velocity arising from primordial perturbations. We find that we are constrained to live within 80 Mpc of the centre of a void--a somewhat weaker constraint than found in previous studies, but nevertheless a strong violation of the Copernican principle. By considering how such a Gpc-scale void would appear on the microwave sky, we also show that there can be a maximum of one of these voids within our Hubble radius. Hence, the constraint on our position corresponds to a fraction of the Hubble volume of order 10^{-8}. Finally, we use the fact that void models only look temporarily similar to a cosmological-constant-dominated universe to argue that these models are not free of temporal fine-tuning.