{"id":"ec47aa71-ec6b-4412-9870-1372afd46563","arxiv_id":"2412.14433","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In simulations, young radio galaxies interacting with a clumpy interstellar medium can mimic double-double restarted sources, and their radio spectral turnover frequency traces ISM density independent of jet power.","lead":"This paper simulates young radio jets from supermassive black holes as they push through the clumpy gas of their host galaxy and out into the surrounding environment, then predicts the radio emission observers would see. The simulations show that young, active sources can be mistaken for restarted ones, and that the spectral turnover frequency depends on the gas density around the source, not the jet power.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'new probe of the ISM' claim is undermined by a degeneracy the paper itself concedes: Section 6.1 states the nu_p-L_S relation 'is likely to also depend on the scales over which the clouds are distributed', so an observed turnover plus source size need not yield a unique ISM density.","rationale":"I read the paper as a careful forward-modeling study whose headline claim is the ISM-density probe. The jet-power independence is well-motivated and demonstrated for the two simulated powers; the paper is honest about idealizations. My concern targets the stronger inference: that turnover frequency can be used to measure ISM density without knowing jet power. That inference requires a one-to-one mapping from (nu_p, L_S) to gas density across plausible ISM geometries. The simulations vary only the density normalisation of a single fixed cloud geometry; the authors explicitly defer 'detailed analysis of alternative cloud distributions to future work' and concede that the nu_p-L_S relation may depend on cloud distribution scales. Because FFA optical depth is an integral over the absorbing column, clumpiness and covering factor are degenerate with mean density. A source with a lower mean density but a more favourable covering factor could produce the same nu_p at the same L_S as a denser, more clumpy ISM. The rotated-environment runs provide some evidence against orientation dependence, but they do not vary the statistical properties of the cloud field. I therefore regard the cloud-distribution degeneracy as the single most load-bearing concern. The reader's weakest assumption (non-relativistic jets and no cooling) could change the amount of cloud clearing and thus the quantitative nu_p-L_S curve, but it would not by itself invalidate the logic of the probe; the missing calibration across cloud geometries does. I agree with the CONDITIONAL verdict: the paper is valuable but the headline claim needs either additional simulations or a semi-analytic calibration across ISM structure parameters before it can be used as a robust ISM probe.","tokens_in":23199,"tokens_out":5590,"duration_ms":47900,"concrete_test":"Hold the volume-averaged ISM density and jet power fixed at the reference values (n_w,0 = 400 cm^-3, Q_jet = 10^44 erg s^-1) and generate two additional cloud fields: one with a different fractal seed and one with clouds extending to 5 kpc (or with a different cloud size distribution). For each, recompute synthetic integrated spectra at matched source sizes using the published method (Section 2.3). If nu_p at L_S = 3 kpc differs by more than the Callingham et al. (2015) fitting uncertainty, the turnover frequency is degenerate with cloud spatial distribution and the 'new probe of the ISM' conclusion requires qualification. A cheaper first step is a semi-analytic FFA calculation along the simulated lines of sight, rescaling cloud clumpiness or radial extent while preserving mean density.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two components: (i) peak frequency nu_p depends on ISM density but not jet power; (ii) turnover in resolved young sources is therefore a new probe of the ISM. Component (i) is supported by two jet-power simulations and by the FFA coefficient being independent of jet properties; component (ii) additionally requires that nu_p(L_S) is a unique function of ISM density across the range of real ISM structures. That uniqueness is not established. The simulated ISM is one fixed realisation of a lognormal fractal cube (correlation lengths 20-250 pc) truncated at r_gal = 2.5 kpc; only the normalisation n_w,0 is varied. FFA optical depth is a path integral over absorbing gas, so at fixed mean density and source size, different cloud covering factors, clumpiness scales, or radial extents will change tau_nu and hence nu_p. Section 6.1 explicitly notes this: 'this relationship is likely to also depend on the scales over which the clouds are distributed', and states that 'the maximum extent of clouds encountered by each jet in our simulations is 2.5 kpc', limiting FFA turnover to source sizes of a few kpc, below the largest observed peaked-spectrum sources. Thus the proposed probe is not yet calibrated: an observed (nu_p, L_S) pair cannot be inverted to a density without independent knowledge of the cloud geometry. The reader's concern about non-relativistic jets and missing cooling is real, but it would shift the quantitative nu_p-L_S relation; the cloud-distribution degeneracy is more load-bearing because it affects the logical status of the headline claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents hydrodynamic simulations of AGN jets propagating through a clumpy, multiphase ISM and into a surrounding cluster or group atmosphere, with synthetic radio observations constructed by combining synchrotron emission from Lagrangian tracer particles with free-free absorption along the line of sight. The authors study how jet–cloud interactions create persistent length and brightness asymmetries, how such asymmetries and the apparent core prominence can lead to misclassification of young sources as restarted or head-tail sources, and how the spectral turnover frequency depends on source size, ISM density, and jet power. They report that the peak frequency is lower for lower ISM density and independent of jet power, and propose that spectral turnover in resolved young sources provides a new probe of the ISM.","tokens_in":23509,"tokens_out":2617,"duration_ms":28454,"significance":"If the central claims hold, the paper would provide a novel observational route to estimating the gas density around young radio sources without needing to know the jet power, and would sharpen the interpretation of double-double morphologies and core prominence as restarted-source indicators. The forward-modeling pipeline is a genuine strength: the simulations self-consistently combine hydrodynamics, synchrotron emissivity with adiabatic and radiative losses, and free-free absorption, and the synthetic observations at different resolutions and redshifts give concrete, falsifiable predictions for survey classification. The qualitative reproduction of the observed inverse peak-frequency–size correlation and the demonstration that jet-power dependence does not appear in the tested cases are useful results. However, the interpretation of turnover as a calibrated ISM probe is not yet supported, because the paper itself concedes that the relation depends on the spatial distribution of the absorbing clouds.","major_comments":[{"comment":"The central claim that spectral turnover provides a new probe of the ISM is not supported as stated. Section 6.1 explicitly concedes that the νp–LS relation 'is likely to also depend on the scales over which the clouds are distributed', and that the maximum cloud extent in the simulations is 2.5 kpc, limiting FFA turnover to source sizes of a few kpc. Because the free-free optical depth in Eq. (9) is a path integral over absorbing gas along each line of sight, at fixed mean density and source size different cloud covering factors, clumpiness scales, or radial extents will change τν and hence νp. Thus an observed (νp, LS) pair cannot be inverted to a unique ISM density without independent knowledge of the cloud geometry. This limitation is load-bearing for the abstract and conclusion (vi), so the claim should either be removed or substantially weakened to state that turnover depends on ISM density and geometry, and that the present simulations demonstrate the density dependence for one fixed cloud distribution.","section":"Section 6.1, abstract, Section 7(vi)"},{"comment":"The jet-power independence of the peak frequency is tested with only two jet powers (1044 and 1043 erg s–1) in a single environment, and the shaded regions in Figure 13 represent fitting uncertainties only, not systematic variations of the environment or jet physics. The absence of explicit jet-power dependence in Eq. (9) makes the result unsurprising at the level of the absorption coefficient, but the simulations are needed to show that jet-induced cloud clearing does not break this dependence for the tested cases. The claim should be framed as applying to the parameter range explored, and the statement in Section 6.2 that 'the turnover frequency and general shape of the spectrum do not change with jet power' should be qualified accordingly.","section":"Section 6.2, Figure 13"},{"comment":"The simulation suite is very small for the breadth of the conclusions: one group environment, two ISM densities, two jet powers, and one lognormal fractal cloud realization (with two 90-degree rotations serving as additional realizations for morphology but not for the spectral analysis). The persistence of asymmetries in poor-group environments is demonstrated for a single group profile and a single host galaxy, and the νp–LS relation is derived from two density normalizations. The paper should state more explicitly that the quantitative results, particularly the proposed density calibration, are proof-of-concept and may shift when the cloud geometry, group profile, or jet velocity is varied.","section":"Section 2.2, Section 4.2, Table 2"}],"minor_comments":[{"comment":"There is a typo: 'making it difficult to to draw conclusions' should read 'making it difficult to draw conclusions'.","section":"Section 5.1"},{"comment":"Several figure labels and captions contain placeholder characters (e.g., '□' in density and surface brightness units), which should be replaced with proper superscripts or symbols in the production version.","section":"Figures 10, 14, 15"},{"comment":"The statement that an ISM density of 150 cm–3 is 'perhaps more representative of the environments surrounding observed sources' is based on only two density normalizations and one cloud geometry; this sentence should be tempered to avoid overinterpreting the comparison with observational data.","section":"Section 6.1"},{"comment":"A brief quantitative statement about the expected effect of radiative cooling on the cloud ablation timescales, rather than only a qualitative reference to Antonuccio-Delogu & Silk (2008), would help the reader assess the robustness of the asymmetry results.","section":"Section 2.2"},{"comment":"The citation to 'Stewart et al. 2024, private communication' for the smoothing-kernel method is not verifiable; if this method is described in a paper in preparation or in a thesis, that should be cited instead.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of PASA and the technical work is solid, but the headline claim of a new ISM probe is stronger than the evidence. I would encourage the editor to ask for a revision that either removes or carefully qualifies the 'new probe' language, and that clearly states the limited parameter coverage of the simulation suite. The manuscript is otherwise publishable after these changes."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this is a solid forward-modeling paper that extends Bicknell et al. (2018) to 60 kpc scales, adds synthetic observations, and produces three genuinely new predictions: environment-dependent persistence of lobe asymmetries, continuously active sources mimicking double-double restarted morphologies, and the insensitivity of the FFA turnover frequency to jet power. The modeling is careful and the authors are honest about their simplifications—they flag missing cooling, non-relativistic jets, and the truncation of clouds at 2.5 kpc. The synthetic radio maps with realistic resolution and sensitivity are a useful step.\n\nThe central claim that the turnover frequency depends on ISM density but not jet power is partly inherited from the FFA coefficient, which has no jet-power dependence. The simulations usefully show that jet-cloud clearing does not break this for the tested parameters, but that is a narrow test: two jet powers, one fractal cloud realization, two ISM densities. The stronger headline claim—that turnover in resolved young sources is a new probe of the ISM—is undermined by a degeneracy the paper itself concedes in Section 6.1: the nu_p-L_S relation likely also depends on the scales over which clouds are distributed. Since FFA optical depth is a path integral over absorbing gas, at fixed mean density and source size, different cloud covering factors, clumpiness scales, or radial extents will change the turnover. The simulated ISM is one fixed lognormal fractal cube truncated at 2.5 kpc; only the normalization is varied. So an observed (nu_p, L_S) pair cannot yet be inverted to a density without independent knowledge of cloud geometry. This is not fatal for the morphology results, but it does mean the 'new probe' is not yet calibrated.\n\nThe missing cooling and non-relativistic jets are worth noting, but as the paper says, the more complex small-scale dynamics are likely missed; that would shift the quantitative relation, not the logical structure. The fitting uncertainties in Figure 13 are also just fitting uncertainties, not model scatter.\n\nCitation pattern looks fine; they build on Bicknell et al. and related work. No public code or data, only 'on reasonable request,' which is a shame for a simulation paper.\n\nWho this is for: radio galaxy evolution people, especially those working on peaked-spectrum sources and restarted-source classification. They will get value from the double-double mimicry and the caution about core prominence. The ISM-density claim needs more simulation variety before it becomes a reliable observing tool.\n\nVerdict: conditional accept. I'd send it to a serious referee. If you work on compact radio sources, cite it; I'd bring it to a reading group focused on observational classification biases.","headline":"Careful simulation study with useful morphology predictions, but the headline 'new probe of the ISM' claim is not yet calibrated because the paper itself concedes the turnover depends on cloud scales, not just density.","tokens_in":24100,"tokens_out":1602,"would_cite":true,"duration_ms":12693,"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":"The paper argues that the turnover frequency of a young radio galaxy's integrated spectrum is fixed by host ISM density and independent of jet power, offering a new way to weigh the gas around these sources.","keywords":["galaxies: jets","radio continuum: galaxies","ISM: jets and outflows","galaxies: active","hydrodynamics","free-free absorption","peaked-spectrum radio sources","spectral turnover"],"falsifier":"Compare two observed resolved peaked-spectrum sources matched in linear size and host ISM density but with jet powers differing by an order of magnitude: the paper predicts identical rest-frame turnover frequencies, so a systematic difference in $\\nu_p$ would falsify the central claim.","tokens_in":22956,"feed_emoji":"📡","tokens_out":7997,"duration_ms":64304,"temperature":0.7,"pith_summary":"The paper uses three-dimensional hydrodynamic simulations of jets launched from an active galactic nucleus to argue that the radio spectrum of a young radio source carries a direct, measurable imprint of the gas it is plowing through. The central claim is that the frequency at which the integrated radio spectrum peaks is set by the density of the host galaxy's multiphase interstellar medium and does not depend on the power of the jet. If this holds, an observer needs only the turnover frequency and the linear size of a resolved young source to estimate the surrounding gas density, with no knowledge of jet power required. The simulations also show that a single continuous young jet can produce double-double morphology and high core prominence, so these commonly used restarted-source signatures do not by themselves prove multiple epochs of activity.","feed_headline":"Spectral peak of young radio galaxies reveals gas, not jet power","feed_subtitle":"Measuring turnover frequency and source size could estimate the host gas density without knowing the jet's power.","key_machinery":"The carrier of the argument is the combination of a three-dimensional hydrodynamic jet simulation with a synthetic radio pipeline: Lagrangian tracer particles track electron packets whose adiabatic and radiative losses set the synchrotron emissivity, and a free-free absorption optical depth is integrated along every line of sight through the simulated gas. The load-bearing identity is the free-free absorption coefficient $\\alpha_\\nu \\propto n_e n_i \\nu^{-2}$, which depends only on the absorbing gas and the observing frequency, so the spectral turnover carries no memory of jet power. The multiphase environment is built from a lognormal distribution of dense clouds in pressure equilibrium with hot diffuse gas inside a double-isothermal galaxy potential, matched to a $\\beta$-profile for the surrounding cluster or group.","core_discovery":"On the paper's own terms, the discovery is that free-free absorption of synchrotron emission by the clumpy multiphase ISM sets the turnover in the radio spectrum of young radio galaxies, and because the absorption coefficient depends on gas density, temperature and frequency but not on any jet property, the turnover frequency at a fixed source size tracks ISM density alone. The simulations reproduce the observed inverse correlation between peak frequency and source size and show that halving the central cloud density shifts the turnover to lower frequencies, while changing jet power by an order of magnitude leaves it unchanged. A second discovery is that the same dense clouds that produce the free-free screen also slow the jet and counterjet unevenly, creating lobe length and brightness asymmetries that persist to tens of kiloparsecs in poor-group environments but wash out in cluster environments.","pith_inferences":["If the calibration holds, turnover frequency could serve as an ISM weighing tool for high-redshift peaked-spectrum samples where direct molecular gas tracers are impractical.","The jet-power independence is unlikely to survive unchanged if relativistic spines or radiative cooling change how jets destroy clouds; rerunning the pipeline with a relativistic jet treatment would test whether the density-only relation persists.","The demonstration that young continuous sources mimic restarted morphology implies that inferred restarted-source fractions in flux-limited samples may be overestimated, and additional spectral curvature metrics may separate the two populations.","The free-free turnover in these runs only persists while the source is inside the 2.5 kpc cloud region, so extending the absorbing cloud distribution or inclining the jet could push free-free turnovers to the larger sizes where observed peaked-spectrum sources live."],"forward_implications":["A measured rest-frame turnover frequency and linear size for a resolved young source yield an estimate of the central ISM density without any assumption about jet power.","Young, continuously active sources can display double-double morphology and core prominence above 0.1, so restarted-source classifications based on these signs alone will include some false positives.","Lengthening asymmetries that persist to tens of kiloparsecs point to a poor-group environment rather than to a restart episode, because flat cluster profiles let the shorter jet catch up.","The simulations reproduce the observed $\\nu_p \\propto L_S^{-0.65}$ trend qualitatively, with lower ISM densities (150 cm$^{-3}$) giving better agreement with observed sources than 400 cm$^{-3}$.","High-resolution and surface-brightness-sensitive observations are needed to avoid misclassifying young sources; at high redshift the connecting bridge can drop below detectability, splitting a single source into unrelated components."],"supporting_citations":[{"why":"Supplies the free-free absorption formalism and the multiphase-ISM simulation setup that this work extends, and is the comparison point for the peak-frequency versus size trend.","marker":"Bicknell et al. (2018)"},{"why":"Provides the observed inverse correlation between turnover frequency and linear source size that the simulations are matched against.","marker":"O’Dea & Baum (1997)"},{"why":"Adds the observational peak-frequency and source-size data used to compare simulated and observed spectra.","marker":"Jeyakumar (2016)"},{"why":"Supplies the semi-analytic synchrotron emissivity calculation used to convert simulated fluid quantities into radio luminosity.","marker":"Turner et al. (2018)"},{"why":"Provides the Lagrangian-particle adiabatic and radiative loss factor used to age the electron population.","marker":"Yates-Jones et al. (2022)"},{"why":"Provides the lognormal fractal cloud distribution, pressure equilibrium setup, and temperature cutoff for free-free absorption.","marker":"Mukherjee et al. (2016)"},{"why":"Provides the spectral model used to fit the turnover frequency and its uncertainty in the simulated spectra.","marker":"Callingham et al. (2015)"},{"why":"Provides the analytic lobe expansion scaling used to explain why asymmetries persist in group but not cluster environments.","marker":"Kaiser & Alexander (1997)"}],"fun_headline_variants":["Turnover in young radio spectra traces gas density, not jet power","Free-free absorption sets radio peak independent of jet power","Spectral cutoff of young jets reveals host ISM density","Peak frequency of young radio galaxies depends on gas, not jets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations treat the jet as a slow, dense flow moving at a tenth of light speed and leave out radiative cooling; if real jets or cooling change how the dense clouds are shredded, the claimed density-only, jet-power-independent turnover could fail.","fun_headline_variants_meta":{"raw":{"variants":["Turnover in young radio spectra traces gas density, not jet power","Free-free absorption sets radio peak independent of jet power","Spectral cutoff of young jets reveals host ISM density","Peak frequency of young radio galaxies depends on gas, not jets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000278,"raw_usage":{"total_tokens":1661,"prompt_tokens":958,"completion_tokens":703,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":633}},"tokens_in":574,"tokens_out":703,"duration_ms":6572,"temperature":1.0,"reasoning_tokens":633,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:15:04.843427+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare two observed resolved peaked-spectrum sources matched in linear size and host ISM density but with jet powers differing by an order of magnitude: the paper predicts identical rest-frame turnover frequencies, so a systematic difference in $\\nu_p$ would falsify the central claim.","supporting_citations":[],"review_version":1}