{"id":"649d2717-8bda-402e-afe2-f6606c22ceb0","arxiv_id":"2607.07462","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":9,"one_line_summary":"2D general-relativistic simulations of Bondi-Hoyle-Lyttleton accretion onto a Kerr black hole in a common envelope yield analytical fits for mass/momentum accretion rates and bremsstrahlung luminosity as functions of spin, Mach number, and density gradient.","lead":"This paper simulates a spinning black hole plowing through the envelope of a giant star, measuring how fast it swallows gas and how the flow changes with spin, speed, and density gradients. It provides fitting formulas for accretion rates and luminosity that could feed into models of how binary stars evolve and produce gravitational-wave sources.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Fitting formulas assume separability of spin and Mach-number dependence, but no simulation varies both simultaneously — the cross-term structure is entirely unconstrained.","rationale":"The reader correctly identified that the fitting formulas lack goodness-of-fit metrics and hold-out validation, and that the 2D geometry introduces biases in volume-integrated quantities. These are legitimate concerns. However, the reader's primary focus on the 2D slab geometry as the 'structurally fragile premise' misses a more fundamental issue: the fitting formulas assume multiplicative separability between a★ and M, but the simulation suite never varies both parameters simultaneously. This means the cross-term structure — which determines whether the formulas are correct in the interior of the claimed parameter space — is entirely unconstrained by data. The 2D geometry issue affects absolute calibration (systematic offsets), which the authors acknowledge. The separability issue affects correctness of the functional form itself, which the authors do not acknowledge. A population synthesis code using these formulas would query regions like (a★=0.5, M=1.4, ερ=0.75), where no simulation exists and the assumed separability is untested. The verdict remains CONDITIONAL — the qualitative morphological results and the individual-slice trends are sound, and the paper is a reasonable first step — but the condition should explicitly note that the fitting formulas are only validated on two orthogonal 2D slices of the 3D parameter space, not on the full volume they claim to describe. The test I propose (a handful of off-slice models) is computationally cheap relative to the 66 already run and would directly determine whether the separability assumption holds.","tokens_in":16322,"tokens_out":2559,"duration_ms":190789,"concrete_test":"Run 4–6 simulations with both a★≠0 and M≠2 (e.g., a★=+0.9375 with M=1.0, 1.5; a★=−0.9375 with M=1.5, 2.0; each at ερ=0.5 and 1.0). Compare the measured Ṁ₀, Ṗ_r, Ṗ_φ, and L_BR against the predictions of Eqs. 7–11. If deviations exceed ~20–30% for any quantity, the separability assumption fails and the formulas cannot be reliably interpolated across the joint parameter space without including cross-terms.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central deliverable is the set of fitting formulas (Eqs. 7–11) that claim to describe accretion rates and luminosity across the full (a★, M, ερ) parameter space. However, the 66-model suite is structured as two non-overlapping slices: Table 1 varies a★ and ερ at fixed M=2 (30 models), and Table 2 varies M and ερ at fixed a★=0 (36 models). No model has both a★≠0 and M≠2. Every fitting formula assumes multiplicative separability between the a★-dependent factor and the M-dependent factor (e.g., Eq. 7: (1+μ₀a²★)×(μ₁+μ₂ερ+μ₃ερ²)×M²; Eqs. 8–9 and 11 have the same structure). This separability is an assumption, not a result: if frame-dragging effects interact with Mach number (physically plausible — the shock-cone geometry and subsonic region size both depend on M, and frame-dragging acts on the post-shock flow), the true dependence would contain irreducible cross-terms like a★×M or a★×M² that the formulas cannot capture. The formulas would then be wrong not just in absolute magnitude (as the 2D geometry issue would cause) but in their qualitative predictions for the joint dependence — which is exactly what a population synthesis code would query. The authors do not acknowledge this limitation. The 2D geometry concern raised by the reader is real but acknowledged by the authors and affects systematic offsets; the separability concern is unacknowledged and affects whether the formulas are correct at all in the parameter-space regions they claim to cover.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This manuscript presents 2D general-relativistic hydrodynamic simulations of Bondi–Hoyle–Lyttleton accretion onto a stellar-mass Kerr black hole embedded in a red supergiant common envelope. The authors use the BHAC code to explore 66 models spanning five black hole spins, six Mach numbers, and six density gradients. The study reports flow morphologies (shock cones, bow shocks, subsonic regions) and provides empirical fitting formulas for the mass accretion rate, radial and angular momentum rates, and bremsstrahlung luminosity as functions of spin, Mach number, and density gradient. The test-fluid approximation is justified via a compactness argument. The paper is positioned as a first GR study of this problem and a stepping stone toward 3D simulations.","tokens_in":16612,"tokens_out":1456,"duration_ms":277593,"significance":"The paper provides a systematic GRHD parameter survey of BHL accretion in a common-envelope context with a rotating black hole, which is a genuinely new combination. The use of a well-tested code (BHAC) and the explicit justification of the test-fluid approximation via the compactness ratio are strengths. The fitting formulas (Eqs. 7–11) are falsifiable in the sense that they make specific quantitative predictions that future 3D simulations can test. The qualitative flow morphologies are consistent with prior BHL studies, lending credibility to the simulation infrastructure. However, the central deliverable—the fitting formulas—rests on assumptions that limit their predictive power, as detailed below.","major_comments":[{"comment":"§3.2, Eqs. (7)–(11): The fitting formulas assume multiplicative separability between the spin-dependent factor and the Mach-number-dependent factor. However, the 66-model suite is structured as two non-overlapping slices: Table 1 varies a★ and ε_ρ at fixed M=2 (30 models), and Table 2 varies M and ε_ρ at fixed a★=0 (36 models). No simulation has both a★≠0 and M≠2. Consequently, any cross-term dependence (e.g., a★×M or a★×M²) is entirely unconstrained by the data. If frame-dragging effects interact with the Mach number—a physically plausible scenario given that both the shock-cone geometry and the subsonic region size depend on M while frame-dragging acts on the post-shock flow—the formulas would be qualitatively wrong in the joint (a★, M) parameter space they claim to cover. This is load-bearing because the formulas are presented as valid across the full parameter space. The authors must","section":null},{"comment":"§2.1, paragraph beginning 'It is worth noting that the 2D simulations captures the essential physical processes': The 2D slab geometry assumes out-of-plane gradients are negligible, and the authors acknowledge this 'may introduce biases in volume-integrated quantities such as the mass and momentum accretion rates and the bremsstrahlung luminosity.' Since the fitting formulas (Eqs. 7–11) are precisely for these volume-integrated quantities, the 2D limitation directly affects the quantitative reliability of the paper's main output. The authors acknowledge this qualitatively in §4, but the fitting formulas are presented without any error bars or quantitative discussion of the expected systematic offset from 3D. At minimum, the manuscript should state explicitly in §3.2–3.3 that the fitting coefficients carry an unquantified systematic uncertainty from the 2D geometry, and the formulas'适用范围","section":null},{"comment":"§3.2, Eqs. (7)–(11): No goodness-of-fit metric (R², χ², RMS residual) is reported for any of the five fitting formulas. Without these metrics, a reader cannot assess whether the formulas provide an adequate description of the simulation data or identify where the largest residuals occur. This is essential for formulas that are intended for use in population synthesis codes. Please add a quantitative measure of fit quality, ideally with a residual plot or table showing the maximum deviation for each formula.","section":null}],"minor_comments":[{"comment":"§2, Eq. (4): The density profile uses ε_ρ in the exponent as ε_ρ(r−r₀)/r_acc, but the text later refers to the gradient parameter as both ε_ρ and ε (e.g., §3.1 uses ε while §2 uses ε_ρ). Please unify the notation.","section":null},{"comment":"§2, Eq. (6): The pressure formula appears to mix a polytropic and ideal-gas EOS. The expression p_in = c²_{s,∞}/(Γ−1) · Γ(Γ−1)^{−c²_{s,∞}/Γ} · ρ_in is unusual; please clarify its derivation or cite the source.","section":null},{"comment":"§3.2, Eq. (7): The formula has a log on the left-hand side but the right-hand side is a product of polynomial factors. Please confirm that the fitting was performed in log-space and state this explicitly.","section":null},{"comment":"Figure 3: The color bar labels use notation like '□1.0' which is unclear. Please use standard colorbar formatting.","section":null},{"comment":"§3.3: The statement 'this supports the assumption in the literature that these two quantities are proportional, i.e., Ṁ₀ ∝ L_BR' is made without a quantitative comparison. A simple scatter plot of L_BR vs. Ṁ₀ across all 66 models would strengthen this claim.","section":null},{"comment":"§4: 'as the stellar black hole approached to the red giant star core' should read 'approaches'.","section":null},{"comment":"The abstract states 'offering first insights into general relativistic hydrodynamics modelling of the secular evolution of the common envelope phase.' The phrase 'secular evolution' may overstate the scope, given that the simulations are local and steady-state; consider softening to 'local accretion dynamics'.","section":null}],"recommendation":"major_revision","confidential_remarks":"The separability issue is the most serious concern. The authors fit formulas with a specific multiplicative structure but have zero data points in the region where both a★ and M vary. This is not a subtle point—it means the formulas are extrapolations, not interpolations, in the joint parameter space. If the authors can add even a handful of models with both a★≠0 and M≠2, this would substantially strengthen the paper. If not, they must at minimum explicitly state that the joint dependence is unconstrained and that the separability is an untested assumption. The 2D limitation is acknowledged but should be more prominently flagged in the sections presenting the formulas."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. All three major comments identify legitimate gaps that we will address in the revised manuscript. Specifically: (1) the fitting formulas assume multiplicative separability between spin and Mach-number dependence, but no simulation probes the joint (a★ ≠ 0, M ≠ 2) parameter space, so cross-terms are unconstrained—we will add explicit caveats and scope restrictions. (2) The 2D geometry introduces unquantified systematic uncertainty in the volume-integrated fitting formulas—we will state this explicitly alongside the formulas. (3) No goodness-of-fit metrics are reported—we will add R² values and maximum residuals for all five formulas. We agree with all three points and will revise accordingly.","responses":[{"response":"The referee is correct. Our 66-model suite consists of two orthogonal slices: Table 1 probes (a★, ε_ρ) at fixed M=2, and Table 2 probes (M, ε_ρ) at fixed a★=0. No simulation simultaneously has a★≠0 and M≠2. The multiplicative separability assumed in Eqs. (7)–(11) is therefore an ansatz, not a result constrained by data in the joint (a★, M) subspace. We agree that cross-terms such as a★×M are physically plausible: frame-dragging acts on the post-shock flow whose geometry and extent depend on M, so an interaction is not implausible. We will revise the manuscript to state this limitation explicitly in §3.2. Specifically, we will: (i) add a paragraph after Eq. (11) noting that the fitting formulas assume separability and that no simulation constrains cross-terms in the joint (a★, M) space; (ii) restrict the claimed validity of the formulas to the two slices actually simulated, noting that extrapolation to the joint parameter space is untested; and (iii) add a corresponding caveat in §4. We note that the M=2 slice does include all five spin values, so the spin dependence at that Mach number is directly constrained, and the a★=0 slice constrains the Mach dependence for a non-rotating hole. The separability assumption interpolates between these two constrained limits. This is a reasonable first-order approximation, but it should be stated as such, and we will do so.","revision_made":"yes","referee_comment":"§3.2, Eqs. (7)–(11): The fitting formulas assume multiplicative separability between the spin-dependent factor and the Mach-number-dependent factor. However, the 66-model suite is structured as two non-overlapping slices: Table 1 varies a★ and ε_ρ at fixed M=2 (30 models), and Table 2 varies M and ε_ρ at fixed a★=0 (36 models). No simulation has both a★≠0 and M≠2. Consequently, any cross-term dependence (e.g., a★×M or a★×M²) is entirely unconstrained by the data. If frame-dragging effects interact with the Mach number—a physically plausible scenario given that both the shock-cone geometry and the subsonic region size depend on M while frame-dragging acts on the post-shock flow—the formulas would be qualitatively wrong in the joint (a★, M) parameter space they claim to cover. This is load-bearing because the formulas are presented as valid across the full parameter space."},{"response":"We agree. The fitting formulas in Eqs. (7)–(11) are for volume-integrated quantities (mass accretion rate, momentum rates, bremsstrahlung luminosity), and the 2D slab geometry introduces a systematic uncertainty that we have acknowledged qualitatively in §4 but not alongside the formulas themselves. This is a valid concern. In the revised manuscript, we will add an explicit statement in §3.2 (immediately after the fitting formulas) and in §3.3 noting that the fitting coefficients carry an unquantified systematic uncertainty from the 2D geometry, and that the formulas should be regarded as first-order approximations requiring calibration against 3D simulations before quantitative use in population synthesis. We note that previous studies (Gracia-Linares & Guzmán 2015; Kim & Most 2024, cited in our manuscript) found that 2D BHL simulations reproduce the main qualitative trends seen in 3D while revealing quantitative differences, which provides some evidence that the systematic offset is moderate rather than order-unity, but we agree this does not substitute for a direct 3D calibration. We will also add a sentence in the abstract noting the 2D limitation more prominently.","revision_made":"yes","referee_comment":"§2.1, paragraph beginning 'It is worth noting that the 2D simulations captures the essential physical processes': The 2D slab geometry assumes out-of-plane gradients are negligible, and the authors acknowledge this 'may introduce biases in volume-integrated quantities such as the mass and momentum accretion rates and the bremsstrahlung luminosity.' Since the fitting formulas (Eqs. 7–11) are precisely for these volume-integrated quantities, the 2D limitation directly affects the quantitative reliability of the paper's main output. The authors acknowledge this qualitatively in §4, but the fitting formulas are presented without any error bars or quantitative discussion of the expected systematic offset from 3D. At minimum, the manuscript should state explicitly in §3.2–3.3 that the fitting coefficients carry an unquantified systematic uncertainty from the 2D geometry, and the formulas'适用范围"},{"response":"This is a fair and straightforward point. We will add goodness-of-fit metrics for all five fitting formulas. Specifically, we will compute and report the R² value and the maximum absolute residual (in dex, since the formulas are in log space) for each of Eqs. (7)–(11). We will present these in a small table in §3.2. We will also add a residual plot (or a panel in Figure 3) showing the deviation of the fitting formula from the simulation data as a function of the relevant parameters, so that readers can identify where the largest residuals occur. This will allow users of the formulas to assess their reliability quantitatively.","revision_made":"yes","referee_comment":"§3.2, Eqs. (7)–(11): No goodness-of-fit metric (R², χ², RMS residual) is reported for any of the five fitting formulas. Without these metrics, a reader cannot assess whether the formulas provide an adequate description of the simulation data or identify where the largest residuals occur. This is essential for formulas that are intended for use in population synthesis codes. Please add a quantitative measure of fit quality, ideally with a residual plot or table showing the maximum deviation for each formula."}],"tokens_in":16303,"tokens_out":1455,"duration_ms":150635,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"The headline result: this is the first general-relativistic hydrodynamic study of Bondi-Hoyle-Lyttleton accretion in a common-envelope setting with a rotating (Kerr) black hole. It extends Cruz-Osorio & Rezzolla (2020), which used Schwarzschild, by adding spin and Mach number variation. The simulations are run with BHAC, a well-tested GRHD code, using HLLE fluxes and second-order TVD reconstruction. The test-fluid approximation is well-justified by the compactness ratio (C_RSG/C_BH ~ 10^-8). The flow morphologies — shock cones, bow shocks, subsonic regions, frame-dragging distortions — are all consistent with prior BHL literature, and the qualitative trends (shock-cone deflection by density gradients, Mach-number dependence of opening angle) are sound. The parameter survey of 66 models is a real effort, and the density-gradient prescription is physically motivated. The authors are honest about the 2D limitation and position the work as a first step toward 3D. That framing is appropriate and the qualitative results have value for the community as a roadmap for what 3D studies should check. The main deliverable is the set of fitting formulas (Eqs. 7–11) for accretion rates and bremsstrahlung luminosity as functions of spin, Mach number, and density gradient. Here is where the central problem lies. The 66-model suite is structured as two non-overlapping slices: Table 1 varies spin and density gradient at fixed Mach number M=2, and Table 2 varies Mach number and density gradient at fixed spin a*=0. No simulation has both nonzero spin and Mach number different from 2. Yet every fitting formula assumes multiplicative separability between the spin-dependent factor and the Mach-dependent factor. If frame-dragging interacts with the shock-cone geometry in a way that depends on Mach number — which is physically plausible, since the subsonic region size and shock-cone angle both depend on M — the true dependence would contain irreducible cross-terms like a*×M that the formulas cannot capture. The formulas would then be qualitatively wrong in the joint parameter space, not just systematically offset. The authors do not acknowledge this. The 2D geometry concern is real but acknowledged and affects systematic offsets; the separability concern is unacknowledged and affects whether the formulas are correct at all in the regions they claim to cover. Additionally, no goodness-of-fit metrics (residuals, R², error bars) are provided for the 4–7 parameter fits, and no hold-out validation is performed. These are minor relative to the separability issue but compound it. The paper is a competent first-pass study with a genuinely new parameter space exploration. The qualitative results are reliable. The fitting formulas should not be used quantitatively in population synthesis codes without 3D validation and at least a few models with simultaneous nonzero spin and non-default Mach number to test separability. The paper deserves a serious referee who should require the authors to explicitly state the separability assumption and its consequences, and to add goodness-of-fit metrics. I lean toward revise-and-resubmit rather than outright rejection — the simulation work is solid and the qualitative findings stand on their own.","headline":"First GRHD study of BHL accretion in a common envelope with a Kerr (spinning) black hole, extending prior Schwarzschild work. The fitting formulas are the main deliverable but have a structural problem the authors don't acknowledge.","tokens_in":17184,"tokens_out":773,"would_cite":false,"duration_ms":104925,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Density gradients dominate black hole spin in common-envelope accretion","keywords":["common envelope evolution","Bondi-Hoyle-Lyttleton accretion","Kerr black hole","general relativistic hydrodynamics","mass accretion rate","bremsstrahlung luminosity","density gradient","frame dragging"],"falsifier":"Run the same simulations in full 3D for a subset of the parameter space and check whether the fitting formulas reproduce the accretion rates and whether density gradients still dominate over spin and Mach number.","tokens_in":16389,"feed_emoji":"","tokens_out":1174,"duration_ms":154608,"temperature":0.7,"pith_summary":"This paper presents the first general-relativistic hydrodynamic simulations of a spinning (Kerr) stellar-mass black hole accreting from the envelope of a red supergiant star, in the Bondi-Hoyle-Lyttleton regime. Across 66 two-dimensional models varying black hole spin, Mach number, and envelope density gradient, the authors find that the density gradient in the stellar envelope has a more significant impact on the mass accretion rate than either the black hole's spin or the flow's Mach number. The authors provide analytical fitting formulas (Eqs. 7-11) with fitted coefficients that encode the dependence of mass accretion rate, radial and angular momentum transfer rates, and bremsstrahlung luminosity on these three parameters. The simulations also show that density gradients deflect the downstream shock cone away from the stellar core, that bow shocks form upstream and create subsonic regions, and that frame-dragging from black hole spin produces measurable but comparatively modest distortions near the event horizon. The fitting formulas are offered as tools for population synthesis models of binary evolution, though the authors note that three-dimensional calibration is needed before direct observational application.","feed_headline":"Density gradients beat black hole spin in envelope accretion","feed_subtitle":"First relativistic simulations of a Kerr black hole in a common envelope show that the star's density structure, not the hole's spin, setsac","key_machinery":"The central objects are the analytical fitting formulas (Eqs. 7-11), which express the logarithm of each normalized accretion quantity as a product of polynomial functions of the dimensionless spin parameter a★, the density gradient parameter ερ, and the Mach number M. The simulations use the BHAC code to solve the general-relativistic hydrodynamics equations in the Valencia formulation on a fixed Kerr spacetime (test-fluid approximation), with an ideal-gas equation of state (γ=5/3) and an exponentially stratified envelope density profile.","core_discovery":"The central finding is that, for a black hole embedded in a common envelope, the local density gradient in the envelope is the dominant driver of variations in mass accretion rate, momentum transfer, and thermal luminosity — more so than the black hole's spin or the Mach number of the relative flow. This is quantified through five fitting formulas that express each quantity as a separable product of functions of spin, density gradient, and Mach number, with fitted coefficients determined from the 66-model parameter survey. A secondary finding is that the shock cone morphology is controlled primarily by the density gradient (which deflects the cone) and the Mach number (which sets the opening","pith_inferences":["If density gradients dominate accretion in 2D, the three-dimensional case may amplify this effect, since real envelopes have gradients along all spatial directions rather than just one — the 2D result could be a lower bound on the influence of stratification.","The test-fluid approximation is well-justified by the compactness ratio, but if the black hole accretes enough mass to grow significantly during the common-envelope phase, the back-reaction on spacetime might become non-negligible over secular timescales, potentially modifying the fitting formulas for late-stage evolution.","The separability of the fitting formulas (each quantity is a product of independent functions of spin, gradient, and Mach number) may break down in regimes where nonlinear coupling between these parameters becomes important, such as very high spin combined with steep gradients — the 66-model grid may not sample these corners adequately."],"forward_implications":["Population synthesis models of binary evolution can use the fitting formulas to estimate black hole growth rates during common-envelope phases without running full hydrodynamic simulations, provided the three-dimensional corrections are calibrated.","The finding that density gradients dominate over spin suggests that simplified common-envelope models that neglect envelope stratification may systematically misestimate accretion rates and momentum transfer.","The proportionality between bremsstrahlung luminosity and mass accretion rate, if it holds in 3D, provides a way to infer accretion rates from thermal signatures of shocked gas, even though the photons themselves are trapped and reprocessed in the optically thick envelope.","The shock-cone deflection angle could serve as a diagnostic of the local density gradient in systems where the cone morphology is observable or inferable."],"fun_headline_variants":["Envelope density gradient dominates over black hole spin in accretion","Spin matters less than envelope structure for black hole accretion rates","Density gradients steer shock cones around spinning black holes","Black hole spin takes back seat to envelope density in accretion dynamics","Five formulas link accretion rates to spin, Mach number, and density gradients"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The two-dimensional slab geometry assumes that density and pressure gradients perpendicular to the simulation plane are negligible. Since the paper's main deliverables — the fitting formulas for mass accretion rate, momentum rates, and luminosity — are volume-integrated quantities, neglecting the third spatial dimension could bias the absolute values and potentially the relative ranking of which parameter matters most.","fun_headline_variants_meta":{"raw":{"variants":["Envelope density gradient dominates over black hole spin in accretion","Spin matters less than envelope structure for black hole accretion rates","Density gradients steer shock cones around spinning black holes","Black hole spin takes back seat to envelope density in accretion dynamics","Five formulas link accretion rates to spin, Mach number, and density gradients"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":623,"prompt_tokens":534,"completion_tokens":89,"prompt_tokens_details":null},"tokens_in":534,"tokens_out":89,"duration_ms":40260,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T10:07:10.109395+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"Run the same simulations in full 3D for a subset of the parameter space and check whether the fitting formulas reproduce the accretion rates and whether density gradients still dominate over spin and Mach number.","supporting_citations":[],"review_version":1}