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 $\sim 1.7\,{\rm{s}}$ with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg² 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}_{\odot }$. 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 $\sim 40\,{\rm{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 $\sim 9$ and $\sim 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.

The Event Horizon Telescope observed the horizon-scale synchrotron emission region around the Galactic center supermassive black hole, Sagittarius A* (Sgr A*), in 2017. These observations revealed a bright, thick ring morphology with a diameter of 51.8 ± 2.3 μas and modest azimuthal brightness asymmetry, consistent with the expected appearance of a black hole with mass M ≈ 4 × 10⁶M_⊙. From these observations, we present the first resolved linear and circular polarimetric images of Sgr A*. The linear polarization images demonstrate that the emission ring is highly polarized, exhibiting a prominent spiral electric vector polarization angle pattern with a peak fractional polarization of ∼40% in the western portion of the ring. The circular polarization images feature a modestly (∼5%–10%) polarized dipole structure along the emission ring, with negative circular polarization in the western region and positive circular polarization in the eastern region, although our methods exhibit stronger disagreement than for linear polarization. We analyze the data using multiple independent imaging and modeling methods, each of which is validated using a standardized suite of synthetic data sets. While the detailed spatial distribution of the linear polarization along the ring remains uncertain owing to the intrinsic variability of the source, the spiraling polarization structure is robust to methodological choices. The degree and orientation of the linear polarization provide stringent constraints for the black hole and its surrounding magnetic fields, which we discuss in an accompanying publication.

In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

We present high-definition observations with the James Webb Space Telescope (JWST) of >1000 Cepheids in a geometric anchor of the distance ladder, NGC 4258, and in five hosts of eight Type Ia supernovae, a far greater sample than previous studies with JWST. These galaxies individually contain the largest samples of Cepheids, an average of >150 each, producing the strongest statistical comparison to those previously measured with the Hubble Space Telescope (HST) in the near-infrared (NIR). They also span the distance range of those used to determine the Hubble constant with HST, allowing us to search for a distance-dependent bias in HST measurements. The superior resolution of JWST negates crowding noise, the largest source of variance in the NIR Cepheid period–luminosity relations (Leavitt laws) measured with HST. Together with the use of two epochs to constrain Cepheid phases and three filters to remove reddening, we reduce the dispersion in the Cepheid P–L relations by a factor of 2.5. We find no significant difference in the mean distance measurements determined from HST and JWST, with a formal difference of −0.01 ± 0.03 mag. This result is independent of zero-points and analysis variants including metallicity dependence, local crowding, choice of filters, and slope of the relations. We can reject the hypothesis of unrecognized crowding of Cepheid photometry from HST that grows with distance as the cause of the "Hubble tension" at 8.2σ, i.e., greater confidence than that of the Hubble tension itself. We conclude that errors in photometric measurements of Cepheids across the distance ladder do not significantly contribute to the tension.

van Dokkum et al. reported the discovery of JWST-ER1, a strong lensing object at redshift z ≈ 2, using data from the James Webb Space Telescope. The lens mass within the Einstein ring is 5.9 times higher than the expected stellar mass from a Chabrier initial mass function, indicating a high dark matter density. In this work, we show that a cold dark matter halo, influenced by gas-driven adiabatic contraction, can account for the observed lens mass. We interpret the measurement of JWST-ER1 in the self-interacting dark matter scenario and show that the cross section per particle mass σ/m ≈ 0.1 cm² g⁻¹ is generally favored. Intriguingly, σ/m ≈ 0.1 cm² g⁻¹ can also be consistent with the strong lensing observations of early-type galaxies at redshift z ≈ 0.2, where adiabatic contraction is not observed overall.

When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 ± 3 μas, which is circular and encompasses a central depression in brightness with a flux ratio ≳10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 ± 0.7) × 10⁹ M_⊙. Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A*), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of λ = 1.3 mm. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of 51.8 ± 2.3 μas (68% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A* are consistent with the expected appearance of a Kerr black hole with mass ∼4 × 10⁶ M_⊙, which is inferred to exist at this location based on previous infrared observations of individual stellar orbits, as well as maser proper-motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination (i > 50°), as well as nonspinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of 10³–10⁵ gravitational radii to event-horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87* shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.

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_⊙ and ${66}_{-18}^{+17}$ M_⊙, 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_⊙. 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_⊙) 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}\,{{\rm{Gpc}}}^{-3}\,{{\rm{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 report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings–Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of 10¹⁴, and this same model is favored over an uncorrelated common power-law spectrum model with Bayes factors of 200–1000, depending on spectral modeling choices. We have built a statistical background distribution for the latter Bayes factors using a method that removes interpulsar correlations from our data set, finding p = 10⁻³ (≈3σ) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of interpulsar correlations yields p = 5 × 10⁻⁵ to 1.9 × 10⁻⁴ (≈3.5σ–4σ). Assuming a fiducial f^−2/3 characteristic strain spectrum, as appropriate for an ensemble of binary supermassive black hole inspirals, the strain amplitude is ${2.4}_{-0.6}^{+0.7}\times {10}^{-15}$ (median + 90% credible interval) at a reference frequency of 1 yr⁻¹. The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings–Downs correlations points to the gravitational-wave origin of this signal.

We present the first Event Horizon Telescope (EHT) images of M87, using observations from April 2017 at 1.3 mm wavelength. These images show a prominent ring with a diameter of ∼40 μas, consistent with the size and shape of the lensed photon orbit encircling the "shadow" of a supermassive black hole. The ring is persistent across four observing nights and shows enhanced brightness in the south. To assess the reliability of these results, we implemented a two-stage imaging procedure. In the first stage, four teams, each blind to the others' work, produced images of M87 using both an established method (CLEAN) and a newer technique (regularized maximum likelihood). This stage allowed us to avoid shared human bias and to assess common features among independent reconstructions. In the second stage, we reconstructed synthetic data from a large survey of imaging parameters and then compared the results with the corresponding ground truth images. This stage allowed us to select parameters objectively to use when reconstructing images of M87. Across all tests in both stages, the ring diameter and asymmetry remained stable, insensitive to the choice of imaging technique. We describe the EHT imaging procedures, the primary image features in M87, and the dependence of these features on imaging assumptions.

Emission lines from Rydberg transitions are detected for the first time from a region close to the surface of Betelgeuse. The H30α line is observed at 231.905 GHz, with an FWHM ∼42 km s⁻¹ and extended wings. A second line at 232.025 GHz (FWHM ∼21 km s⁻¹ ), is modeled as a combination of Rydberg transitions of abundant low first ionization potential metals. Both H30α and the Rydberg combined line X30α are fitted by Voigt profiles, and collisional broadening with electrons may be partly responsible for the Lorentzian contribution, indicating electron densities of a few 10⁸ cm⁻³. X30α is located in a relatively smooth ring at a projected radius of 0.9× the optical photospheric radius R_⋆, whereas H30α is more clumpy, reaching a peak at ∼1.4 R_⋆. We use a semiempirical thermodynamic atmospheric model of Betelgeuse to compute the 232 GHz (1.29 mm) continuum and line profiles making simple assumptions. Photoionized abundant metals dominate the electron density, and the predicted surface of continuum optical depth unity at 232 GHz occurs at ∼1.3 R_⋆, in good agreement with observations. Assuming a Saha–Boltzmann distribution for the level populations of Mg, Si, and Fe, the model predicts that the X30α emission arises in a region of radially increasing temperature and turbulence. Inclusion of ionized C and non-LTE effects could modify the integrated fluxes and location of emission. These simulations confirm the identity of the Rydberg transition lines observed toward Betelgeuse and reveal that such diagnostics can improve future atmospheric models.

We present the first observations of in situ electron acceleration at corotating interaction region (CIR) shocks near 1 au, utilizing measurements from Wind and Magnetospheric Multiscale (MMS) mission in the interplanetary medium. As the forward (reverse) shock of the 2018 January CIR (the 2020 February CIR) moves from Wind at [206, 92, −7]R_E ([257, 25, 3]R_E) to MMS1 at [24, 2, 7]R_E ([25, 3, 0.5]R_E), the shock's thickness becomes 8 (3) times thinner, but the convective electric field E_drift gets weaker (stronger) along the shock; both the upstream and shocked suprathermal electrons exhibit a flatter flux energy spectrum, while the electron shock acceleration becomes less (more) significant. For the shocked suprathermal electrons with significant flux enhancement, the flux ratio across the shock appears to peak in the direction perpendicular to the magnetic field. Therefore, the CIR shock acceleration of solar wind suprathermal electrons at 1 au exhibits an efficiency increasing with the E_drift strength. These results also suggest that such acceleration through the interplanetary medium can contribute to the formation of solar wind suprathermal electrons.

The unprecedented medium-resolution (R_λ ∼ 1500–3500) near- and mid-infrared (1–18 μm) spectrum provided by JWST for the young (140 ± 20 Myr) low-mass (12–20 M_Jup) L–T transition (L7) companion VHS 1256 b gives access to a catalog of molecular absorptions. In this study, we present a comprehensive analysis of this data set utilizing a forward-modeling approach applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: T_eff, log(g), [M/H], C/O, γ, f_sed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS 1256 b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a T_eff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS 1256 b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models.

The study of temperate sub-Neptunes is the new frontier in exoplanetary science. A major development in the past year has been the first detection of carbon-bearing molecules in the atmosphere of a temperate sub-Neptune, K2-18 b, a possible Hycean world, with the James Webb Space Telescope (JWST). The JWST is poised to characterize the atmospheres of several other such planets, with important implications for planetary processes in the temperate regime. Meanwhile, ground-based high-resolution spectroscopy has been highly successful in detecting chemical signatures of giant exoplanets, though low-mass planets have remained elusive. In the present work, we report the atmospheric reconnaissance of a temperate sub-Neptune, TOI-732 c, using ground-based high-resolution transmission spectroscopy. The long orbital period and the low systemic velocity result in a low planetary radial velocity during transit, making this system a valuable test bed for high-resolution spectroscopy of temperate sub-Neptunes. We observe high-resolution time-series spectroscopy in the H and K bands during the planetary transit with the IGRINS instrument (R ∼ 45,000) on Gemini-South. Using observations from a single transit, we find marginal evidence (2.2σ) for the presence of methane (CH₄) in the atmosphere and no evidence for ammonia (NH₃) despite its strong detectability for a cloud-free H₂-rich atmosphere. We assess our findings using injection tests with different atmospheric scenarios and find them to be consistent with a high CH₄/NH₃ ratio and/or the presence of high-altitude clouds. Our results demonstrate the capability of Gemini-S/IGRINS for atmospheric characterization of temperate sub-Neptunes and the complementarity between space- and ground-based facilities in this planetary regime.

We model the spectral formation occurring in the binary X-ray pulsar (XRP) RX J0209.6−7427 during the 2019 super-Eddington outburst. Using a theoretical model previously developed by the authors, we are able to produce spectra that closely resemble the phase-averaged X-ray spectra observed using the Nuclear Spectroscopic Telescope Array and Insight-HXMT during low- and high-luminosity states of the outburst, respectively. The theoretical model simulates the accretion of fully ionized gas in a dipole magnetic field and includes a complete description of the radiation hydrodynamics, matter distribution, and spectral formation. Type II X-ray outbursts provide an opportunity to study accretion over a large range of luminosities for the same neutron star. The analysis performed here represents the first time both the outburst low and high states of an accretion-powered XRP are modeled using a physics-based model rather than standard phenomenological fitting with arbitrary mathematical functions. We find that the outer polar cap radius remains constant and the column is more fully filled with increasing luminosity, Comptonized bremsstrahlung dominates the formation of the phase-averaged X-ray spectrum, and a negative correlation exists between cyclotron centroid energy and luminosity, as expected. The super-Eddington nature of the outburst is rendered possible owing to the low scattering cross section for photons propagating parallel to the magnetic field. We also find that emission through the column top dominates in both the low and high states, implying that the pulse profiles should have a roughly sinusoidal shape, which agrees with observed properties of ultraluminous XRPs.