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).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

Significant observational effort has been devoted to determining volumetric type Ia supernova rates at high redshifts, leading to clues about the nature of Ia supernova progenitors and constraints on the iron production in the universe. A complementary approach is to investigate type Ia supernova rates in individual, more nearby, galaxies. The popular A+B model for the specific supernova rate, while reliable for a wide range of galaxy properties, diverges for large specific star formation rates. Applying it outside its range of validity could lead to the prediction of excessive type Ia supernova rates. Moreover, the A+B model it is not directly derived from a delay time distribution. We here introduce a new model which is explicitly motivated by a simple delay time distribution composed of a prompt and a delayed component. The model is in remarkably good agreement with current observational constraints. It yields a prompt fraction of $f_{p}= 0.11^{+0.10}_{-0.06}$ in agreement with results based on volumetric rates of type Ia supernovae at high redshift (Rodney et al. 2014). The model is tested against realistic star formation rates from the Illustris-1 simulation and is found to be self consistent in the asymptotic limits. An analytic function that encapsulates the features of the new model is shown to be in excellent agreement with the data. In terms of goodness of fit, the new model is strongly preferred over the A+B model. At $\log{({\rm sSFR})} \gtrsim -9$ there are no constraints from observations. Observations in this regime will further constrain the delay time distribution of type Ia supernovae at short delay times.

Following the success of type Ia supernovae in constraining cosmologies at lower redshift $(z\lesssim2)$, effort has been spent determining if a similarly useful standardisable candle can be found at higher redshift. {In this work we determine the largest possible magnitude discrepancy between a constant dark energy $\Lambda$CDM cosmology and a cosmology in which the equation of state $w(z)$ of dark energy is a function of redshift for high redshift standard candles $(z\gtrsim2)$}. We discuss a number of popular parametrisations of $w(z)$ with two free parameters, $w_z$CDM cosmologies, including the Chevallier-Polarski-Linder and generalisation thereof, $n$CPL, as well as the Jassal-Bagla-Padmanabhan parametrisation. For each of these parametrisations we calculate and find extrema of $\Delta \mu$, the difference between the distance modulus of a $w_z$CDM cosmology and a fiducial $\Lambda$CDM cosmology as a function of redshift, given 68\% likelihood constraints on the parameters $P=(\Omega_{m,0}, w_0, w_a)$. The parameters are constrained using cosmic microwave background, baryon acoustic oscillations, and type Ia supernovae data using CosmoMC. We find that none of the tested cosmologies can deviate more than 0.05 mag from the fiducial $\Lambda$CDM cosmology at high redshift, implying that high redshift standard candles will not aid in discerning between a $w_z$CDM cosmology and the fiducial $\Lambda$CDM cosmology. Conversely, this implies that if high redshift standard candles are found to be in disagreement with $\Lambda$CDM at high redshift, then this is a problem not only for $\Lambda$CDM but for the entire family of $w_z$CDM cosmologies.

In this paper we investigate how observational effects could possibly bias cosmological inferences from peculiar velocity measurements. Specifically, we look at how bulk flow measurements are compared with theoretical predictions. Usually bulk flow calculations try to approximate the flow that would occur in a sphere around the observer. Using the Horizon Run 2 simulation we show that the traditional methods for bulk flow estimation can overestimate the magnitude of the bulk flow for two reasons: when the survey geometry is not spherical (the data do not cover the whole sky), and when the observations undersample the velocity distributions. Our results may explain why several bulk flow measurements found bulk flow velocities that seem larger than those expected in standard {\Lambda}CDM cosmologies. We recommend a different approach when comparing bulk flows to cosmological models, in which the theoretical prediction for each bulk flow measurement is calculated specifically for the geometry and sampling rate of that survey. This means that bulk flow values will not be comparable between surveys, but instead they are comparable with cosmological models, which is the more important measure.