{"id":"1151d0a0-eb93-479a-88d2-4b615ec87f52","arxiv_id":"2507.23161","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A thin-shell model shows that radio-emitting outflows from tidal disruption events expand with a t^(2/3) law under black hole gravity, allowing SMBH masses to be estimated from observed shell radius and velocity.","lead":"This paper models the radio-emitting gas shell produced when a star is torn apart by a supermassive black hole, and shows the shell's expansion speed and size encode the black hole's mass. The authors propose monitoring such events at millimeter wavelengths to weigh black holes within months, an independent route to masses usually measured by other means.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed universality of K=1/2 in the mass estimator fails at ambient density slope s=2: the governing ODEs then admit a one-parameter family of t^(2/3) solutions with K set by initial conditions, so M_bh is not determinable from R and V alone.","rationale":"The paper's central claim is that a gravity-dominated, adiabatic thin shell follows R∝t^(2/3) with a known coefficient, giving an independent SMBH mass estimate. The weak point is not the finite-thickness K formula per se: re-deriving the coefficient balance from Eqs. (19)–(20) reproduces Eq. (23) and the thin-shell limit K=1/2 for s≠2. The genuine problem is the s=2 degeneracy. At s=2 the ambient mass grows linearly with R, and all terms in the momentum and energy equations share the same time dependence, so the coefficient comparison no longer fixes the ratio G M_bh/R*V*^2. The exact power-law family R=C t^(2/3) exists for arbitrary C, with the thermal energy adjusting to match C and M_bh; the constant C is selected by initial conditions. This contradicts the paper's advertised independence of the ambient density slope and undermines the mass estimator for an isothermal ambient, a realistic galactic-center profile. The reader's weakest_assumption focused on geometry/clumpiness; the s=2 degeneracy is a separate, more algebraic threat to Eq. (23). The verdict should remain CONDITIONAL: the paper can be made correct by restricting the universality claim to s≠2, quantifying the s≈2 convergence time, and testing the estimator against numerical integrations, but as written the central claim is overbroad.","tokens_in":10140,"tokens_out":28998,"duration_ms":320554,"concrete_test":"Integrate Eqs. (19)–(20) numerically for s=2 with fixed M_bh and two different initial conditions (e.g., initial shell velocity 0.8 and 1.2 times the local escape speed, with the same initial radius and zero or equal thermal energy). At t/t0=100, extract K_i = G M_bh/(R V^2) in each run. If K_1 and K_2 differ by more than ~20%, the universal mass estimator fails for isothermal ambients; repeat for s=1.9 to check whether the convergence to the claimed K is fast enough to be relevant at observed ages.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In §III.A, substituting R=R*(t/t*)^m into Eqs. (19)–(20) and equating coefficients gives m=2/3 and, after eliminating the thermal-energy amplitude, Eq. (23) with the stated K. The weakly secured point is the coefficient comparison for the ambient density slope s=2 (an isothermal profile, within the paper's range s<5/2). For s=2, M_am ∝ R and the momentum balance is R^2 Rddot + R Rdot^2 + G M_bh = α ET/A, while the energy equation reduces to a first integral. Substituting a power law R=C t^(2/3) shows that the coefficient equation is identically satisfied for any C: with x=G M_bh/C^3, the equation becomes x+2/9 = x+2/9. Thus K=(9/4)x is not determined by the dynamics; the ODE system has an exact one-parameter family R=C t^(2/3), ET=(A/α)[(2/9)C^3+G M_bh], and initial conditions select C (hence K). The value K→1/2 is recovered only by taking a limit from s≠2 or by imposing that the ejecta is exactly marginally bound, neither of which follows from the stated equations. Consequently the central claim that SMBH mass can be inferred as M_bh=K R* V*^2/G 'irrespective of the ambient density profile' is false at s=2, and for s near 2 the observable K may retain a long-lived dependence on initial conditions, making the estimator unreliable in exactly the regime the paper targets.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter presents a one-zone model for the adiabatic expansion of a thin, quasi-spherical radio-emitting shell in a non-jetted tidal disruption event (TDE), driven by a disk wind and decelerated by the gravity of the central supermassive black hole (SMBH). The authors combine momentum conservation (Eq. 14) and energy conservation (Eq. 18) for a shell that sweeps up a power-law ambient medium ρ_am ∝ r^{-s}, and seek asymptotic power-law solutions R ∝ t^m. Their central result is that when SMBH gravity dominates or is comparable to thermal pressure, m = 2/3 regardless of s, and the SMBH mass is given by M_bh = K R_* V_*^2 / G (Eq. 23), with K → 1/2 in the thin-shell limit ξ/χ ≪ 1. They also recover the Sedov-Taylor solution when gravity is negligible (R ∝ t^{2/(5-s)}) and a momentum-driven snowplow solution (R ∝ t^{1/(4-s)}). The paper applies the mass estimator to six non-jetted TDEs, introduces a discriminant D_c = M_•/M_c to decide which regime applies, and argues that the estimator overestimates M_• for AT2019dsg because that event lies in the Sedov-Taylor regime.","tokens_in":10414,"tokens_out":24325,"duration_ms":256633,"significance":"If valid, the proposed mass estimator would provide a genuinely independent, radio-only route to SMBH masses in non-jetted TDEs, complementing reverberation mapping and scaling relations, and it would make a concrete, falsifiable prediction for high-frequency radio monitoring campaigns. The analytic structure is a strength: the governing equations are physically transparent, the Sedov-Taylor and snowplow limits are recovered, and the m = 2/3 solution has a clear dimensional interpretation. The authors also cite a public Supplemental Material archive for the detailed algebra, which is good scholarly practice. However, the central claim of universality fails at the special ambient slope s = 2, where the coefficient K is not fixed by the dynamics, and the observational validation contains a selection loop through D_c. These issues do not invalidate the model for s ≠ 2, but they require correction before the Letter can be accepted.","major_comments":[{"comment":"The coefficient matching that fixes m = 2/3 and K fails at ambient density slope s = 2, which lies inside the stated range s < 5/2. For s = 2, M_am ∝ R, and substituting R = C t^{2/3} into Eqs. (19)-(20) yields a one-parameter family of exact solutions: with x ≡ G M_bh / C^3, the energy equation determines E_T = (A/α)(2 C^3 / 9 + G M_bh), and the momentum equation is then satisfied identically for any x. Thus K = (9/4)x is not determined by R_* and V_* alone. The value K → 1/2 in the thin-shell limit is recovered only as the s → 2 limit of the non-degenerate formula, not as a property of the s = 2 solutions, where K can take any positive value set by initial conditions. The abstract's claim that the scaling and the mass estimator hold 'irrespective of the ambient density profile' is therefore false at s = 2, and Eq. (23) is not a mass estimator in that regime. The paper should either exclude s = 2 explicitly, supply the additional condition that selects K, or state the resulting uncertainty for s near 2.","section":"III.A, Eq. (23)"},{"comment":"The statement that D_c 'depends only on the outflow radius and expansion speed' is not consistent with Eq. (25), which contains M_• explicitly. Consequently, deciding whether the m = 2/3 solution applies (D_c > 1) requires an independent estimate of M_•, the very quantity Eq. (23) is meant to provide. The paper explains the order-of-magnitude overestimate for AT2019dsg by assigning it to the Sedov-Taylor regime with D_c ≪ 1, but that classification uses the conventional masses (X-ray and M-σ) that Eq. (23) is supposed to replace; if M_• were instead the value predicted by Eq. (23), D_c would be greater than unity and the same data would be classified as gravity-dominated. This selection loop weakens the observational validation. It should be addressed, for example, by computing D_c from the Eq. (23) mass and checking self-consistency, or by treating the conventional mass as a prior in a Bayesian comparison.","section":"IV, Eq. (25)"}],"minor_comments":[{"comment":"Equation (24) invokes M_am,* , E_T,* , and T_p,* that are defined only in the Supplemental Material; please provide one-line definitions in the main text so that Eq. (24) is self-contained.","section":"III.B, Eq. (24)"},{"comment":"The abstract's prediction that 10-100 GHz monitoring can yield masses 'within months of disruption' is not derived in the main text; connect it to the normalization time t_* in Eq. (21) and to the expected R_*, V_* for events with D_c > 1.","section":"Abstract and V"},{"comment":"The parameters η and δ are introduced in Eqs. (6)-(7) but play no role in the asymptotic analysis or in the mass estimator; state explicitly that they only affect the early-time behavior, or remove them from the main-text definitions.","section":"II.A, Eqs. (6)-(7)"},{"comment":"Reference [25] appears both as the general 'Cendes and et al.' citation and as the detailed 'Cendes, Alexander, Berger et al.' entry; consolidate the bibliography.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the s = 2 degeneracy is a genuine mathematical defect in the coefficient matching and should be the main focus of the revision. The D_c-based validation also needs restructuring to avoid the circular use of conventional masses. I did not have access to the Zenodo Supplemental Material; the revision should ensure that the supplement's derivation does not silently divide by (2 - s) or otherwise hide the exceptional case. The paper is otherwise a useful contribution for the s ≠ 2 regime."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. The paper builds an explicit thin-shell model for radio-emitting TDE outflows and shows that in the gravity-dominated regime asymptotic solutions scale as R ∝ t^(2/3), giving M_bh = K R V^2/G. That derivation is clean for s < 5/2, s ≠ 2, and the D_c criterion (equation 25) is a genuinely useful, scale-free way to decide between the gravitational and Sedov-Taylor regimes. The authors are also honest about the earlier t^(2/3) solutions in Sakashita (1974) and Pen (1994); the new content is the explicit shell equations, the regime classification, and the application to radio TDEs.\n\nThe main problem is in Section III.A. For s=2, the coefficient matching that fixes m and K fails. Substituting R = C t^(2/3) into equations (19)-(20), the momentum equation gives (2/9)C^3 + G M_bh = (χ(γ-1)/ξ) E_T / A, and the energy equation is then identically satisfied by the same relation. It contains no new information. Thus C, and hence K = (9/4) G M_bh / C^3, is selected by initial conditions, not by the dynamics. The claim that K → 1/2 \"regardless of s\" in the thin-shell limit is not correct at s=2; taking ξ/χ small does not remove the degeneracy. Since s=2 is within the stated range s < 5/2 and corresponds to an isothermal ambient profile, this is a real gap in the central mass estimator. For s near but not equal to 2, the memory of initial conditions can persist, making the estimator unreliable unless this case is excluded or an additional physical constraint is supplied.\n\nThe observational validation is thin: four of the six events are firmly in the Sedov-Taylor regime, no uncertainties are given for R* and V*, and the closest event to the gravity-dominated boundary (eRASStJ2344) matches an independent mass only at the 40% level. That is not enough to calibrate the method.\n\nBottom line: this is a motivated analytic framework with a useful diagnostic, but the central mass formula has a genuine degeneracy at s=2. I would send it to a serious referee because the derivation is explicit and the flaw is fixable by a caveat or by removing s=2 from the claimed universality. It should not be accepted as is. Good reading-group material for a discussion of coefficient matching and degeneracies in self-similar solutions.","headline":"Useful analytic shell framework for radio TDEs, but the mass estimator degenerates for ambient density index s=2 and needs a caveat before the claimed universality holds.","tokens_in":11085,"tokens_out":9181,"would_cite":false,"duration_ms":101604,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Radio-emitting shells around tidally disrupted stars can weigh their central supermassive black holes.","keywords":["tidal disruption events","supermassive black hole masses","radio-emitting outflows","thin shell dynamics","Sedov-Taylor solution","disk winds","t^(2/3) scaling","non-jetted TDEs"],"falsifier":"For a non-jetted TDE whose black hole mass is known independently, measure the radio shell radius and expansion velocity at three or more epochs. The model requires a clean $R(t)\\propto t^{2/3}$ phase when $D_c>1$, so an observed power-law index differing from $2/3$, or an inferred $M_{\\rm bh}=K R_* V_*^2/G$ disagreeing with the independent mass by more than the stated uncertainties, would falsify the central claim.","tokens_in":9816,"feed_emoji":"🕳️","tokens_out":7188,"duration_ms":76897,"temperature":0.7,"pith_summary":"The paper tries to establish that in non-jetted tidal disruption events, the radio-emitting shell driven outward by disk winds is decelerated strongly enough by the black hole's gravity that its late-time expansion follows a universal $t^{2/3}$ scaling, independent of the ambient density profile. If correct, measuring the shell radius and expansion velocity from radio spectra gives the supermassive black hole mass directly through $M_{\\rm bh}=K R_* V_*^2/G$, with $K\\simeq1/2$ for a thin shell. This would matter because it supplies an independent, purely radio-based mass estimator at a time when most non-jetted TDEs lack jet-related constraints. The paper also recovers the classical Sedov-Taylor expansion when gravity is negligible, and introduces a criterion $D_c$ that tells which regime applies to a given event.","feed_headline":"Radio shells weigh black holes via a universal 2/3 power law","feed_subtitle":"When gravity dominates, radio size and speed alone give the black hole mass in months.","key_machinery":"The central machinery is a one-dimensional thin-shell model combining momentum conservation (equation 14), including ram-pressure injection, gravitational force, and thermal pressure, with thermal-energy evolution (equation 18). The load-bearing step is taking the asymptotic limit $R\\gg r_0$, where $M(t)\\approx M_{\\rm am}(t)$, inserting the power-law ansatz $R(t)=R_*(t/t_*)^m$, and equating the time exponents: this fixes $m=2/3$ when gravity dominates and fixes the coefficient $K$ in $M_{\\rm bh}=K R_* V_*^2/G$. The companion object is the critical mass $M_c$ (equation 24) and its dimensionless form $D_c=M/M_c$, a scale-free diagnostic that selects whether the gravity-dominated, Sedov-Taylor, or transition solution applies.","core_discovery":"The governing equations are momentum and energy conservation for a quasi-spherical thin shell (equations 14 and 18) with mass from the disk wind, the ambient medium, and an initial shell. At radii $R\\gg r_0$ the shell mass is dominated by swept-up ambient matter, $M\\approx M_{\\rm am}$, and assuming a power-law radius $R(t)=R_*(t/t_*)^m$, coefficient matching yields exactly $m=2/3$ when gravity is comparable to or stronger than thermal pressure, independent of the ambient density slope $s$. The same matching gives equation (23), $M_{\\rm bh}=K R_* V_*^2/G$, with $K$ approaching $1/2$ in the thin-shell limit $\\xi/\\chi\\to 0$. In the opposite limit, $M_{\\rm bh}\\ll M_c$, the solution returns to the density-dependent Sedov-Taylor index $2/(5-s)$, and when both gravity and thermal pressure are negligible it reduces to the momentum-driven snowplow solution $t^{1/(4-s)}$. The paper asserts that the gravity-dominated $m=2/3$ solution is the correct description for events with the dimensionless ratio $D_c\\gtrsim 1$.","pith_inferences":["The $D_c$ criterion uses only radius and speed, so it could be exported to other transient accretion outflows, such as X-ray binary winds and AGN outflows, to decide when gravity-dominated dynamics matter; the paper gestures at this but does not develop it.","If the $m=2/3$ scaling is confirmed, the same dimensional relation implies a link between the radio-measured expansion rate and the tidal radius, which could be used to constrain the spin-dependent Hills mass for the most massive events.","The model assumes spherical symmetry, so radio imaging of a candidate TDE that resolves the shell morphology would be a direct check: an anisotropic or clumpy outflow would break the $t^{2/3}$ prediction."],"forward_implications":["Non-jetted TDEs with $D_c>1$ give black hole masses purely from radio observables, with no reliance on X-ray scaling relations or reverberation mapping.","The $t^{2/3}$ law is independent of the ambient density slope $s$, so finding the same expansion index across events in different environments would support the gravity-dominated picture.","Because the synchrotron peak frequency drops with source size, the model motivates 10-100 GHz follow-up with ALMA and ngVLA to catch the small radii where $D_c>1$ for lower-mass black holes, yielding masses within months of disruption.","Application to the current sample explains why the formula overestimates AT2019dsg's mass by an order of magnitude ($D_c\\simeq0.06$, deep Sedov-Taylor regime) while it is only about 40% above the independent estimate for eRASSt J2344 ($D_c\\simeq0.84$)."],"supporting_citations":[{"why":"Supplies the previous momentum-conservation shell model that this paper extends by adding energy conservation, thermal pressure, and SMBH gravity.","marker":"[10]"},{"why":"Provides the $n=5/3$ mass fallback rate that sets the disk-wind injection profile used in the shell equations.","marker":"[2]"},{"why":"The Sedov solution is the gravity-free baseline whose density-dependent index $2/(5-s)$ the gravity-dominated $2/3$ result is contrasted with.","marker":"[18]"},{"why":"Together with [18], defines the classical Sedov-Taylor regime recovered when black hole gravity is negligible.","marker":"[19]"},{"why":"Provides the method for inferring shell radius and expansion velocity from radio spectra, the two observables entering the mass formula.","marker":"[27]"},{"why":"Supplies the AT2019dsg radio measurements used for the concrete mass estimate in the discussion.","marker":"[25]"},{"why":"Gives the X-ray derived mass for AT2019dsg that the new formula overestimates, motivating the regime criterion.","marker":"[28]"},{"why":"Provides the $M$-$\\sigma$ mass estimate for AT2019dsg used in comparing the radio-shell mass with conventional estimators.","marker":"[14]"},{"why":"Defines the Hills mass limit relevant to the claim that the method is best suited to massive spinning black holes that can still disrupt stars.","marker":"[1]"}],"fun_headline_variants":["Black hole masses from radio shell expansion","Weighing black holes with radio shells in months","Universal 2/3 law ties radio shells to black hole mass","Radio observations reveal black hole mass via shell dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the outflow is a smooth, roughly spherical, adiabatic thin shell whose swept-up surrounding gas dominates its mass at late times; if the outflow is clumpy, collimated, or strongly radiative, both the $t^{2/3}$ law and the mass formula stop applying.","fun_headline_variants_meta":{"raw":{"variants":["Black hole masses from radio shell expansion","Weighing black holes with radio shells in months","Universal 2/3 law ties radio shells to black hole mass","Radio observations reveal black hole mass via shell dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000484,"raw_usage":{"total_tokens":2426,"prompt_tokens":1015,"completion_tokens":1411,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":1349}},"tokens_in":631,"tokens_out":1411,"duration_ms":12664,"temperature":1.0,"reasoning_tokens":1349,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:02:37.910724+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a non-jetted TDE whose black hole mass is known independently, measure the radio shell radius and expansion velocity at three or more epochs. The model requires a clean $R(t)\\propto t^{2/3}$ phase when $D_c>1$, so an observed power-law index differing from $2/3$, or an inferred $M_{\\rm bh}=K R_* V_*^2/G$ disagreeing with the independent mass by more than the stated uncertainties, would falsify the central claim.","supporting_citations":[{"cited_title":"Krolik, T","cited_arxiv_id":null,"evidence_quote":"Supplies the previous momentum-conservation shell model that this paper extends by adding energy conservation, thermal pressure, and SMBH gravity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the $n=5/3$ mass fallback rate that sets the disk-wind injection profile used in the shell equations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Sedov solution is the gravity-free baseline whose density-dependent index $2/(5-s)$ the gravity-dominated $2/3$ result is contrasted with."},{"cited_title":"Taylor, Proceedings of the Royal Society of London Series A 201, 159 (1950)","cited_arxiv_id":null,"evidence_quote":"Together with [18], defines the classical Sedov-Taylor regime recovered when black hole gravity is negligible."},{"cited_title":"A General Class of Self-similar Self-gravitating Fluids","cited_arxiv_id":"astro-ph/9402039","evidence_quote":"Supplies the AT2019dsg radio measurements used for the concrete mass estimate in the discussion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the X-ray derived mass for AT2019dsg that the new formula overestimates, motivating the regime criterion."},{"cited_title":"Cendes and et al., The Astrophysical Journal 919, 127 (2021)","cited_arxiv_id":null,"evidence_quote":"Provides the $M$-$\\sigma$ mass estimate for AT2019dsg used in comparing the radio-shell mass with conventional estimators."},{"cited_title":"We introduce the initial shell mass as ∆m = 4πδρ ej,0r3 0, (7) where δ is a parameter to decide the amount of the initial shell mass","cited_arxiv_id":null,"evidence_quote":"Defines the Hills mass limit relevant to the claim that the method is best suited to massive spinning black holes that can still disrupt stars."}],"review_version":1}