We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

The observed tension ($\sim 9\%$ difference) between the local distance ladder measurement of the Hubble constant, $H_0$, and its value inferred from the cosmic microwave background (CMB) could hint at new, exotic, cosmological physics. We test the impact of the assumption about the expansion history of the universe ($0.01<z<2.3$) on the local distance ladder estimate of $H_0$. In the fiducial analysis, the Hubble flow Type Ia supernova (SN~Ia) sample is truncated to $z < 0.15$ and the deceleration parameter ($q_0$) fixed to -0.55. We create realistic simulations of the calibrator and Pantheon samples and account for a full systematics covariance between these two sets. We fit several physically motivated dark energy models and derive combined constraints from calibrator and Pantheon SNe~Ia and simultaneously infer $H_0$ and dark energy properties. We find that the assumption on the dark energy model does not significantly change the local distance ladder value of $H_0$, with a maximum difference ($\Delta H_0$) between the inferred value for different models of 0.47 km$^{-1}$ s$^{-1}$ Mpc $^{-1}$, i.e. a 0.6$\%$ shift in $H_0$, significantly smaller than the observed tension. Additional freedom in the dark energy models does not increase the error in the inferred value of $H_0$. Including systematics covariance between the calibrators, low redshift SNe, and high redshift SNe can induce small shifts in the inferred value for $H_0$. The SN~Ia systematics in this study contribute $\lesssim 0.8 \%$ to the total uncertainty on $H_0$.

The Wide Field InfraRed Survey Telescope (WFIRST) was the highest ranked large space-based mission of the 2010 New Worlds, New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid-2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral-field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated, we find the most successful to be WFC-focused. However, further work in constraining systematics is required to fully optimize the use of the IFC. Even without improvements to other cosmological probes, the WFIRST SN survey has the potential to increase the FoM by more than an order of magnitude from the current values. Although the survey strategies presented here have not been fully optimized, these initial investigations are an important step in the development of the final hardware design and implementation of the WFIRST mission.

We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, $\Lambda$CDM and $w$CDM, by using a combined analysis of galaxy clustering and weak gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: 1) nonzero curvature $\Omega_k$, 2) number of relativistic species $N_{\rm eff}$ different from the standard value of 3.046, 3) time-varying equation-of-state of dark energy described by the parameters $w_0$ and $w_a$ (alternatively quoted by the values at the pivot redshift, $w_p$, and $w_a$), and 4) modified gravity described by the parameters $\mu_0$ and $\Sigma_0$ that modify the metric potentials. We also consider external information from Planck CMB measurements; BAO measurements from SDSS, 6dF, and BOSS; RSD measurements from BOSS; and SNIa information from the Pantheon compilation. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find $\Omega_k=0.0020^{+0.0037}_{-0.0032}$ at the 68% confidence level, and $N_{\rm eff}<3.28\, (3.55)$ at 68% (95%) confidence. For the time-varying equation-of-state, we find the pivot value $(w_p, w_a)=(-0.91^{+0.19}_{-0.23}, -0.57^{+0.93}_{-1.11})$ at pivot redshift $z_p=0.27$ from DES alone, and $(w_p, w_a)=(-1.01^{+0.04}_{-0.04}, -0.28^{+0.37}_{-0.48})$ at $z_p=0.20$ from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be $\Sigma_0= 0.43^{+0.28}_{-0.29}$ from DES Y1 alone, and $(\Sigma_0, \mu_0)=(0.06^{+0.08}_{-0.07}, -0.11^{+0.42}_{-0.46})$ from DES Y1 combined with external data, consistent with predictions from GR.

We present the analysis underpinning the measurement of cosmological parameters from 207 spectroscopically classified type Ia supernovae (SNe Ia) from the first three years of the Dark Energy Survey Supernova Program (DES-SN), spanning a redshift range of 0.017<$z$<0.849. We combine the DES-SN sample with an external sample of 122 low-redshift ($z$<0.1) SNe Ia, resulting in a "DES-SN3YR" sample of 329 SNe Ia. Our cosmological analyses are blinded: after combining our DES-SN3YR distances with constraints from the Cosmic Microwave Background (CMB; Planck Collaboration 2016), our uncertainties in the measurement of the dark energy equation-of-state parameter, $w$, are .042 (stat) and .059 (stat+syst) at 68% confidence. We provide a detailed systematic uncertainty budget, which has nearly equal contributions from photometric calibration, astrophysical bias corrections, and instrumental bias corrections. We also include several new sources of systematic uncertainty. While our sample is <1/3 the size of the Pantheon sample, our constraints on $w$ are only larger by 1.4$\times$, showing the impact of the DES SN Ia light curve quality. We find that the traditional stretch and color standardization parameters of the DES SNe Ia are in agreement with earlier SN Ia samples such as Pan-STARRS1 and the Supernova Legacy Survey. However, we find smaller intrinsic scatter about the Hubble diagram (0.077 mag). Interestingly, we find no evidence for a Hubble residual step ( 0.007 $\pm$ 0.018 mag) as a function of host galaxy mass for the DES subset, in 2.4$\sigma$ tension with previous measurements. We also present novel validation methods of our sample using simulated SNe Ia inserted in DECam images and using large catalog-level simulations to test for biases in our analysis pipelines.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

We investigate the bounds on the sum of neutrino masses in a cosmic-acceleration scenario where the equation of state $w(z)$ of dark energy (DE) is constructed in a model-independent way, using a basis of principal components (PCs) that are allowed to cross the phantom barrier $w(z)=-1$. We find that the additional freedom provided to $w(z)$ means the DE can undo changes in the background expansion induced by massive neutrinos at low redshifts. This has two significant consequences: (1) it leads to a substantial increase in the upper bound for the sum of the neutrino masses ($M_{\nu} < 0.33 - 0.55$ eV at the 95\% C.L. depending on the data sets and number of PCs included) compared to studies that choose a specific parametrization for $w(z)$; and (2) it causes $\sim1\sigma$ deviations from $\Lambda$CDM in the luminosity distance and the Hubble expansion rate at higher redshifts ($z \gtrsim 2$), where the contribution of DE is subdominant and there is little constraining data. The second point consequently means that there are also observable deviations in the shear power spectrum and in the matter power spectrum at low redshift, since the clustering of matter throughout cosmic time depends on the expansion rate. This provides a compelling case to pursue high-$z$ BAO and SN measurements as a way of disentangling the effects of neutrinos and dark energy. Finally, we find that the additional freedom given to the dark energy component has the effect of lowering $S_8$ with respect to $\Lambda$CDM.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

The largest structures in the cosmic web probe the dynamical nature of dark energy through their integrated Sachs-Wolfe imprints. In the strength of the signal, typical cosmic voids have shown good consistency with expectation $A_{\rm ISW}=\Delta T^{\rm data} / \Delta T^{\rm theory}=1$, given the substantial cosmic variance. Discordantly, large-scale hills in the gravitational potential, or supervoids, have shown excess signals. In this study, we mapped out 87 new supervoids in the total 5000 deg$^2$ footprint of the Dark Energy Survey at $0.2<z<0.9$ to probe these anomalous claims. We found an excess imprinted profile with $ A_{\rm ISW}\approx4.1\pm2.0$ amplitude. The combination with independent BOSS data reveals an ISW imprint of supervoids at the $3.3\sigma$ significance level with an enhanced $A_{\rm ISW}\approx5.2\pm1.6$ amplitude. The tension with $\Lambda$CDM predictions is equivalent to $2.6\sigma$ and remains unexplained.

Features during inflation and reionization leave corresponding features in the temperature and polarization power spectra that could potentially explain anomalies in the Planck 2015 data but require a joint analysis to disentangle. We study the interplay between these two effects using a model-independent parametrization of the inflationary power spectrum and the ionization history. Preference for a sharp suppression of large scale power is driven by a feature in the temperature power spectrum at multipoles $\ell \sim 20$, whereas preference for a component of high redshift ionization is driven by a sharp excess of polarization power at $\ell \sim 10$ when compared with the lowest multipoles. Marginalizing inflationary freedom does not weaken the preference for $z \gtrsim 10$ ionization, whereas marginalizing reionization freedom slightly enhances the preference for an inflationary feature but can also mask its direct signature in polarization. The inflation and reionization interpretation of these features makes predictions for the polarization spectrum which can be tested in future precision measurements especially at $10\lesssim \ell \lesssim 40$.

We develop an integral form for the bispectrum in general single-field inflation whose domain of validity includes models of inflation where the background evolution is not constrained to be slowly varying everywhere. Our integral form preserves the squeezed-limit consistency relation, allows for fast evaluation of the bispectrum for all triangle configurations expediting the efficient comparison of slow-roll violating models with data, and provides complete and compact slow-roll expressions correct to first order in slow-roll parameters. Motivated by the recent Planck results, we consider as an example a sharp step in the warped-brane tension of DBI inflation and provide analytic solutions for the peak of the resulting bispectrum. For the step in the warp that reproduces the oscillations in the power spectrum favored by the Planck data, the corresponding equilateral bispectrum is both extremely large and highly scale dependent. The bispectrum serves as a means of distinguishing such a model from alternative scenarios that generate otherwise indistinguishable power spectra, such as a step in the potential in canonical single-field inflation.

Hubble Space Telescope observations from December 11 2015 detected the expected fifth counter image of SN Refsdal at $z = 1.49$. In this letter, we compare the time delay predictions from numerous models with the measured value derived by \citet{Kelly2016a} from very early data in the light curve of the SN Refsdal, and find a best value for $H_0 = 64^{+9}_{-11}~\mathrm{km~s^{-1}~Mpc^{-1}}$ (68\% CL), in excellent agreement with predictions from CMB and recent weak lensing data + BAO + BBN (from the DES Collaboration). This is the first constraint on $H_0$ derived from time delays between multiple lensed SN images, and the first with a galaxy cluster lens, so subject to systematic effects different from other time-delay $H_0$ estimates. Additional time delay measurements from new multiply-imaged SNe will allow derivation of competitive constraints on $H_0$.

Through likelihood analyses of both current and future data that constrain both the expansion history of the universe and the clustering of matter fluctuations, we provide falsifiable predictions for three broad classes of models that explain the accelerated expansions of the universe: $\Lambda$CDM, the quintessence scenario and a more general class of smooth dark energy models that can cross the phantom barrier $w(z)=-1$. Our predictions are model independent in the sense that we do not rely on a specific parametrization, but we instead use a principal component (PC) basis function constructed a priori from a noise model of supernovae and Cosmic Microwave Background observations. For the supernovae measurements, we consider two type of surveys: the current JLA and the upcoming WFIRST surveys. We show that WFIRST will be able to improve growth predictions in curved models significantly. The remaining degeneracy between spatial curvature and $w(z)$ could be overcome with improved measurements of $\sigma_8 \Omega_m^{1/2}$, a combination that controls the amplitude of the growth of structure. We also point out that a PC-based Figure of Merit reveals that the usual two-parameter description of $w(z)$ does not exhaust the information that can be extracted from current data (JLA) or future data (WFIRST).

We reexamine interactions between the dark sectors of cosmology, with a focus on robust constraints that can be obtained using only mildly nonlinear scales. While it is well known that couplings between dark matter and dark energy can be constrained to the percent level when including the full range of scales probed by future optical surveys, calibrating matter power spectrum emulators to all possible choices of potentials and couplings requires many computationally expensive n-body simulations. Here we show that lensing and clustering of galaxies in combination with the Cosmic Microwave Background (CMB) is capable of probing the dark sector coupling to the few percent level for a given class of models, using only linear and quasi-linear Fourier modes. These scales can, in principle, be described by semi-analytical techniques such as the effective field theory of large-scale structure.

We present the methodology for and detail the implementation of the Dark Energy Survey (DES) 3x2pt DES Year 1 (Y1) analysis, which combines configuration-space two-point statistics from three different cosmological probes: cosmic shear, galaxy-galaxy lensing, and galaxy clustering, using data from the first year of DES observations. We have developed two independent modeling pipelines and describe the code validation process. We derive expressions for analytical real-space multi-probe covariances, and describe their validation with numerical simulations. We stress-test the inference pipelines in simulated likelihood analyses that vary 6-7 cosmology parameters plus 20 nuisance parameters and precisely resemble the analysis to be presented in the DES 3x2pt analysis paper, using a variety of simulated input data vectors with varying assumptions. We find that any disagreement between pipelines leads to changes in assigned likelihood $\Delta \chi^2 \le 0.045$ with respect to the statistical error of the DES Y1 data vector. We also find that angular binning and survey mask do not impact our analytic covariance at a significant level. We determine lower bounds on scales used for analysis of galaxy clustering (8 Mpc$~h^{-1}$) and galaxy-galaxy lensing (12 Mpc$~h^{-1}$) such that the impact of modeling uncertainties in the non-linear regime is well below statistical errors, and show that our analysis choices are robust against a variety of systematics. These tests demonstrate that we have a robust analysis pipeline that yields unbiased cosmological parameter inferences for the flagship 3x2pt DES Y1 analysis. We emphasize that the level of independent code development and subsequent code comparison as demonstrated in this paper is necessary to produce credible constraints from increasingly complex multi-probe analyses of current data.

We explore the relationship between features in the Planck 2015 temperature and polarization data, shifts in the cosmological parameters, and features from inflation. Residuals in the temperature data at low multipole $\ell$, which are responsible for the high $H_0\approx 70$ km s$^{-1}$Mpc$^{-1}$ and low $\sigma_8\Omega_m^{1/2}$ values from $\ell<1000$ in power-law $\Lambda$CDM models, are better fit to inflationary features with a $1.9\sigma$ preference for running of the running of the tilt or a stronger $99\%$ CL local significance preference for a sharp drop in power around $k=0.004$ Mpc$^{-1}$ in generalized slow roll and a lower $H_0\approx 67$ km s$^{-1}$Mpc$^{-1}$. The same in-phase acoustic residuals at $\ell>1000$ that drive the global $H_0$ constraints and appear as a lensing anomaly also favor running parameters which allow even lower $H_0$, but not once lensing reconstruction is considered. Polarization spectra are intrinsically highly sensitive to these parameter shifts, and even more so in the Planck 2015 TE data due to an outlier at $\ell \approx 165$, which disfavors the best fit $H_0$ $\Lambda$CDM solution by more than $2\sigma$, and high $H_0$ value at almost $3\sigma$. Current polarization data also slightly enhance the significance of a sharp suppression of large-scale power but leave room for large improvements in the future with cosmic variance limited $E$-mode measurements.

We conduct an analysis of the Planck 2015 data that is complete in reionization observables from the large angle polarization $E$-mode spectrum in the redshift range $6 < z < 30$. Based on 5 principal components, all of which are constrained by the data, this single analysis can be used to infer constraints on any model for reionization in the same range; we develop an effective likelihood approach for applying these constraints to models. By allowing for an arbitrary ionization history, this technique tests the robustness of inferences on the total optical depth from the usual step-like transition assumption, which is important for the interpretation of many other cosmological parameters such as the dark energy and neutrino mass. The Planck 2015 data not only allow a high redshift $z>15$ component to the optical depth but prefer it at the $2\sigma$ level. This preference is associated with excess power in the multipole range $10 \lesssim \ell \lesssim 20$ and may indicate high redshift ionization sources or unaccounted for systematics and foregrounds in the 2015 data.

Standard analyses of the reionization history of the universe from Planck cosmic microwave background (CMB) polarization measurements consider only the overall optical depth to electron scattering ($\tau$), and further assume a step-like reionization history. However, the polarization data contain information beyond the overall optical depth, and the assumption of a step-like function may miss high redshift contributions to the optical depth and lead to biased $\tau$ constraints. Accounting for its full reionization information content, we reconsider the interpretation of Planck 2015 Low Frequency Instrument (LFI) polarization data using simple, yet physically-motivated reionization models. We show that these measurements still, in fact, allow a non-negligible contribution from metal-free (Pop-III) stars forming in mini-halos of mass $M \sim 10^5-10^6 M_\odot$ at $z \gtrsim 15$, provided this mode of star formation is fairly inefficient. Our best fit model includes an early, self-regulated phase of Pop-III star formation in which the reionization history has a gradual, plateau feature. In this model, $\sim$20\% of the volume of the universe is ionized by $z \sim 20$, yet it nevertheless provides a good match to the Planck LFI measurements. Although preferred when the full information content of the data is incorporated, this model would spuriously be disfavored in the standard analysis. This preference is driven mostly by excess power from E-mode polarization at multipoles of $10 \lesssim \ell \lesssim 20$, which may reflect remaining systematic errors in the data, a statistical fluctuation, or signatures of the first stars. Measurements from the Planck High Frequency Instrument (HFI) should be able to confirm or refute this hint and future cosmic-variance limited E-mode polarization surveys can provide substantially more information on these signatures

We develop methods to calculate the curvature power spectrum in models where features in the inflaton potential nonlinearly excite modes and generate high frequency features in the spectrum. The first nontrivial effect of excitations generating further excitations arises at third order in deviations from slow roll. If these further excitations are contemporaneous, the series can be resummed, showing the exponential sensitivity of the curvature spectrum to potential features. More generally, this exponential approximation provides a power spectrum template which nonlinearly obeys relations between excitation coefficients and whose parameters may be appropriately adjusted. For a large sharp step in the potential, it greatly improves the analytic power spectrum template and its dependence on potential parameters. For axionic oscillations in the potential, it corrects the mapping between the potential and the amplitude, phase and zero point of the curvature oscillations, which might otherwise cause erroneous inferences in for example the tensor-scalar ratio, formally even when that amplitude is $10^3$ times larger than the slow roll power spectrum. It also estimates when terms that produce double frequency oscillations that are usually omitted when analyzing data should be included. These techniques should allow future studies of high frequency features in the CMB and large scale structure to extend to higher amplitude and/or higher precision.

We conduct a model-independent analysis of temporal features during inflation in the large-scale CMB temperature power spectrum allowing for the possibility of non-negligible tensor contributions. Of 20 principal components of the inflationary history, the suppression of power at low multipoles beginning with a glitch at multipoles $\ell \sim 20-40$ implies deviations in 2-3 of them with 2-3$\sigma$ deviations in each, with larger values reflecting cases where tensors are allowed.If tensors are absent, the corresponding $E$-mode polarization features follow a similar pattern but are predicted to be up to twice as large. They offer the opportunity to soon double the significance of inflationary features or eliminate them as an explanation of temperature features. The tensor degeneracy with features in the temperature power spectrum is broken not only by $B$ but also by $E$-polarization. A precision measurement of $E$-mode polarization at multipoles from $\ell\sim 20-60$ can potentially provide an independent constraint on tensors that is less subject to dust foreground uncertainties.

The recent BICEP2 B-mode polarization determination of an inflationary tensor-scalar ratio $r=0.2^{+0.07}_{-0.05}$ is in tension with simple scale-free models of inflation due to a lack of a corresponding low multipole excess in the temperature power spectrum which places a limit of $r_{0.002}<0.11$ (95% CL) on such models. Single-field inflationary models that reconcile these two observations, even those where the tilt runs substantially, introduce a scale into the scalar power spectrum. To cancel the tensor excess, and simultaneously explain the excess already present in $\Lambda$CDM, ideally the model should introduce this scale as a relatively sharp transition in the tensor-scalar ratio around the horizon at recombination. We consider models which generate such a step in this quantity and find that they can improve the joint fit to the temperature and polarization data by up to $2\Delta \ln{\cal L} \approx -14$ without changing cosmological parameters. Precision E-mode polarization measurements should be able to test this explanation.

We extend and improve the modeling and analysis of large-amplitude, sharp inflationary steps for second order corrections required by the precision of the Planck CMB power spectrum and for arbitrary Dirac-Born-Infeld sound speed. With two parameters, the amplitude and frequency of the resulting oscillations, step models improve the fit by $\Delta \chi^2 = -11.4$. Evidence for oscillations damping before the Planck beam scale is weak: damping only improves the fit to $\Delta \chi^2 = -14.0$ for one extra parameter, if step and cosmological parameters are jointly fit, in contrast to analyses which fix the latter. Likewise, further including the sound speed as a parameter only marginally improves the fit to $\Delta \chi^2 = -15.2$ but has interesting implications for the lowest multipole temperature and polarization anisotropy. Since chance features in the noise can mimic these oscillatory features, we discuss tests from polarization power spectra, lensing reconstruction and squeezed and equilateral bispectra that should soon verify or falsify their primordial origin.

In Dirac-Born-Infeld inflation, changes in the sound speed that transiently break the slow roll approximation lead to features in the power spectrum. We develop and test the generalized slow roll approximation for calculating such effects and show that it can be extended to treat order unity features. As in slow-roll, model independent constraints on the potential of canonical inflation can be directly reinterpreted in the DBI context through this approximation. In particular, a sharp horizon scale step in the warped brane tension can explain oscillatory features in the WMAP7 CMB power spectrum as well as features in the potential. Differences appear only as a small suppression of power on horizon scales and larger.

Several authors have argued that self-consistent $f(R)$ gravity models distinct from $\Lambda $CDM are almost ruled out. Confronting such claims, we present a particular two-parameter $f(R)$ model that: (a) is cosmologically viable and distinguishable from $\Lambda $CDM; (b) is compatible with the existence of relativistic stars; (c) is free of singularities of the Ricci scalar during the cosmological evolution and (d) allows the addition of high curvature corrections that could be relevant for inflation.