Elisabetta Liuzzo

In 2026, the Event Horizon Telescope (EHT) carried out the first time-lapse monitoring campaign comprising a series of 18 observations of the 6.1 6 − 0.57 + 0.71 × 1 0 9 M ⊙ black hole at the center of M87 with a cadence of approximately twice per week over a few months, matched to the expected dynamical timescale. Here we present a first-look image from the campaign, reconstructed from a single four-minute scan obtained on 2026 April 3. The resulting image displays the now-familiar asymmetric ring morphology first revealed by the 2017 EHT observations and subsequently confirmed by observations obtained in 2018 and 2021. While representing only a very small subset of the full program, this image demonstrates a significant advance in the operational capabilities of the EHT. Through improvements in array operations, data transport and processing, the collaboration has reduced the interval between observation and imaging to only a few months, establishing the EHT as a facility capable of routinely undertaking time-domain studies of the dynamic universe at extremely high resolutions of ∼20 μ as.

The 2021 Event Horizon Telescope observations resolve the innermost jet region of the blazar 3C 279 with unprecedented detail. The reconstructed images consistently reveal a compact core elongated nearly orthogonal to the large-scale jet axis. This rarely observed morphology recurs across multiple epochs and from 22–230 GHz and is therefore intrinsic rather than an imaging artifact. Geometric model fitting identifies several components with apparent speeds up to β _app ∼ 10, requiring bulk Lorentz factors Γ ≳ 10.3 and constraining viewing angles to extremely small values (≲1°). Rest-frame brightness temperatures are systematically low ( T _B ∼ 10 ⁹ –10 ¹⁰ K), consistent with optically thin emission at 230 GHz. These results suggest that the jet bends toward the observer on sub-parsec scales, producing strong relativistic beaming. Possible drivers of the observed jet bending and temporal evolution include the jet’s interaction with the interstellar medium, kink or Kelvin-Helmholtz instabilities, magnetic reconnection near the horizon, or binary-induced precession. However, the current temporal coverage of VLBI data remains insufficient to distinguish between these mechanisms. Continued multifrequency VLBI monitoring will be essential to constraining the dynamics and geometry of the jet base in 3C 279.

Context . The Atacama Large Millimeter/submillimeter Array (ALMA) is the most sensitive interferometric array at millimeter and submillimeter wavelengths. Through the ALMA Phasing System (APS), it can participate in global very long baseline interferometry (VLBI) arrays, significantly enhancing their sensitivity and resolution. However, processing and analyzing the ALMA data obtained in APS mode during VLBI observations remains a complex task, requiring specialized expertise and time-consuming calibration and imaging procedures. Aims . In this paper, we present “VAPOLA”—the first online, multi-epoch, multiband repository of high-level data products from ALMA observations of active galactic nuclei (AGNs) and Sgr A* during global VLBI campaigns. Methods . VAPOLA is built around an automated pipeline that processes fully calibrated ALMA data from the quality assurance level 2 (QA2) process and generates science-ready products with minimal user intervention. Results . The repository includes fully calibrated interferometric visibilities, full-Stokes images across individual and combined spectral windows, polarimetric and spectral index maps, and tabulated polarimetric parameters from visibility-domain polarization fitting. By offering ready-to-use data through a user-friendly web portal, VAPOLA enables nonexpert users to perform advanced science analyses without needing an in-depth knowledge of ALMA polarization calibration and imaging procedures. Conclusions . This resource facilitates a broad range of scientific investigations, including the characterization of magnetic field properties in accretion flows and relativistic jets, the structure and kinematics of dusty and molecular tori in AGNs, and absorption studies of the interstellar medium toward the Galactic center. In addition, the dataset provides source-integrated parameters and calibration metadata essential for refining VLBI calibration and imaging workflows as well as for placing robust observational constraints on theoretical models of supermassive black holes and their environments.

Context. Sagittarius A* (Sgr A ^* ), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. Aims. We investigated the variability and polarization properties of Sgr A ^* using ALMA observations during the 2018 Event Horizon Telescope campaign. Methods. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Results. We find low variability in total intensity ( σ / μ < 10%), but significantly higher variability in linear and circular polarization (∼30% and ∼50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between −2 and −3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Conclusions. Our results offer new constraints on the physical mechanisms driving variability in Sgr A ^* , and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

Event Horizon Telescope (EHT) images of the supermassive black hole M87* depict an asymmetric ring of emission. General relativistic magnetohydrodynamic (GRMHD) models of M87* and its accretion disk predict that the amplitude and location of the ring’s peak brightness asymmetry should fluctuate due to turbulence in the source plasma. We compare the observed distribution of brightness asymmetry amplitudes to the simulated distribution in GRMHD models, across varying black hole spin a _* . We show that, for strongly magnetized (MAD) models, three epochs of EHT data marginally disfavor ∣ a _* ∣ ≲ 0.2. This is consistent with the Blandford–Znajek model for M87’s jet, which predicts that M87* should have nonzero spin. We show quantitatively how future observations could improve spin constraints and discuss how improved spin constraints could distinguish between differing jet-launching mechanisms and black hole growth scenarios.

We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the significantly enhanced ( u,v ) coverage in the 2021 Event Horizon Telescope (EHT) observations. The integration of the 12−m Kitt Peak Telescope (USA) and NOEMA (France) stations into the array introduces two critical intermediate-length baselines to SMT (USA) and IRAM 30−m (Spain), providing sensitivity to emission structures at spatial scales of ∼250 μas and ∼2500 μas (∼ 0.02 pc and ∼ 0.02 pc). Without these new baselines, previous EHT observations of the source in 2017 and 2018 lacked the capability to constrain emission on large scales, where a “missing flux” of order ∼1 Jy is expected to reside. To probe these scales, we analyzed closure phases–robust against station-based gain calibration errors–and model the jet base emission using a simple Gaussian component offset from the compact ring emission at spatial separations > 100 μas. Our analysis revealed a Gaussian feature centered at (ΔRA ≈ 320 μas, ΔDec. ≈ 60 μ as), projected separation of ≈ 5500 AU, with an estimated flux density of only ∼60 mJy, implying that most of the missing flux identified in previous EHT studies had to originate from different, larger scales. Brighter emission at the relevant spatial scales is firmly ruled out, and the data do not favor more complex models. This component aligns with the inferred position of the large-scale jet and is therefore physically consistent with the emission of the jet base. While our findings point to detectable jet base emission at 230 GHz, the limited coverage provided by only two intermediate baselines limits our ability to robustly reconstruct its morphology. Consequently, we treated the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate baseline coverage will be essential to constrain the structure and nature of this component with higher precision.

In Very-Long Baseline Interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show with an expansion about co-located stations, how these trivial closure phases become non-zero with brightness distribution on smaller scales than their short baseline would suggest. When applied to sources that are made up of a bright compact and large-scale diffuse component, the trivial closure phases directly measure the centroid relative to the compact source and higher-order image moments. We present a technique to measure these image moments with minimal model assumptions and validate it on synthetic Event Horizon Telescope (EHT) data. We then apply this technique to 2017 and 2018 EHT observations of M87* and find a weak preference for extended emission in the direction of the large-scale jet. We also apply it to 2021 EHT data and measure the source centroid about 1 mas northwest of the compact ring, consistent with the jet observed at lower frequencies.

We report three epochs of polarized images of M87* at 230 GHz using data from the Event Horizon Telescope (EHT) taken in 2017, 2018, and 2021. The baseline coverage of the 2021 observations is significantly improved through the addition of two new EHT stations: the 12 m Kitt Peak Telescope and the Northern Extended Millimetre Array (NOEMA). All observations result in images dominated by a bright, asymmetric ring with a persistent diameter of 43.9 ± 0.6 μas, consistent with expectations for lensed synchrotron emission encircling the apparent shadow of a supermassive black hole. We find that the total intensity and linear polarization of M87* vary significantly across the three epochs. Specifically, the azimuthal brightness distribution of the total intensity images varies from year to year, as expected for a stochastic accretion flow. However, despite a gamma-ray flare erupting in M87 quasi-contemporaneously to the 2018 observations, the 2018 and 2021 images look remarkably similar. The resolved linear polarization fractions in 2018 and 2021 peak at ∼5%, compared to ∼15% in 2017. The spiral polarization pattern on the ring also varies from year to year, including a change in the electric vector position angle helicity in 2021 that could reflect changes in the magnetized accretion flow or an external Faraday screen. The improved 2021 coverage also provides the first EHT constraints on jet emission outside the ring, on scales of ≲1 mas. Overall, these observations provide strong proof of the reliability of the EHT images and probe the dynamic properties of the horizon-scale accretion flow surrounding M87*.