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.

Sagittarius 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. We investigated the variability and polarization properties of using ALMA observations during the 2018 Event Horizon Telescope campaign. 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 on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. Chandra 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. Our results offer new constraints on the physical mechanisms driving variability in and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes.

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.

Recent multiwavelength observations of M87* revealed a high-energy γ -ray flare without a millimeter counterpart. We present a theoretical study of potential millimeter flares in M87*, using general relativistic magnetohydrodynamical simulations with varying black hole (BH) spins and magnetic field configurations. The emergence of a millimeter flare is strongly influenced by both spin and magnetic structure, with limited sensitivity to the electron distribution. We model the intensity light curve with a damped random walk and compare the characteristic timescale ( τ ) with recent Submillimeter Array observations, finding that simulated τ exceeds observed values by over an order of magnitude. For BH spin a = +0.5, we identify a distinct flare followed by an order-of-magnitude flux drop. All Stokes parameters vary near the flare, including a sign reversal in the electric vector position angle. Most β _m modes remain stable, but the EB -correlation phase is sensitive to both the flare peak and decay. We examine polarimetric signatures in photon subrings, focusing on modes n _s = (0, 1). The n _s = 0 signal closely matches the full image, while n _s = 1 exhibits distinct behavior, highlighting the potential of space very long baseline interferometry to isolate subring features. Magnetic reconnection near the disk plane triggers a flux eruption that injects magnetized plasma into the jet funnel, enhancing outflow speeds and temporarily disrupting jet collimation before gradual recollimation. These results suggest transient variability near flares encodes key information about BH spin and magnetic structure, offering new probes of active galactic nuclei.

Context. The archetypal blazar 3C 279 has a prominent relativistic jet and strong broadband variability across the electromagnetic spectrum. In April 2017, the Event Horizon Telescope (EHT) observed 3C 279 with an unprecedented angular resolution of about 20 μas, accompanied by one of the most extensive quasi-simultaneous multiwavelength (MWL) campaigns ever conducted, spanning from radio to TeV γ -ray energies. Aims. Taking advantage of this comprehensive MWL dataset, we investigated the physical processes governing 3C 279, with a particular focus on the formation, collimation, and acceleration of its relativistic jet and on the origin of its high-energy emission, including the underlying particle-acceleration mechanisms. Methods. We analyzed individual observations and multiband light curves. We also constructed a new quasi-simultaneous spectral energy distribution covering frequencies from the radio band to very high-energy (VHE) γ rays. We further performed a phenomenological modeling using the turbulent extreme multi-zone (TEMZ) model to constrain the fundamental physical properties of the source. Results. The EHT observations reveal a clear flux increase in the innermost core between April 5 and 11, 2017. Over a broader time span, radio observations at longer wavelengths reveal concurrent enhancements in the core flux and polarization around mid-April, coinciding with the ejection of a superluminal knot moving at (25 ± 2) c . Record UV–optical flares with strong polarization variability occurred in late March, followed by high-energy γ -ray activity that declined before the end of the EHT observing period. During this time, the source remained in a low X-ray state and exhibited no detectable VHE emission. Conclusions. The results of the TEMZ modeling indicate that the broadband spectrum and variability of 3C 279 might be explained with a jet scenario in which turbulent plasma cells are compressed by a stationary conical shock. Nonetheless, alternative interpretations, such as magnetic reconnection or a moving shock-in-jet event, remain possible. This coordinated MWL campaign advances our understanding of the origin of the jet and γ -ray emission in the blazar 3C 279, and it also provides a comprehensive publicly available dataset that will serve as a valuable reference for future studies.

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 divided by a(*)divided by less than or similar to 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.

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.

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage.

The 2021 Event Horizon Telescope observations resolve the innermost jet region of the blazar 3C279 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 10c, requiring bulk Lorentz factors greater than 10.3 and constraining viewing angles to extremely small values (smaller than one degree). Rest-frame brightness temperatures are systematically low (between 10^9 and 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 3C279.

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 a 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 how these trivial closure phases become nonzero with a 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, which is consistent with the jet observed at lower frequencies.

Recent Event Horizon Telescope results have demonstrated unique and transformative science in gravitational physics and black hole astrophysics enabled by event-horizon-scale imaging of supermassive black holes (SMBHs). Nevertheless, the angular resolution of current ground-based very long baseline interferometry (VLBI) arrays limits such studies to only two sources, precluding systematic investigations of horizon-scale emission across a nearby SMBH population. The proposed Black Hole Explorer (BHEX), a millimeter/submillimeter space VLBI mission, would overcome this limitation by delivering substantially higher angular resolution. Here, we present a series of simulated observations to assess a population of nearby horizon-scale targets accessible with BHEX. Based on a recently developed SMBH number density model, we find that BHEX could infer black hole masses for 70-90 sources from size measurements, constrain magnetic field structures through linear polarization imaging for 20-30 sources, and resolve black hole shadows for 20-25 sources. Targeted observations of 50 nearby SMBHs are expected to yield measurements for 30 source sizes and 10 shadows and linear-polarization patterns. These projections are supported by detailed imaging simulations of general relativistic magnetohydrodynamic (GRMHD) models for eleven nearby SMBHs. Together, our results highlight BHEX as a powerful facility for revealing the demographics of SMBH properties across diverse accretion states, radio loudness, host galaxy environments, and viewing geometries.

Event horizons are a defining feature of black holes. Consequently, there have been many efforts to probe their existence in astrophysical black hole candidates, spanning ten orders of magnitude in mass. Nevertheless, horizons remain an obstacle to unifying general relativity and quantum mechanics, most sharply presented by the information paradox. This has motivated a proliferation of horizonless alternatives (black hole foils) that avoid event horizons and are therefore benign. We show that for typical accreting astronomical targets, largely independent of a foil's underlying microphysics, a horizonless compact surface will generically be ensconced within an optically thick, scattering dominated baryonic settling layer that efficiently reprocesses the kinetic energy of infalling matter into observable thermal emission. The emergent photosphere luminosity is driven toward the accretion-powered equilibrium value and is only weakly sensitive to the foil redshift. These atmospheres are convectively stable and naturally imply that the emitting photosphere forms at modest redshift even when the surface redshift is extreme. Moreover, local gas-surface interaction provides a microphysical lower bound on the effective base temperature, insulating the atmosphere from arbitrarily cold foils. The unknown properties of the foil enter only through local boundary conditions controlling baryon processing and thermal coupling at the surface, making the solutions broadly applicable to horizonless alternatives that do not invoke significant additional nonlocal interactions. Thus, under minimal assumptions (GR exterior and local surface interactions), horizonless foils are generically observationally exposed: the absence of a thermal photosphere directly constrains or rules out broad classes of such models.

The Event Horizon Telescope (EHT) enables the exploration of black hole accretion flows at event-horizon scales. Fitting ray-traced physical models to EHT observations requires the generation of synthetic images, a task that is computationally demanding. This study leverages ALINet , a generative machine learning model, to efficiently produce radiatively inefficient accretion flow (RIAF) images as a function of the specified physical parameters. ALINet has previously been shown to be able to interpolate black hole images and their associated physical parameters after training on a computationally tractable set of library images. We utilize this model to estimate the uncertainty introduced by a number of anticipated unmodeled physical effects, including interstellar scattering and intrinsic source variability. We then use this to calibrate physical parameter estimates and their associated uncertainties from RIAF model fits to mock EHT data via a library of general relativistic magnetohydrodynamics models.

Context. The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN. Aims. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Konigl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion. Methods. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132-1654, J0006-0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Results. Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factor delta proportional to r(<= 0.5). This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radius B proportional to r(-3), hinting at jet acceleration or efficient magnetic energy dissipation.

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 mu 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 similar to 5%, compared to similar 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 less than or similar to 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*.

The Event Horizon Telescope (EHT) Collaboration recently published the first images of the supermassive black holes in the cores of the Messier 87 and Milky Way galaxies. These observations have provided a new means to study supermassive black holes and probe physical processes occurring in the strong-field regime. We review the prospects of future observations and theoretical studies of supermassive black hole systems. Current ground-based very-long-baseline interferometry (VLBI) arrays like the EHT and proposed future extensions like the next-generation Event Horizon Telescope will greatly enhance the capabilities of black-hole imaging interferometry. These enhancements will open up several previously inaccessible avenues of investigation, thereby providing important new insights into the properties of supermassive black holes and their environments. This review describes the current state of knowledge for five key science cases, summarising the unique challenges and opportunities for fundamental physics investigations that future mm/sub-mm VLBI developments will enable.

Recent multi-wavelength observations of M87* \citep{2024A&A...692A.140A} revealed a high-energy $\gamma$-ray flare without a corresponding millimeter counterpart. We present a theoretical polarimetric study to evaluate the presence and nature of a potential millimeter flare in M87*, using a suite of general relativistic magnetohydrodynamical simulations with varying black hole (BH) spins and magnetic field configurations. We find that the emergence of a millimeter flare is strongly influenced by both spin and magnetic structure, with limited sensitivity to the electron distribution (thermal vs. non-thermal). We model the intensity light curve with a damped random walk (DRW) and compare the characteristic timescale ($\tau$) with recent SMA observations, finding that the simulated $\tau$ exceeds observed values by over an order of magnitude. In a flaring case with BH spin a=+0.5, we identify a distinct millimeter flare followed by an order-of-magnitude flux drop. All Stokes parameters show variability near the flare, including a sign reversal in the electric vector position angle. While most $\beta_m$ modes remain stable, the $EB$-correlation phase is highly sensitive to both the flare peak and decay. We examine polarimetric signatures in photon sub-rings, focusing on modes ns=0 and ns=1. The ns=0 signal closely matches the full image, while ns=1 reveals distinct behaviors, highlighting the potential of space VLBI to isolate sub-ring features. Finally, we analyze the magnetic and velocity field evolution during the flare, finding that magnetic reconnection weakens during the flux decay, and the clockwise velocity flow transitions into an outflow-dominated regime. These results suggest that transient radio variability near flares encodes key information about black hole spin and magnetic field structure, offering a novel probe into the physics of active galactic nuclei.

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of tau = 0.08(-0.02)(+0.03) with a position angle of xi = 50.1(-7.6)(+6.2) degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.

Aims. We investigated the polarization and Faraday properties of Messier 87 (M87) and seven other radio-loud active galactic nuclei (AGNs) at λ 0.87 mm (345 GHz) using the Atacama Large Millimeter/submillimeter Array (ALMA). Our goal was to characterize the linear polarization (LP) fractions, measure Faraday rotation measures (RMs), and examine the magnetic field structures in the emission regions of these AGNs. Methods. We conducted full-polarization observations as part of the ALMA Band 7 very long baseline interferometry (VLBI) commissioning during the April 2021 Event Horizon Telescope (EHT) campaign. We analyzed the LP fractions and RMs to assess the nature of Faraday screens and magnetic fields in the submillimeter emission regions. Results. We find LP fractions between 1% and 17% and RMs exceeding 10 5 rad m −2 , which are 1–2 orders of magnitude higher than typically observed at longer wavelengths ( λ >3 mm). This suggests denser Faraday screens or stronger magnetic fields. Additionally, we present the first submillimeter polarized images of the M87 jet and the observed AGNs, revealing RM gradients and sign reversals in the M87 jet indicative of a kiloparsec-scale helical magnetic field structure. Conclusions. Our results provide essential constraints for calibrating, analyzing, and interpreting VLBI data from the EHT at 345 GHz, representing a critical step toward submillimeter VLBI imaging.

The Event Horizon Telescope (EHT) Collaboration recently published the first images of the supermassive black holes in the cores of the Messier 87 and Milky Way galaxies. These observations have provided a new means to study supermassive black holes and probe physical processes occurring in the strong-field regime. We review the prospects of future observations and theoretical studies of supermassive black hole systems with the next-generation Event Horizon Telescope (ngEHT), which will greatly enhance the capabilities of the existing EHT array. These enhancements will open up several previously inaccessible avenues of investigation, thereby providing important new insights into the properties of supermassive black holes and their environments. This review describes the current state of knowledge for five key science cases, summarising the unique challenges and opportunities for fundamental physics investigations that the ngEHT will enable.