Showing 1–50 of 129 results for author: Wielgus, M

Full-polarization millimeter wavelength variability of Sagittarius A* during the 2018 EHT campaign

Authors: Ezequiel Albentosa-Ruiz, Jasmin E. Washington, Nicola Marchili, Iván Martí-Vidal, Ciriaco Goddi, Maciek Wielgus, Alejandro Mus, Angelo Ricarte, Daniel P. Marrone, León D. S. Salas, Yuhei Iwata, Douglas F. Carlos, Alexandra J. Tetarenko, Kotaro Moriyama, Vedant Dhruv, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball , et al. (250 additional authors not shown)

Abstract: Sagittarius A* (Srg 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 Srg A* using ALMA observations during the 2018 Event Horizon Telescope campaign. We analyzed high-cadence full-polarization light curves from ALMA at… ▽ More Sagittarius A* (Srg 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 Srg A* 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 Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. 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 Srg A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes. △ Less

Submitted 11 April, 2026; originally announced April 2026.

Journal ref: Astronomy & Astrophysics 708 (2026) A179

Radiative GRMHD simulations of puffy accretion discs: Numerical versus analytical models of sub-Eddington accretion

Authors: Debora Lančová, Maciek Wielgus, Marek Abramowicz, Agata Różańska, Włodek Kluźniak, Jiří Horák, David Abarca, Aleksander Sądowski, Gabriel Török

Abstract: A widely accepted picture of an accretion flow in the luminous soft spectral state of X-ray binary systems is a geometrically thin disc structure much like the classic analytic solution of Shakura \& Sunyaev. Although the analytic models are troubled by instabilities and miss important aspects of physics, such as magnetic fields, they are successfully used as a framework for interpreting observati… ▽ More A widely accepted picture of an accretion flow in the luminous soft spectral state of X-ray binary systems is a geometrically thin disc structure much like the classic analytic solution of Shakura \& Sunyaev. Although the analytic models are troubled by instabilities and miss important aspects of physics, such as magnetic fields, they are successfully used as a framework for interpreting observational data. Here, we compare the results of general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of optically thick, mildly sub-Eddington accretion on a stellar-mass black hole (the puffy disc) with established analytic and semi-analytic accretion models in the same regime. From the simulations, we find that the accretion flow is stabilised by the magnetic field, with a puffed-up, optically thick region resembling a warm corona surrounding a denser and cooler disc core. However, the stratified vertical structure of the disc significantly influences the observational picture of such a system. We analyse the inner disc structure, flow properties, effective viscosity, and inner edge position, and compare them to the predictions of standard models. We find that the simulated discs share some similarities with the models; however, they differ in several important aspects, most notably: the photosphere is geometrically thick, the inner edge is located closer to the central black hole than the analytic models assume, the surface density is significantly lower than analytically predicted, and the effective viscosity parameter is not constant but rises steeply in the innermost region. △ Less

Submitted 18 March, 2026; originally announced March 2026.

Comments: 12 pages, 8 figures, submitted to MNRAS

Universal behaviour of $α$-viscosity in black hole accretion discs

Authors: Marek A. Abramowicz, Axel Brandenburg, Jiří Horák, Debora Lančová, John C. Miller, Ewa Szuszkiewicz, Maciek Wielgus

Abstract: The Shakura-Sunyaev $α$-viscosity coefficient, defined as the ratio of total stress to total pressure, $α= \mathbb{T}/p$, played an important role in the development of the accretion disc theory in the early 1970s. The origin of turbulence that causes the stress $\mathbb{T}$ was unknown at that time. Shakura and Sunyaev assumed $α=$ const. Today we know that this was not quite realistic - the mode… ▽ More The Shakura-Sunyaev $α$-viscosity coefficient, defined as the ratio of total stress to total pressure, $α= \mathbb{T}/p$, played an important role in the development of the accretion disc theory in the early 1970s. The origin of turbulence that causes the stress $\mathbb{T}$ was unknown at that time. Shakura and Sunyaev assumed $α=$ const. Today we know that this was not quite realistic - the modern general relativistic magneto-hydrodynamic simulations (GRMHD) of black hole accretion discs revealed that $α$ changes by about an order of magnitude within the disc, being smaller far away from the black hole and larger in the plunging region close in. It was found that the behaviour of $α$ reflects some underlying, fundamental properties of the stress $\mathbb{T}$ itself. In particular, as argued by several authors, the stress must be zero at the black hole horizon. We notice that the stress calculated in GRMHD simulations by different authors, including us, has a maximum rather close to the location of the circular photon orbit. We propose a formula that accurately describes this universal behaviour of $α$ in terms of the "gyration radius'', a physical characteristic of rotation well known in Newtonian dynamics and in the black hole case uniquely defined by the Kerr space-time geometry. Analytic and semi-analytic models of black hole accretion discs provide an invaluable insight into fundamental physics, and the GRMHD simulations do not aspire to replace them. Rather, simulations could help to improve analytic models by making them more realistic. For example, our $α$-formula, deduced from the GRMHD simulations, may be handy in the construction of improved versions of thin and slim disc models. △ Less

Submitted 11 March, 2026; originally announced March 2026.

Comments: 8 pages, 3 figures, submitted to Astronomy&Astrophysics

Report number: NORDITA-2026-028

Locating the missing large-scale emission in the jet of M87* with short EHT baselines

Authors: Boris Georgiev, Paul Tiede, Sebastiano D. von Fellenberg, Michael Janssen, Iniyan Natarajan, Lindy Blackburn, Jongho Park, Erandi Chavez, Andrew T. West, Kotaro Moriyama, Jun Yi Koay, Hendrik Müller, Dhanya G. Nair, Avery E. Broderick, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko , et al. (258 additional authors not shown)

Abstract: 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,… ▽ More 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. △ Less

Submitted 19 January, 2026; originally announced January 2026.

Comments: 15 pages, 8 figures, Submitted to Astronomy & Astrophysics

Ring Asymmetry and Spin in M87*

Authors: Vadim Bernshteyn, Nicholas S. Conroy, Michi Bauböck, Paul Tiede, Abhishek V. Joshi, Ben S. Prather, Charles F. Gammie, the Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Bidisha Bandyopadhyay, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (241 additional authors not shown)

Abstract: 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 amp… ▽ More 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_{*}| \lesssim 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. △ Less

Submitted 1 April, 2026; v1 submitted 1 January, 2026; originally announced January 2026.

Comments: Published by ApJ. 10 figures, 3 tables

Journal ref: ApJ 1000 231 (2026)

Correlations of Simulated Black-Hole Movies Reveal Extreme-Lensing Signatures

Authors: Barbora Bezděková, Shahar Hadar, George Wong, Maciek Wielgus

Abstract: A black hole's gravitational pull can deflect light rays to an arbitrary degree. As a result, any source fluctuation near the black hole creates multiple lagged images on an observer's screen. For optically thin stochastic emission, these light echoes give rise to correlations of brightness fluctuations across time-dependent images (movies). The correlation pattern disentangles source-specific cha… ▽ More A black hole's gravitational pull can deflect light rays to an arbitrary degree. As a result, any source fluctuation near the black hole creates multiple lagged images on an observer's screen. For optically thin stochastic emission, these light echoes give rise to correlations of brightness fluctuations across time-dependent images (movies). The correlation pattern disentangles source-specific characteristics from universal features dictated by general relativity. This picture has motivated a proposal to use the two-point image correlation function as a probe of extreme gravitational lensing in upcoming black-hole imaging campaigns. In this work, we test the feasibility of this method by computing the two-point correlation function of brightness fluctuations in a black-hole movie of state-of-the-art realism. The movie is generated by ray tracing a general relativistic magnetohydrodynamic simulation, which can then be blurred to any angular resolution. At an effective resolution expected to be achieved by next-generation terrestrial very-long-baseline interferometric arrays, the lensing signatures appear in neither time-averaged images nor light-curve autocorrelations. However, we demonstrate that they are clearly visible in the more fine-grained two-point image correlation function. Our positive findings motivate a more comprehensive investigation into the instrument specifications and inference techniques needed to resolve extreme lensing effects through correlations. △ Less

Submitted 10 December, 2025; originally announced December 2025.

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

Authors: Saurabh, Hendrik Müller, Sebastiano D. von Fellenberg, Paul Tiede, Michael Janssen, Lindy Blackburn, Avery E. Broderick, Erandi Chavez, Boris Georgiev, Thomas P. Krichbaum, Kotaro Moriyama, Dhanya G. Nair, Iniyan Natarajan, Jongho Park, Andrew Thomas West, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach , et al. (260 additional authors not shown)

Abstract: We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and… ▽ More We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and $\sim2500~μ$as (0.02 pc and 0.2 pc). Without these baselines, earlier EHT observations lacked the capability to constrain emission on large scales, where a "missing flux" of order $\sim1$ Jy is expected. To probe these scales, we analyzed closure phases, robust against station-based gain errors, and modeled the jet base emission using a simple Gaussian offset from the compact ring emission at separations $>100~μ$as. Our analysis reveals a Gaussian feature centered at ($Δ$RA $\approx320~μ$as, $Δ$Dec $\approx60~μ$as), a projected separation of $\approx5500$ AU, with a flux density of only $\sim60$ mJy, implying that most of the missing flux in previous studies must arise from larger scales. Brighter emission at these scales is ruled out, and the data do not favor more complex models. This component aligns with the inferred direction of the large-scale jet and is consistent with emission from the jet base. While our findings indicate detectable jet base emission at 230 GHz, coverage from only two intermediate baselines limits reconstruction of its morphology. We therefore treat 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. △ Less

Submitted 1 December, 2025; originally announced December 2025.

Comments: 19 pages, 18 figures. Abstract shortened with respect to the manuscript. Accepted in Astronomy & Astrophysics

Journal ref: A&A, 706, A27 (2026)

Scattering of massive particles from black holes and naked singularities

Authors: Angelos Karakonstantakis, Włodek Kluźniak, Maciek Wielgus

Abstract: We performed a numerical study of the dynamics of massive particles orbiting black holes and naked singularities in the Reissner-Nordström geometry. We modeled a stream of particles with a constant angular momentum and with a range of energies. We then solved the geodesic equation of motion and compared the trajectories around black holes and naked singularities by tuning the charge parameter of t… ▽ More We performed a numerical study of the dynamics of massive particles orbiting black holes and naked singularities in the Reissner-Nordström geometry. We modeled a stream of particles with a constant angular momentum and with a range of energies. We then solved the geodesic equation of motion and compared the trajectories around black holes and naked singularities by tuning the charge parameter of the metric. The setup {allows us to explore the orbital dynamics relevant for} astrophysical scenarios such as tidal disruption events{, particularly for deep encounters}. We discussed differences and similarities in the orbital dynamics and deflection angles. We found that particles reflected by a black hole follow a stream-like family of orbits within a narrow range of deflection angles, whereas in the case of naked singularities particles are scattered in all directions on the plane of motion. We explained this behavior as an interplay between the presence of a centrifugal barrier at the location of the unstable circular orbit and an absorbing event horizon in the case of a black hole or a {repulsive core} in the case of a naked singularity. These qualitative differences are expected to impact the observable signatures of tidal disruption events. △ Less

Submitted 13 April, 2026; v1 submitted 28 September, 2025; originally announced September 2025.

Comments: Accepted to MNRAS; 13 pages, 12 figures

Revisiting 3C 279 jet morphology with space VLBI at 26 microarcsecond resolution

Authors: Teresa Toscano, José L. Gómez, Guang-Yao Zhao, Rocco Lico, Antonio Fuentes, Tuomas Savolainen, Jan Röder, Maciek Wielgus, Alexander B. Pushkarev, Efthalia Traianou, Ai-Ling Zeng, Leonid I. Gurvits, Yuri Y. Kovalev, Felix M. Pötzl, Mikhail Lisakov

Abstract: We present observations of the blazar 3C 279 at 22 GHz using the space VLBI mission RadioAstron on 2018 January 15. Images in both total intensity and fractional polarization are reconstructed using RML method implemented in the eht-imaging library. The electric vector position angles are found to be mostly aligned with the general jet direction, suggesting a predominantly toroidal magnetic field,… ▽ More We present observations of the blazar 3C 279 at 22 GHz using the space VLBI mission RadioAstron on 2018 January 15. Images in both total intensity and fractional polarization are reconstructed using RML method implemented in the eht-imaging library. The electric vector position angles are found to be mostly aligned with the general jet direction, suggesting a predominantly toroidal magnetic field, in agreement with the presence of a helical magnetic field. Ground-space fringes were detected up to a projected baseline length of $\sim 8$G$λ$, achieving the angular resolution of around 26$μ$as. The fine-scale structure of the relativistic jet is found in our study extending to a projected distance of $\sim 180$ parsec from the radio core. However, the filamentary structure reported by previous RadioAstron observations of 2014 is not detected in our current study. We discuss potential causes for this phenomenon, together with a comparison using public 43 GHz data from the BEAM-ME program, showing a significant drop in the jet's total intensity. The optically thick core is observed with a brightness temperature of $ 1.6 \times 10^{12}$ K, consistent with equipartition between the energy densities of the relativistic particles and the magnetic field. This yields an estimated magnetic field strength of 0.2 G. △ Less

Submitted 26 September, 2025; originally announced September 2025.

Journal ref: A&A 704, A225 (2025)

Semi-analytic studies of accretion disk and magnetic field geometry in M87*

Authors: Saurabh, Maciek Wielgus, Arman Tursunov, Andrei P. Lobanov, Razieh Emami

Abstract: Context: Magnetic fields play a pivotal role in dynamics of black hole accretion flows and formation of relativistic jets. Observations by the Event Horizon Telescope (EHT) provided unprecedented insights into accretion structures near black holes. Interpreting these observations requires a theoretical framework linking polarized emission to underlying system properties and magnetic field geometri… ▽ More Context: Magnetic fields play a pivotal role in dynamics of black hole accretion flows and formation of relativistic jets. Observations by the Event Horizon Telescope (EHT) provided unprecedented insights into accretion structures near black holes. Interpreting these observations requires a theoretical framework linking polarized emission to underlying system properties and magnetic field geometries. Aims: We investigate how system properties, particularly magnetic field geometry in the event horizon scale region, influence the structure of the observable synchrotron emission in M87*. Specifically, we aim to quantify the sensitivity of observables used by the EHT to black hole spin, plasma dynamics, accretion disk thickness, and magnetic field geometry. Methods: We adopt a semi-analytic radiatively inefficient accretion flow model in Kerr spacetime. We vary magnetic field geometry, black hole spin, accretion disk dynamics, and geometric thickness of the disk. We perform general relativistic ray tracing with a full polarized radiative transfer to obtain synthetic images of M87*. We extract EHT observables, such as disk diameter, asymmetry, and polarimetric metrics from synthetic models. We also consider a number of general relativistic magnetohydrodynamics simulations and compare them with the semi-analytical models. Results: We see limited impact of the disk thickness on observables. On the other hand, toroidal field dominated and poloidal field dominated magnetic configurations can be distinguished reliably. The flow dynamics, in particular presence of radial inflow, also significantly impacts the EHT observables. Conclusions: The M87* system is most consistent with a poloidal magnetic field dominated flow with partially radial inflow. While the spin remain elusive, moderate or large positive values are preferred. △ Less

Submitted 28 January, 2026; v1 submitted 15 August, 2025; originally announced August 2025.

Comments: 12 pages and 8 figures

Journal ref: A&A 705, A166 (2026)

Revealing a ribbon-like jet in OJ 287 with RadioAstron

Authors: E. Traianou, J. L. Gomez, I. Cho, A. Chael, A. Fuentes, I. Myserlis, M. Wielgus, G. -Y. Zhao, R. Lico, K. Moriyama, L. Dey, G. Bruni, R. Dahale, T. Toscano, L. I. Gurvits, M. M. Lisakov, Y. Y. Kovalev, A. P. Lobanov, A. B. Pushkarev, K. V. Sokolovsky

Abstract: We present space-based very long baseline interferometry observations of the BL Lac type object OJ 287 taken with RadioAstron at 22 GHz on April 25, 2016, in conjunction with a ground array comprising 27 radio telescopes. We detect ground-space fringes at projected baselines extending up to 4.6 Earth diameters, which allowed us to image the jet in OJ 287 with an angular resolution of ~47 μas. Appl… ▽ More We present space-based very long baseline interferometry observations of the BL Lac type object OJ 287 taken with RadioAstron at 22 GHz on April 25, 2016, in conjunction with a ground array comprising 27 radio telescopes. We detect ground-space fringes at projected baselines extending up to 4.6 Earth diameters, which allowed us to image the jet in OJ 287 with an angular resolution of ~47 μas. Applying an advanced regularized maximum likelihood imaging method, we resolved the innermost jet structure with a complex morphology at a resolution of ~15 μas (~0.1 pc projected distance). For the first time, due to a favorable geometrical position of the jet in tandem with high data quality, we detect multiple sharp bends that form a "ribbon-like" jet structure that extends down to 1 mas. Two-dimensional Gaussian model-fitting reveals regions of the jet with brightness temperatures of more than 10^13 K, indicative of strong Doppler boosting. Polarimetric imaging reveals that the electric vector position angles are predominantly perpendicular to the innermost jet direction, implying a dominant poloidal magnetic field component near the central engine. Complementary multi-epoch Very Long Baseline Array observations at 43 GHz provide a multifrequency view of the jet evolution. Ridgeline analysis of the 43 GHz data shows significant variations in the jet position angle from 2014 to 2017, behavior consistent with a rotating helical jet structure. Finally, we confirm the emergence of a new jet component (B15 or K), which may be associated with the source's first TeV flare, and offer new observational constraints relevant to models involving a supermassive black hole binary. △ Less

Submitted 3 August, 2025; originally announced August 2025.

Comments: 12 pages, 6 figures, 2 tables. Submitted to Astronomy & Astrophysics

Journal ref: A&A 700, A16 (2025)

The infrared jet of M87 observed with JWST

Authors: Jan Röder, Maciek Wielgus, Joseph B. Jensen, Gagandeep S. Anand, R. Brent Tully

Abstract: We present the first JWST+NIRCam images of the giant elliptical active galaxy M87 and its jet at 0.90, 1.50, 2.77 and 3.56 $μ$m. We analysed the large-scale jet structure, identifying prominent components, and determined the near-infrared spectral index. The data were calibrated using the standard JWST pipeline. We subtracted a constant background level and a smooth model of the galaxy surface bri… ▽ More We present the first JWST+NIRCam images of the giant elliptical active galaxy M87 and its jet at 0.90, 1.50, 2.77 and 3.56 $μ$m. We analysed the large-scale jet structure, identifying prominent components, and determined the near-infrared spectral index. The data were calibrated using the standard JWST pipeline. We subtracted a constant background level and a smooth model of the galaxy surface brightness to isolate the jet. The total image fluxes measured in the NIRCam filters follow the infrared bump pattern seen near 1.6 $μ$m in the spectral energy distribution of M87, caused by the surrounding stellar population in the galaxy. The residual jet images broadly agree with the radio to optical synchrotron power law $S_λ\proptoλ^α$ with $α=$0.7-1.0. We identified the most upstream knot L at a distance of (320$\pm$50) mas from the core. The component HST-1, at (950$\pm$50) mas from the core, is transversely resolved, and both the individual images and the spectral index map clearly indicate its double-component substructure with two elements of similar size and flux density, with centroids separated by (150$\pm$20) mas and with a significantly larger spectral index $α$ observed in the downstream component ($α_{\rm do}=0.3$) than in the upstream one ($α_{\rm up}= -0.15$). We also observe the counter-jet component located about 24 arcsec away from the nucleus. △ Less

Submitted 8 September, 2025; v1 submitted 24 July, 2025; originally announced July 2025.

Journal ref: A&A 701, L12 (2025)

Deep learning inference with the Event Horizon Telescope III. Zingularity results from the 2017 observations and predictions for future array expansions

Authors: M. Janssen, C. -k. Chan, J. Davelaar, M. Wielgus

Abstract: (abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting supermassive black hole systems. The considered models are ray-traced GRMHD simulations of Sgr A* and M87*. In this work, we infer the best-fitting accretion and… ▽ More (abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting supermassive black hole systems. The considered models are ray-traced GRMHD simulations of Sgr A* and M87*. In this work, we infer the best-fitting accretion and black hole parameters from 2017 EHT data and predict improvements that will come with future upgrades of the array. Compared to previous EHT analyses, we considered a substantially larger synthetic data library and the most complete set of information from the observational data. We made use of the Bayesian nature of the trained neural networks and apply bootstrapping of known systematics in the observational data to obtain parameter posteriors. Within a wide GRMHD parameter space, we find M87* to be best described by a spin between 0.5 and 0.94 with a retrograde MAD accretion flow and strong synchrotron emission from the jet. Sgr A* has a high spin of $\sim$ 0.8 $-$ 0.9 and a prograde accretion flow beyond the standard MAD/SANE models with a comparatively weak jet emission, seen at a $\sim$ 20$^\circ$ $-$ 40$^\circ$ inclination and $\sim$ 106$^\circ$ $-$ 137$^\circ$ position angle. While previous EHT analyses could rule out specific regions in the model parameter space considered here, we are able to obtain narrow parameter posteriors with our Zingularity framework without being impacted by the unknown foreground Faraday screens and data calibration biases. We further demonstrate that the AMT extension to the EHT will reduce parameter inference errors by a factor of three for non-Kerr models, enabling more robust tests of general relativity. It will be instructive to produce new GRMHD models with the inferred interpolated parameters to study their accretion rate plus jet power. △ Less

Submitted 16 June, 2025; originally announced June 2025.

Comments: 9 pages, 6 figures, A&A submitted 16/01/2025, accepted 01/04/2025, published 07/06/2025

Journal ref: A&A, 698, A62, 2025

Deep learning inference with the Event Horizon Telescope II. The Zingularity framework for Bayesian artificial neural networks

Authors: M. Janssen, C. -k. Chan, J. Davelaar, M. Wielgus

Abstract: (abridged) In this second paper in our publication series, we present the open-source Zingularity framework for parameter inference with deep Bayesian artificial neural networks. We carried out out supervised learning with synthetic millimeter very long baseline interferometry observations of the EHT. Our ground-truth models are based on GRMHD simulations of Sgr A* and M87* on horizon scales. We i… ▽ More (abridged) In this second paper in our publication series, we present the open-source Zingularity framework for parameter inference with deep Bayesian artificial neural networks. We carried out out supervised learning with synthetic millimeter very long baseline interferometry observations of the EHT. Our ground-truth models are based on GRMHD simulations of Sgr A* and M87* on horizon scales. We investigated how well Zingularity neural networks are able to infer key model parameters from EHT observations, such as the black hole spin and the magnetic state of the accretion disk, when uncertainties in the data are accurately taken into account. Zingularity makes use of the TensorFlow Probability library and is able to handle large amounts of data with a combination of the efficient TFRecord data format plus the Horovod framework. Our approach is the first analysis of EHT data with Bayesian neural networks, where an unprecedented training data size, under consideration of a closely modeled EHT signal path, and the full information content of the observational data are used. Zingularity infers parameters based on salient features in the data and is containerized. Through parameter surveys and dedicated validation tests, we identified neural network architectures, that are robust against internal stochastic processes and unaffected by noise in the observational and model data. We give examples of how different data properties affect the network training. We show how the Bayesian nature of our networks gives trustworthy uncertainties and uncovers failure modes for uncharacterizable data. It is easy to achieve low validation errors during training on synthetic data with neural networks, particularly when the forward modeling is too simplified. Through careful studies, we demonstrate that our trained networks can generalize well so that reliable results can be obtained from observational data. △ Less

Submitted 16 June, 2025; originally announced June 2025.

Comments: 19 pages, 6 figures, A&A submitted 16/01/2025, accepted 30/03/2025, published 07/06/2025

Journal ref: A&A, 698, A61, 2025

Deep learning inference with the Event Horizon Telescope I. Calibration improvements and a comprehensive synthetic data library

Authors: M. Janssen, C. -k. Chan, J. Davelaar, I. Natarajan, H. Olivares, B. Ripperda, J. Röder, M. Rynge, M. Wielgus

Abstract: (abridged) In a series of publications, we describe a comprehensive comparison of Event Horizon Telescope (EHT) data with theoretical models of Sgr A* and M87*. Here, we report on improvements made to our observational data reduction pipeline and present the generation of observables derived from the EHT models. We make use of ray-traced GRMHD simulations that are based on different black hole spa… ▽ More (abridged) In a series of publications, we describe a comprehensive comparison of Event Horizon Telescope (EHT) data with theoretical models of Sgr A* and M87*. Here, we report on improvements made to our observational data reduction pipeline and present the generation of observables derived from the EHT models. We make use of ray-traced GRMHD simulations that are based on different black hole spacetime metrics and accretion physics parameters. These broad classes of models provide a good representation of the primary targets observed by the EHT. To generate realistic synthetic data from our models, we took the signal path as well as the calibration process, and thereby the aforementioned improvements, into account. We could thus produce synthetic visibilities akin to calibrated EHT data and identify salient features for the discrimination of model parameters. We have produced a library consisting of an unparalleled 962,000 synthetic Sgr A* and M87* datasets. In terms of baseline coverage and noise properties, the library encompasses 2017 EHT measurements as well as future observations with an extended telescope array. We differentiate between robust visibility data products related to model features and data products that are strongly affected by data corruption effects. Parameter inference is mostly limited by intrinsic model variability, which highlights the importance of long-term monitoring observations with the EHT. In later papers in this series, we will show how a Bayesian neural network trained on our synthetic data is capable of dealing with the model variability and extracting physical parameters from EHT observations. With our calibration improvements, our newly reduced EHT datasets have a considerably better quality compared to previously analyzed data. △ Less

Submitted 16 June, 2025; originally announced June 2025.

Comments: 15 pages, 8 figures, A&A submitted 16/01/2025, accepted 08/04/2025, published 07/06/2025

Journal ref: A&A, 698, A60, 2025

Gravity versus astrophysics in black hole images and photon rings: Equatorial emissions and spherically symmetric space-times

Authors: I. Urso, F. H. Vincent, M. Wielgus, T. Paumard, G. Perrin

Abstract: The EHT collaboration released in 2019 the first horizon-scale images of a black hole accretion flow, opening a novel route for plasma physics comprehension and gravitational tests. Although the present unresolved images deeply depend on the astrophysical properties of the accreted matter, GR predicts that they contain highly lensed observables, the "photon rings", embodying the effects of strong-… ▽ More The EHT collaboration released in 2019 the first horizon-scale images of a black hole accretion flow, opening a novel route for plasma physics comprehension and gravitational tests. Although the present unresolved images deeply depend on the astrophysical properties of the accreted matter, GR predicts that they contain highly lensed observables, the "photon rings", embodying the effects of strong-field gravity. Focusing on the supermassive black hole M87* and adopting a geometrically thin, equatorial disc as a phenomenological configuration for the accreting matter, our goal is to study the degeneracy of spacetime curvature and of physically-motivated emission processes on EHT-like images observed at 230 and 345 GHz. In a parametric framework, we simulate adaptively ray-traced images using GYOTO in various spherically-symmetric spacetime geometries, for a comprehensive class of disc velocities and realistic synchrotron emission profiles. We then extract the width and the peak position of 1D intensity cross sections on the direct image and the first photon ring. We show that, among the investigated quantities, the most appropriate observables to probe the geometry are the peak positions of the first photon ring. Small geometric deviations can be unequivocally detected, regardless of the motion of the disc, ranging from a Keplerian to a radially infalling one, if the black hole mass-to-distance estimate is accurate up to around 2%; the current uncertainty of 11% being sufficient just to access extreme deviations. The equatorial set-up of this paper, favoured by present EHT observations of M87*, is adapted to model future measurements at higher observing frequencies, where absorption effects are negligible, and with higher resolution, indispensable to resolve the photon rings. Additional work is needed to investigate if our conclusions hold for more realistic disc configurations. △ Less

Submitted 29 August, 2025; v1 submitted 16 June, 2025; originally announced June 2025.

Comments: Published in A&A, 16+8 pages, 9+10 figures, 3 tables

Journal ref: A&A 700, A193 (2025)

Origin of the ring ellipticity in the black hole images of M87*

Authors: Rohan Dahale, Ilje Cho, Kotaro Moriyama, Kaj Wiik, Paul Tiede, José L. Gómez, Chi-kwan Chan, Roman Gold, Vadim Y. Bernshteyn, Marianna Foschi, Britton Jeter, Hung-Yi Pu, Boris Georgiev, Abhishek V. Joshi, Alejandro Cruz-Osorio, Iniyan Natarajan, Avery E. Broderick, León D. S. Salas, Koushik Chatterjee, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua , et al. (251 additional authors not shown)

Abstract: 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) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity usi… ▽ More 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) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of $τ= 0.08_{-0.02}^{+0.03}$ with a position angle $ξ= 50.1_{-7.6}^{+6.2}$ degrees. To interpret this measurement, we compare against General Relativistic Magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including thermal or non-thermal 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 non-thermal 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. △ Less

Submitted 15 May, 2025; originally announced May 2025.

Comments: 18 pages, 13 figures

Journal ref: A&A 699, A279 (2025)

Novel Polarimetric Analysis of Near Horizon Flaring Episodes in M87* in Millimeter Wavelength

Authors: Razieh Emami, Matthew Liska, Koushik Chatterjee, Geoffrey C. Bower, Wystan Benbow, Douglas Finkbeiner, Maciek Wielgus, Lars Hernquist, Randall Smith, Grant Tremblay, Angelo Ricarte, James F. Steiner, Avery E. Broderick, Saurabh, Jordy Davelaar, Josh Grindlay, Mark Vogelsberger, Chi-Kwan Chan

Abstract: Recent multi-wavelength observations of M87* \citep{2024A&A...692A.140A} revealed a high-energy $γ$-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 f… ▽ More Recent multi-wavelength observations of M87* \citep{2024A&A...692A.140A} revealed a high-energy $γ$-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 ($τ$) with recent SMA observations, finding that the simulated $τ$ 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 $β_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. △ Less

Submitted 14 April, 2025; v1 submitted 6 April, 2025; originally announced April 2025.

Comments: 24 pages, 23 figures

Dynamical Inference from Polarized Light Curves of Sagittarius A*

Authors: Angelo Ricarte, Nicholas S. Conroy, Maciek Wielgus, Daniel Palumbo, Razieh Emami, Chi-kwan Chan

Abstract: Polarimetric light curves of Sagittarius A* (Sgr A*) sometimes exhibit loops in the Stokes $Q$ and $U$ plane over time, often interpreted as orbiting hotspot motion. In this work, we apply the differential geometry of planar curves to develop a new technique for estimating polarimetric rotation rates. Applying this technique to 230 GHz light curves of Sgr A*, we find evidence of clockwise motion n… ▽ More Polarimetric light curves of Sagittarius A* (Sgr A*) sometimes exhibit loops in the Stokes $Q$ and $U$ plane over time, often interpreted as orbiting hotspot motion. In this work, we apply the differential geometry of planar curves to develop a new technique for estimating polarimetric rotation rates. Applying this technique to 230 GHz light curves of Sgr A*, we find evidence of clockwise motion not only during a post-flare period on 2017 April 11th, as previously discovered, but also during the quiescent days imaged by the Event Horizon Telescope (EHT). The data exhibit a clockwise fraction of $0.65 \pm 0.09$ and an overall $Q-U$ rotation rate of $-2.6 \pm 0.6 \ \mathrm{deg}\,t_g^{-1}$. We analyze a library of General Relativistic Magnetohydrodynamic (GRMHD) simulations and find that face-on, clockwise-rotating models with strong magnetic fields are most likely to be consistent with the observations. These results are consistent with EHT and GRAVITY Collaboration studies, and indirectly support an interpretation in which the polarized image of Sgr A* has been rotated by an external Faraday screen. This technique offers a novel probe of event horizon scale dynamics that complements dynamical reconstructions. △ Less

Submitted 11 July, 2025; v1 submitted 1 April, 2025; originally announced April 2025.

Comments: Accepted by ApJ. Note that some values shifted during the referee process

A helical magnetic field in quasar NRAO150 revealed by Faraday rotation

Authors: J. D. Livingston, A. S. Nikonov, S. A. Dzib, L. C. Debbrecht, Y. Y. Kovalev, M. M. Lisakov, N. R. MacDonald, G. F. Paraschos, J. Röder, M. Wielgus

Abstract: Active Galactic Nuclei (AGN) are some of the most luminous and extreme environments in the Universe. The central engines of AGN, believed to be super-massive black-holes, are fed by accretion discs threaded by magnetic fields within a dense magneto-ionic medium. We report our findings from polarimetric Very-long-baseline Interferometry (VLBI) observations of quasar NRAO150 taken in October 2022 us… ▽ More Active Galactic Nuclei (AGN) are some of the most luminous and extreme environments in the Universe. The central engines of AGN, believed to be super-massive black-holes, are fed by accretion discs threaded by magnetic fields within a dense magneto-ionic medium. We report our findings from polarimetric Very-long-baseline Interferometry (VLBI) observations of quasar NRAO150 taken in October 2022 using a combined network of the Very Long Baseline Array (VLBA) and Effelsberg 100-m Radio Telescope. These observations are the first co-temporal multi-frequency polarimetric VLBI observations of NRAO150 at frequencies above 15GHz. We use the new VLBI polarization calibration procedure, GPCAL, with polarization observations of frequencies of 12GHz, 15GHz, 24GHz, and 43GHz of NRAO150. From these observations, we measure Faraday rotation. Using our measurement of Faraday rotation, we also derive the intrinsic electric vector position angle (EVPA0) for the source. As a complementary measurement we determine the behavior of polarization as a function of observed frequency. The polarization from NRAO150 only comes from the core region, with a peak polarization intensity occurring at 24GHz. Across the core region of NRAO150 we see clear gradients in Faraday rotation and EVPA0 values that are aligned with the direction of the jet curving around the core region. We find that for the majority of the polarized region the polarization fraction is greater at higher frequencies, with intrinsic polarization fractions in the core 3%. The Faraday rotation gradients and circular patterns in EVPA0 are strong evidence for a helical/toroidal magnetic field, and the presence of low intrinsic polarization fractions indicate that the polarized emission and hence the helical/toroidal magnetic field, occur within the innermost jet. △ Less

Submitted 5 March, 2025; originally announced March 2025.

Comments: 9 pages, 9 figures

Journal ref: A&A 695, A260 (2025)

The putative center in NGC 1052

Authors: Anne-Kathrin Baczko, Matthias Kadler, Eduardo Ros, Christian M. Fromm, Maciek Wielgus, Manel Perucho, Thomas P. Krichbaum, Mislav Baloković, Lindy Blackburn, Chi-kwan Chan, Sara Issaoun, Michael Janssen, Luca Ricci, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, David Ball, Bidisha Bandyopadhyay, John Barrett , et al. (262 additional authors not shown)

Abstract: Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combine… ▽ More Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 250 RS (Schwarzschild radii) perpendicular to the jet axes. This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 300 RS and 10000 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 μas from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at about 130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 26000 Gauss, which is consistent with previous measurements. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength. △ Less

Submitted 15 January, 2025; originally announced January 2025.

Comments: 22 pages, 10 figures, published in A&A

Journal ref: A&A, 692, A205 (2024)

A multi-frequency study of sub-parsec jets with the Event Horizon Telescope

Authors: Jan Röder, Maciek Wielgus, Andrei P. Lobanov, Thomas P. Krichbaum, Dhanya G. Nair, Sang-Sung Lee, Eduardo Ros, Vincent L. Fish, Lindy Blackburn, Chi-kwan Chan, Sara Issaoun, Michael Janssen, Michael D. Johnson, Sheperd S. Doeleman, Geoffrey C. Bower, Geoffrey B. Crew, Remo P. J. Tilanus, Tuomas Savolainen, C. M. Violette Impellizzeri, Antxon Alberdi, Anne-Kathrin Baczko, José L. Gómez, Ru-Sen Lu, Georgios F. Paraschos, Efthalia Traianou , et al. (265 additional authors not shown)

Abstract: 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: the main science targets, Sgr A* and M87 along with various calibrators. We… ▽ More 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: the main science targets, Sgr A* and M87 along with various calibrators. 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-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. △ Less

Submitted 9 January, 2025; originally announced January 2025.

Journal ref: A&A 695, A233 (2025)

Pseudo-Newtonian simulation of a thin accretion disk around a Reissner-Nordström naked singularity

Authors: Miljenko Čemeljić, Włodek Kluźniak, Ruchi Mishra, Maciek Wielgus

Abstract: We present the first numerical simulations of a thin accretion disk around a Reissner-Nordström (RN) naked singularity (a charged point mass). The gravity of the RN naked singularity is modeled with a pseudo-Newtonian potential that reproduces exactly the radial dependence of the RN Keplerian orbital frequency; in particular, orbital angular velocity vanishes at the zero gravity radius and has a m… ▽ More We present the first numerical simulations of a thin accretion disk around a Reissner-Nordström (RN) naked singularity (a charged point mass). The gravity of the RN naked singularity is modeled with a pseudo-Newtonian potential that reproduces exactly the radial dependence of the RN Keplerian orbital frequency; in particular, orbital angular velocity vanishes at the zero gravity radius and has a maximum at 4/3 of that radius. Angular momentum is transported outwards by viscous stresses only outside the location of this maximum. Nonetheless, even at that radius, accretion proceeds at higher latitudes, the disk having thickened there owing to excess pressure. The accretion stops at a certain distance away from the singularity, with the material accumulating in a toroidal structure close to the zero-gravity sphere. The shape of the structure obtained in our simulations is reminiscent of fluid figures of equilibrium analytically derived in full general relativity for the RN singularity. The presence of a rotating ring, such as the one found in our simulations, could be an observational signature of a naked singularity. For charge to mass ratios close to but larger than unity, the inner edge of the quasi-toroidal inner accretion structure would be located well within the Schwarzschild marginally stable orbit (ISCO), and the maximum orbital frequency in thin accretion disks would be much higher than the Schwarzschild ISCO frequency. △ Less

Submitted 6 January, 2025; originally announced January 2025.

Comments: 7 pages, 4 figures, submitted to ApJ

Demographics of black holes at $<$100 R$_{\rm g}$ scales: accretion flows, jets, and shadows

Authors: Dhanya G. Nair, Neil M. Nagar, Venkatessh Ramakrishnan, Maciek Wielgus, Vicente Arratia, Thomas P. Krichbaum, Xinyue A. Zhang, Angelo Ricarte, Silpa S., Joaquín Hernández-Yévenes, Nicole M. Ford, Bidisha Bandyopadhyay, Mark Gurwell, Roman Burridge, Dominic W. Pesce, Sheperd S. Doeleman, Jae-Young Kim, Daewon Kim, Michael Janssen, Sebastiano D. von Fellenberg, Christian M. Fromm, Deokhyeong Lee, Heino Falcke, Jan Wagner, Geoffrey C. Bower , et al. (65 additional authors not shown)

Abstract: Using the Event Horizon Telescope (EHT), the gravitationally lensed rings around the supermassive black holes (SMBHs) in Messier 87 (M87) and Sagittarius A* (Sgr A*) have now been successfully imaged at a resolution under 10 gravitational radii (R$_{\rm g}$ $ = \rm{GM/c^2}$). To expand studies beyond M87 and Sgr A*, we have constructed the Event Horizon and Environs (ETHER) sample, a comprehensive… ▽ More Using the Event Horizon Telescope (EHT), the gravitationally lensed rings around the supermassive black holes (SMBHs) in Messier 87 (M87) and Sagittarius A* (Sgr A*) have now been successfully imaged at a resolution under 10 gravitational radii (R$_{\rm g}$ $ = \rm{GM/c^2}$). To expand studies beyond M87 and Sgr A*, we have constructed the Event Horizon and Environs (ETHER) sample, a comprehensive database encompassing approximately 3.15 million SMBH mass estimates, $\sim$ 20,000 Very-Long Baseline Interferometry (VLBI) radio flux densities, and $\sim$ 36,000 hard X-ray flux densities. This database is designed to identify and optimize target selection for the EHT and its upgrades on the ground and in space. We have identified a Gold Sample (GS) of nearby low-luminosity Active Galactic Nuclei (AGNs) within it that are ideal for studying jet bases and potentially imaging black hole shadows. We observed 27 of these AGNs using the EHT from 2022 to 2024, providing an opportunity to resolve and image accretion flows and jets at resolutions of $\leq$ 100 R$_{\rm g}$. Only a few SMBHs have sufficiently high enough flux density to be imaged at scales of $\leq$ 50 R$_{\rm g}$ with the present EHT. Among these are M87, Sgr A*, NGC4594 (Sombrero/M104), NGC4261, and NGC4374 (Messier 84/M84). Of these, NGC4261, Sombrero, and M84 have been observed and/or are scheduled for deep imaging with EHT+ALMA from 2023 to 2025. Sombrero, NGC4261, M84, NGC4278, and NGC5232 are clearly detected in our EHT+ALMA observations in 2022, indicating that the 230 GHz flux density from the accretion flows is significantly high. Ongoing imaging of the ETHER GS will enable measurements of black hole mass and spin, help constrain General Relativity, and enrich our understanding of jet launching and accretion inflows across a broad multi-parameter space, including black hole mass, spin, accretion rate, and orientation. △ Less

Submitted 28 December, 2024; originally announced December 2024.

Comments: 9 pages, 6 figures, 1 table, published in Proceedings of the 16th EVN Symposium, Ed. E. Ros, P. Benke, S.A. Dzib, I. Rottmann, & J.A. Zensus, Bonn: Max-Planck-Institut für Radioastronomie, 2024, pages 75-84, https://cloud.mpifr-bonn.mpg.de/index.php/s/BkX2CC2Xjn2aKR4

Helical magnetic field structure in 3C 273. A Faraday rotation analysis using multi-frequency polarimetric VLBA data

Authors: Teresa Toscano, Sol N. Molina, José L. Gómez, Ai-Ling Zeng, Rohan Dahale, Ilje Cho, Kotaro Moriyama, Maciek Wielgus, Antonio Fuentes, Marianna Foschi, Efthalia Traianou, Jan Röder, Ioannis Myserlis, Emmanouil Angelakis, Anton Zensus

Abstract: We present a study on rotation measure (RM) of the quasar 3C 273. This analysis aims to discern the magnetic field structure and its temporal evolution. The quasar 3C 273 is one of the most studied active galactic nuclei due to its high brightness, strong polarization, and proximity, which enables resolving the transverse structure of its jet in detail. We used polarized data from 2014, collected… ▽ More We present a study on rotation measure (RM) of the quasar 3C 273. This analysis aims to discern the magnetic field structure and its temporal evolution. The quasar 3C 273 is one of the most studied active galactic nuclei due to its high brightness, strong polarization, and proximity, which enables resolving the transverse structure of its jet in detail. We used polarized data from 2014, collected at six frequencies (5, 8, 15, 22, 43, 86 GHz) with the Very Long Baseline Array, to produce total and linear polarization intensity images, as well as RM maps. Our analysis reveals a well-defined transverse RM gradient across the jet, indicating a helical, ordered magnetic field that threads the jet and likely contributes to its collimation. Furthermore, we identified temporal variations in the RM magnitude when compared with prior observations. These temporal variations show that the environment around the jet is dynamic, with changes in the density and magnetic field strength of the sheath that are possibly caused by interactions with the surrounding medium. △ Less

Submitted 24 December, 2024; originally announced December 2024.

Journal ref: A&A 698, A210 (2025)

Kilogauss magnetic field and jet dynamics in the quasar NRAO 530

Authors: Mikhail Lisakov, Svetlana Jorstad, Maciek Wielgus, Evgeniya V. Kravchenko, Aleksei S. Nikonov, Ilje Cho, Sara Issaoun, Juan-Carlos Algaba, Thomas P. Krichbaum, Uwe Bach, Eduardo Ros, Helge Rottmann, Salvador S'anchez, Jan Wagner, Anton Zensus

Abstract: The advancement of the Event Horizon Telescope has enabled the study of relativistic jets in active galactic nuclei down to sub-parsec linear scales even at high redshift. Quasi-simultaneous multifrequency observations provide insights into the physical conditions in compact regions and allow testing accretion theories. Initially we aimed at measuring the magnetic field strength close to the centr… ▽ More The advancement of the Event Horizon Telescope has enabled the study of relativistic jets in active galactic nuclei down to sub-parsec linear scales even at high redshift. Quasi-simultaneous multifrequency observations provide insights into the physical conditions in compact regions and allow testing accretion theories. Initially we aimed at measuring the magnetic field strength close to the central supermassive black hole in NRAO 530 (1730-130) by studying frequency-dependent opacity of the jet matter, Faraday rotation and the spectral index in the mm-radio bands. NRAO 530 was observed quasi-simultaneously at 15, 22, 43, 86, and 227 GHz at four different very long baseline interferometer (VLBI) networks. By the means of imaging and model-fitting, we aligned the images, taken at different frequencies. We explored opacity along the jet and distribution of the linearly polarized emission in it. Our findings reveal that the jet of NRAO 530 at 86 and 227 GHz is transparent down to its origin, with 70 mJy emission detected at 227 GHz potentially originating from the accretion disk. The magnetic field strength near the black hole, estimated at $5r_\mathrm{g}$, is $3\times10^3-3\times10^4$ G (depending on the central black hole mass). These values represent some of the highest magnetic field strengths reported for active galaxies. We also report the first ever VLBI measurement of the Faraday rotation at 43-227 GHz, which reveals rotation measure values as high as -48000 rad/m2 consistent with higher particle density and stronger magnetic fields at the jet's outset. The complex shape of the jet in NRAO 530 is in line with the expected behavior of a precessing jet, with a period estimated to be around $6\pm4$~years. △ Less

Submitted 5 November, 2024; originally announced November 2024.

Comments: 18 pages, accepted to A&A

Journal ref: A&A 693, A9 (2025)

First Very Long Baseline Interferometry Detections at 870μm

Authors: Alexander W. Raymond, Sheperd S. Doeleman, Keiichi Asada, Lindy Blackburn, Geoffrey C. Bower, Michael Bremer, Dominique Broguiere, Ming-Tang Chen, Geoffrey B. Crew, Sven Dornbusch, Vincent L. Fish, Roberto García, Olivier Gentaz, Ciriaco Goddi, Chih-Chiang Han, Michael H. Hecht, Yau-De Huang, Michael Janssen, Garrett K. Keating, Jun Yi Koay, Thomas P. Krichbaum, Wen-Ping Lo, Satoki Matsushita, Lynn D. Matthews, James M. Moran , et al. (254 additional authors not shown)

Abstract: The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescop… ▽ More The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$λ$ corresponding to an angular resolution, or fringe spacing, of 19$μ$as. The Allan deviation of the visibility phase at 870$μ$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$μ$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Comments: Corresponding author: S. Doeleman

Journal ref: The Astronomical Journal, Volume 168, Issue 3, id.130, 19 pp. 2024

Mid-Range Science Objectives for the Event Horizon Telescope

Authors: The Event Horizon Telescope Collaboration

Abstract: The first images of the black holes in Sagittarius A* and M87* have created a wide range of new scientific opportunities in gravitational physics, compact objects, and relativistic astrophysics. We discuss here the scientific opportunities that arise from the rich data sets that have already been obtained and the new data sets that will be obtained, exploiting a wide range of technical advances, i… ▽ More The first images of the black holes in Sagittarius A* and M87* have created a wide range of new scientific opportunities in gravitational physics, compact objects, and relativistic astrophysics. We discuss here the scientific opportunities that arise from the rich data sets that have already been obtained and the new data sets that will be obtained, exploiting a wide range of technical advances, including observational agility, receiver upgrades, and the addition of new stations. This document provides a 5-year framework for Event Horizon Telescope (EHT) science structured around four fundamental questions that are used to prioritize the analysis of existing data, guide technical upgrades, and determine the optimal use of future observational opportunities with EHT, ALMA, and multi-wavelength facilities. Through enhancements over this period, the EHT will create the first movie of M87* connecting black hole and jet physics, provide detailed studies of the structure and dynamics of Sgr A*, characterize the magnetospheres of both systems through polarimetric imaging, and explore the spacetime properties of black holes with greater precision and range. △ Less

Submitted 3 October, 2024; originally announced October 2024.

Comments: 22 pages, 4 figures

Selective Dynamical Imaging of Interferometric Data

Authors: Joseph Farah, Peter Galison, Kazunori Akiyama, Katherine L. Bouman, Geoffrey C. Bower, Andrew Chael, Antonio Fuentes, José L. Gómez, Mareki Honma, Michael D. Johnson, Yutaro Kofuji, Daniel P. Marrone, Kotaro Moriyama, Ramesh Narayan, Dominic W. Pesce, Paul Tiede, Maciek Wielgus, Guang-Yao Zhao, The Event Horizon Telescope Collaboration

Abstract: Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT's $(u, v)$-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the… ▽ More Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT's $(u, v)$-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the EHT can contain regions of time over the course of a single observation that facilitate dynamical imaging. These optimal time regions typically have projected baseline distributions that are approximately angularly isotropic and radially homogeneous. We derive a metric of coverage quality based on baseline isotropy and density that is capable of ranking array configurations by their ability to produce accurate dynamical reconstructions. We compare this metric to existing metrics in the literature and investigate their utility by performing dynamical reconstructions on synthetic data from simulated EHT observations of sources with simple orbital variability. We then use these results to make recommendations for imaging the 2017 EHT Sgr A* data set. △ Less

Submitted 12 September, 2024; originally announced September 2024.

Comments: 20 pages, 14 figures, published in ApJL (2022), uploaded for open access

Journal ref: ApJL, Volume 930, L18, 2022

Imaging a ring-like structure and the extended jet of M87 at 86 GHz

Authors: Jong-Seo Kim, Hendrik Mueller, Aleksei S. Nikonov, Ru-Sen Lu, Jakob Knollmueller, Torsten A. Ensslin, Maciek Wielgus, Andrei P. Lobanov

Abstract: The galaxy M87 is one of the prime targets for high resolution radio imaging to investigate the supermassive black hole, accretion flow, and relativistic jet. However, it remains challenging to observe them jointly. In 2018, GMVA+ALMA observations at 86 GHz enabled the simultaneous reconstruction of a ring structure and the extended jet emission. In order to analyze the ring and jet of M87, conven… ▽ More The galaxy M87 is one of the prime targets for high resolution radio imaging to investigate the supermassive black hole, accretion flow, and relativistic jet. However, it remains challenging to observe them jointly. In 2018, GMVA+ALMA observations at 86 GHz enabled the simultaneous reconstruction of a ring structure and the extended jet emission. In order to analyze the ring and jet of M87, conventional 'CLEAN' algorithms were mainly employed alongside the RML method 'SMILI' in previous work. To test the robustness of the reconstructed structures of M87 GMVA+ALMA observations at 86 GHz, we estimate the ring diameter, width, and the extended jet emission with the possible central spine by two different novel imaging algorithms: 'resolve' and 'DoG-HiT'. Overall, reconstructions are consistent with the 'CLEAN' and 'SMILI' images. The ring structure of M87 is resolved at a higher resolution and the posterior distribution of M87 ring features is explored. The 'resolve' images show that the ring diameter is 60.9 +- 2.2 muas and its width is 16.0 +- 0.9 muas. The ring diameter and the ring width measured from the 'DoG-HiT' image are 61.0 muas and 20.6 muas, respectively. The ring diameter is therefore in agreement with the estimation in the previous work (64 +4 -8 muas). Two bright spots in the ring are reconstructed by four independent imaging methods. Therefore, the substructure in the ring most likely results from the data. A consistent limb-brightened jet structure is reconstructed by 'resolve' and 'DoG-HiT', albeit with a less pronounced central spine. Modern data-driven imaging methods confirm the ring and jet structure in M87, and complement traditional VLBI methods with novel perspectives on evaluating the significance of the recovered features. They confirm the result of the previous report. △ Less

Submitted 29 May, 2025; v1 submitted 31 August, 2024; originally announced September 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 696, A169 (2025)

Hot-spots around Sagittarius A*: Joint fits to astrometry and polarimetry

Authors: A. I. Yfantis, M. Wielgus, M. A. Mościbrodzka

Abstract: Observations of Sagittarius A* (Sgr A*) in near-infrared (NIR) show irregular flaring activity. Flares coincide with astrometric rotation of brightness centroid and with looping patterns in fractional linear polarization. These signatures can be explained with a model of a bright hot-spot, transiently orbiting the black hole. We extend the capabilities of the existing algorithms to perform paramet… ▽ More Observations of Sagittarius A* (Sgr A*) in near-infrared (NIR) show irregular flaring activity. Flares coincide with astrometric rotation of brightness centroid and with looping patterns in fractional linear polarization. These signatures can be explained with a model of a bright hot-spot, transiently orbiting the black hole. We extend the capabilities of the existing algorithms to perform parameter estimation and model comparison in the Bayesian framework using NIR observations from the GRAVITY instrument, and simultaneously fitting to the astrometric and polarimetric data. Using the numerical radiative transfer code ipole, we defined several geometric models describing a hot-spot orbiting Sgr A*, threaded with magnetic field, and emitting synchrotron radiation. We then explored the posterior space of our models in the Bayesian framework with a nested sampling code dynesty. We have used 11 models to sharpen our understanding of the importance of various aspects of the orbital model, such as non-Keplerian motion, hot-spot size, and off-equatorial orbit. All considered models converge to realizations that fit the data well, but the equatorial super-Keplerian model is favored by the currently available NIR dataset. We have inferred an inclination of ~ 155 deg, which corroborates previous estimates, a preferred period of ~ 70 minutes, and an orbital radius of ~ 12 gravitational radii with the orbital velocity of ~ 1.3 times the Keplerian value. A hot-spot of a diameter smaller than 5 gravitational radii is favored. Black hole spin is not constrained well. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Journal ref: A&A 691, A327 (2024)

Circular Polarization of Simulated Images of Black Holes

Authors: Abhishek V. Joshi, Ben S. Prather, Chi-kwan Chan, Maciek Wielgus, Charles F. Gammie

Abstract: Models of the resolved Event Horizon Telescope (EHT) sources Sgr A* and M87* are constrained by observations at multiple wavelengths, resolutions, polarizations, and time cadences. In this paper we compare unresolved circular polarization (CP) measurements to a library of models, where each model is characterized by a distribution of CP over time. In the library we vary the spin of the black hole,… ▽ More Models of the resolved Event Horizon Telescope (EHT) sources Sgr A* and M87* are constrained by observations at multiple wavelengths, resolutions, polarizations, and time cadences. In this paper we compare unresolved circular polarization (CP) measurements to a library of models, where each model is characterized by a distribution of CP over time. In the library we vary the spin of the black hole, the magnetic field strength at the horizon (i.e. both SANE and MAD models), the observer inclination, a parameter for the maximum ion-electron temperature ratio assuming a thermal plasma, and the direction of the magnetic field dipole moment. We find that ALMA observations of Sgr A* are inconsistent with all edge-on ($i = 90^\circ$) models. Restricting attention to the magnetically arrested disk (MAD) models favored by earlier EHT studies of Sgr A*, we find that only models with magnetic dipole moment pointing away from the observer are consistent with ALMA data. We also note that in 26 of the 27 passing MAD models the accretion flow rotates clockwise on the sky. We provide a table of the mean and standard deviation of the CP distributions for all model parameters along with their trends. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 33 pages, 17 figures, 2 tables. Accepted for publication in ApJ

The Black Hole Explorer: Motivation and Vision

Authors: Michael D. Johnson, Kazunori Akiyama, Rebecca Baturin, Bryan Bilyeu, Lindy Blackburn, Don Boroson, Alejandro Cardenas-Avendano, Andrew Chael, Chi-kwan Chan, Dominic Chang, Peter Cheimets, Cathy Chou, Sheperd S. Doeleman, Joseph Farah, Peter Galison, Ronald Gamble, Charles F. Gammie, Zachary Gelles, Jose L. Gomez, Samuel E. Gralla, Paul Grimes, Leonid I. Gurvits, Shahar Hadar, Kari Haworth, Kazuhiro Hada , et al. (43 additional authors not shown)

Abstract: We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery… ▽ More We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Proceedings for SPIE Astronomical Telescopes and Instrumentation

Journal ref: Proceedings Volume 13092, Space Telescopes and Instrumentation 2024: Optical, Infrared, and Millimeter Wave; 130922D (2024)

First VLBI detection of Fornax A

Authors: G. F. Paraschos, M. Wielgus, P. Benke, V. Mpisketzis, F. Rösch, K. Dasyra, E. Ros, M. Kadler, R. Ojha, P. G. Edwards, L. Hyland, J. F. H. Quick, S. Weston

Abstract: Radio galaxies harbouring jetted active galactic nuclei are a frequent target of very-long-baseline interferometry (VLBI) because they play an essential role in exploring how jets form and propagate. Hence, only few have not been detected with VLBI yet; Fornax A is one of the most famous examples. Here we present the first detection of the compact core region of Fornax A with VLBI. At 8.4 GHz the… ▽ More Radio galaxies harbouring jetted active galactic nuclei are a frequent target of very-long-baseline interferometry (VLBI) because they play an essential role in exploring how jets form and propagate. Hence, only few have not been detected with VLBI yet; Fornax A is one of the most famous examples. Here we present the first detection of the compact core region of Fornax A with VLBI. At 8.4 GHz the faint core is consistent with an unresolved point source. We constrained its flux density to be $S_0 = 47.5-62.3\,\textrm{mJy}$ and its diameter to be $D^\textrm{min}_0 \leq 70\,μ\textrm{as}$. The high values of the measured brightness temperature ($T_\textrm{B} \gtrsim 10^{11}\,\textrm{K}$) imply that the observed radiation is of non-thermal origin, likely associated with the synchrotron emission from the active galactic nucleus. We also investigated the possibility of a second radio source being present within the field of view. Adding a second Gaussian component to the geometrical model-fit does not significantly improve the quality of the fit and we, therefore, conclude that our detection corresponds to the compact core of Fornax A. Analysis of the non-trivial closure phases provides evidence for the detection of more extended flux density, on the angular scale of $\sim4000\,μ\textrm{as}$. Finally, the fractional circular polarisation of the core is consistent with zero, with a conservative upper limit being $m_\textrm{circ} \leq 4\%$. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 4 pages, 2 figures, accepted for publication in A&A

Journal ref: A&A 687, L6 (2024)

Broadband Multi-wavelength Properties of M87 during the 2018 EHT Campaign including a Very High Energy Flaring Episode

Authors: J. C. Algaba, M. Balokovic, S. Chandra, W. Y. Cheong, Y. Z. Cui, F. D'Ammando, A. D. Falcone, N. M. Ford, M. Giroletti, C. Goddi, M. A. Gurwell, K. Hada, D. Haggard, S. Jorstad, A. Kaur, T. Kawashima, S. Kerby, J. Y. Kim, M. Kino, E. V. Kravchenko, S. S. Lee, R. S. Lu, S. Markoff, J. Michail, J. Neilsen , et al. (721 additional authors not shown)

Abstract: The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physi… ▽ More The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling. △ Less

Submitted 5 December, 2024; v1 submitted 24 April, 2024; originally announced April 2024.

Comments: 46 pages, 23 figures, accepted by Astronomy & Astrophysics on August. 29, 2024

Journal ref: A&A 692, A140 (2024)

Hotspots and Photon Rings in Spherically-Symmetric Spacetimes

Authors: Prashant Kocherlakota, Luciano Rezzolla, Rittick Roy, Maciek Wielgus

Abstract: Future black hole (BH) imaging observations are expected to resolve finer features corresponding to higher-order images of hotspots and of the horizon-scale accretion flow. In spherical spacetimes, the image order is determined by the number of half-loops executed by the photons that form it. Consecutive-order images arrive approximately after a delay time of $\approxπ$ times the BH shadow radius.… ▽ More Future black hole (BH) imaging observations are expected to resolve finer features corresponding to higher-order images of hotspots and of the horizon-scale accretion flow. In spherical spacetimes, the image order is determined by the number of half-loops executed by the photons that form it. Consecutive-order images arrive approximately after a delay time of $\approxπ$ times the BH shadow radius. The fractional diameters, widths, and flux-densities of consecutive-order images are exponentially demagnified by the lensing Lyapunov exponent, a characteristic of the spacetime. The appearance of a simple point-sized hotspot when located at fixed spatial locations or in motion on circular orbits is investigated. The exact time delay between the appearance of its zeroth and first-order images agrees with our analytic estimate, which accounts for the observer inclination, with $\lesssim 20\%$ error for hotspots located about $\lesssim 5M$ from a Schwarzschild BH of mass $M$. Since M87$^\star$ and Sgr A$^\star$ host geometrically-thick accretion flows, we also explore the variation in the diameters and widths of their first-order images with disk scale-height. Using a simple conical torus model, for realistic morphologies, we estimate the first-order image diameter to deviate from that of the shadow by $\lesssim 30\%$ and its width to be $\lesssim 1.3M$. Finally, the error in recovering the Schwarzschild lensing exponent ($π$), when using the diameters or the widths of the first and second-order images is estimated to be $\lesssim 20\%$. It will soon become possible to robustly learn more about the spacetime geometry of astrophysical BHs from such measurements. △ Less

Submitted 24 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: Discussion on degeneracies added, matches accepted version. For implications for gravity tests, see part II of this series at arXiv:2307.16841

Particle motion around luminous neutron stars: effects of deviation from Schwarzschild spacetime

Authors: Ronaldo S. S. Vieira, Maciek Wielgus

Abstract: We study trajectories of test particles around a luminous, static, spherically symmetric neutron star, under the combined influence of gravity and radiation. In general relativity, for Schwarzschild spacetime, an equilibrium sphere (the Eddington Capture Sphere) is formed for near-Eddington luminosities. We generalize these results to a broad class of static, spherical spacetimes. We also study th… ▽ More We study trajectories of test particles around a luminous, static, spherically symmetric neutron star, under the combined influence of gravity and radiation. In general relativity, for Schwarzschild spacetime, an equilibrium sphere (the Eddington Capture Sphere) is formed for near-Eddington luminosities. We generalize these results to a broad class of static, spherical spacetimes. We also study the dynamics of particles in a strong radiation field in spherical spacetimes. The results are illustrated for two cases, Reissner-Nordström spacetime of a charged spherical object in general relativity and Kehagias-Sfetsos spacetime, arising from the Horava-Lifshitz gravity theory. Our findings apply to neutron stars under gravitational field equations different from the vacuum Einstein field equations of general relativity, such as in modified theories of gravity, the only requirement being that test particles follow geodesics in the absence of the radiation field. The effects that we describe are, in principle, measurable through observations of X-ray bursts of neutron stars. Hence, detailed future studies could use such observations to test gravity theories in the strong-field regime, provided that the impact of the spacetime geometry can be disentangled from the astrophysical uncertainties. △ Less

Submitted 5 September, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 8 pages, 4 figures

Journal ref: Physical Review D 110, 064010 (2024)

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Energy dissipation in astrophysical simulations: results of the Orszag-Tang test problem

Authors: Fatemeh Kayanikhoo, Miljenko Cemeljic, Maciek Wielgus, Wlodek Kluzniak

Abstract: The magnetic field through the magnetic reconnection process affects the dynamics and structure of astrophysical systems. Numerical simulations are the tools to study the evolution of these systems. However, the resolution, dimensions, resistivity, and turbulence of the system are some important parameters to take into account in the simulations. In this paper, we investigate the evolution of magn… ▽ More The magnetic field through the magnetic reconnection process affects the dynamics and structure of astrophysical systems. Numerical simulations are the tools to study the evolution of these systems. However, the resolution, dimensions, resistivity, and turbulence of the system are some important parameters to take into account in the simulations. In this paper, we investigate the evolution of magnetic energy in astrophysical simulations by performing a standard test problem for MHD codes, Orszag-Tang. We estimate the numerical dissipation in the simulations using state-of-the-art numerical simulation code in astrophysics, PLUTO. The estimated numerical resistivity in 2D simulations corresponds to the Lundquist number $\approx 10^{4}$ in the resolution of $512\times512$ grid cells. It is also shown that the plasmoid unstable reconnection layer can be resolved with sufficient resolutions. Our analysis demonstrates that in non-relativistic magnetohydrodynamics simulations, magnetic and kinetic energies undergo conversion into internal energy, resulting in plasma heating. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 8 pages, 3 figures, Accepted in Proceedings of RAGtime 23-25

MSC Class: 85-08 ACM Class: J.2

Fundamental Physics Opportunities with the Next-Generation Event Horizon Telescope

Authors: Dimitry Ayzenberg, Lindy Blackburn, Richard Brito, Silke Britzen, Avery E. Broderick, Raúl Carballo-Rubio, Vitor Cardoso, Andrew Chael, Koushik Chatterjee, Yifan Chen, Pedro V. P. Cunha, Hooman Davoudiasl, Peter B. Denton, Sheperd S. Doeleman, Astrid Eichhorn, Marshall Eubanks, Yun Fang, Arianna Foschi, Christian M. Fromm, Peter Galison, Sushant G. Ghosh, Roman Gold, Leonid I. Gurvits, Shahar Hadar, Aaron Held , et al. (23 additional authors not shown)

Abstract: 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 supermass… ▽ More 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. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: To be submitted to journal. Comments are welcome

First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization

Authors: The Event Horizon Telescope Collaboration

Abstract: Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here, we present the first analysis of circular polarization using EHT da… ▽ More Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here, we present the first analysis of circular polarization using EHT data, acquired in 2017, which can potentially provide additional insights into the magnetic fields and plasma composition near the black hole. Interferometric closure quantities provide convincing evidence for the presence of circularly polarized emission on event-horizon scales. We produce images of the circular polarization using both traditional and newly developed methods. All methods find a moderate level of resolved circular polarization across the image ($\langle|v|\rangle < 3.7\%$), consistent with the low image-integrated circular polarization fraction measured by the ALMA array ($|v_{\rm int}| < 1\%$). Despite this broad agreement, the methods show substantial variation in the morphology of the circularly polarized emission, indicating that our conclusions are strongly dependent upon the imaging assumptions because of the limited baseline coverage, uncertain telescope gain calibration, and weakly polarized signal. We include this upper limit in an updated comparison to general relativistic magnetohydrodynamic (GRMHD) simulation models. This analysis reinforces the previously reported preference for magnetically arrested accretion flow models. We find that most simulations naturally produce a low level of circular polarization consistent with our upper limit, and that Faraday conversion is likely the dominant production mechanism for circular polarization at 230 GHz in M87*. △ Less

Submitted 18 November, 2023; originally announced November 2023.

Comments: 47 pages, 31 figures, published in The Astrophysical Journal Letters on Nov. 8, 2023. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 957, Number 2, id L20, 42pp., 2023 Nov. 8

Polarimetric Geometric Modeling for mm-VLBI Observations of Black Holes

Authors: Freek Roelofs, Michael D. Johnson, Andrew Chael, Michael Janssen, Maciek Wielgus, Avery E. Broderick, The Event Horizon Telescope Collaboration

Abstract: The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, foc… ▽ More The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, focusing on approaches that fit data products with differing degrees of invariance to broad classes of calibration errors. We establish a fitting procedure using a polarimetric "m-ring" model to approximate the image structure near a black hole. By fitting this model to synthetic EHT data from general relativistic magnetohydrodynamic models, we show that the linear and circular polarization structure can be successfully approximated with relatively few model parameters. We then fit this model to EHT observations of M87* taken in 2017. In total intensity and linear polarization, the m-ring fits are consistent with previous results from imaging methods. In circular polarization, the m-ring fits indicate the presence of event-horizon-scale circular polarization structure, with a persistent dipolar asymmetry and orientation across several days. The same structure was recovered independently of observing band, used data products, and model assumptions. Despite this broad agreement, imaging methods do not produce similarly consistent results. Our circular polarization results, which imposed additional assumptions on the source structure, should thus be interpreted with some caution. Polarimetric geometric modeling provides a useful and powerful method to constrain the properties of horizon-scale polarized emission, particularly for sparse arrays like the EHT. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Comments: 34 pages, 15 figures, published in The Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 2023, Volume 957, Issue 2, id.L21, 29 pp

First Sagittarius A* Event Horizon Telescope Results. VI: Testing the Black Hole Metric

Authors: The Event Horizon Telescope Collaboration

Abstract: Astrophysical black holes are expected to be described by the Kerr metric. This is the only stationary, vacuum, axisymmetric metric, without electromagnetic charge, that satisfies Einstein's equations and does not have pathologies outside of the event horizon. We present new constraints on potential deviations from the Kerr prediction based on 2017 EHT observations of Sagittarius A* (Sgr A*). We c… ▽ More Astrophysical black holes are expected to be described by the Kerr metric. This is the only stationary, vacuum, axisymmetric metric, without electromagnetic charge, that satisfies Einstein's equations and does not have pathologies outside of the event horizon. We present new constraints on potential deviations from the Kerr prediction based on 2017 EHT observations of Sagittarius A* (Sgr A*). We calibrate the relationship between the geometrically defined black hole shadow and the observed size of the ring-like images using a library that includes both Kerr and non-Kerr simulations. We use the exquisite prior constraints on the mass-to-distance ratio for Sgr A* to show that the observed image size is within $\sim$ 10$\%$ of the Kerr predictions. We use these bounds to constrain metrics that are parametrically different from Kerr as well as the charges of several known spacetimes. To consider alternatives to the presence of an event horizon we explore the possibility that Sgr A* is a compact object with a surface that either absorbs and thermally re-emits incident radiation or partially reflects it. Using the observed image size and the broadband spectrum of Sgr A*, we conclude that a thermal surface can be ruled out and a fully reflective one is unlikely. We compare our results to the broader landscape of gravitational tests. Together with the bounds found for stellar mass black holes and the M87 black hole, our observations provide further support that the external spacetimes of all black holes are described by the Kerr metric, independent of their mass. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 55 pages, 27 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L17, 44pp., 2022 May 12

First Sagittarius A* Event Horizon Telescope Results. III: Imaging of the Galactic Center Supermassive Black Hole

Authors: The Event Horizon Telescope Collaboration

Abstract: We present the first event-horizon-scale images and spatiotemporal analysis of Sgr A* taken with the Event Horizon Telescope in 2017 April at a wavelength of 1.3 mm. Imaging of Sgr A* has been conducted through surveys over a wide range of imaging assumptions using the classical CLEAN algorithm, regularized maximum likelihood methods, and a Bayesian posterior sampling method. Different prescriptio… ▽ More We present the first event-horizon-scale images and spatiotemporal analysis of Sgr A* taken with the Event Horizon Telescope in 2017 April at a wavelength of 1.3 mm. Imaging of Sgr A* has been conducted through surveys over a wide range of imaging assumptions using the classical CLEAN algorithm, regularized maximum likelihood methods, and a Bayesian posterior sampling method. Different prescriptions have been used to account for scattering effects by the interstellar medium towards the Galactic Center. Mitigation of the rapid intra-day variability that characterizes Sgr A* has been carried out through the addition of a "variability noise budget" in the observed visibilities, facilitating the reconstruction of static full-track images. Our static reconstructions of Sgr A* can be clustered into four representative morphologies that correspond to ring images with three different azimuthal brightness distributions, and a small cluster that contains diverse non-ring morphologies. Based on our extensive analysis of the effects of sparse $(u,v)$-coverage, source variability and interstellar scattering, as well as studies of simulated visibility data, we conclude that the Event Horizon Telescope Sgr A* data show compelling evidence for an image that is dominated by a bright ring of emission with a ring diameter of $\sim$ 50 $μ$as, consistent with the expected "shadow" of a $4\times10^6 M_\odot$ black hole in the Galactic Center located at a distance of 8 kpc. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 81 pages, 45 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L14, 64pp., 2022 May 12

First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole

Authors: The Event Horizon Telescope Collaboration

Abstract: In this paper, we provide a first physical interpretation for the Event Horizon Telescope (EHT)'s 2017 observations of Sgr A*. Our main approach is to compare resolved EHT data at 230 GHz and unresolved non-EHT observations from radio to X-ray wavelengths to predictions from a library of models based on time-dependent general relativistic magnetohydrodynamics (GRMHD) simulations, including aligned… ▽ More In this paper, we provide a first physical interpretation for the Event Horizon Telescope (EHT)'s 2017 observations of Sgr A*. Our main approach is to compare resolved EHT data at 230 GHz and unresolved non-EHT observations from radio to X-ray wavelengths to predictions from a library of models based on time-dependent general relativistic magnetohydrodynamics (GRMHD) simulations, including aligned, tilted, and stellar wind-fed simulations; radiative transfer is performed assuming both thermal and non-thermal electron distribution functions. We test the models against 11 constraints drawn from EHT 230 GHz data and observations at 86 GHz, 2.2 $μ$m, and in the X-ray. All models fail at least one constraint. Light curve variability provides a particularly severe constraint, failing nearly all strongly magnetized (MAD) models and a large fraction of weakly magnetized (SANE) models. A number of models fail only the variability constraints. We identify a promising cluster of these models, which are MAD and have inclination $i \le$ 30$^\circ$. They have accretion rate $(5.2$-$9.5)\times10^{-9}M_\odot$yr$^{-1}$, bolometric luminosity $(6.8$--$9.2)\times10^{35}$ erg s$^{-1}$, and outflow power $(1.3$--$4.8)\times10^{38}$ erg s$^{-1}$. We also find that: all models with $i \ge$ 70$^\circ$ fail at least two constraints, as do all models with equal ion and electron temperature; exploratory, non-thermal model sets tend to have higher 2.2 $μ$m flux density; the population of cold electrons is limited by X-ray constraints due to the risk of bremsstrahlung overproduction. Finally we discuss physical and numerical limitations of the models, highlighting the possible importance of kinetic effects and duration of the simulations. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 62 pages, 35 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L16, 49pp., 2022 May 12

First Sagittarius A* Event Horizon Telescope Results. IV. Variability, Morphology, and Black Hole Mass

Authors: The Event Horizon Telescope Collaboration

Abstract: In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fra… ▽ More In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fraction of the correlated flux density, reaching $\sim$100\% on some baselines. Through an exploration of simple geometric source models, we demonstrate that ring-like morphologies provide better fits to the Sgr A* data than do other morphologies with comparable complexity. We develop two strategies for fitting static geometric ring models to the time-variable Sgr A* data; one strategy fits models to short segments of data over which the source is static and averages these independent fits, while the other fits models to the full dataset using a parametric model for the structural variability power spectrum around the average source structure. Both geometric modeling and image-domain feature extraction techniques determine the ring diameter to be $51.8 \pm 2.3$ $μ$as (68\% credible intervals), with the ring thickness constrained to have an FWHM between $\sim$30\% and 50\% of the ring diameter. To bring the diameter measurements to a common physical scale, we calibrate them using synthetic data generated from GRMHD simulations. This calibration constrains the angular size of the gravitational radius to be $4.8_{-0.7}^{+1.4}$ \mathrm{μas}, which we combine with an independent distance measurement from maser parallaxes to determine the mass of Sgr A* to be $4.0_{-0.6}^{+1.1} \times 10^6$ M$_{\odot}$. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 65 pages, 35 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L15, 52pp., 2022 May 12

First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way

Authors: The Event Horizon Telescope Collaboration

Abstract: 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\,{\rm mm}$. The EHT data resolve a compact emission region with intrahour variability. A variety of… ▽ More 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\,{\rm 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 \pm 2.3$\,\uas (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 ${\sim}4 \times 10^6\,{\rm M}_\odot$, 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^\circ$), as well as non-spinning 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 galaxy, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of $10^3-10^5$ 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. △ Less

Submitted 16 November, 2023; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 24 pages, 6 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L12, 21pp., 2022 May 12

First Sagittarius A* Event Horizon Telescope Results. II. EHT and Multi-wavelength Observations, Data Processing, and Calibration

Authors: The Event Horizon Telescope Collaboration

Abstract: We present Event Horizon Telescope (EHT) 1.3 mm measurements of the radio source located at the position of the supermassive black hole Sagittarius A* (Sgr A*), collected during the 2017 April 5--11 campaign. The observations were carried out with eight facilities at six locations across the globe. Novel calibration methods are employed to account for Sgr A*'s flux variability. The majority of the… ▽ More We present Event Horizon Telescope (EHT) 1.3 mm measurements of the radio source located at the position of the supermassive black hole Sagittarius A* (Sgr A*), collected during the 2017 April 5--11 campaign. The observations were carried out with eight facilities at six locations across the globe. Novel calibration methods are employed to account for Sgr A*'s flux variability. The majority of the 1.3 mm emission arises from horizon scales, where intrinsic structural source variability is detected on timescales of minutes to hours. The effects of interstellar scattering on the image and its variability are found to be subdominant to intrinsic source structure. The calibrated visibility amplitudes, particularly the locations of the visibility minima, are broadly consistent with a blurred ring with a diameter of $\sim$50 $μ$as, as determined in later works in this series. Contemporaneous multi-wavelength monitoring of Sgr A* was performed at 22, 43, and 86 GHz and at near infrared and X-ray wavelengths. Several X-ray flares from Sgr A* are detected by Chandra, one at low significance jointly with Swift on 2017 April 7 and the other at higher significance jointly with NuSTAR on 2017 April 11. The brighter April 11 flare is not observed simultaneously by the EHT but is followed by a significant increase in millimeter flux variability immediately after the X-ray outburst, indicating a likely connection in the emission physics near the event horizon. We compare Sgr A*'s broadband flux during the EHT campaign to its historical spectral energy distribution and find both the quiescent and flare emission are consistent with its long-term behaviour. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 30 pages, 13 figures, published in The Astrophysical Journal Letters on May 12, 2022. See the published paper for the full authors list

Journal ref: The Astrophysical Journal Letters, Volume 930, Number 2, id L13, 31pp., 2022 May 12

Fitting the light curves of Sagittarius A* with a hot-spot model

Authors: A. I. Yfantis, M. A. Mościbrodzka, M. Wielgus, J. T. Vos, A. Jimenez-Rosales

Abstract: Sagittarius A* exhibits frequent flaring activity across the electromagnetic spectrum. Signatures of an orbiting hot spot have been identified in the polarized millimeter wavelength light curves observed with ALMA in 2017 immediately after an X-ray flare. The nature of these hot spots remains uncertain. We expanded existing theoretical hot-spot models created to describe the Sgr A* polarized emiss… ▽ More Sagittarius A* exhibits frequent flaring activity across the electromagnetic spectrum. Signatures of an orbiting hot spot have been identified in the polarized millimeter wavelength light curves observed with ALMA in 2017 immediately after an X-ray flare. The nature of these hot spots remains uncertain. We expanded existing theoretical hot-spot models created to describe the Sgr A* polarized emission at millimeter wavelengths. We sampled the posterior space, identifying best-fitting parameters and characterizing uncertainties. Using the numerical radiative transfer code ipole, we defined a semi-analytical model describing a ball of plasma orbiting Sgr A*, threaded with a magnetic field and emitting synchrotron radiation. We then explored the posterior space in the Bayesian framework of dynesty. We fit the static background emission separately, using a radiatively inefficient accretion flow model. We considered eight models with a varying level of complexity, distinguished by choices regarding dynamically important cooling, non-Keplerian motion, and magnetic field polarity. All models converge to realizations that fit the data, but one model without cooling, non-Keplerian motion, and magnetic field pointing toward us improved the fit significantly and also matched the observed circular polarization. Our models represent observational data well and allow testing various effects in a systematic manner. From our analysis, we have inferred an inclination of $155-160$ deg, which corroborates previous estimates, a preferred period of 90 minutes, and an orbital radius of $9-12$ gravitational radii. Our non-Keplerian models indicate a preference for an orbital velocity of $0.6-0.9$ times the Keplerian value. Last, all our models agree on a high dimensionless spin value ($a_{*}>0.8$), but the impact of spin on the corresponding light curves is subdominant with respect to other parameters. △ Less

Submitted 5 March, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Orbital Polarimetric Tomography of a Flare Near the Sagittarius A* Supermassive Black Hole

Authors: Aviad Levis, Andrew A. Chael, Katherine L. Bouman, Maciek Wielgus, Pratul P. Srinivasan

Abstract: The interaction between the supermassive black hole at the center of the Milky Way, Sagittarius A*, and its accretion disk occasionally produces high-energy flares seen in X-ray, infrared, and radio. One proposed mechanism that produces flares is the formation of compact, bright regions that appear within the accretion disk and close to the event horizon. Understanding these flares provides a wind… ▽ More The interaction between the supermassive black hole at the center of the Milky Way, Sagittarius A*, and its accretion disk occasionally produces high-energy flares seen in X-ray, infrared, and radio. One proposed mechanism that produces flares is the formation of compact, bright regions that appear within the accretion disk and close to the event horizon. Understanding these flares provides a window into accretion processes. Although sophisticated simulations predict the formation of these flares, their structure has yet to be recovered by observations. Here we show the first three-dimensional (3D) reconstruction of an emission flare recovered from ALMA light curves observed on April 11, 2017. Our recovery shows compact, bright regions at a distance of roughly six times the event horizon. Moreover, it suggests a clockwise rotation in a low-inclination orbital plane, consistent with prior studies by GRAVITY and EHT. To recover this emission structure, we solve an ill-posed tomography problem by integrating a neural 3D representation with a gravitational model for black holes. Although the recovery is subject to, and sometimes sensitive to, the model assumptions, under physically motivated choices, our results are stable, and our approach is successful on simulated data. △ Less

Submitted 16 April, 2024; v1 submitted 11 October, 2023; originally announced October 2023.