{"id":"579db616-3ec5-4d01-a8f3-4c7e1f50f8b0","arxiv_id":"2502.03026","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A method to measure ring diameter, width, and asymmetry from the zeros of the interferometric visibility function is derived and tested on M87*, yielding r≈23 μas and B=0.48.","lead":"The paper derives analytical visibility functions for asymmetric Gaussian ring models and extracts the ring diameter, width, and asymmetry from zero-crossings of the visibility amplitude. Applied to Event Horizon Telescope data for M87*, the method reports a ring radius of about 23 microarcseconds and an asymmetry parameter of 0.48.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The M87 B estimate ignores the model's free position angle ϕ0; the two fixed-axis fits cannot separate B from orientation, so B=0.48 may be an artifact.","rationale":"The paper's analytical framework is internally consistent for the idealized model: the visibility expressions (Eqs. 7, 8, 11, 30) are plausible, and the simulated tests demonstrate that zeros of J0 track the ring radius for thin rings with reasonable accuracy. However, the application to M87* is where the central claim becomes load-bearing, and there the fitting procedure omits a key model parameter. The model's asymmetry is defined by both an amplitude B and a position angle ϕ0 (Eq. 2), and these two parameters enter the visibility cuts in a degenerate way unless ϕ0 is known or fitted. The paper instead fixes the two cuts to the coordinate axes without justification, so the reported B=0.48 is not identifiable from the chosen data products. This is not merely an external disagreement about realism; it is an internal incompleteness in the inference step. The reader's weakest assumption—that the image model may not match M87*—is related but broader; our concern is a specific, testable flaw within the model as applied. If the proposed test shows that ϕ0 is unconstrained or that B changes dramatically with ϕ0, the M87 result should be downgraded. If the test shows that the fit is stable and B≈0.48 regardless of ϕ0, the concern would be resolved. Because the paper does not currently provide the information needed to settle this, the conditional status is appropriate, and the verdict need not change.","tokens_in":13621,"tokens_out":14668,"duration_ms":131539,"concrete_test":"Re-fit the digitized EHT visibility points of Fig. 8 using Eq. (7) with ϕ0 free (scanning 0 to π) and compare the best-fit (r, B) and chi-square with the published fixed-axis fit. If the best-fit ϕ0 is not consistent with the assumed alignment, or if B changes significantly, then the published B=0.48 is not the intrinsic asymmetry. A stronger synthetic test: generate visibility from a Gaussian ring with known r=21 μas, B=0.48, and ϕ0=30°, sample it at the actual M87* (u,v) coverage, and apply the same two-cut, fixed-axis fitting; if the recovered B deviates from 0.48, the M87 application is invalid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central M87 result—r≈23 μas and B=0.48—is obtained by fitting the thin-ring expression (Eq. 7) along two perpendicular cuts fixed to the left-right and up-down axes of the (u,v) plane. In Eq. (7), the visibility amplitude is |V(u,ϕu)| = 1/2 sqrt((2−B)^2 J0^2 + B^2 J1^2 cos^2(ϕu−ϕ0)), where ϕ0 is the position angle of the brightness maximum. The difference between the two cuts therefore depends jointly on B and ϕ0. Section III.A rotates the model so that the maximum is on the X-axis before taking sections, but Section V does not perform or report such an alignment for M87*; it simply chooses left-right and up-down axes and never fits or quotes ϕ0. If the true bright-spot direction is not aligned with the chosen axes, the left-right cut is not the pure J0 term (so the fitted radius is biased), and the up-down cut yields only a projection of the asymmetry—not the intrinsic ring-asymmetry parameter B. Since the paper claims B=0.48 as a first measurement, this degeneracy directly undermines the headline application. A separate, smaller issue is that finite ring width biases zero-derived radii by ~5–7% at w=10 μas (their own Table I), but the orientation degeneracy is the more fundamental problem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical model of an asymmetric Gaussian ring, derives closed-form visibility expressions for thin, thick, parabolic, and Gaussian-profile rings (Eqs. 5–11 and 30 in the Appendix), and proposes to estimate the ring radius from the positions of visibility-function zeros and the ring width from the first envelope zero. The method is tested on simulated rings of the same model family and on two-ring configurations. In the final section, the thin-ring formula (7) is fitted to two perpendicular cuts of the M87* EHT visibility amplitudes, yielding a ring radius r ≈ 23 μas and an asymmetry parameter B = 0.48.","tokens_in":13916,"tokens_out":13968,"duration_ms":131480,"significance":"If made robust, the paper would provide a simple and transparent frequency-domain route to ring parameters from VLBI visibilities, and the explicit Fourier transforms and the zero/envelope relations are useful additions to the methods toolkit. The analytical derivation is a genuine strength: the visibility formulas are obtained directly from a stated brightness model, and the zero-locating approach is standard. However, the headline M87* inference currently depends on an unstated alignment assumption and lacks uncertainty propagation, so the claimed B = 0.48 is not yet a reliable measurement.","major_comments":[{"comment":"The extraction of B = 0.48 from two fixed perpendicular cuts is not invertible as presented. For the thin-ring model, |V(u,φu)| = 1/2 sqrt((2−B)^2 J0^2(2πur0) + B^2 J1^2(2πur0) cos^2(φu−φ0)). Two cuts at φu = 0 and φu = π/2 therefore constrain two projections of the three parameters (r0, B, φ0), and B is determined only after fixing φ0. Section III.A explicitly states that the model is rotated so that the brightness maximum lies on the negative X-axis before sections are taken, but Section V does not report or perform such an alignment for M87*; it simply selects left–right and up–down axes. Unless φ0 is fitted jointly or fixed with a justified position angle, the “left–right” cut is generally not a pure J0 section, the zero-derived radius is biased, and the “up–down” cut gives only a projection of the asymmetry rather than the intrinsic parameter B. A concrete remedy is to fit Eq. (7) to the full visibility amplitudes with (r0, B, φ0) free and to compare that fit with the two-axis result.","section":"Section V, Eq. (7), Section III.A"},{"comment":"The M87* estimates are based on visibility points digitized from published figures without uncertainty estimates, and the fitted parameters are quoted without error bars. The statement that the asymmetry parameter was “estimated for the first time” therefore cannot be evaluated: there is no variance from the digitization, no systematic from the choice of the two cuts, and no goodness-of-fit. The paper should report parameter uncertainties and, ideally, fit the model to the calibrated visibility data; at minimum it should provide a Monte Carlo sensitivity test that randomizes the digitized points and the section orientations.","section":"Section V and Fig. 8"},{"comment":"The zero-location radius estimator carries a finite-width bias that is visible in the paper's own Table I: for a true ring width w = 10 μas, the first three zero-based radii are 21.0–21.4 μas instead of 20.0 μas, i.e. a 5–7% bias. Section V applies the infinitely thin formula (7) to M87* and quotes r ≈ 23 μas without applying a finite-width correction or adding this systematic error. Since the fit does not constrain the M87* width, this bias should be propagated into the quoted radius before comparing with EHT's r ≈ 21 μas estimate.","section":"Section V and Table I"},{"comment":"The validation is performed only on simulated images generated from the same Gaussian asymmetric-ring family used to construct the estimators; the method has not been tested on GRMHD simulation images or on models with ellipticity, jet emission, or n > 1 angular profiles. This is a correctness-risk for the M87* application, which interprets B = 0.48 as a physical asymmetry. A concrete test would be to run the estimator on a set of GRMHD images with EHT-like (u,v) coverage and report the resulting scatter and bias.","section":"Sections II–IV"}],"minor_comments":[{"comment":"The text says “for a thin ring (w ≫ r0)”; this should be w ≪ r0.","section":"Section VI"},{"comment":"The caption states that both the red and black dots correspond to φu = π; the second direction should presumably be φu = π/2.","section":"Fig. 4 caption"},{"comment":"The three rows labeled “r0, μas” are not identified with their zero orders; adding column headers such as u0, u1, u2 would make the bias trend clearer.","section":"Table I"},{"comment":"The use of “plotdigitizer” and “curvefit” is informal; the paper should specify the digitization procedure, the number of points, the fitting algorithm, and the initial values used.","section":"Section V"},{"comment":"There are several typographical errors, including “obtaoined”, “t M87*”, and inconsistent use of commas in numerical values; these should be corrected in a careful revision.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of astro-ph.IM and the analytical core is sound, so I would not reject it. The main weakness is the M87* orientation degeneracy in Section V; this is fixable by fitting the position angle and by reporting uncertainties. I also recommend that the authors check the EHT literature carefully before repeating the claim that the asymmetry parameter is estimated for the first time, since that claim may be stronger than the current evidence supports."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has a solid analytic core: the Gaussian-ring visibility model (Eq. 30) is new, the Fourier derivations look correct, and the recipe of using Bessel zeros for radius and envelope zeros for width is clearly laid out. It is honest about the prior work of Johnson et al. and Tiede et al., and the application to real EHT data is a genuine attempt. The figures are helpful, and I spot-checked the thin-ring and thick-ring expressions without finding an error.\n\nThe soft spot is load-bearing, though. In Section V they fit the M87 visibility along two fixed axes (left-right and up-down) and report B = 0.48. But the model amplitude depends on cos²(ϕu − ϕ0). Section III.A carefully rotates the ring so the brightness maximum sits on an axis before taking sections; Section V does not do that, and it never fits or reports ϕ0. So the left-right cut is not necessarily the pure J0 term, the up-down cut is only a projection of the asymmetry, and the quoted B is degenerate with the orientation. That is a real flaw, not a cosmetic one. The fit also uses digitized points with no uncertainties and hand-picked directions, so no error bars follow. The simulation tests are self-consistency checks on rings from the same model family, which flatters the method.\n\nThe orientation degeneracy is the main thing to fix. A joint fit of the two perpendicular cuts with ϕ0 as a free parameter, plus some sensitivity checks, would probably resolve it. As it stands, the M87 asymmetry claim is not reliable, though the analytic model itself remains a useful tool for visibility-domain work.\n\nThis paper is for VLBI people who want a quick analytic estimate of ring parameters, especially for future space interferometry. It deserves a serious referee because the derivations are correct and the fix is tractable. I would send it to review but with an explicit warning about the position-angle degeneracy and the lack of uncertainty estimates.","headline":"Solid visibility math, but the M87 asymmetry measurement is not supported because the fit ignores the model's own position angle.","tokens_in":656,"tokens_out":565,"would_cite":false,"duration_ms":52741,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The radius, width, and angular asymmetry of a black-hole ring can be recovered from the zeros and envelope of its VLBI visibility function.","keywords":["black hole shadow","visibility function","VLBI","ring parameters","Event Horizon Telescope","M87*","Bessel function zeros","asymmetry parameter"],"falsifier":"Take a VLBI dataset synthesized from a simulated black-hole image with known radius, width, and asymmetry but with ellipticity or non-Gaussian radial brightness; if the zero-based recipe returns a radius or B that differs from the known input by more than the scatter in the paper's own Table I, the central circular-Gaussian assumption fails.","tokens_in":13464,"feed_emoji":"🔭","tokens_out":9571,"duration_ms":97029,"temperature":0.7,"pith_summary":"This paper proposes a direct analytical route from VLBI visibility measurements to the basic parameters of a ring-like black-hole image: radius, width, and angular asymmetry. The key idea is that each parameter leaves a distinctive fingerprint in the visibility function—the ring radius sets the positions of the visibility minima, the ring width sets the first zero of the oscillation envelope, and the asymmetry makes the visibility complex and different in two perpendicular directions. The authors test the recipe on synthetic rings of various widths and on two-ring configurations, then apply it to the Event Horizon Telescope's M87* data, obtaining a ring radius r≈23 μas and an asymmetry parameter B=0.48. If the method holds, it gives a fast, imaging-free way to extract black-hole shadow parameters from current and future interferometric data, and a step toward estimating black-hole spin.","feed_headline":"M87* ring radius and asymmetry are read from visibility zeros","feed_subtitle":"Zero positions fix the shadow radius; complex visibility fixes B=0.48; width stays unresolved.","key_machinery":"The load-bearing object is the analytic visibility function of a Gaussian asymmetric ring, obtained by substituting the brightness model $I(r,\\phi_r)=I_r(r)I_\\phi(\\phi_r)$ into the Fourier-transform definition of the visibility and integrating over the polar angle. That angular integration collapses the brightness asymmetry $I_\\phi(\\phi_r)=(1-B\\sin^2((\\phi_r-\\phi_0)/2))^n$ into the compact expression $V(u,\\phi_u)=\\pi\\int[(2-B)J_0(2\\pi u r)-iB J_1(2\\pi u r)\\cos(\\phi_u-\\phi_0)]I_r(r)r\\,dr$. The Bessel zeros $j_{0,n}$ then mark the visibility minima that fix the radius; the approximate large-baseline form shows a low-frequency envelope $\\sin(\\pi u w)$ whose first zero gives the width $w\\approx 1/u$; and the imaginary part proportional to $B$ carries the asymmetry. This reduction turns a two-dimensional imaging problem into a one-dimensional curve-fitting problem.","core_discovery":"On its own terms, the paper's discovery is that a black-hole shadow's ring can be parameterized without image reconstruction: the radius sits in the zeros of the visibility function, the width sits in the zeros of its envelope, and the asymmetry sits in the complex part of the visibility. The authors show this for an asymmetric Gaussian ring by reducing the two-dimensional Fourier integral to a one-dimensional combination of $J_0$ and $J_1$ Bessel functions, give exact expressions for thin, thick, and parabolic rings, and verify the zero-based recipe on synthetic rings and on two-ring configurations. Applied to the EHT observations of M87*, the recipe yields $r\\approx 23\\,\\mu\\text{as}$ and $B=0.48$, with the width left undetermined because the relevant envelope zero falls outside the available baseline coverage.","pith_inferences":["Editorial inference: the $B=0.48$ value should be read as the asymmetry of the best-fitting circular Gaussian model, not as a direct spin measurement; connecting $B$ to spin and inclination requires comparison with general-relativistic radiative-transfer simulations.","Editorial inference: the recipe's reliance on a single circular radius means it will be biased by ellipticity or off-center structure; checking whether visibility zero positions vary with position angle would be a cheap test of that assumption on real data.","Editorial inference: the same zero-and-envelope logic could be applied to closure amplitudes instead of complex visibilities, which would make the extraction less sensitive to phase-calibration errors in VLBI.","Editorial inference: if future space-VLBI baselines reach the first envelope zero for M87*, the currently unresolved ring width could be measured and the radius estimate cross-checked, potentially tightening black-hole mass estimates."],"forward_implications":["The shadow radius of a black hole can be estimated from the minima of the visibility amplitude alone, without reconstructing an image; in the paper's synthetic tests the first three minima put the radius of rings of width 0.1–5 μas at 20.0–20.2 μas against the input 20 μas.","The ring width can be read from the first zero of the visibility envelope whenever the baseline coverage reaches $u\\approx 1/w$; the paper's estimates for widths of 5.0, 7.5, and 10.0 μas are 3.41, 6.91, and 10.5 μas, with accuracy improving for wider rings.","For a bright thin ring superimposed on a thicker ring, both radii come from Bessel zeros and the thick ring's width is recovered to better than 7% by fitting the approximate Gaussian-ring visibility.","Applied to the EHT M87* observations, the method gives $r\\approx 23$ μas and, the paper reports, the first estimate of the asymmetry parameter $B=0.48$; the ring width could not be determined because the envelope zero lies beyond the available baselines.","Because the visibility is complex only when the ring is asymmetric, measuring two perpendicular cuts through the visibility plane is a direct way to detect and quantify angular brightness asymmetry."],"supporting_citations":[{"why":"Supplies the M87* EHT visibility and image data used for the application, including the ring diameter and width values the estimate is compared with.","marker":"[1]"},{"why":"Supplies the space-interferometer geometry whose (u,v) coverage underlies the numerical visibility observations for the two-ring cases.","marker":"[13]"},{"why":"Supplies the earlier analytic visibility-function treatment of an infinitely thin black-hole photon ring that this paper extends to asymmetric and thick rings.","marker":"[14]"},{"why":"Supplies the angular brightness-asymmetry model used to parametrize the ring's non-uniform brightness.","marker":"[15]"},{"why":"Supplies the standard definition of the complex visibility function as the Fourier transform of sky brightness, used as the starting integral.","marker":"[16]"},{"why":"Supplies the exact thick-homogeneous-ring and parabolic-ring visibility expressions used as comparison models in the width-estimation analysis.","marker":"[17]"},{"why":"Supplies the method for locating the brightness maxima of successive photon rings for non-zero spin, used to set up the two-thin-ring configurations.","marker":"[18]"}],"fun_headline_variants":["Black hole ring radius from visibility zeros","M87* shadow radius and asymmetry from visibility zeros","Visibility zeros reveal black hole ring size and shape","New method extracts black hole ring parameters without imaging","Ring radius from visibility zeros: M87* application"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the true image is a perfectly circular, centered ring whose radial profile is Gaussian and whose angular asymmetry has the fixed form of Eq. (2); if M87*'s actual brightness is elliptical, has jet emission, or falls off non-Gaussianly, the fitted radius and B will be biased.","fun_headline_variants_meta":{"raw":{"variants":["Black hole ring radius from visibility zeros","M87* shadow radius and asymmetry from visibility zeros","Visibility zeros reveal black hole ring size and shape","New method extracts black hole ring parameters without imaging","Ring radius from visibility zeros: M87* application"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1213,"prompt_tokens":828,"completion_tokens":385,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":444,"completion_tokens_details":{"reasoning_tokens":314}},"tokens_in":444,"tokens_out":385,"duration_ms":4019,"temperature":1.0,"reasoning_tokens":314,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T10:08:41.948902+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a VLBI dataset synthesized from a simulated black-hole image with known radius, width, and asymmetry but with ellipticity or non-Gaussian radial brightness; if the zero-based recipe returns a radius or B that differs from the known input by more than the scatter in the paper's own Table I, the central circular-Gaussian assumption fails.","supporting_citations":[{"cited_title":"On Optimal Geometry for Space Interferometers","cited_arxiv_id":"2305.19072","evidence_quote":"Supplies the space-interferometer geometry whose (u,v) coverage underlies the numerical visibility observations for the two-ring cases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the standard definition of the complex visibility function as the Fourier transform of sky brightness, used as the starting integral."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the exact thick-homogeneous-ring and parabolic-ring visibility expressions used as comparison models in the width-estimation analysis."},{"cited_title":"Flares and their echoes can help distinguish photon rings from black holes with space-Earth very long baseline interferometry","cited_arxiv_id":"2203.00577","evidence_quote":"Supplies the method for locating the brightness maxima of successive photon rings for non-zero spin, used to set up the two-thin-ring configurations."}],"review_version":1}