{"id":"984dd6fe-3908-4ca5-a873-4823025c7365","arxiv_id":"2412.09679","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In a single simulated Milky Way-mass halo, increasing Seyfert jet accretion efficiency pushes star formation to larger radii, reduces stellar mass, and can nearly quench star formation.","lead":"This paper uses high-resolution cosmological simulations to test how jets from supermassive black holes in Seyfert galaxies change the host galaxy's shape, gas, and star formation. It finds that stronger jets push star formation outward and can nearly shut it off, which helps explain observed differences in Seyfert galaxies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'only epsilon varied' control is compromised by a systematic jet-orientation trend: the simulated jet angle relative to the disk rises from 21° to 52° with epsilon, and jet-disk coupling is known to depend strongly on angle.","rationale":"Read in good faith: this is a well-structured controlled experiment on a single zoom-in halo, and the internal trends in SFR, gas distribution, and morphology are qualitatively clear. The strongest claim about outward star formation migration and near-quenching at epsilon50 is credible as a simulation result. My concern is about the attribution of that trend to epsilon alone. Table 4 shows that the final jet angle increases monotonically with epsilon, and because the SMBH spin is explicitly evolved from accreted gas angular momentum, the jet orientation is an emergent, uncontrolled variable that is known to modulate feedback. This does not make the simulation wrong, but it means the controlled experiment has a second, unisolated driver, and the causal language in the abstract and Section 4 goes beyond what the design can establish. The reader's stated weakest assumption (unresolved jet head) is a different limitation; I partly agree because both are injection-model concerns, but the orientation confound is more directly testable from the existing suite. A fixed-axis rerun would settle the issue. The CONDITIONAL verdict remains appropriate: accept the qualitative result conditional on an orientation-controlled test and broader halo sampling.","tokens_in":29347,"tokens_out":15403,"duration_ms":182649,"concrete_test":"Rerun the epsilon50 model with the SMBH spin axis artificially held fixed at the epsilon5 orientation (20.7° with respect to the initial inner stellar disk) for the entire evolution, keeping all other inputs identical, and compare the resulting SFR radial profile (Fig. 5) and B/T (Table 3) with the original epsilon50 and epsilon5 runs. If the orientation-locked epsilon50 remains similar to the free-spin epsilon50, the trend is robust to orientation; if it moves toward epsilon5, the claimed epsilon-only causal chain is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that increasing the SMBH accretion efficiency epsilon—and hence jet power—pushes star formation to larger radii and reduces B/T. But epsilon is not the only quantity that differs between models: the SMBH spin axis, which sets the jet direction, is self-consistently evolved from the angular momentum of accreted gas, and Table 4 shows the final jet angle relative to the inner stellar disk rising monotonically with epsilon: 20.7°, 37.0°, and 51.5° for epsilon5, epsilon15, and epsilon50. The same paper cites Mukherjee et al. (2018) and Talbot et al. (2022) for the strong dependence of jet feedback on the jet angle relative to the disk: jets directed into the disk launch slower, colder outflows and can trigger star formation, while out-of-plane jets couple differently to the CGM. Because the high-epsilon runs are also the runs whose jets point increasingly out of the disk plane, the outward shift of the SFR peak (Fig. 5), the central cavity, and the B/T decline (Table 3) may be a consequence of the orientation sequence rather than of the intended epsilon/energy sequence. The paper reports the angle but does not control for it or analyze it as a covariate, so the controlled-experiment interpretation is not yet established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents cosmological zoom-in simulations of a single MW-mass dark matter halo (log M_vir/M_sun ~ 11.8) at z=0, comparing four models: a baseline with supernova feedback only (epsilon0) and three models with Seyfert-type jet feedback at increasing SMBH accretion efficiency (epsilon = 4.5%, 15%, 50%). The central claim is that higher accretion efficiency, and hence higher jet power, pushes star formation to larger radii, depletes central gas, reduces the bulge-to-total stellar mass ratio, lowers the SFR by up to three orders of magnitude, and eventually quenches star formation in the highest-efficiency model. The paper also tracks jet cocoons expanding into the CGM/IGM, examines baryon fractions, and compares the galaxies to observed Seyferts and scaling relations, explicitly noting the lack of a statistical observational sample.","tokens_in":29668,"tokens_out":6609,"duration_ms":62685,"significance":"If the claimed trend holds, the paper is a useful contribution because coordinated cosmological simulations of low-luminosity Seyfert jets in MW-mass halos are rare, and the identical-initial-conditions setup is a genuine strength. The authors are also transparent about key limitations: they state that the jet head is unresolved, that only one halo is used, and that the observational comparison lacks statistical analysis. The qualitative trends, particularly the outward shift of star formation with jet strength, are interesting and likely to motivate follow-up work. However, the controlled-experiment interpretation is presently weakened by a systematic covariance between the jet orientation and epsilon, and the quantitative morphological claims lack uncertainty estimates.","major_comments":[{"comment":"The claim that the models constitute a controlled experiment varying only the accretion efficiency is compromised by the systematic variation of the jet angle with respect to the inner stellar disk. Table 4 reports the jet angle increasing monotonically with epsilon: 20.7, 37.0, and 51.5 degrees for epsilon5, epsilon15, and epsilon50. The introduction itself cites Mukherjee et al. (2018) and Talbot et al. (2022) for the strong dependence of jet feedback on jet angle, noting that disk-directed jets launch slower/colder outflows and can trigger star formation while out-of-plane jets couple differently. Since the high-epsilon runs are also the most out-of-plane, the observed radial SFR peak shift (Fig. 5), central cavity (Figs. 6 and 9), and B/T decline (Table 3) may be caused in part by the orientation sequence rather than solely by the intended epsilon/energy sequence. The manuscript reports the angle but does not analyze it as a covariate, fix the jet direction across runs, or run constant-angle controls. This should be addressed before the controlled-experiment interpretation is accepted.","section":"Section 3.4, Table 4; Section 4"},{"comment":"The unresolved jet head is load-bearing for the feedback mechanism. The authors write: 'Our numerical simulations do not resolve the jet's head (the hot spot), where the jet energy is decollimated.' The injection prescription launches perfectly collimated, zero-opening-angle particles at 3e4 km/s and 1e10 K, and the cocoon-driven feedback and the energy deposition shown in Figure 14 depend on how these particles decelerate and merge with the ambient gas. If the real decollimation and dissipation at the hot spot differ from this subgrid treatment, the outward shift of star formation and the quenching trend could be artifacts of the injection model. A resolution study or a comparison with an alternative injection geometry (e.g., different injection radius, opening angle, or energy partition) would materially strengthen the robustness of the central claim.","section":"Section 3.7; Section 2.3"},{"comment":"The Sersic decomposition is reported without any uncertainties or fit-quality metrics. The B/D and B/T values for epsilon0 and epsilon5 are close (1.20 vs 1.12 and 0.55 vs 0.53, respectively), while epsilon15 and epsilon50 differ (0.28 vs 0.38), and the bulge Sersic index is non-monotonic (1.12, 1.28, 0.71, 1.00). Without error bars on the fitted parameters, a statement of the number of free parameters, or residuals/chi-square diagnostics from the fits, the claimed systematic morphological sequence is not quantitatively established. Please add uncertainty estimates or at least a fit-quality table.","section":"Section 3.2, Table 3, Figure 15, Appendix A"},{"comment":"The paper is based on a single dark matter halo (Table 1: log M_vir/M_sun = 11.8) and explicitly states in Section 2.1 that 'This paper focuses on the results using just one of the halos.' Nevertheless, the title and conclusions generalize to 'Jetted Seyfert Galaxies at z = 0,' and Section 4.1 claims 'close agreement' with observed nearby Seyferts while acknowledging that 'statistical analysis of these properties is currently absent.' A single, noise-free realization cannot support general conclusions about the population. The authors should either add additional halos or reframe the claims as a case study and qualify the observational comparison accordingly.","section":"Section 2.1, Table 1; Section 4.1"}],"minor_comments":[{"comment":"The abstract says 'in a controlled experiment, we vary only the efficiency of the SMBH accretion,' but the epsilon0 model has no SMBH. Please clarify whether the controlled sequence refers only to the three AGN models or includes the SMBH-free baseline.","section":"Abstract; Section 2.3"},{"comment":"There is a typo: 'The amount of younger stars decreases with increasing effciency' should read 'efficiency.' Also, the sentence 'The ϵ50 model shows a very low signal and only between 20–25 kpc' is grammatically incomplete and should be revised.","section":"Section 3.3"},{"comment":"Equation (A1) defines the Sersic plus double-exponential disk fit but does not define the Sersic b_n parameter, the fitting algorithm, the radial binning, or the number of free parameters. Please provide these details so the fits can be reproduced.","section":"Appendix A, Eq. (A1)"},{"comment":"The statement 'B/T decreases ∼1/2 along the ϵ sequence' is ambiguous. Since the values go from ~0.55 to ~0.28, it should say 'decreases by a factor of about 2' or 'decreases to about half its initial value.'","section":"Section 4, bullet list"},{"comment":"Section 2.2 states that metals 'can be transported by mechanical feedback from SN and AGN,' but Section 3.5 says 'metals are distributed to the gas through SN feedback only.' Please reconcile these statements, as the AGN jet particles should carry metals from the accreted gas.","section":"Section 2.2; Section 3.5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about several limitations, which is commendable, but the aggregate effect is that the central 'controlled experiment' interpretation is not yet established. The most serious issue is the jet-angle covariance with epsilon (Table 4); because the authors themselves cite the strong dependence of jet feedback on orientation, this is a genuine confound rather than a speculative one. The single-halo and no-statistics issues further limit the strength of the conclusions. I do not think the paper should be rejected, as the qualitative trends are plausible and the simulation setup is valuable, but the revision needs either additional runs/analysis or a carefully qualified reframing."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this if you care about low-luminosity AGN feedback. It's a clean simulation experiment: same halo, same seed, only the SMBH accretion efficiency epsilon is varied, and the authors follow four models (including SN-only) to z=0. The central claims hold up as a qualitative result: higher epsilon pushes star formation outward, lowers the SFR by up to three orders of magnitude, reduces B/T, and drives expanding cocoons out to ~2 Mpc. I found the radial SFR shift in Figure 5 persuasive, and the morphology trend is monotonic and consistent. The paper is also honest about what it cannot do: it states plainly in Section 3.7 that the jet head is unresolved, and the abstract concedes that the observational comparison lacks statistics. That matters, and the authors deserve credit for saying it.\n\nThe new piece is the controlled variation of accretion efficiency in a cosmological zoom-in with collimated jet feedback, and the specific trends in SFR radius, double-exponential disk structure, and cocoon propagation. The methods are established (GIZMO, Torrey/Su spawning), and the parameters are mostly from prior literature, so there is no circular fitting to the observations being compared.\n\nSoft spots, in order of importance. First, it is one halo, one realization. The morphological fits and SFR averages carry no error bars, so the monotonic trends rest on four points without scatter estimates. Second, the unresolved jet head means the injection of 10^10 K, 3x10^4 km/s particles is a prescription; if hot-spot decollimation differs, the cocoon pressure and the radial gas displacement could change. Third, the stress-test point about jet orientation is fair: Table 4 shows the jet angle relative to the inner disk rising from 20.7 to 51.5 degrees along the epsilon sequence. That is a real covariate, and jet-disk coupling is angle-dependent, as the paper itself cites in Mukherjee et al. (2018) and Talbot et al. (2022). The authors report the angle but do not analyze it as a control. I don't think this kills the paper, because the angle is not an imposed parameter—it is the self-consistent evolution of the SMBH spin—and the direction of the bias is not obvious (out-of-plane jets may couple differently, but they still drive the cocoon feedback). Still, 'vary only epsilon' is an overstatement as written; a fixed-spin control or a time-dependent angle analysis would strengthen the claim. Finally, data and code are not released, and the observational comparison is qualitative.\n\nWho is this for: anyone working on AGN feedback in disk galaxies will get value from it. It deserves a serious referee. My recommendation would be major revision: ask for an angle-oriented robustness check, at least a discussion of the covariate, error/uncertainty estimates for the morphological fits, and ideally a second halo. If the authors address the covariate, I would take the central trend as solid.","headline":"A competent single-halo simulation study of Seyfert jet feedback showing a plausible efficiency-driven radial shift of star formation, but the claimed controlled experiment is partly muddied by a jet-orientation covariate and unresolved hot-spot physics.","tokens_in":30238,"tokens_out":5727,"would_cite":true,"duration_ms":54782,"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":"This paper shows, in controlled cosmological simulations of Milky Way-mass halos, that stronger Seyfert jets push star formation to larger radii, shrink the bulge, and nearly quench the galaxy at the highest accretion efficiency.","keywords":["Seyfert galaxies","AGN jet feedback","cosmological zoom-in simulations","star formation quenching","galactic morphology","bulge-to-total ratio","circumgalactic medium","jet cocoons"],"falsifier":"Run the same halo and epsilon sequence with the jet head resolved at sub-100 pc scales or with a different subgrid decollimation prescription; if the outward migration of the star formation peak and the central cavity disappear, the claim is an injection artifact. Observationally, a matched sample of nearby disk Seyferts with radio-detected jets, compared to non-jetted Seyferts of equal stellar mass, should show the predicted central star formation suppression and gas cavities scaling with jet power.","tokens_in":29150,"feed_emoji":"🌌","tokens_out":7913,"duration_ms":68208,"temperature":0.7,"pith_summary":"The paper asks whether the weak, collimated jets of Seyfert galaxies can shape their host galaxies as much as the powerful jets of quasars. Using cosmological zoom-in simulations of identical dark matter halos, the authors vary only the efficiency of gas accretion onto a seeded ~$10^{6}$ solar-mass black hole, producing an 'epsilon sequence' from no AGN to 50% efficiency. They find that stronger jets heat and push the interstellar gas outward, moving the peak of star formation from ~1 kpc to tens of kpc, lowering the bulge-to-total ratio from ~0.55 to ~0.28, and cutting the star formation rate by three orders of magnitude. At the highest efficiency the galaxy is essentially quenched after z ~ 1, with a central gas cavity and an extended outer gaseous ring. If correct, this means low-luminosity jetted AGN are not passive bystanders: they can determine morphology, gas content, and quiescence of Milky Way-class disk galaxies.","feed_headline":"Seyfert jets push star formation out of galactic cores","feed_subtitle":"Stronger jets shrink bulges, cut star formation 1,000-fold, and nearly quench simulated Milky Way-mass galaxies.","key_machinery":"The load-bearing mechanism is the collimated bipolar jet built from hyper-refined gas particles spawned along the SMBH spin axis at 3e4 km/s and 1e10 K, with mechanical luminosity L_jet = 1/2 eta Mdot $v^{2}$. Their energy is released when they decelerate and merge with interstellar gas; this drives an overpressured cocoon that expands perpendicular to the jet and sweeps up ambient gas. The only varied parameter is the accretion efficiency epsilon, which scales down the gravitational-torque accretion rate and hence the jet power. The radial shift of star formation is traced through Sersic decomposition of the face-on stellar surface density into a bulge plus double-exponential disk.","core_discovery":"The central discovery is a monotonic response of a galaxy to the mechanical power of its Seyfert jet in a controlled simulation. With all initial conditions, feedback recipes, and halo properties fixed, raising the accretion efficiency epsilon from 0 to 0.5 moves the star formation peak from roughly 1 kpc to 5, 10, and finally 20-25 kpc, reduces the stellar mass by a factor of ~2.5, lowers B/T from ~0.55 to ~0.28, and decreases the SFR from 2 solar masses per year to 2e-3 solar masses per year. The jets deposit most of their energy into the ISM when the injected particles merge with ambient gas, and the resulting overpressured cocoons expand to ~750 kpc, ~2 Mpc, and ~2.4 Mpc, respectively, enriching and heating the CGM. The paper presents the z=0 endpoint of these simulations and argues that the same physics explains observed central cavities, displaced gas rings, and suppressed central star formation in nearby Seyferts.","pith_inferences":["If the radial shift of star formation is real, nearby disk Seyferts with more powerful radio jets should show systematically suppressed central star formation and central gas cavities in spatially resolved maps, and the suppression should track jet power rather than bulge mass.","The near-quenching at $\\epsilon=0.5$ suggests a jet-driven path into the green valley for Milky Way-mass disks, one that does not require a major merger or a starburst.","Because the black holes are seeded only at $z\\sim3.7$ and half the stars form before that, earlier seeding would likely move the quenching epoch earlier; the authors flag this for Paper II, but it is a testable prediction of the same mechanism.","Megaparsec-scale cocoons imply that Seyfert jets can pre-heat gas in the cosmic web and suppress gas accretion onto neighboring halos, an environmental effect the paper does not quantify."],"forward_implications":["Stronger jets reduce the stellar mass within the halo and move the modeled galaxies closer to the observed $M_\\star$–$M_{\\mathrm{halo}}$ relation, indicating that jet feedback can counteract the over-cooling problem in Milky Way-mass halos.","The peak of star formation moves from $\\sim 1$ kpc to $\\sim 5$, $\\sim 10$, and $\\sim 20$–$25$ kpc as $\\epsilon$ increases, so the outer stellar disk grows in mass and the face-on profile is best described as a double-exponential disk.","At $\\epsilon = 0.5$ the central gas is evacuated into an outer ring and the star formation rate falls to $\\sim 2\\times 10^{-3}\\,M_\\odot\\,\\mathrm{yr}^{-1}$, making the galaxy essentially quiescent after $z \\sim 1$ while still retaining a stellar disk.","Jet cocoons expand to $\\sim 0.75$–$2.4$ Mpc, heating and enriching the circumgalactic medium, and the baryon fraction within the virial radius falls from $17\\%$ to $10\\%$ along the $\\epsilon$ sequence.","The simulated galaxies sit close to the observed $M_\\star$–$M_{\\mathrm{halo}}$, star-forming main sequence, and $M_{\\mathrm{bulge}}$–$\\sigma$ relations, so the jet models behave like real Seyfert hosts at $z=0$."],"supporting_citations":[{"why":"Supplies the hydrodynamics code and the meshless finite-mass solver used for the zoom-in simulations.","marker":"Hopkins (2015)"},{"why":"Introduces the accretion efficiency parameter epsilon that scales the gravitational-torque accretion rate and is the sole varied parameter.","marker":"Anglés-Alcázar et al. (2017)"},{"why":"Provides the gravitational-torque accretion formula used to compute the SMBH growth rate.","marker":"Hopkins & Quataert (2011)"},{"why":"Supplies the hyper-refined particle spawning technique used to inject jet particles into the simulation.","marker":"Torrey et al. (2020)"},{"why":"Provides the modified jet launching and merging prescription that the feedback model is built on.","marker":"Su et al. (2021)"},{"why":"Gives the stellar mass-halo mass relation used to locate the modeled galaxies and assess over-cooling.","marker":"Behroozi et al. (2019)"},{"why":"Earlier radio-mode feedback study in Milky Way-mass halos that the morphology and star formation trends are compared against.","marker":"Irodotou et al. (2022)"},{"why":"Previous collimated jet feedback simulations in Milky Way-mass galaxies showing star formation reduction and morphological changes that motivate this work.","marker":"Byrne et al. (2024)"},{"why":"Demonstrates that jet feedback can match galaxy scaling relations across a wide mass range, providing the comparison baseline for this study.","marker":"Wellons et al. (2023)"},{"why":"Supplies the M_bulge-sigma relation used to verify that the modeled Seyferts sit near observed galaxies.","marker":"Kormendy & Ho (2013)"}],"fun_headline_variants":["Seyfert jets push star formation to galaxy outskirts","Jet feedback moves star formation outward in Seyferts","Stronger jets shrink bulges, extend star-forming disks","Simulated Seyfert jets quench cores, shift star formation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulation does not resolve the jet's working surface (the hot spot) where the jet decollimates, so the entire cocoon-driven feedback picture assumes that injecting $10^{10}$ K, 3e4 km/s hydrodynamic particles faithfully represents real Seyfert jets; if the unresolved hot-spot physics dissipates energy differently, the predicted outward shift of star formation and the quenching trend could be artifacts of the injection scheme.","fun_headline_variants_meta":{"raw":{"variants":["Seyfert jets push star formation to galaxy outskirts","Jet feedback moves star formation outward in Seyferts","Stronger jets shrink bulges, extend star-forming disks","Simulated Seyfert jets quench cores, shift star formation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00061,"raw_usage":{"total_tokens":2916,"prompt_tokens":1097,"completion_tokens":1819,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":713,"completion_tokens_details":{"reasoning_tokens":1753}},"tokens_in":713,"tokens_out":1819,"duration_ms":15009,"temperature":1.0,"reasoning_tokens":1753,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:50:30.491018+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same halo and epsilon sequence with the jet head resolved at sub-100 pc scales or with a different subgrid decollimation prescription; if the outward migration of the star formation peak and the central cavity disappear, the claim is an injection artifact. Observationally, a matched sample of nearby disk Seyferts with radio-detected jets, compared to non-jetted Seyferts of equal stellar mass, should show the predicted central star formation suppression and gas cavities scaling with jet power.","supporting_citations":[],"review_version":1}