On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

We develop a new formalism to study nonlinear evolution in the growth of large-scale structure, by following the dynamics of gravitational clustering as it builds up in time. This approach is conveniently represented by Feynman diagrams constructed in terms of three objects: the initial conditions (e.g. perturbation spectrum), the vertex (describing non-linearities) and the propagator (describing linear evolution). We show that loop corrections to the linear power spectrum organize themselves into two classes of diagrams: one corresponding to mode-coupling effects, the other to a renormalization of the propagator. Resummation of the latter gives rise to a quantity that measures the memory of perturbations to initial conditions as a function of scale. As a result of this, we show that a well-defined (renormalized) perturbation theory follows, in the sense that each term in the remaining mode-coupling series dominates at some characteristic scale and is subdominant otherwise. This is unlike standard perturbation theory, where different loop corrections can become of the same magnitude in the nonlinear regime. In companion papers we compare the resummation of the propagator with numerical simulations, and apply these results to the calculation of the nonlinear power spectrum. Remarkably, the expressions in renormalized perturbation theory can be written in a way that closely resembles the halo model.

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 Λ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.026 −0.020 and Ωm = 0.267 +0.030 −0.017 for ΛCDM; for wCDM, we find S8 = 0.782 +0.036 −0.024 , Ωm = 0.284 +0.033 −0.030 , and w = −0.82 +0.

We study the nonlinear evolution of the baryon acoustic oscillations (BAO) in the dark matter power spectrum and correlation function using renormalized perturbation theory (RPT). In a previous paper we showed that RPT successfully predicts the damping of acoustic oscillations; here we extend our calculation to the enhancement of power due to mode-coupling. We show that modecoupling generates additional oscillations that are out of phase with those in the linear spectrum, leading to shifts in the scales of oscillation nodes defined with respect to a smooth spectrum. When Fourier transformed, these out of phase oscillations induce percent-level shifts in the acoustic peak of the two-point correlation function. We present predictions for these shifts as a function of redshift; these should be considered as a robust lower limit to the more realistic case that includes in addition redshift distortions and galaxy bias. We show that these nonlinear effects occur at very large scales, leading to a breakdown of linear theory at scales much larger than commonly thought. We discuss why virialized halo profiles are not responsible for these effects, which can be understood from basic physics of gravitational instability. Our results are in excellent agreement with numerical simulations, and can be used as a starting point for modeling BAO in future observations. To meet this end, we suggest a simple physically motivated model to correct for the shifts caused by mode-coupling.