REVIEW 3 major objections 5 minor 55 references
Enhanced imaging of M87*: Simulations with the EHT and extended-KVN
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Adding the four-telescope eKVN array to the EHT would roughly halve residual noise in recovered M87* jet images, simulations show.
desk verdict Solid simulation study: eKVN short baselines improve EHT jet reconstruction for M87*, but the headline gain is tuned to the single ground-truth model and needs robustness tests before being quoted. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is short-baseline uv coverage—the set of baseline lengths and orientations the array samples. eKVN is a four 21 m dish array in South Korea with baselines of roughly 130–500 km; at 230 GHz those provide 74–370 Mλ spacings, sensitivity to angular scales of about 2.1–0.5 mas, and at 86 GHz they supply the array's shortest spacings. The EHT alone has only two baselines below 1 Gλ, leaving a gap near 0.7 Gλ. Image reconstruction proceeds by regularized maximum likelihood, first using gain-independent closure quantities, then self-calibration, then visibility amplitudes; multi-frequency synthesis fits a spectral-index map with smoothness regularizers. The short baselines constrain the Fourier components that carry the diffuse jet, which is why their absence shows up as residual noise rather than as a missing ring.
What would settle it
Rerun the same synthetic-observation pipeline on a different GRMHD snapshot or a different weather realization: if $\rho_{\rm NX}^{\rm jet}$ does not rise from about $0.92$ to $0.98$ when eKVN is added, the gain is model-specific. A data-side check is to flag the KYS station in real April 2024 EHT observations of M87* and measure whether jet-region residual noise increases by the roughly factor of two the simulations predict.
Extended reading notes
Core claim
On the paper's own terms, the discovery is that the EHT's short-baseline deficit—not angular resolution or long-baseline sensitivity—is what limits jet recovery, and eKVN removes that deficit. With eKVN added to the full 2022-era EHT array, the root-mean-square difference between the ground-truth and reconstructed image becomes a half for the jet region, $\rho_{\rm NX}^{\rm jet}$ improves from $0.92$ to $0.98$, and the log-scale correlation improves markedly. Jackknife tests that flag one or two EHT stations show the gain is largest exactly where the EHT alone is weakest: when Chile, Europe, or the US-mainland stations are removed, eKVN keeps the jet visible while EHT-only jet similarity drops to $0.76$–$0.78$. In multi-frequency synthesis at 86+230 GHz, eKVN removes diffuse residual emission outside the true source structure, so the spectral-index map's high-fidelity region traces the jet rather than noise, and the ring ellipticity is recovered more stably.
Load-bearing premise
The load-bearing premise is that the single general-relativistic magnetohydrodynamic (GRMHD) simulation used as ground truth faithfully represents the real M87* source, and that rescaling the total flux to that known model before self-calibration is a fair stand-in for real data; if the real jet morphology, flux scale, or weather differs, the reported improvement could shrink.
Editorial extensions
If this is right
- EHT+eKVN observations of M87* at 230 GHz should recover the jet alongside the black-hole shadow with about half the residual noise of the EHT alone.
- Losing one or two EHT stations—for example in Chile, Europe, or the US mainland—no longer removes the short-baseline coverage needed for the jet, since eKVN restores most of the lost fidelity in jackknife tests.
- Simultaneous 86 and 230 GHz imaging with eKVN improves both single-frequency images and gives a spectral-index map whose high-fidelity region follows the jet, whereas EHT-only synthesis leaves diffuse residual emission.
- At 86 GHz the eKVN baselines are the shortest in the array, giving a better constraint on the compact flux that separates the bright ring from the faint jet.
- For monitoring campaigns and Sgr A* dynamic imaging, the added snapshot uv coverage should improve movie fidelity.
Reading between the lines
- If the gain scales with the number of baselines below about 0.7 Gλ, then any future short-spacing addition to the EHT or ngEHT—not only the Korean sites—should produce comparable jet-recovery improvements; this is testable by repeating the simulation with a different short-baseline station.
- Because all metrics are computed against one ground-truth model, the factor-of-two noise reduction is best read as a demonstration of mechanism; rerunning with other jet morphologies would show whether the benefit is universal.
- The 86-plus-230 GHz result suggests a cheaper path to jet science: the partial-participation case with only two eKVN stations at 230 GHz and all four at 86 GHz already delivers most of the multi-frequency gain, so simultaneous-band capability may matter more than full 230 GHz coverage.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents an end-to-end imaging simulation for M87* with the 2022-era EHT and with the addition of the four telescopes of the extended KVN. Synthetic observations at 86 and 230 GHz are generated with ngehtsim from a Chael et al. (2019) GRMHD image, reconstructed with eht-imaging RML, and evaluated via normalized cross-correlation against the ground truth, with separate ring/jet masks, jackknifed station loss, and 86+230 GHz multi-frequency synthesis with spectral-index recovery. The central claim is that eKVN's six short baselines below 1 Gλ substantially improve jet-structure recovery, reduce residual noise, and make the array more robust to station loss.
Significance. If the reported gains are genuine, the paper provides quantitative support for a concrete, near-term array addition to the EHT: eKVN operates at 230 GHz, is already partially equipped, and helps fill a known short-baseline gap. The simulation pipeline uses public, community-standard tools (ngehtsim, eht-imaging), a realistic 2024 array schedule, and transparent uv-coverage and jackknife diagnostics. The qualitative conclusion that short baselines help large-scale jet and spectral-index recovery is physically plausible and is well supported by the uv-coverage comparison. However, all headline numbers are in-sample estimates derived from a single ground-truth model and evaluation choices that use that model's known flux and morphology; the magnitudes, and possibly even the sign, of the improvement for real observations are not yet established.
major comments (3)
- [Section 2.2 and Appendix B] The quantitative comparison is not blind. Section 2.2 states that the fiducial RML hyperparameters are chosen by maximizing ρNX to the groundtruth image, and Appendix B chooses the 8 μas blurring kernel as the one that maximizes ρNX between the reconstructed and groundtruth images. Sections 3.1 and 3.2 then report ρNX and r.m.s. differences computed against the same blurred ground truth. Consequently, the 0.92→0.98 ρjet_NX gain and the factor-of-two r.m.s. reduction are optimized with respect to the test image; they quantify best-case in-sample fidelity, not predictive image fidelity for an unknown M87* realization. The spectral-index regularizers (l2=10, TV=50) are selected the same way. The qualitative uv-coverage argument survives, but the quantitative headline numbers need a blinded evaluation, for example by fixing hyperparameters a priori or selecting them on a training model and then applying them to a held-out model or snapshot.
- [Section 2.2 and Section 3.2] The reconstruction uses information that real EHT+eKVN observations will not have. Section 2.2 rescales the total flux before self-calibration using the known groundtruth flux density, and Section 3.2 estimates the residual noise σs outside a source region that is 'guided by the groundtruth image.' Because the jet is faint and its absolute flux is weakly constrained by closure-only imaging, the reported reduction of residual noise and the spectral-index reliability partly reflect this injected information. The authors should rerun at least one comparison without groundtruth-based flux rescaling, or with the flux estimated from the data itself (for example from short-baseline or zero-spacing constraints), and restate the conclusions accordingly.
- [Sections 2.1, 3.1, and 3.2] All quantitative claims rest on one ground-truth image (Chael et al. 2019), one observing date (April 6, 2024), one ngehtsim weather/SEFD draw, and one noise realization. No uncertainty estimates accompany the reported ρNX and r.m.s. values, and no test shows that the eKVN gain persists for different jet morphologies, source flux scales, or weather conditions. Since the paper's Discussion emphasizes stability and monitoring, it should include at least a small ensemble with multiple GRMHD snapshots or another M87* model and independent weather realizations, reporting the distributions of the fidelity metrics and of the ring parameters.
minor comments (5)
- [Equation (1)] The sign of the regularizer term appears to be wrong: in a minimization of J(I), positive penalty terms should be added, not subtracted, unless S_R is defined with an unconventional sign.
- [Figure 1 caption] The caption contains a typo: 'telesctopes' should be 'telescopes.'
- [Section 4] The count of 'six more baselines shorter than 1 Gλ' is used to support the main argument, but it appears only in the Discussion; stating this count in Section 2.1 alongside the uv-coverage discussion would make the argument easier to follow.
- [Appendix A] Figure A.1 is described as 'Same as Figure 4 (top)', but it shows only 230 GHz images and adds the KYS+KPC configuration; the caption should state exactly which panels are reproduced.
- [General] The paper would benefit from a data and software availability statement specifying ngehtsim and eht-imaging versions and whether the synthetic data and imaging scripts are publicly available.
Circularity Check
Headline fidelity metrics are tuned against the same ground truth used to measure them; the directional eKVN improvement remains independent.
-
fitted input called prediction
[Section 2.2 (Imaging) and Appendix B (Effective angular resolution); results reported in Section 3.1]
"As a result, 81 images were obtained and the fiducial image was selected based on the normalized cross correlation (ρNX) to the groundtruth image ... By this way, we found that a beam size of 8µas yields the highest ρNX corresponding to ∼ 1/3 of the nominal resolution in our data."
The headline fidelity metrics (ρNX and the r.m.s. difference computed with the same 8 µas blur) are measured against the same groundtruth image that was used to choose the RML hyperparameters and the blurring kernel. The kernel is explicitly chosen to maximize ρNX, so the reported ρNX and r.m.s. values are partly constructed by the evaluation target rather than being blind measurements of recovery quality. This does not force the EHT vs EHT+eKVN difference, since the identical choices are applied to both arrays, but the absolute 'recovery' numbers are optimistic and are not out-of-sample predictions.
-
fitted input called prediction
[Section 2.2 (Imaging) and Section 3.2 (86−230 GHz synthesis)]
"After the survey, the fiducial parameters for α were selected as l2 = 10 and TV = 50 based on the highest ρNX at both frequencies ... Note that σs was estimated outside the source structure that was guided by the groundtruth image."
The spectral-index regularizer hyperparameters and the residual-noise masking region are selected using the same groundtruth image against which the spectral-index fidelity is then claimed ('the spectral index distribution is better constrained'). The selection criterion is exactly the reported metric (ρNX), and the mask is groundtruth-guided, so the spectral-index comparison is not a blind prediction. As with the 230 GHz case, the differential EHT vs EHT+eKVN claim is not mathematically forced, but the absolute fidelity of the α map is partly determined by groundtruth inputs rather than independently predicted.
full rationale
This is a simulation and evaluation paper rather than a derivation, so the circularity burden falls on the evaluation design. The central claim—that adding eKVN short baselines improves jet recovery—is not equivalent to any input: it is a differential comparison between two arrays using the same groundtruth, the same data-generation pipeline, and the same imaging parameters. The groundtruth model (Chael et al. 2019) and the software tools (ngehtsim, eht-imaging) are external resources, and no load-bearing uniqueness theorem is imported from the authors' prior work. However, several evaluation choices are tuned against the same groundtruth that later serves as the yardstick: the RML hyperparameters and the 8 µas blur kernel are selected to maximize ρNX, the spectral-index regularizers are chosen by the same ρNX criterion, and total flux is rescaled to the known groundtruth flux before self-calibration. These make the absolute fidelity numbers partly constructed by the target and optimistic relative to real observations. Because the identical choices are applied to both arrays, the directional claim that eKVN improves jet reconstruction retains independent content. Score 4 reflects partial circularity in the evaluation metrics, not a reduction of the headline claim to its inputs.
Assumptions & free parameters
free parameters (4)
- RML regularizer hyperparameters =
l1=0.1, MEM=1, log(TSV)=0.1, compact2=1
- Spectral index regularizer hyperparameters =
l2=10, TV=50
- Total flux normalization =
groundtruth total flux
- Effective resolution blurring kernel =
8 uas Gaussian
assumptions (5)
- domain assumption The Chael et al. (2019) GRMHD simulation image is an adequate stand-in for the true M87* structure over the observed scales.
- domain assumption ngehtsim's heuristic detection criteria, atmospheric models, and SEFD estimates reproduce realistic EHT/eKVN observing conditions.
- standard math The Fourier relation between sky brightness and interferometric visibility (van Cittert-Zernike theorem) underlies image reconstruction.
- domain assumption The same fiducial RML hyperparameters are suitable for both EHT-only and EHT+eKVN data.
- domain assumption Total flux density of the source is known a priori and can be used to rescale images.
Cite this review
Pith. "Pith review of Enhanced imaging of M87*: Simulations with the EHT and extended-KVN." pith.science (2026). https://pith.science/paper/WWJFGGF2
@misc{pith2026250111822,
author = {Pith},
title = {Pith review of: Enhanced imaging of M87*: Simulations with the EHT and extended-KVN},
year = {2026},
howpublished = {\url{https://pith.science/paper/WWJFGGF2}},
note = {Machine review of arXiv:2501.11822}
}
read the original abstract
The Event Horizon Telescope (EHT) has successfully revealed the shadow of the supermassive black hole, M87*, with an unprecedented angular resolution of approximately 20 uas at 230 GHz. However, because of limited short baseline lengths, the EHT has been constrained in its ability to recover larger scale jet structures. The extended Korean VLBI Network (eKVN) is committed to joining the EHT from 2024 that can improve short baseline coverage. This study evaluates the impact of the participation of eKVN in the EHT on the recovery of the M87* jet. Synthetic data, derived from a simulated M87* model, were observed using both the EHT and the combined EHT+eKVN arrays, followed by image reconstructions from both configurations. The results indicate that the inclusion of eKVN significantly improves the recovery of jet structures by reducing residual noise. Furthermore, jackknife tests, in which one or two EHT telescopes were omitted - simulating potential data loss due to poor weather - demonstrate that eKVN effectively compensates for these missing telescopes, particularly in short baseline coverage. Multi-frequency synthesis imaging at 86-230 GHz shows that the EHT+eKVN array enhances the recovered spectral index distribution compared to the EHT alone and improves image reconstruction at each frequency over single-frequency imaging. As the EHT continues to expand its array configuration and observing capabilities to probe black hole physics more in depth, the integration of eKVN into the EHT will significantly enhance the stability of observational results and improve image fidelity. This advancement will be particularly valuable for future regular monitoring observations, where consistent data quality is essential.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
- [3]
- [4]
- [5]
- [6]
- [7]
-
[8]
Peterson, B. M. 1997, An Introduction to Active Galactic Nuclei, Cambridge: Cambridge University Press
work page 1997
Show all 55 references
-
[9]
2007, A&A, 463, 551
Sch\" o del, R., et al. 2007, A&A, 463, 551
2007
-
[10]
1995, A&A, 300, 707
Timmer, J., & K\" o nig, M. 1995, A&A, 300, 707
1995
-
[11]
2007, MNRAS, 375, 764
Trippe, S., et al. 2007, MNRAS, 375, 764
2007
-
[12]
2013, JKAS, 46, 133
Trippe, S. 2013, JKAS, 46, 133
2013
-
[13]
2022, Galaxies, 10, 113
Akiyama , K., Algaba , J.-C., An , T., et al. 2022, Galaxies, 10, 113
2022
-
[14]
2012, , 745, L28
Asada , K., & Nakamura , M. 2012, , 745, L28
2012
-
[15]
2017, arXiv e-prints, arXiv:1705.04776
Asada , K., Kino , M., Honma , M., et al. 2017, arXiv e-prints, arXiv:1705.04776
2017 arXiv
-
[16]
2023, arXiv e-prints, arXiv:2312.02130
Ayzenberg , D., Blackburn , L., Brito , R., et al. 2023, arXiv e-prints, arXiv:2312.02130
2023 arXiv
-
[17]
D., & Payne , D
Blandford , R. D., & Payne , D. G. 1982, , 199, 883
1982
-
[18]
D., & Znajek , R
Blandford , R. D., & Znajek , R. L. 1977, , 179, 433
1977
-
[19]
L., Johnson , M
Bouman , K. L., Johnson , M. D., Dalca , A. V., et al. 2017, arXiv e-prints, arXiv:1711.01357
2017 arXiv
-
[20]
D., Falcke , H., Law , C
Brinkerink , C. D., Falcke , H., Law , C. J., et al. 2015, , 576, A41
2015
-
[21]
W., et al
Chael , A., Issaoun , S., Pesce , D. W., et al. 2023, , 945, 40
2023
-
[22]
Chael , A., Narayan , R., & Johnson , M. D. 2019, , 486, 2873
2019
-
[23]
A., Johnson , M
Chael , A. A., Johnson , M. D., Bouman , K. L., et al. 2018, , 857, 23
2018
-
[24]
A., Johnson , M
Chael , A. A., Johnson , M. D., Narayan , R., et al. 2016, , 829, 11
2016
-
[25]
2022, , 926, 108
Cho , I., Zhao , G.-Y., Kawashima , T., et al. 2022, , 926, 108
2022
-
[26]
2023, , 621, 711
Cui , Y., Hada , K., Kawashima , T., et al. 2023, , 621, 711
2023
- [27]
-
[28]
S., Barrett , J., Blackburn , L., et al
Doeleman , S. S., Barrett , J., Blackburn , L., et al. 2023, Galaxies, 11, 107
2023
-
[29]
2019a, , 875, L1
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019a, , 875, L1
-
[30]
2019b, , 875, L2
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019b, , 875, L2
-
[31]
2019c, , 875, L3
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019c, , 875, L3
-
[32]
2019d, , 875, L4
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019d, , 875, L4
-
[33]
2019e, , 875, L5
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019e, , 875, L5
-
[34]
2019f, , 875, L6
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2019f, , 875, L6
-
[35]
2022 a , , 930, L14
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2022 a , , 930, L14
2022
-
[36]
2022 b , , 930, L12
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2022 b , , 930, L12
2022
-
[37]
2024, , 681, A79
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2024, , 681, A79
2024
-
[38]
2017, Journal of Infrared, Millimeter, and Terahertz Waves, 38, 1487
Han , S.-T., Lee , J.-W., Lee , B., et al. 2017, Journal of Infrared, Millimeter, and Terahertz Waves, 38, 1487
2017
-
[39]
W., Roelofs , F., et al
Issaoun , S., Pesce , D. W., Roelofs , F., et al. 2023, Galaxies, 11, 28
2023
-
[40]
2023, Galaxies, 11, 3
Jiang , W., Zhao , G.-Y., Shen , Z.-Q., et al. 2023, Galaxies, 11, 3
2023
-
[41]
D., Bouman , K
Johnson , M. D., Bouman , K. L., Blackburn , L., et al. 2017, , 850, 172
2017
-
[42]
D., Akiyama , K., Blackburn , L., et al
Johnson , M. D., Akiyama , K., Blackburn , L., et al. 2023, Galaxies, 11, 61
2023
-
[43]
A., et al
Kim , J.-Y., Lee , S.-S., Hodgson , J. A., et al. 2018, , 610, L5
2018
-
[44]
P., et al
Lu , R.-S., Asada , K., Krichbaum , T. P., et al. 2023, , 616, 686
2023
-
[45]
Mo \'s cibrodzka , M., Falcke , H., Shiokawa , H., & Gammie , C. F. 2014, , 570, A7
2014
-
[46]
2019, , 887, 147
Park , J., Hada , K., Kino , M., et al. 2019, , 887, 147
2019
-
[47]
W., Blackburn , L., Chaves , R., et al
Pesce , D. W., Blackburn , L., Chaves , R., et al. 2024 a , arXiv e-prints, arXiv:2404.01482
2024 arXiv
-
[48]
W., Blackburn , L., Chaves , R., et al
Pesce , D. W., Blackburn , L., Chaves , R., et al. 2024 b , ngEHT simulation tools , v1.0.0, Zenodo, 10.5281/zenodo.10722363
2024 doi
-
[49]
2011, , 141, 114
Rioja , M., & Dodson , R. 2011, , 141, 114
2011
-
[50]
J., Dodson , R., & Asaki , Y
Rioja , M. J., Dodson , R., & Asaki , Y. 2023, Galaxies, 11, 16
2023
-
[51]
2023, Galaxies, 11, 12
Roelofs , F., Blackburn , L., Lindahl , G., et al. 2023, Galaxies, 11, 12
2023
-
[52]
2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol
Shin , N., Park , Y.-S., Kang , H., et al. 2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 11453, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, ed. J. Zmuidzinas & J.-R. Gao , 114532Z, 10.1...
2020 doi
-
[53]
E., & Palumbo , D
Tiede , P., Broderick , A. E., & Palumbo , D. C. M. 2022, , 925, 122
2022
-
[54]
2022, , 930, L19
Wielgus , M., Marchili , N., Mart \' -Vidal , I., et al. 2022, , 930, L19
2022
-
[55]
Salpeter, E. E. 1955, The Luminosity Function and Stellar Evolution, ApJ, 121, 161
1955
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.