{"id":"6f59ff8e-e26d-4ef7-9716-80d54d4b92d7","arxiv_id":"2412.17725","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Fast AGN outflows with embedded star formation can drive tens-of-km/s radial oscillations in spherical bulges that decay after the AGN shuts off.","lead":"AGN outflows that form stars inside them can remove gas from a galaxy's bulge and set the bulge's stars gently sloshing, with radial motions of about 10 km/s. If real, such fossil oscillations would let astronomers spot galaxies that hosted powerful black-hole outflows long after the black hole went quiet.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted few-10 km/s velocities depend linearly on the ad hoc fg=0.1 mass-loss budget of Eq 37; realistic lower gas fractions would shrink them to ~1 km/s, undercutting the observational claim.","rationale":"The paper's central mechanism is a linear response of a spherical stellar bulge to a time-dependent external potential dominated by the removal of gas mass (Φ_e,bul). The amplitude of that potential is set by Eq 37, which fixes the total lost mass as fg M_bul with fg=0.1. The reader's weakest_assumption identifies exactly this as load-bearing, and I agree: the predicted few-10 km/s velocities scale essentially linearly with fg, and the chosen value is not observationally calibrated for the massive bulges the paper targets. A lower, more realistic gas fraction (e.g., 1%) would bring the predicted signal down to ~1 km/s, making the claimed observational signature vanish. The assumption of uniform removal with the original potential shape is also unphysical, since the outflow propagates from the center, but the fg dependence alone is sufficient to make the headline amplitude fragile. I do not see a more fatal internal inconsistency: the mathematics of the linear response is standard, the sign and order of magnitude of the mass-loss effect are plausible, and the figures are coherent with the stated equations except for minor wording issues (the conclusion's claim that radial velocity is 'not so dependent' on outflow velocity contradicts Fig 4, and the decay time in Fig 2 appears longer than τ_AGN). These are secondary. Because the reader already returned CONDITIONAL, this concern reinforces that verdict rather than changing it.","tokens_in":16974,"tokens_out":11350,"duration_ms":110630,"concrete_test":"Rerun the fiducial run (M_bul=2×10^10 Msun, ⌀M_out=500 Msun/yr, σ0=200 km/s, Vout(Rin)=900 km/s) with fg=0.01 in Eq 37, keeping all other settings identical, and record max |̅v_R| over 0≤τ≤2. If the peak is below ~2 km/s, the 'few 10 km/s' claim is not robust to plausible gas fractions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central amplitude claim is carried by the mass-loss term Φ_e,bul in Eq 36, whose final strength is fixed by Eq 37 to be -G fg M_bul/R_bul times a spatial factor, with fg=0.1. This is not a derived quantity: the paper chooses fg=0.1 because it gives a 2×10^9 Msun gas reservoir, but massive early-type bulges typically have much lower cold-gas fractions (≲1–5%). Since the response velocity is approximately linear in the total removed mass, reducing fg to 0.01 would reduce the quoted few-10 km/s signal to order 1 km/s, below the stated observability threshold. The assumption is also spatial: Eq 36 removes gas uniformly at every radius from t=0, despite the outflow sweeping outward from Rin; a delayed or centrally concentrated removal would alter the response and could change its sign. The outflow potential and injection terms are smaller contributors, so the headline number stands or falls with this ad hoc budget.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a linear-response calculation for a non-rotating, spherically symmetric isothermal bulge perturbed by a fast, massive AGN outflow whose cold gas forms stars. The perturbation consists of three pieces: the potential of the escaping outflow, the potential change from removing the bulge's gas, and the combined effects of stellar injection and dynamical friction. The authors solve the linearized Boltzmann equation in angle-action variables with a response-matrix formalism and compute the induced mean radial velocity of the bulge. Their central claim, stated in the abstract and Section 4, is that an episodic outflow with mass rate ~500 M_sun/yr acting on a 1e10 M_sun bulge produces radial velocities of a few tens of km/s, and that after the AGN switches off the radial velocity decays toward zero on a timescale ~tau_AGN, leaving a potentially observable kinematic fossil of past SMBH activity.","tokens_in":17228,"tokens_out":7017,"duration_ms":70697,"significance":"The idea that bulge kinematics may retain a signature of past AGN episodes is attractive, and the paper is genuinely forward-modeling: no parameter is fitted to reproduce a target radial velocity, and the response-matrix machinery is imported from established stellar-dynamics references. If the quantitative prediction were robust, it would motivate new spectroscopic searches for radial breathing motions in nearby bulges. However, the headline amplitude currently rests on an unsupported gas-removal budget and on perturbation amplitudes that are not small compared with sigma0^2. With conservative gas-fraction values typical of massive early-type bulges, the predicted signal drops to roughly the 1 km/s level, below the observability claimed in the paper. The significance is therefore conditional on a revised, better-justified amplitude estimate.","major_comments":[{"comment":"The quantitative claim is set by the assumption that one AGN episode removes the entire gas content of the bulge, Mdot*tau_AGN = f_g M_bul with f_g = 0.1. This is not derived from an observational gas budget or a self-consistent supply model, and the predicted radial velocity scales essentially linearly with the removed mass. Typical cold-gas fractions in massive early-type bulges are generally at or below a few percent, so f_g = 0.1 likely overestimates the signal by a factor of 2-10; for f_g = 0.01 the quoted few-times-10 km/s velocities would fall to order 1 km/s, below the observational threshold asserted in Section 5. The authors should either justify f_g with referenced gas-fraction measurements for the relevant bulge population or present the signal explicitly as a function of f_g and of the radial distribution of the removed gas.","section":"Section 2.4, Eqs. (36)-(37)"},{"comment":"The linearization underlying Eq. (10) requires the perturbing potential to be small compared with the unperturbed potential scale sigma0^2. With Eq. (37) and M_bul = 2 sigma0^2 R_bul/G, the final amplitude of Phi_e,bul from Eq. (36) is about 0.66 sigma0^2 at R = 0.1 R_bul and about 0.34 sigma0^2 at R = 0.5 R_bul; Figure 2 likewise shows the response potential reaching values of order 0.2-0.4 sigma0^2. These amplitudes are not a small perturbation, so the linear response calculation and the numerical velocities in Figures 3-4 are being used outside their strict validity regime. The authors should either restrict the quoted results to parameter combinations where |Phi_1|/sigma0^2 is demonstrably small or extend the calculation beyond first order, and they should quantify how the radial velocity changes within the valid regime.","section":"Sections 2.2-2.3, Eqs. (10) and (36)"},{"comment":"The temporal and spatial treatment of the mass-loss term is too idealized for the quoted parameters. The outflow is switched on with a step function H(t), and Phi_e,bul removes gas uniformly from all radii with a linear time ramp starting at t = 0. For Mdot = 800 M_sun/yr, Eq. (37) gives tau_AGN = 2.5e6 yr, while the text states the outflow crossing time is ~1e6 yr, so the two timescales are not widely separated. A centrally concentrated or radially propagating removal history would change both the amplitude and the sign of the induced response. The authors should test the sensitivity of v_R to a delayed or propagating mass-loss profile rather than assuming instantaneous, spatially uniform gas removal.","section":"Section 2.4, Eqs. (12) and (36)"},{"comment":"The numerical results are presented without convergence tests. The response matrix is sensitive to the number and shape of the potential-density basis functions, the radial grid, the temporal step Delta_tau, and the number of Fourier modes n in Eqs. (56)-(57); the paper does not state these choices or demonstrate that the results in Figures 2-4 have converged. Given that the central numbers depend on the response-matrix solution, the authors should report convergence checks or at least specify the numerical resolution and show that the radial velocities are stable under refinement.","section":"Section 3, after Eq. (67)"}],"minor_comments":[{"comment":"The normalization of f0 appears inconsistent with Eq. (30) when the potential of Eq. (31) is used; integrating f0 over velocity gives a factor exp(2) relative to the stated rho0 unless the constant is adjusted. The authors should verify the prefactor.","section":"Section 2.4, Eq. (29)"},{"comment":"The statements that the radial velocity 'exceeds a few 10 km/s' and is 'not so dependent on the outflow velocity' should be quantified with the actual ranges shown in Figures 3-4, including the dependence on radius and on time after quenching.","section":"Section 5"},{"comment":"Several figure captions and axis labels lack units or parameter definitions, such as the number of basis functions in Figure 1 and the time normalization in the right panel of Figure 2. The reader should be able to reproduce each panel from the caption alone.","section":"Figures 1-4"},{"comment":"The citation to Dootson & Magorrian (2022) is given as an arXiv preprint; if the paper has appeared in a journal, the published reference should be used.","section":"References"},{"comment":"The sentence 'Still, we find non-zero radial velocity of bulges will be driven by the episodic outflows' is grammatically awkward and should be rewritten for clarity.","section":"Abstract and Section 1"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the manuscript is within the scope of the journal and the formalism is a reasonable application of existing response-matrix methods. My main concern is quantitative robustness: the headline amplitude depends on an unvalidated f_g = 0.1 gas-removal budget and on perturbation amplitudes that violate linearity. Both issues are addressable in revision, for example by adding a parameter scan over f_g, a linearity check, and a propagating mass-loss model, so I do not recommend rejection. I saw no ethical or citation-pattern concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read of Peng et al. on bulge oscillations from AGN outflows. The core idea is genuinely new: use the time-dependent response formalism (Murali; Dootson & Magorrian) to compute how a non-rotating spherical bulge responds to episodic fast outflows carrying star formation, then predict radial velocities that persist after the AGN quenches, decaying on tau_AGN. That is a clean, falsifiable diagnostic for past AGN activity, and they do not fit anything to the target — it is a forward calculation. The formalism is imported properly and the appendices show real care with the isothermal sphere boundary conditions and the Fokker-Planck collision term. Credit where due: this is a serious linear-response exercise, not a toy.\n\nNow the soft spots. The quantitative claim, a few 10 km/s, is carried almost entirely by the mass-loss term, Eq 36, whose amplitude is fixed by Eq 37 with fg=0.1. That is not derived; it is chosen. Observed cold gas fractions in massive early-type bulges are typically well below 10%, often 1% or less. The response velocity is roughly linear in the removed mass, so a realistic fg of 0.01 gives signals around 1 km/s — below the stated observability threshold. That is the load-bearing assumption, and the paper would be far more convincing if it tested the sensitivity to fg and to the spatial distribution of the removed gas.\n\nAlso, the perturbation amplitude reaches roughly 20% of sigma0^2 by the end of the episode, which sits uncomfortably with the linear assumption. And there is a small editorial inconsistency: Section 4 says larger Vout gives smaller radial velocity, while the conclusion says the velocity is 'not so dependent' on Vout. That is not fatal — the physical dependence is through the density scaling — but it should be fixed.\n\nThe numerical implementation is not reproducible from the text alone: no code, no convergence tests, and the basis selection is brushed off as a detail. For a paper whose headline is a number, a bit more transparency would help.\n\nBottom line: the mechanism is plausible and the formalism is sound, but the amplitude claim is fragile. I would send it to a serious referee — the idea deserves to be worked out, and the referee can push for sensitivity analysis and a realistic gas budget. I would not build my own work on the few-10 km/s figure until that is done. Bring it to reading group if you want a discussion of how to make linear response predictions falsifiable.","headline":"A credible forward-model argument that episodic AGN outflows can leave fossil radial motions in bulges, but the headline few-10 km/s amplitude leans on an ad hoc fg=0.1 mass-loss budget and marginal linearity.","tokens_in":17744,"tokens_out":2241,"would_cite":false,"duration_ms":20898,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Episodic fast outflows from active galactic nuclei, with stars forming inside them, can drive coherent radial breathing in non-rotating spherical bulges at speeds of a few tens of km/s, leaving a kinematic fossil of past black hole…","keywords":["AGN outflows","galaxy bulges","bulge oscillations","dynamical friction","stellar kinematics","star formation in outflows","supermassive black holes","linear perturbation theory"],"falsifier":"Look for the predicted coherent radial velocity pattern in a bulge whose AGN is known to have been active within the last roughly one AGN lifetime: if a $10^{10}$ solar-mass bulge with an outflow of about 500 solar masses per year shows no few-tens-of-km/s expansion or contraction signal that decays after the AGN switches off, or if the measured gas-loss fraction over an episode falls well below f_g = 0.1, the central prediction is contradicted.","tokens_in":16765,"feed_emoji":"🌌","tokens_out":5905,"duration_ms":49240,"temperature":0.7,"pith_summary":"The paper argues that when a supermassive black hole goes through an active episode, the fast outflow it launches—which can contain newly formed stars—acts as a gravitational and dynamical perturbation on the surrounding bulge. Using linearized kinetic theory of stellar systems, the authors show that a spherical, non-rotating bulge responds to the outflow by developing a net radial velocity: it first shrinks slightly, then expands as gas is lost, reaching a few tens of km/s for a $10^{10}$ solar-mass bulge with a 500 solar-mass-per-year outflow. Once the AGN switches off, this radial motion does not stop instantly; it decays on roughly the AGN lifetime, meaning bulges could retain a measurable kinematic memory of past outflows. The significance is that bulge kinematics, not just gas or emission lines, could serve as a fossil diagnostic of a galaxy's black hole feeding history.","feed_headline":"Black hole winds can make galaxy bulges breathe at tens of km/s","feed_subtitle":"If the model holds, a galaxy's slow expansion and contraction records past supermassive black hole activity.","key_machinery":"The machinery is linearized kinetic perturbation theory of a collisionless stellar system in angle-action coordinates, following the response-kernel method of Dootson & Magorrian (2022) and Murali (1999). The central object is the self-gravity response potential, whose time evolution is obtained by solving the linearized Boltzmann equation with a source term from injected outflowing stars and a dynamical-friction term from their passage through the bulge; the final observable is the angle-averaged radial velocity of bulge stars. The external perturbation is the sum of two pieces: the potential of the escaping outflow, determined by mass conservation, and the potential change from gas mass lost by the bulge, normalized so that the total mass lost in one AGN episode equals the bulge's gas content. Angle-action variables allow each orbit's response to be computed by a discrete Fourier transform along the radial angle, and the time-convolution integral yields the response coefficients from which the radial velocity is assembled.","core_discovery":"The central claim is that episodic fast AGN outflows with star formation inside them act as an external potential perturbation plus a dynamical-friction source on an otherwise equilibrium spherical bulge, and that this combination drives a nonzero bulk radial velocity. The response potential oscillates and damps after the outflow ends, and the mean radial velocity follows the outflow on/off duty cycle, reaching a few 10 km/s for a $10^{10}$ solar-mass bulge with an outflow rate of 500 solar masses per year, and tending to zero within roughly the AGN lifetime after quenching. The paper identifies two competing perturbation channels: the potential of the escaping outflow itself, which pulls the bulge inward at early times, and the mass loss from the bulge's gas reservoir, which pushes the bulge outward as it expands toward a new equilibrium. Larger outflow rates produce faster radial motion, while larger outflow velocities produce weaker motion, because the outflowing density is lower and fewer outflowing stars are captured by the bulge.","pith_inferences":["If the predicted signal is real, bulge radial velocities could be used as an independent clock of AGN duty cycles in galaxies whose current AGN is weak or off, complementing ionization echoes and absorption-line proximity effects.","The same formalism could be inverted: observed residual radial motion in a quiescent bulge could be used to estimate the mass outflow rate and lifetime of the last supermassive-black-hole episode.","Because the model assumes spherical symmetry and no rotation, rotating bulges and disk contamination would mix the breathing signal with rotational kinematics, likely requiring spatially resolved stellar velocity maps to isolate it.","The authors note that slower outflows, with velocities comparable to the bulge velocity dispersion, enter a regime where the perturbation approximation breaks down and mixing efficiency approaches unity; if such outflows dominate, the observable signatures could be stronger and more dissipative than the fast-outflow case treated here."],"forward_implications":["A bulge that recently hosted a fast AGN outflow should show a coherent radial velocity pattern of order 10 km/s that tracks the AGN on/off state, not random stellar motions.","After the AGN quenches, the radial velocity decays to zero on a timescale comparable to the AGN lifetime, so the kinematic signal is a short-lived fossil of the last active episode.","Larger mass outflow rates produce faster bulge radial motion, while faster outflow velocities produce slower motion, because the outflowing density and the capture probability both drop.","The model implies that the Milky Way's bulge may currently be contracting from past activity of its central black hole, a signature testable with stellar-survey and integral-field observations of nearby bulges.","Because the same outflow episodes also shape the bulge, the predicted residual radial motion offers a direct kinematic counterpart to feedback energy traced by galaxy binding energy and stellar mass relations."],"supporting_citations":[{"why":"Supplies the isothermal sphere model, the Fokker-Planck collision operator, the dynamical-friction formula, and the angle-action coordinates used throughout.","marker":"Binney & Tremaine 2008"},{"why":"Provides the time-dependent response-kernel method, including the time-convolution integral used to evolve the bulge's self-gravity response.","marker":"Dootson & Magorrian 2022"},{"why":"Earlier time-dependent response calculation that the numerical scheme extends by adding source and dynamical-friction terms.","marker":"Murali 1999"},{"why":"Provides the scenario of star formation in cold gas of AGN outflows, the basis for treating outflowing stars as injected into the bulge.","marker":"Zubovas & King 2014"},{"why":"Offers observational evidence for star formation inside AGN outflows that motivates the injection term in the Boltzmann equation.","marker":"Maiolino et al. 2017"},{"why":"Supplies the AGN lifetime range of 1-30 Myr used to compare the bulge oscillation period with the outflow duration.","marker":"Khrykin et al. 2021"},{"why":"Links AGN winds to large-scale outflows and the two-phase interstellar medium that hosts the star formation assumed in the model.","marker":"King & Pounds 2015"},{"why":"Provides the bulge mass and M-sigma relation used to normalize the model and check its internal consistency.","marker":"Kormendy & Ho 2013"}],"fun_headline_variants":["AGN outflows drive galaxy bulges to breathe at tens of km/s","AGN winds can make bulges oscillate, leaving a velocity signature","Outflow-driven star formation makes bulges pulsate, then fade","Black hole winds induce bulge breathing, revealing past activity","AGN outflows push and pull bulges at tens of km/s"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predicted velocities rest on the assumption that over one AGN episode the outflow removes the full gas content of the bulge, f_g M_bul with f_g = 0.1, and that this gas is removed uniformly following the original potential shape; if the expelled fraction is lower or the removal is spatially concentrated, the computed tens-of-km/s velocities shrink or change sign.","fun_headline_variants_meta":{"raw":{"variants":["AGN outflows drive galaxy bulges to breathe at tens of km/s","AGN winds can make bulges oscillate, leaving a velocity signature","Outflow-driven star formation makes bulges pulsate, then fade","Black hole winds induce bulge breathing, revealing past activity","AGN outflows push and pull bulges at tens of km/s"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000879,"raw_usage":{"total_tokens":3791,"prompt_tokens":930,"completion_tokens":2861,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":2769}},"tokens_in":546,"tokens_out":2861,"duration_ms":19960,"temperature":1.0,"reasoning_tokens":2769,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:15:13.268528+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Look for the predicted coherent radial velocity pattern in a bulge whose AGN is known to have been active within the last roughly one AGN lifetime: if a $10^{10}$ solar-mass bulge with an outflow of about 500 solar masses per year shows no few-tens-of-km/s expansion or contraction signal that decays after the AGN switches off, or if the measured gas-loss fraction over an episode falls well below f_g = 0.1, the central prediction is contradicted.","supporting_citations":[{"cited_title":"1999, ApJ, 519, 580, doi: 10.1086/307408","cited_arxiv_id":null,"evidence_quote":"Earlier time-dependent response calculation that the numerical scheme extends by adding source and dynamical-friction terms."}],"review_version":1}