1 of 44

Event Horizon Telescope observations of Sgr A*

and the future of black hole imaging

Freek Roelofs

Center for Astrophysics | Harvard & Smithsonian

2 of 44

The EHT Collaboration

300+ members, 60 institutes, 20 countries in Europe, Asia, Africa, North and South America.

Slide credit: E. Traianou, MPIfR

EHT Collaboration Meeting

Granada, Spain

June 2022

3 of 44

The Event Horizon Telescope in 2017

Image Credits: ALMA/ESO, Sven Dornbusch, Junhan Kim, Helge Rottmann, David Sanchez, Daniel Michalik, Jonathan Weintroub, William Montgomerie, Tom Lowe, Serge Brunier

SMT, Arizona

JCMT, Hawaii

APEX, Chile

IRAM 30m Spain

LMT, Mexico

Photos: ALMA, Sven Dornbusch, Junhan Kim, Helge Rottmann, David Sanchez, Daniel Michalik, Jonathan Weintroub, William Montgomerie, Tom Lowe, Serge Brunier

Slide: Sara Issaoun

ALMA, Chile

SMA, Hawaii

SPT, South Pole

6 different locations

6 single-dish telescopes

2 phased arrays

4 of 44

Our Multi-Wavelength Partners

VERITAS

MAGIC

Fermi

Image credits: NSF/VERITAS, Juan Cortina, Vikas Chander, NASA, NASA/JPL-Caltech, NASA/CXC/SAO, NASA, ESO, P. Kranzler & A. Phelps, NRAO/AUI/NSF, HyeRyung, NAOJ, MPIfR/N. Tacken

Slide: Sara Issaoun

NuSTAR

Chandra

Swift

HESS

GBT

VLA

VERA

KVN

VLBA

VLT

GMVA

ALMA

5 of 44

The 2017 observations

6 of 44

Horizon-scale targets: M87* and Sagittarius A*

M87*

  • Mass
  • Distance
  • Inclination
  • Spin
  • Astrophysical model

Sagittarius A*

  • Mass
  • Distance
  • Inclination
  • Spin
  • Astrophysical model
  • Variability
  • Scattering

43 GHz with the Very Long Baseline Array/ Walker et al. 2018

22 GHz with the Very Large Array/ NRAO

Slide credit: S. Issaoun

7 of 44

First results: M87

8 of 44

Multi-wavelength view of M87

EHT 2021 Summer Collaboration Meeting

9 of 44

The supermassive black hole Sagittarius A*

20 as

Gravity Collaboration+ 2018

Closest supermassive black hole

  • Mass: 4.1x106 solar masses
  • Distance: 8.1 kpc

(Gravity Collaboration+ 2019, Do+ 2019)

Image Credits:

X-ray: NASA/CXC/UCLA/Z. Li et al

Radio 22 GHz: NRAO/VLA

S-stars: UCLA Galactic Center Group (Keck), Genzel et al. (2010), Yuan et al. (2003)

S2: Gravity Collaboration+ 2018, ESO/Gravity

Slide Credit: S. Issaoun

JVLA, 1.3cm

10 of 44

The supermassive black hole Sagittarius A*

20 as

Gravity Collaboration+ 2018

What does Sgr A* look like?

Expected size of the shadow of Sgr A*:

~50 𝜇as ~ 5 Schwarzschild radii

(Falcke+2000, Doeleman+2008, Fish+2011,

Johnson+2015, Fish+2016, Lu+ 2018)

What is the orientation of the black hole?

Is it spinning?

Long-standing debate: what emission process dominates in the radio (disk versus jet)?

Image Credits:

X-ray: NASA/CXC/UCLA/Z. Li et al

Radio 22 GHz: NRAO/VLA

S-stars: UCLA Galactic Center Group (Keck), Genzel et al. (2010), Yuan et al. (2003)

S2: Gravity Collaboration+ 2018, ESO/Gravity

Slide Credit: S. Issaoun

JVLA, 1.3cm

11 of 44

The SED of Sgr A* is highly variable, simultaneous coverage is crucial for modeling

Slide credit: D. Haggard, McGill / EHT MWL Science WG

The supermassive black hole Sagittarius A*

12 of 44

Our 2017 Multi-Wavelength Effort on Sagittarius A*

Figure credit: M. Johnson, J. Farah, SAO/EHT MWL Science WG.

13 of 44

First results: Sgr A*

14 of 44

A Tale of Two Black Holes

55 million light-years away from us

6.5 billion solar masses

27 thousand light-years away from us

4 million solar masses

15 of 44

A Tale of Two Black Holes

Visualization by L. Medeiros, IAS/ xkcd

16 of 44

From recordings to image

17 of 44

What is imaging?

Aperture synthesis: Earth rotation helps us fill our virtual mirror and combine data on multiple temporal and spatial scales

M87 becomes visible to stations as Earth rotates

Our coverage, or “virtual mirror” fills up

Our data, with the help of scan-averaging, give high SNR information on multiple spatial scales that improves our image

Animation credit: D. Palumbo, Harvard/ M. Wielgus, MPIfR

18 of 44

The imaging process

Data calibration

Blind imaging in independent teams

Pipeline building per independent software

Image validation

Synthetic Models

Consistency with calibrator analysis

Goodness of fit to data

Structure robustness

Step 1

Step 2

Step 3

Feedback on data quality

Slide credit: S. Issaoun

19 of 44

Event-horizon scale variability in Sgr A* data

Visibility amplitudes

Closure phases

Q-metric: Roelofs et al. (2017)

Broderick et al. (2022)

Georgiev et al. (2022)

GRMHD simulations

Interstellar scattering

The detected variability is most statistically significant on baselines between ~2.5 and 6 Gλ.

20 of 44

Imaging a variable source

Visualization by C. Fromm, U. Wurzburg/ L. Rezzola U. Frankfurt

21 of 44

Variability mitigation and geometric modeling

  1. “Snapshot” approach fits to short segments of data over which the source is static, and then averages these independent fits�

  • “Full-track” approach fits the average source structure alongside a parameterized model for the variability about that average source structure�

The detected variability is most statistically significant on baselines between ~2.5 and 6 Gλ.

Ring

Angular Profile

(Fourier Modes)

22 of 44

Modeling and imaging results

The detected variability is most statistically significant on baselines between ~2.5 and 6 Gλ.

ring-like model

non-ring-like model

We quantify various morphological parameters using both geometric modeling and image-domain feature extraction techniques

  • Ring diameter: 51.8 ± 2.3 μas
  • Ring thickness: FWHM between ~30-50% of the ring diameter
  • other morphological quantities – magnitude and orientation of the asymmetry, depth of the central brightness depression – are poorly constrained and depend on the measurement method

23 of 44

Slide credit: H. Shiokawa

Simulating a black hole and its environment: GRMHD

24 of 44

Ray tracing

25 of 44

GRMHD + ray tracing parameters

Mass

Distance

Astrophysics

  • Electron temperature distribution
  • Magnetic field strength (MAD/SANE)

Inclination

Spin

Scattering

M. Moscibrodzka, Radboud U.

Moscibrodzka & Gammie 2018

26 of 44

Sgr A* GRMHD library

Animation Credit: B. Prather / UIUC

27 of 44

Gravitational radius and black hole mass

Calibrating size measurements using synthetic data generated from GRMHD simulations:

  • θg = G M / (c2 D) = 4.8 (+1.4,-0.7) μas
  • Uncertainty arises from
    • Model flexibility needed to capture structural variability in the source
    • Broad morphological diversity of the GRMHD calibration suite (reflecting the unknown inclination of Sgr A*)

Combining the gravitational radius constraint with an independent distance measurement from maser parallaxes (Reid et al. 2019), we determine the mass of Sgr A* to be 4.0 (+1.1,-0.6) x 106 M

Both the gravitational radius and mass measurements are consistent with the more precise constraints obtained from stellar orbits (Do et al. 2019, Gravity Collaboration et al. 2019, 2020)

28 of 44

Testing General Relativity

Even though tests with gravitational waves and black hole images span black hole masses that are different by 8 orders of magnitude, they are all consistent with the GR predictions that all black holes are described by the same metric, independent of their mass.

29 of 44

Future Arrays: Synthetic Data

Plus:

→ Telescope locations

→ Observation schedule

→ Calibration pipeline

Aperture efficiency

Absorption, emission, turbulence, and delays in the atmosphere

→resulting in lower signal-to-noise ratios, phase corruptions, and amplitude systematics

Instrumental (thermal) noise, gains, and leakage

Antenna pointing offsets

→resulting in systematic amplitude gain offsets

Antenna diameter

Source model

Adapted from diagram by S. Issaoun

Data generation codes:

  • eht-imaging (Chael+ 2016, 2018)
  • SYMBA (Roelofs+ 2020, Janssen+ 2019, Natarajan+ 2022)

30 of 44

Array expansion

Davelaar et al. (2019)

Roelofs, Janssen et al. (2020)

Gamsberg, Namibia

31 of 44

32 of 44

GRMHD and hotspot reconstructions

Noemi La Bella et al., in prep.

Western

Eastern + Africa

33 of 44

The next-generation Event Horizon Telescope

230 GHz

Spans ~2 weeks

Observes ~7 nights/year

EHT Sites

34 of 44

The next-generation Event Horizon Telescope

86+230+345 GHz

Spans 3+ months

Observes 60+ nights/year

EHT Sites

ngEHT Sites (Phase 1)

35 of 44

The next-generation Event Horizon Telescope

86+230+345 GHz

Spans 3+ months

Observes 60+ nights/year

EHT Sites

ngEHT Sites (Phase 2)

36 of 44

ngEHT threshold science goals

M. Johnson

P. Tiede

37 of 44

ngEHT objective science goals

M. Johnson

38 of 44

ngEHT Analysis Challenges: process overview

  • Start from a “ground truth” source model
    • Can be tailored to a specific science case
    • Can be multi-frequency
    • Static or time-variable�
  • Simulate ngEHT observations of the source model
    • Using different arrays/design choices
    • Adding in realistic observation and calibration effects�
  • Participants analyze the synthetic datasets, trying to maximize the science output
    • Image and movie reconstruction
    • Model fitting and parameter estimation
    • Any other analysis�
  • Submissions are evaluated using quantitative metrics

C. Fromm

Roelofs et al. 2022 (in prep.)

39 of 44

The next-generation EHT: Reference Array

EHT2022 + 10 additional sites (incl. HAY, OVRO)

ngEHT reference array 1 (A. Raymond)

Performs well on uv-coverage metrics

6 m dishes, 8 GHz bandwidth, median April weather

Roelofs et al. 2022 (in prep.)

40 of 44

ngEHT Analysis Challenges: M87

Roelofs et al. 2022 (in prep.)

86 GHz

230 GHz

Model by K. Chatterjee and R. Emami

Reconstructions by P. Arras and J. Knollmueller using Resolve Bayesian imaging algorithm

Weekly observations of M87 for 5 months

41 of 44

ngEHT Analysis Challenges: Sgr A*

Roelofs et al. 2022 (in prep.)

Model by P. Tiede

Reconstructions by A. Fuentes

Intra-hour variability of Sgr A*

10-minute scans with 10-minute gaps

42 of 44

Space VLBI: Event Horizon Imager

Two or three satellites in circular orbits with slightly different radius

Dense spiral-shaped uv-coverage

Frequency up to 690 GHz

4 μas nominal resolution

Roelofs et al. (2019)

43 of 44

Image reconstruction and parameter estimation

Davelaar et al. (2019)

Roelofs et al. (2019)

Sub-percent precision mass measurements → GR tests

Fitting EHT 2017 M87 GRMHD library to simulated data

Roelofs et al. (2021)

Constraints on black hole spin and other GRMHD parameters

44 of 44

Summary

  • We have imaged the black hole shadow of M87 and Sgr A*
  • Expansions of the array are currently happening
  • The ngEHT will a be transformative improvement, opening new windows in angular resolution, sensitivity, frequency coverage, and time coverage
  • On the longer term, a Space VLBI mission will take high-precision black hole imaging even further