{"id":"d10c5b09-748f-45bd-8aca-b15a557f8bab","arxiv_id":"2507.01273","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Simulated tidal disruption outflows crashing into surrounding gas produce radio emission matching several observed TDE radio flares.","lead":"This paper simulates the collision between outflows from a tidal disruption event and surrounding gas, and finds that the resulting shocks produce radio emission resembling what telescopes see. It supports the idea that radio flares from tidal disruption events come from outflow-gas collisions rather than jets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim depends on an untested smooth, spherical CNM shell (Sec. 2.1): mass chosen to make collisions significant, slope from AT2019dsg; clumpy or thinner CNM may break fitted spectra, so sufficiency is not yet established.","rationale":"The reader's weakest_assumption pinpoints the CNM modeling, and my reading agrees: this is the most load-bearing input because the collision strength, shock properties, and the synthetic peak-frequency evolution all scale with the assumed density profile. The mass 0.1 M_sun is selected to guarantee a significant interaction, and the slope n=-1.7 comes from AT2019dsg, which is itself one of the two events with good spectral fits. This creates a circularity risk that is not addressed by any test with realistic CNM structure. I also considered whether the strongest concern might instead be the per-event choice of epsilon_e and epsilon_B in Fig. 5 (the bottom-row fits use values different from the 0.1/0.01 adopted in the conclusions), but the two good fits use the same microphysics, so the residual issue is the environmental assumption. Independent support is present: the paper provides open simulation data, a resolution study showing peak frequency and flux converged to within a factor of two (Appendix A), and a benchmark of the ray-tracing scheme (Appendix B). These strengthen the numerical execution, but they do not test the physical CNM input. A clumpy-CNM simulation is the most direct way to settle whether the smooth shell is essential; if the fits survive, the mechanism's sufficiency is substantially more credible, and the conditional verdict can be upgraded.","tokens_in":22462,"tokens_out":23678,"duration_ms":246869,"concrete_test":"Run the beta=5, e=0.95 radio simulation with a clumpy CNM generated from the same mean profile (n=-1.7, M=0.1 M_sun) but lognormal density fluctuations with sigma_ln_rho ~ 1.5, then ray-trace spectra with epsilon_e=0.1, epsilon_B=0.01, p=2.5 and compare nu_peak(t) and the spectral fits to AT2020vwl and AT2019dsg (Fig. 5, top row). If peak flux or nu_peak shifts by more than ~0.3 dex relative to the smooth-shell run, or the fitted decay slope changes by more than ~0.2, the smooth-spherical CNM assumption is load-bearing for the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim of sufficiency rests on the CNM model in Sec. 2.1: a smooth, static, spherically symmetric shell with rho = rho0 (r/r0)^n, n = -1.7, and M_CNM = 0.1 M_sun. The mass is explicitly chosen 'to ensure significant collisions', and the slope is adopted from AT2019dsg, one of the events later fitted in Fig. 5. Eq. (21) and Sec. 3.6 show that the synthetic peak frequency scales with the CNM density, and Sec. 4 concedes that a steeper/thinner or clumpy CNM could explain the observed peak-flux variations that the model fails to reproduce. Because no simulation with a clumpy or truncated CNM is presented, the good fits to AT2020vwl and AT2019dsg may be an artifact of this smooth-shell assumption rather than robust evidence for the mechanism; the degeneracy between M_CNM and beta noted in Sec. 3.3 further means that shock properties alone cannot discriminate the model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses general relativistic smoothed particle hydrodynamics (Phantom) simulations of tidal disruption events (TDEs) with penetration factors beta = 1 and 5 on eccentric orbits (e = 0.95) and injects the resulting super-Eddington outflows into a static, spherically symmetric circumnuclear medium (CNM) shell. The authors identify shocked material via an entropy threshold, extract shock radius, velocity, and energy, compare these with equipartition-model inferences for six radio TDEs, and ray-trace synchrotron emission from the shocked particles to produce synthetic images and spectra. The central claim is that the outflow-CNM collision produces shocks as early as ~10 days with radius ~1e17 cm, velocity ~0.15c, and energy ~1e51 erg, and that the synthetic spectra show a continuously decaying peak frequency that matches prompt radio TDE observations, supporting the hypothesis that such collisions, rather than jets, can explain at least some prompt radio flares from TDEs.","tokens_in":22744,"tokens_out":5736,"duration_ms":60969,"significance":"If the central result holds, this is a significant step: it provides a self-consistent hydrodynamic demonstration that quasi-spherical super-Eddington TDE outflows colliding with a CNM shell can produce shocks with properties consistent with prompt radio observations, and that the synchrotron peak frequency decays in time as observed. The paper includes a resolution study (Appendix A), a benchmark of the ray-tracing scheme against analytic synchrotron spectra (Appendix B), and public data releases on Zenodo, which strengthen confidence in the numerical machinery. The significance is tempered by the ad hoc construction of the CNM model (Section 2.1), the calibration of epsilon_e and epsilon_B to the same observations (Sections 2.6.3 and 4), and the admitted discrepancies in peak-flux evolution and spectral shape for several events (Section 3.5).","major_comments":[{"comment":"The CNM density profile is the main external input and is not independently constrained: n = -1.7 is adopted from AT2019dsg (Cendes et al. 2021) and M_CNM = 0.1 M_sun is chosen to 'ensure significant collisions' (Section 2.1). Because the model is subsequently compared with AT2019dsg in Figure 5 (top right), the agreement for that event is at least partly by construction. The paper varies M_CNM only for beta = 5 (0.01 and 1 M_sun; Figures 3 and 6) and does not vary n or test a clumpy or truncated CNM, although Section 3.6 and Section 4 argue that the peak frequency and flux evolution depend sensitively on the CNM structure. Without a systematic exploration or an independent constraint on the CNM, the sufficiency claim in the abstract is not established.","section":"§2.1, §3.5"},{"comment":"The spectral comparison is made after selecting epsilon_e and epsilon_B from a grid (Section 2.6.3) and choosing a viewing direction per event (Figure 5 caption) to best match the observed spectra. The reported 'best fit' values (epsilon_e ~ 0.1, epsilon_B ~ 0.01, Section 4) are therefore calibrated to the same data that the model is claimed to reproduce; no goodness-of-fit or uncertainty is provided, and three of the seven events (AT2020opy, AT2019azh, CNSS J0019) are not fitted at all. This weakens the inference that the synthetic spectra independently support the model over the jet hypothesis.","section":"§2.6.3, §3.5, §4"},{"comment":"The abstract states that the synthetic spectra 'match prompt radio TDE observations,' but Figure 5 shows that the spectral shapes are not reproduced for eRASSt J2344 and ASAASN-14li (bottom row) and that the peak-flux evolution differs from observations for all events (Section 3.5). The paper's own Section 4 lists several possible explanations for these discrepancies, including a steeper or clumpy CNM and incomplete cooling. The evidence supports a more limited claim: the model reproduces the decaying peak-frequency trend and, for two of four fitted events, the early-time spectral shape at order-of-magnitude level.","section":"Abstract, §3.5, §4"}],"minor_comments":[{"comment":"There is a typo in the sentence beginning 'consistent with Spectrum 2 shown in Figure B.4. s The scaling of νa ...' where an unattached 's' appears between the two sentences.","section":"§3.6"},{"comment":"The notation q is used both for the exponent in the scaling relation and as a variable in the expression q ≡ 2/(p+4)n − (p+2)/(p+4)m + 2/(p+4); please clarify in the text that q is the resulting exponent and avoid reusing q later without definition.","section":"Eq. (21) and following text"},{"comment":"Equation (3) is only valid for n > -3; since the paper uses n = -1.7 this is fine, but the general expression should state this restriction to avoid a formally negative density for steeper profiles.","section":"Eq. (3)"},{"comment":"The distinction between the temperature used for interpreting internal energy (Section 2.2) and the electron energy distribution assumed in the synchrotron calculation (Section 2.6.1) is clear in the text, but it would help to restate in Section 2.6 that the large post-shock temperatures shown in Figure 2 are not used directly as electron temperatures in the radiative transfer.","section":"§2.2 and §2.6.1"},{"comment":"The resolution study is presented only at t = 0.3 yr; the statement that the synthetic images are 'not critically resolution dependent' should be qualified as applying to the early-time epoch shown in Figures A.1 and A.2, since the text notes that the flux below ~20 GHz is not yet converged.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good candidate for a Letter if the authors either add a minimal robustness test of the CNM assumptions (e.g., varying n or adding a clumpy component) or moderate the sufficiency claim to 'consistent with'. The circularity concern regarding AT2019dsg and the calibration of epsilon_e/epsilon_B should be addressed explicitly in the text. I would not recommend rejection; the numerical machinery and the basic shock properties are solid and the paper is clearly of interest to the TDE radio community."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the method: self-consistent GR-SPH TDE outflows injected into a CNM shell and followed with Phantom, then ray-traced for synchrotron images and spectra. That is a real step beyond the prior analytical estimates, and the main hydrodynamic numbers come directly from the simulations: shock radius ~1e17 cm, velocity ~0.15c, total energy ~1e51 erg, with about 10% of the outflow energy available in the shocked region. Matching those order-of-magnitude quantities to observed TDE radio properties is a solid result, and the paper is honest that this supports the outflow-CNM hypothesis rather than proving it. The resolution study and the benchmark spectra in the appendices are genuine checks, and the synthetic images showing an aspherical but reflection-symmetric emitter are a nice, concrete prediction.\n\nThe soft spots are real but mostly not disqualifying. The CNM is a smooth, static, spherical shell with n=-1.7 taken from AT2019dsg and M_CNM=0.1 Msun chosen to make collisions significant, and then the model is compared to AT2019dsg's spectra. That is partly circular for that one event, and the paper says so in spirit by noting the degeneracy between beta and M_CNM and the sensitivity of peak frequency to CNM density. The spectral fitting also chooses eps_e=0.1, eps_B=0.01, p=2.5 to best match the same data, so the good fits to AT2020vwl and AT2019dsg are demonstrations of consistency, not predictions. And it is openly acknowledged that three of seven events cannot be fitted and that the peak flux barely evolves in the models while observations vary strongly. The listed explanations for that---clumpy CNM, steeper density profiles, cooling, time-varying eps_B---are plausible, but none is simulated here. The sufficiency claim is therefore conditional on a smooth-shell CNM, which is exactly the weakest link.\n\nAll that said, the paper's central hydrodynamic result holds up: if you put a super-Eddington outflow into a dense-enough CNM shell, you naturally get prompt radio shocks in the right energy and radius range. That is a useful and previously missing data point. The paper is well-cited, the simulations are archived on Zenodo, and the parameter choices and limitations are laid out clearly rather than hidden.\n\nWho should read this: anyone working on TDE radio emission or outflow-CNM interaction models. It deserves a serious referee. My recommendation: send it to review, and push the authors to clarify which claims are simulation results versus fitted consistency checks, and to add at least one clumpy or truncated-CNM test to see whether the spectral fits survive. But if you need a quick answer now: the paper is worth engaging with, and I would cite it for the shock-energy result while treating the spectral matches with caution.","headline":"First simulations of TDE outflow-CNM collisions show the mechanism can produce the right prompt-radio energies and peak-frequency decays, but the smooth-shell CNM setup and fitted microphysics mean it is still a proof-of-mechanism, not a predictive model.","tokens_in":23319,"tokens_out":1429,"would_cite":true,"duration_ms":115264,"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":"Prompt radio flares from tidal disruption events can be produced when the event's quasi-spherical super-Eddington outflow slams into surrounding circumnuclear gas, with no jet required.","keywords":["tidal disruption events","radio transients","hydrodynamical simulations","synchrotron emissions","ray tracing","synthetic observations","supermassive black holes"],"falsifier":"A very long baseline interferometry image of a prompt radio TDE that resolves the emitting region as a narrow, one-sided jet feature rather than a flattened shell with mirror symmetry about the stellar orbital plane would contradict the model; alternatively, a TDE with an independently measured CNM density whose radio peak frequency does not decay as a power law with time would break the $\\nu_{\\rm peak}\\simeq\\nu_{\\rm a}$ scaling.","tokens_in":22243,"feed_emoji":"📡","tokens_out":8957,"duration_ms":100959,"temperature":0.7,"pith_summary":"Tidal disruption events—stars torn apart by supermassive black holes—often produce radio flares within weeks, but the origin of this radio emission has been debated: jets or outflows? This paper argues that the prompt radio flare requires no jet at all. It simulates the collision between the quasi-spherical, super-Eddington outflow produced during a TDE and a surrounding cloud of circumnuclear gas, and finds that shocks form within about ten days with radius, speed, and energy that match the values inferred from actual radio observations. Ray-traced synthetic spectra show a continuously decaying peak frequency, matching the prompt radio data, and the emitting region has a flattened, doughnut-like shape symmetric about the original stellar orbital plane. The authors conclude that outflow–gas collisions, not jets, can be the dominant source of synchrotron radio emission from TDEs.","feed_headline":"Outflow crashing into gas explains tidal disruption radio flares","feed_subtitle":"Simulations match observed shock size, speed, and decaying radio spectra, with no jet required.","key_machinery":"The central mechanism is the forward shock at the interface between the injected super-Eddington outflow and an initially static, spherically symmetric circumnuclear shell. Shocked gas is identified by a >30% entropy rise over its background value, and only that gas is treated as a synchrotron emitter. The paper's new physical identification is that the spectral peak is the self-absorption frequency, $\\nu_{\\rm peak}\\simeq\\nu_{\\rm a}$, with $\\nu_{\\rm a}\\propto \\rho^{2/(p+4)} B^{(p+2)/(p+4)} R^{2/(p+4)}$; combined with the radial expansion $R\\propto t^{\\delta}$ this yields the peak-frequency decay $\\nu_{\\rm peak}\\propto t^{q\\delta}$, reproducing the simulated $t^{-0.78}$ and connecting to Sedov–Taylor and snow-plow expansion phases.","core_discovery":"The paper claims that the prompt synchrotron radio flare of a tidal disruption event can be powered by the collision between the quasi-spherical, super-Eddington outflow produced during the disruption and a pre-existing circumnuclear gas cloud, with no jet required. In the simulations a strong forward shock appears as early as $\\sim 10$ days after disruption, and by three years the shocked shell has radius $\\approx 10^{17}$ cm, velocity $\\sim 0.15c$, and total energy $\\sim 10^{51}$ erg, in line with values inferred from observations. The synthetic spectra, produced by ray tracing the shocked gas, peak at the self-absorption frequency rather than at the synchrotron characteristic or cooling frequencies, and this peak frequency decays continuously, roughly as $t^{-0.78}$, matching the prompt-radio observations. The emitting region is aspherical with a flattened, ring-like morphology and reflection symmetry about the original stellar orbital plane, even though the outflow itself is quasi-spherical. The authors conclude that outflow–circumnuclear-medium collisions are a sufficient and likely dominant source of radio synchrotron emission from TDEs.","pith_inferences":["Beyond the paper: the $\\nu_{\\rm peak}\\simeq\\nu_{\\rm a}$ scaling predicts that combining a measured CNM density profile (from X-ray or free-free absorption) with a measured radio peak-frequency decay slope would pin down the expansion index $\\delta$ and the magnetic-field decay index in a single event.","Beyond the paper: if the CNM is clumpy, the smooth-shell simulations underestimate peak-flux variability; radio light curves with re-brightening episodes would be a natural test of clump-driven bow shocks.","Beyond the paper: the same external shocks should accelerate protons, but with only about 10% of the outflow energy in the shock, neutrino and gamma-ray predictions scaled to the full outflow energy are likely too optimistic.","Beyond the paper: the apparent asymmetry in the $\\beta=5$ outflow that the authors flag as possibly numerical could be tested by repeating the run with different noise seeds or resolution; if physical, it would break the exact mirror symmetry of the radio image."],"forward_implications":["Prompt radio TDE flares can be produced without jets: a quasi-spherical super-Eddington outflow colliding with a modest circumnuclear shell is sufficient.","The radio-emitting region should appear as a flattened, mirror-symmetric shell rather than a round blob, so the viewing angle relative to the star's orbital plane shapes the observed image.","Only about 10% of the outflow kinetic energy ends up in the shocked, radio-emitting gas, so radio-inferred energy budgets should not be equated with the total outflow energy.","The decaying peak frequency of prompt radio TDEs can be read as the self-absorption frequency of the shocked shell, giving a physical scaling rather than an empirical fit.","Because TDE outflows are common while jets are rare, outflow-driven radio emission may dominate the radio TDE population whenever enough circumnuclear gas is present."],"supporting_citations":[{"why":"Supplies the two eccentric-disruption simulations whose outflows are injected into the circumnuclear shell.","marker":"Hu et al. 2024"},{"why":"Establishes that these super-Eddington outflows arise from self-collision and accretion and can explain the optical TDE emission.","marker":"Price et al. 2024"},{"why":"Provides the AT2019dsg data and the CNM density power-law index n=-1.7 used in the shell setup, plus the epsilon_B comparison.","marker":"Cendes et al. 2021"},{"why":"Provides the prompt-radio TDE observations (including AT2020vwl) that the shock properties and synthetic spectra are compared against.","marker":"Goodwin et al. 2023a"},{"why":"Supplies the equipartition model used to convert observed spectral peaks into the radii, velocities, and energies shown in Figure 3.","marker":"Barniol Duran et al. 2013"},{"why":"Supplies the slow-cooling synchrotron spectrum used to ray-trace the synthetic images and spectra.","marker":"Sari et al. 1998"},{"why":"Provides the ASAASN-14li radio observations that the synthetic spectra are matched to.","marker":"Alexander et al. 2016"},{"why":"Provides the analytic snow-plow phase solution used to interpret the peak-frequency decay slopes.","marker":"Hayasaki & Yamazaki 2023"}],"fun_headline_variants":["No jet needed: outflow collision powers TDE radio flare","Prompt TDE radio flare explained by outflow–gas shock","TDE radio flare emerges from outflow–cloud collision","Shock within 10 days powers TDE radio flare","Simulations show TDE radio flare from outflow collision"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumed circumnuclear medium: a smooth, initially static, spherically symmetric shell with density $\\rho=\\rho_0(r/r_0)^{-1.7}$ and total mass $0.1\\,M_\\odot$; if the real gas is clumpy, thinner, or has a different density slope, the predicted shock properties and radio spectra could shift substantially.","fun_headline_variants_meta":{"raw":{"variants":["No jet needed: outflow collision powers TDE radio flare","Prompt TDE radio flare explained by outflow–gas shock","TDE radio flare emerges from outflow–cloud collision","Shock within 10 days powers TDE radio flare","Simulations show TDE radio flare from outflow collision"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0008,"raw_usage":{"total_tokens":3555,"prompt_tokens":1018,"completion_tokens":2537,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":2458}},"tokens_in":634,"tokens_out":2537,"duration_ms":21029,"temperature":1.0,"reasoning_tokens":2458,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:56:25.297348+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A very long baseline interferometry image of a prompt radio TDE that resolves the emitting region as a narrow, one-sided jet feature rather than a flattened shell with mirror symmetry about the stellar orbital plane would contradict the model; alternatively, a TDE with an independently measured CNM density whose radio peak frequency does not decay as a power law with time would break the $\\nu_{\\rm peak}\\simeq\\nu_{\\rm a}$ scaling.","supporting_citations":[],"review_version":1}