We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg2 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 ΛCDM and wCDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for ΛCDM) or 7 (for wCDM) cosmological parameters including the neutrino mass density and including the 457×457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions and from their combination obtain S8≡σ8(Ωm/0.3)0.5=0.773-0.020+0.026 and Ωm=0.267-0.017+0.030 for ΛCDM; for wCDM, we find S8=0.782-0.024+0.036, Ωm=0.284-0.030+0.033, and w=-0.82-0.20+0.21 at 68% C.L. The precision of these DES Y1 constraints 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 S8 and Ωm are lower than the central values from Planck for both ΛCDM and wCDM, the Bayes factor indicates that the DES Y1 and Planck data sets are consistent with each other in the context of ΛCDM. Combining DES Y1 with Planck, baryonic acoustic oscillation measurements from SDSS, 6dF, and BOSS and type Ia supernovae from the Joint Lightcurve Analysis data set, we derive very tight constraints on cosmological parameters: S8=0.802±0.012 and Ωm=0.298±0.007 in ΛCDM and w=-1.00-0.04+0.05 in wCDM. Upcoming Dark Energy Survey analyses will provide more stringent tests of the ΛCDM model and extensions such as a time-varying equation of state of dark energy or modified gravity.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single-epoch images, co-added images, co-added source catalogs, and associated products and services assembled over the first 3 yr of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (2013 August to 2016 February) by the Dark Energy Camera mounted on the 4 m Blanco telescope at the Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ∼5000 deg2 of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.″90 FWHM, a photometric precision of <1% in all bands, and an astrometric precision of 151 . The median co-added catalog depth for a 1.″95 diameter aperture at signal-to-noise ratio (S/N) = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44. DES DR1 includes nearly 400 million distinct astronomical objects detected in ∼10,000 co-add tiles of size 0.534 deg2 produced from ∼39,000 individual exposures. Benchmark galaxy and stellar samples contain ∼310 million and ∼80 million objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive co-add image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

The gravitational-wave signal GW190521 is consistent with a binary black hole (BBH) merger source at redshift 0.8 with unusually high component masses, 85-14+21 M o˙ and 66-18+17 M o˙, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range 65-120 M o˙. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger (142-16+28 M o˙) classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular BBH coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of the coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be 0.13-0.11+0.30 Gpc-3 yr-1. We discuss the astrophysical implications of GW190521 for stellar collapse and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescences, or via hierarchical mergers of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

We present a search for subsolar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational-wave candidates consistent with subsolar mass ultracompact binaries with at least one component between 0.2  M_{⊙}-1.0  M_{⊙}. We use the null result to constrain the binary merger rate of (0.2  M_{⊙}, 0.2  M_{⊙}) binaries to be less than 3.7×10^{5}  Gpc^{-3} yr^{-1} and the binary merger rate of (1.0  M_{⊙}, 1.0  M_{⊙}) binaries to be less than 5.2×10^{3}  Gpc^{-3} yr^{-1}. Subsolar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with subsolar masses. Assuming a particular primordial black hole (PBH) formation model, we constrain a population of merging 0.2  M_{⊙} black holes to account for less than 16% of the dark matter density and a population of merging 1.0  M_{⊙} black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict subsolar mass binaries.