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 make use of a new subgrid model of turbulent mixing to accurately follow the cosmological evolution of the first stars, the mixing of their supernova (SN) ejecta, and the impact on the chemical composition of the Galactic Halo. Using the cosmological adaptive mesh refinement code ramses, we implement a model for the pollution of pristine gas as described in Pan et al. Tracking the metallicity of Pop III stars with metallicities below a critical value allows us to account for the fraction of Z < Z crit stars formed even in regions in which the gas's average metallicity is well above Z crit . We demonstrate that such partially-mixed regions account for 0.5 to 0.7 of all Pop III stars formed up to z = 5. Additionally, we track the creation and transport of "primordial metals" (PM) generated by Pop III SNe. These neutron-capture deficient metals are taken up by second-generation stars and likely lead to unique abundance signatures characteristic of carbon-enhanced, metal-poor (CEMP-no) stars. As an illustrative example, we associate primordial metals with abundance ratios used by Keller et al.

Direct observational searches for Population III (Pop III) stars at high redshift are faced with the question of how to select the most promising targets for spectroscopic follow-up. To help answer this, we use a large-scale cosmological simulation, augmented with a new subgrid model that tracks the fraction of pristine gas, to follow the evolution of high-redshift galaxies and the Pop III stars they contain. We generate rest-frame ultraviolet (UV) luminosity functions for our galaxies and find that they are consistent with current z ≥ 7 observations. Throughout the redshift range 7 ≤ z ≤ 15 we identify "Pop III-bright" galaxies as those with at least 75% of their flux coming from Pop III stars. While less than 1% of galaxies brighter than m UV,AB = 31.4 mag are Pop III-bright in the range 7 ≤ z ≤ 8, roughly 17% of such galaxies are Pop III-bright at z = 9, immediately before reionization occurs in our simulation. Moving to z = 10, m UV,AB = 31.4 mag corresponds to larger, more luminous galaxies and the Pop III-bright fraction falls off to 5%. Finally, at the highest redshifts, a large fraction (29% at z = 14 and 41% at z = 15) of all galaxies are Pop III-bright regardless of magnitude. While m UV,AB = 31.4 mag galaxies are extremely rare during this epoch, we find that 13% of galaxies at z = 14 are Pop III-bright with m UV,AB ≤ 33 mag, an intrisic magnitude within reach of the James Webb Space Telescope using lensing. Thus, we predict that the best redshift to search for luminous Pop III-bright galaxies is just before reionization, while lensing surveys for fainter galaxies should push to the highest redshifts possible.

We use combined South Pole Telescope (SPT)+Planck temperature maps to analyze the circumgalactic medium (CGM) encompassing 138,235 massive, quiescent 0.5 ≤ z ≤ 1.5 galaxies selected from data from the Dark Energy Survey (DES) and Wide-Field Infrared Survey Explorer (WISE). Images centered on these galaxies were cut from the 1.85 arcmin resolution maps with frequency bands at 95, 150, and 220 GHz. The images were stacked, filtered, and fit with a graybody dust model to isolate the thermal Sunyaev–Zel’dovich (tSZ) signal, which is proportional to the total energy contained in the CGM of the galaxies. We separated these M ⋆ = 1010.9 M ⊙–1012 M ⊙ galaxies into 0.1 dex stellar mass bins, detecting tSZ per bin up to 5.6σ and a total signal-to-noise ratio of 10.1σ. We also detect dust with an overall signal-to-noise ratio of 9.8σ, which overwhelms the tSZ at 150 GHz more than in other lower-redshift studies. We corrected for the 0.16 dex uncertainty in the stellar mass measurements by parameter fitting for an unconvolved power-law energy-mass relation, E therm = E therm , peak M ⋆ / M ⋆ , peak α , with the peak stellar mass distribution of our selected galaxies defined as M ⋆,peak = 2.3 × 1011 M ⊙. This yields an E therm , peak = 5.98 − 1.00 + 1.02 × 10 60 erg and α = 3.77 − 0.74 + 0.60 . These are consistent with z ≈ 0 observations and within the limits of moderate models of active galactic nucleus feedback. We also computed the radial profile of our full sample, which is similar to that recently measured at lower-redshift by Schaan et al.

We study the observational consequences of several unknown properties of Population III (Pop III) stars using large-scale cosmological simulations that include a subgrid model to track the unresolved mixing of pollutants. Varying the value of the critical metallicity that marks the boundary between Pop III and Population II (Pop II) star formation across 2 dex has a negligible effect on the fraction of Pop III stars formed and the subsequent fraction of Pop III flux from high-redshift galaxies. However, adopting a log normal initial mass function (IMF) for Pop III stars, in place of a baseline Salpeter IMF, results in a Pop III star formation rate density (SFRD) that is 1/4 of the baseline rate. The flux from high-redshift galaxies modeled with this IMF is highly bimodal, resulting in a tiny fraction of z ≤ 8 galaxies with more than 75% of their flux coming from Pop III stars. However, at z = 9, right before reionization in our simulations, ≈ 20% of galaxies are Pop III-bright with m UV ≤ 31.4 mag and at least 75% of their flux generated by Pop III stars . Additionally, the log normal Pop III IMF results in a population of carbon enhanced, metal poor stars in reasonable agreement with MW halo observations. Our analysis supports the conclusion that the Pop III IMF was dominated by stars in the 20-120 M range that generate SN with carbon-enhanced ejecta.