{"id":"812cfb11-9898-4298-be98-a87466ca8256","arxiv_id":"2411.10864","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Five 3D relativistic MHD simulations with varied jet power and triaxial ambient media all grow one-sided lobes of 500 to 700 kiloparsecs, implying giant radio galaxies can form through multiple channels.","lead":"Giant radio galaxies are jets from supermassive black holes that stretch over a million light-years, and their rarity is puzzling. This paper simulates five different jet and environment combinations and finds megaparsec-scale structures can form in all of them, suggesting they may be more common than surveys suggest.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Resolution convergence tested for only one of five runs; low-power and edge jets are most vulnerable to the 2-cell-per-diameter setup, so the generic GRG-formation claim is not yet numerically secured.","rationale":"The reader's weakest_assumption is that the resolution study was performed for only one configuration and that the other four runs may be systematically affected. This is the most load-bearing concern because the central claim—generic GRG-scale formation across all five setups—depends on each run's lobe length exceeding ~500 kpc within a plausible dynamical age. The paper's own Appendix A demonstrates that resolution changes the very quantities used to support the claim: one-sided length growth, internal energy, and jet-beam stability. The fact that the tested case is the one least likely to suffer from jet decollimation makes the untested low-power runs the critical missing evidence. I considered other potential weaknesses—such as doubling the one-sided length to infer total extent, idealized ambient profiles, and absence of radiative cooling—but none of these undermines the central claim as directly as numerical fidelity of the individual runs. The doubling assumption is justified by the central symmetry of the triaxial profile; the idealized environment is an acknowledged modeling choice; and the lack of cooling affects observed morphology and spectral ages more than the dynamical feasibility claimed here. The resolution concern is also actionable: the authors already performed a 1440×600×600 run, so extending this to the other configurations is computationally feasible and would settle the issue. The paper is otherwise honest, well-structured, and appropriately hedged about the speculative phase transition; my non-finding on the other points and agreement on this one leave the reader's CONDITIONAL verdict unchanged.","tokens_in":31057,"tokens_out":4693,"duration_ms":52193,"concrete_test":"Repeat the resolution study for GRG_lp_min and GRG_hp_edge with 4 cells per jet diameter (grids 1440×720×720 for both, matching the domain sizes in Table 1), integrating to the same physical end times as the reference runs (~170 Myr for lp_min, ~50 Myr for hp_edge). Compare one-sided lobe length versus time (using tr1≥10−7), lobe expansion speed versus length, and fractional pressure change versus length. If the 4-cell runs match the reference lengths to within 20% at each age and the break around 350 kpc persists, the resolution concern is mitigated. If lengths come up >20% short or the break shifts by >20%, the quantitative age and transition claims need revision, and the generic-formation conclusion should be downgraded to conditional.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central conclusion—that all five RMHD configurations reach one-sided lobe lengths of 500–700 kpc, implying GRG-scale growth is dynamically generic—rests on runs resolving the jet with only 2 grid cells per jet diameter. The convergence study in Appendix A covers only GRG_hp_min, the high-power minor-axis case. Even for that case, the higher-resolution run shows slower cocoon expansion, stronger jet-beam instabilities, and altered internal energy evolution (Figs. A.1–A.3); only the cocoon volume converges at late times. The paper explicitly notes that low-power jets decollimate and are more susceptible to instabilities, precisely the regime where an under-resolved narrow beam can artificially suppress or enhance growth. The lp_min and lp_maj runs are therefore more likely to be affected than the tested hp_min case. If their converged lobe lengths at 166–196 Myr are substantially shorter than 500–600 kpc, the claim that low-power FR I-like jets and major-axis propagation can form GRGs in plausible times weakens. The proposed 350 kpc 'phase transition' in expansion speed and pressure (Fig. 9) is also extracted from these under-resolved runs; a resolution-dependent break would compromise that signature. Appendix B's argument that the jet rapidly spreads over several cells after injection is qualitative and does not replace a convergence test for the decollimation and bending regimes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Giri et al. present three-dimensional relativistic MHD simulations of one-sided AGN jets in triaxial beta-profile group environments, using the PLUTO code. They run five configurations: low and high jet power, propagation along the minor and major axes, and propagation at the edge of the environment. They report that all five runs produce one-sided lobe lengths of about 500 to 700 kpc within dynamical ages of 39 to 196 Myr, which they double to total extents of 1.0 to 1.4 Mpc. The lobes are overpressured; the paper identifies a potential phase transition near a one-sided length of 350 kpc, self-similar expansion for high-power jets, and X-shaped morphologies for major-axis cases. The central claim is that GRG-scale growth is dynamically generic across these configurations, implying that observed GRG rarity is a selection or duty-cycle effect rather than a hard dynamical limit.","tokens_in":31390,"tokens_out":5963,"duration_ms":53812,"significance":"If the results hold, the paper provides a useful counterpoint to the prevailing underdense-filament picture, showing that dense group-center environments, major-axis propagation, and low-power FR I-like jets can still produce Mpc-scale structures. The predicted overpressure of active GRG lobes, the possible SRG/GRG transition, and the formation of GRG-XRGs are falsifiable with current and upcoming radio and X-ray observations. The simulation setup is described in enough detail to be reproducible, and the comparison to Kaiser-Alexander theory provides a useful benchmark. However, the significance is currently capped by the unresolved numerical fidelity issue, since the generic claim rests on results from runs with only two grid cells per jet diameter.","major_comments":[{"comment":"The production runs in Table 1 use only two grid cells per jet diameter, and the resolution study in Appendix A is limited to the single case GRG_hp_min. In that case, doubling resolution to four cells per diameter (High_Res extended) produces slower lobe-length growth and stronger jet-beam bending (Figs. A.2, A.3), so the reference-resolution run is not demonstrated to converge for the quantity that anchors the central claim. Since low-power and edge jets are explicitly described as more susceptible to decollimation and instabilities (§3.1.1, §3.1.5), the unavailability of convergence tests for GRG_lp_min, GRG_lp_maj, GRG_hp_maj, and GRG_hp_edge leaves the generic statement in §1 and the abstract unsecured. This is a load-bearing numerical fidelity issue rather than a cosmetic one.","section":"§2.2, Table 1, Appendix A"},{"comment":"The proposed phase transition near a one-sided length of 350 kpc (shaded region in Fig. 9) is derived from the same five under-resolved runs. The higher-resolution hp_min run in Appendix A shows that internal energy evolution and the expansion rate change with resolution (Fig. A.3), so a break in lobe speed and pressure derivative at 350 kpc could be a numerical artifact rather than a physical signature. The paper appropriately hedges this as a 'potential' transition, but the current evidence does not yet justify a sharp SRG/GRG transition claim; a convergence study targeting this observable is needed.","section":"§3.3.1, Fig. 9"},{"comment":"The conversion from one-sided lobe length to total source extent by doubling (e.g., 'total extent is expected to reach 1.2 Mpc' in §3.1.1) is an assumption that is not tested in the paper. The ambient medium is triaxial and rotated by 10 degrees, so the counter-jet would propagate through a different ambient column; in principle the two sides may evolve differently. Because the GRG classification in Sections 3.1.1–3.1.5 depends on the total extent exceeding 700 kpc or 1 Mpc, this assumption should be stated explicitly and ideally validated with a two-sided run or an explicit symmetry argument.","section":"§2.2, §3.1.1–3.1.5"}],"minor_comments":[{"comment":"The dynamical age of GRG_lp_min is given as approximately 166 Myr in §3.1.1 but 157 Myr in Table 2; this inconsistency should be corrected.","section":"§3.1.1, Table 2"},{"comment":"Similarly, GRG_hp_min is quoted at nearly 68 Myr in §3.1.2 but 59 Myr in Table 2, and GRG_hp_edge at 49 Myr in §3.1.5 but 39 Myr in Table 2.","section":"§3.1.2, Table 2"},{"comment":"The caption contains a typo: 'Deatil' should be 'Detail'.","section":"Fig. 4 caption"},{"comment":"The notation ∇P/P for a finite-difference fractional pressure change is unconventional; suggest ΔP/P or a clearly defined finite difference operator.","section":"§3.3.1, Eq. (13)"},{"comment":"The text refers to a 'tri-axial spheroid' shape; since the profile is a triaxial ellipsoid, this should be rephrased.","section":"§2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is relevant to A&A and the GRG community. The main concern is numerical convergence; I am not requesting rejection, but the central claim requires either additional convergence tests for the other runs or a substantial softening of the generic conclusion. The age inconsistencies and the doubling assumption should also be fixed. No issues with citation pattern or novelty disclosure were noted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what you should know: this is the first 3D relativistic MHD campaign explicitly aimed at megaparsec-scale GRG formation. Five runs—low/high power, minor-axis, major-axis, and edge propagation in a triaxial King-profile medium—all produce one-sided lobes of 500–700 kpc at dynamical ages of 39–196 Myr, which the authors double to total extents of 1.0–1.4 Mpc. The claim that GRG-scale growth is dynamically feasible across these configurations is new.\n\nWhat the paper does well: the setup is physically motivated and carefully documented. The resolution study in Appendix A is a good-faith effort; they openly admit 2 cells per jet diameter is low. The comparison to Kaiser-Alexander theory is fair, and the deviations for low-power jets are sensibly attributed to decollimation breaking self-similarity. The overpressured lobes and the X-shaped morphology from backflow along the major axis are plausible and tied to observations. The two most speculative claims are properly hedged: the 350-kpc 'phase transition' is called potential, and they call for more parameter-space coverage. No circularity: the runs are self-contained and the theory comparison is external.\n\nSoft spots, in proportion. The resolution issue is the real one. Only hp_min is convergence-tested; at 4 and 10 cells per jet diameter, the higher-resolution run shows slower lobe growth, stronger jet-beam instabilities, and altered internal energy—only cocoon volume converges at late times. The lp_min and lp_maj runs are exactly the decollimating, bending regime where under-resolution is most likely to matter, and they are not tested. Appendix B's argument that the jet spreads over several cells after injection is qualitative and does not replace a convergence test. That said, the one test that exists shows slower growth at higher resolution, so the quoted one-sided lengths are probably upper bounds rather than artificially inflated lower limits. The central qualitative picture likely survives, but the ages and the phase-transition location should be treated as approximate until the remaining runs are tested.\n\nMinor items: the total extent is taken as twice the one-sided lobe length, which overstates sources with asymmetric lobes; the 350-kpc transition is identified by eye from five runs; and there is no data-availability statement for the run configurations. None of these are load-bearing, but they are easy fixes.\n\nWho this is for: jet simulators, radio-galaxy population folks, and anyone testing GRG formation theories. It deserves a serious referee. I would send it to review, asking for at least one additional resolution test on a low-power run (or a softening of the quantitative claims), a data-availability statement, and a clearer caveat on the total-extent doubling.","headline":"First 3D RMHD campaign targeting Mpc-scale GRGs; the qualitative formation story holds up, but the quantitative ages and the 350-kpc phase transition rest on a resolution test that covers only one of five runs.","tokens_in":31936,"tokens_out":4953,"would_cite":true,"duration_ms":50219,"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":"Simulations show megaparsec radio jets can grow in dense environments, not just cosmic voids.","keywords":["giant radio galaxies","relativistic magnetohydrodynamics","jet propagation","triaxial ambient medium","X-shaped radio galaxies","lobe dynamics","AGN jets"],"falsifier":"Re-run all five configurations with at least four grid cells per jet diameter for the full domain and compare one-sided lobe length, internal pressure, and the roughly 350-kpc transition as functions of age; the central claim fails if any configuration no longer reaches about 500 kpc one-sided length or if the transition in expansion speed disappears in most runs.","tokens_in":30841,"feed_emoji":"📡","tokens_out":5508,"duration_ms":53990,"temperature":0.7,"pith_summary":"This paper uses relativistic magnetohydrodynamical simulations to ask whether megaparsec-scale giant radio galaxies (GRGs) can form through ordinary jet propagation in a range of realistic environments, or whether they require a special mechanism such as an underdense cosmic-web filament. The authors build five jet–ambient configurations—low- and high-power jets along the minor axis, along the major axis, and at the edge of a triaxial galaxy-group atmosphere—and find that every configuration develops a one-sided lobe of roughly 500–700 kpc within 39–196 Myr, which they double to total extents of 1.0–1.4 Mpc. If this result holds, it means megaparsec growth is dynamically generic rather than restricted to low-density environments, so the observed rarity of GRGs would be more likely a selection, lifetime, or duty-cycle effect. It would also mean GRGs can be born from low-power FR I-like jets and even from brightest cluster galaxies, and that X-shaped and mini-winged giant sources arise naturally from back-flow in triaxial atmospheres. A predicted SRG-to-GRG transition around 350 kpc (one-sided) in lobe expansion speed and pressure evolution gives observers a concrete length scale to test.","feed_headline":"Five jet simulations all grow roughly megaparsec radio galaxies","feed_subtitle":"Low- and high-power jets reach ~1 Mpc in groups and cluster centers within 200 Myr; rarity may be selection, not physics.","key_machinery":"The machinery is a set of five relativistic magnetohydrodynamic jet simulations in a triaxial King $\\beta$-profile atmosphere (core density $10^{-3}\\,\\mathrm{cm^{-3}}$, temperature 1.6 keV, roughly a poor galaxy group or warm-hot intergalactic medium filament), with jets injected carrying a toroidal magnetic field, magnetization $\\sigma=0.01$, and bulk Lorentz factors $\\Gamma=3$ (low power, $Q_j \\simeq 2.3\\times10^{44}\\,\\mathrm{erg\\,s^{-1}}$) or $\\Gamma=5$ (high power, $Q_j \\simeq 7.2\\times10^{44}\\,\\mathrm{erg\\,s^{-1}}$). The jet propagates along the minor axis, the major axis, or the edge of the triaxial atmosphere, producing different degrees of jet frustration. The paper tracks cocoon morphology through passive tracers, lobe shape through axial ratio, and thermodynamical state through lobe pressure, expansion speed, magnetic field, and total energy; the Kaiser–Alexander self-similar length–age relation is the theoretical benchmark for the simulated growth curves.","core_discovery":"The central claim is that giant radio galaxy formation is a generic dynamical outcome of relativistic jets propagating through triaxial galaxy-group atmospheres, not a special-case phenomenon. In all five RMHD runs the one-sided lobe exceeds 500 kpc: GRG_lp_min reaches 600 kpc after roughly 166 Myr, GRG_hp_min passes 700 kpc at about 68 Myr, GRG_lp_maj reaches about 500 kpc after 196 Myr, GRG_hp_maj reaches about 700 kpc at 137 Myr, and GRG_hp_edge reaches about 650 kpc after only 49 Myr, so the authors report total extents of 1.0 to 1.4 Mpc. Along the way they identify a dynamical phase change: after a one-sided length of about 350 kpc, lobe expansion speed and the fractional change in lobe pressure shift behavior in four of five runs, marking a possible transition from smaller radio galaxies to giants. All simulated active lobes remain overpressured relative to the ambient medium, by factors of roughly 1.9 to 7.9, and are confined by bow shocks, and the cocoon magnetic field converges to about 0.15 $\\mu$G regardless of evolutionary path.","pith_inferences":["Inference: the authors' claim implies GRG rarity is dominated by duty cycle and detection or selection effects; a direct test would measure the GRG fraction among sources above a fixed radio power in a complete survey, split by environment density.","Inference: the 350-kpc transition may track the jet escaping the core radius of the atmosphere, so it should shift when the core radius or the $\\beta$ slope changes—a parameter-space prediction the paper does not itself make.","Inference: the convergence of the dynamical magnetic field near 0.15 $\\mu$G across very different histories suggests that equipartition-based field estimates for GRG lobes may be systematically high, which would push spectral ages upward relative to dynamical ages.","Inference: extending the runs to cessation of jet activity could show whether the overpressured lobes relax on timescales that distinguish active from relic GRGs; the pressure-jump diagnostics in the paper set up exactly that comparison."],"forward_implications":["If giant phases are generic, deep radio surveys should find many more GRGs than current catalogues, including in dense group centers and around brightest cluster galaxies.","Low-power, FR I-like jets that decollimate into fat lobes can still reach megaparsec total extents, so FR I GRGs do not require a separate formation channel.","The roughly 350 kpc one-sided transition gives a specific length scale for comparing lobe expansion speeds and pressure profiles between smaller radio galaxies and giants in future observations.","X-shaped and mini-winged giant sources are natural products of back-flow diverted along the minor axis of a triaxial atmosphere, so their incidence should correlate with environmental asymmetry.","Edge-of-environment jets can cover about 1 Mpc in roughly 50 Myr, explaining how multi-megaparsec sources could grow within plausible source ages."],"supporting_citations":[{"why":"Supplies the PLUTO code used to solve the RMHD equations in all simulations.","marker":"Mignone et al. 2007"},{"why":"Provides the analytical length–age relation used as the theoretical benchmark for simulated lobe growth.","marker":"Kaiser & Alexander 1997"},{"why":"Gives the underdense, pressure-matched jet prescription and the back-flow model for wing formation.","marker":"Rossi et al. 2017"},{"why":"Defines the passive-tracer technique used to identify the cocoon and measure its length.","marker":"Mukherjee et al. 2020"},{"why":"Motivates the mapping between jet kinetic power and expected radio power and the environmental scalings used in the setup.","marker":"Hardcastle 2018"},{"why":"Supplies the injection and ambient-medium prescription adapted for the present GRG-scale runs.","marker":"Giri et al. 2023"}],"fun_headline_variants":["Simulations show giant radio galaxies may be universal","MHD runs grow megaparsec jets in every scenario","All five jet models reach giant radio galaxy size","Giant radio galaxies: common result of jet propagation","Jet simulations: megaparsec radio galaxies not rare"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that two grid cells across the jet diameter are enough: the resolution check is run for only one of the five configurations, so if that resolution systematically alters lobe length, pressure, or the 350-kpc transition in the other runs, the universal-giant-phase conclusion would weaken.","fun_headline_variants_meta":{"raw":{"variants":["Simulations show giant radio galaxies may be universal","MHD runs grow megaparsec jets in every scenario","All five jet models reach giant radio galaxy size","Giant radio galaxies: common result of jet propagation","Jet simulations: megaparsec radio galaxies not rare"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00066,"raw_usage":{"total_tokens":3108,"prompt_tokens":1127,"completion_tokens":1981,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":743,"completion_tokens_details":{"reasoning_tokens":1906}},"tokens_in":743,"tokens_out":1981,"duration_ms":14301,"temperature":1.0,"reasoning_tokens":1906,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:12:45.379278+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run all five configurations with at least four grid cells per jet diameter for the full domain and compare one-sided lobe length, internal pressure, and the roughly 350-kpc transition as functions of age; the central claim fails if any configuration no longer reaches about 500 kpc one-sided length or if the transition in expansion speed disappears in most runs.","supporting_citations":[],"review_version":1}