{"id":"1a0e26fe-9714-4633-b53c-faf8b0538d41","arxiv_id":"2411.18545","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"GRMHD simulations of magnetized advective accretion onto a 20 solar mass black hole produce outflow power in the ULX range, supporting the magnetized hard-state ULX model.","lead":"This paper uses GRMHD simulations to test whether highly magnetized gas falling onto a 20-solar-mass black hole can power ultraluminous X-ray sources (ULXs). The simulated outflows carry power comparable to ULX luminosities, with magnetic fields near 10^7 Gauss, and a rotating-observer frame separates spin-powered and disk-wind components.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The simulated power is mechanical only; without radiative cooling the link to observed ULX X-ray luminosity is unestablished, and the claimed ULX-range normalization inherits the arbitrary Mdot_phy=0.05 Mdot_Edd, MBH=20 Msun choices.","rationale":"Read in good faith: this is a short proceedings paper presenting standard GRMHD simulations, and the cross-code comparison with Raha et al. (2024) is a genuine robustness check. The central issue is not internal inconsistency in the simulation itself but the strength of the astrophysical inference. The most load-bearing gap is the mechanical-to-radiative conversion: without radiative cooling, the computed outflow power is not the observed ULX X-ray luminosity. The second gap is the arbitrary physical normalization, which controls whether the dimensionless output actually lands in the ULX band. Both are identified in the reader's weakest_assumption, and I agree with that assessment. The paper's own Section 3.2 labels the power as only a precursor or upper bound, so the appropriate verdict remains conditional rather than reject. If a radiative post-processing test places L_X in the ULX range and the normalization sensitivity is modest, the central claim would be supported; if not, the abstract and conclusion overclaim.","tokens_in":6086,"tokens_out":6449,"duration_ms":64137,"concrete_test":"Take the final ~5000 rg/c of the MAD and SANE runs and post-process the time-averaged snapshots with a radiative transfer or cooling module (or rerun with a cooling function for synchrotron, Bremsstrahlung, and Compton emission) to compute the emergent 0.1-10 keV luminosity L_X, using the same Mdot_phy=0.05 Mdot_Edd, MBH=20 Msun normalization. If L_X falls below 3e39 erg/s, the mechanical power alone does not establish ULX-like luminosity. In the same post-processing, recompute P with Mdot_phy=0.01 Mdot_Edd and MBH=10 Msun; if the power leaves the ULX band, the claimed ULX-range result is normalization-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim ('high outflow power, well within the observed ULX luminosity range', Section 6) requires two independent conversions to hold. First, the code-unit energy flux in Eqs. (5)-(6) is a mechanical energy flux (thermal, kinetic, and Poynting) from an ideal GRMHD simulation, not a radiated luminosity. Section 3.2 explicitly states that radiative cooling is absent, so the computed power is 'only the precursor to ULX luminosities or the upper bound.' A mechanical outflow of 10^40 erg/s can be radiatively inefficient (optically thin, adiabatic expansion, or particle-dominated) and emit well below 3e39 erg/s in X-rays; hard-state ULX observations are X-ray luminosities, not kinetic powers. Second, the physical normalization P = [(Mdot-Edot)/Mdot(r')] Mdot_phy c^2 uses hand-chosen Mdot_phy = 0.05 Mdot_Edd and MBH = 20 Msun. Because the power scales linearly with both, a lower accretion rate or smaller black hole drops the same dimensionless profiles below the ULX band. Neither conversion is tested against a specific observed hard-state ULX, so the inference from simulation to ULX luminosity is not yet established. The paper partially acknowledges this in Section 3.2, but the abstract and Section 6 present the ULX-range power as a finding without that caveat.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the BHAC general-relativistic MHD code to simulate axisymmetric, magnetized advective accretion flows around a Kerr black hole with spin a=0.9375, starting from SANE and MAD initial conditions with initial plasma beta=100. It computes mass accretion rates, inward energy fluxes, and net outflow power profiles, normalizes these to physical power using Mdot_phy=0.05 Mdot_Edd and MBH=20 Msun, and reports that the resulting powers lie within the ULX range 3e39-3e41 erg/s. It further reports magnetic field strengths of order 1e7 G, compares the MAD run with an independent HARMPI simulation, and presents ZAMO-frame power profiles attributed to Blandford-Znajek and Blandford-Payne mechanisms. The central claim is that highly magnetized advective accretion flows can explain hard-state ULX luminosities without invoking intermediate-mass black holes or modified Eddington limits.","tokens_in":6431,"tokens_out":5853,"duration_ms":49080,"significance":"If fully established, the result would provide numerical support for the Mondal-Mukhopadhyay (MM19) steady-state model of hard-state ULXs and would strengthen the case that stellar-mass black holes with highly magnetized advective disks can power ULX-like outputs. The paper has several genuine strengths: it uses standard GRMHD methods and publicly available code, it performs a cross-code robustness check with HARMPI, it presents multiple time-averaging definitions of power, and the 1e7 G field strength is an emergent outcome rather than an input. The dimensionless power profiles are also a useful quantity for future comparisons. However, the step from simulated mechanical power to observed ULX X-ray luminosity is not established, and the physical normalization relies on hand-picked values; these issues are load-bearing for the paper's main claim.","major_comments":[{"comment":"The translation from code units to physical power uses hand-picked values Mdot_phy=0.05 Mdot_Edd and MBH=20 Msun. The reported power scales linearly with both quantities, so a lower accretion rate or a smaller black hole mass places the same dimensionless profiles below the ULX band (3e39 erg/s). The paper should present the dimensionless outflow efficiency P/(Mdot c^2) as the primary simulation output and either display the physical power as a function of the assumed Mdot_phy and MBH or clearly state that the ULX-range normalization is an assumption, not a prediction of the simulation.","section":"§3.2, Eq. (6), Fig. 4"},{"comment":"The simulations do not include radiative cooling, and Section 3.2 explicitly states that the computed outflow power is 'only the precursor to ULX luminosities or the upper bound.' Yet the Abstract claims the systems 'produce high luminosities like ULXs' and Section 6 claims 'high outflow power, well within the observed ULX luminosity range.' The simulated quantity is mechanical power (thermal, kinetic, and Poynting), which can be radiatively inefficient and need not emerge as X-ray luminosity. No specific hard-state ULX is compared with the model, and no radiative efficiency or emission mechanism is supplied. Please qualify the abstract and conclusion or add an explicit conversion from mechanical power to observable X-ray luminosity.","section":"§3.2 and §6 (also Abstract)"},{"comment":"The ZAMO-frame energy flux is not defined unambiguously. The contraction T^mu_nu u^nu e^mu has mismatched indices, and the stated radial vector e^mu=(0,1/grr,0,0) is not a unit vector in the Kerr metric unless grr is defined carefully; the unit radial vector should be e^r=1/sqrt(g_rr) or equivalently sqrt(g^rr). In addition, the jet region criterion '-T^mu_nu > 0' is not a well-defined scalar condition; the relevant component and observer frame should be specified. These issues affect the claimed decomposition into Blandford-Znajek and Blandford-Payne contributions.","section":"§5.1, Eq. (7)"}],"minor_comments":[{"comment":"Please fix typographical issues: 'around107 G' in the Abstract should be 'around 10^7 G', and 'Blanford-Znajek' in Section 1 should be 'Blandford-Znajek'.","section":"Abstract, §1"},{"comment":"The setup text says 'MBM is the mass of the black hole'; this should read 'MBH'. Also, 'Fishbone Moncrief (FM) tours setup' should be 'torus setup'.","section":"§2"},{"comment":"The caption for Fig. 4 states that 'all definitions' of time averaging are shown, but the text does not specify which line type corresponds to which of the three definitions in the list. Please identify each curve explicitly.","section":"§3.2, Fig. 4"},{"comment":"The steady-flow radii req=10 rg and req=20 rg are attributed to an in-preparation paper [13]. Since Fig. 3 displays the red dashed lines that define these radii, the values are directly supported by the present paper; citing an unpublished work for them is unnecessary and makes verification harder.","section":"§3.1, Ref. [13]"},{"comment":"The sentence 'The jet region is considered to be the part of the simulation domain in which -T^mu_nu > 0 [17]' is unclear as written; please specify whether this is a particular component such as -T^r_t, and in which basis or observer frame the condition is applied.","section":"§5.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings contribution and its scope is modest, but the central claim as stated in the abstract and conclusion overreaches the simulation evidence. The authors themselves acknowledge in Section 3.2 that radiative cooling is absent, which is the key limitation; the abstract and conclusion should be brought in line with that caveat. The dependence of the ULX-range claim on the arbitrary Mdot_phy and MBH choices should also be made transparent. The cross-code check and the emergent field-strength result are valuable, so I see the work as salvageable through revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll get straight to it. The paper does two useful things. It runs BHAC GRMHD simulations of magnetized advective accretion flows around a stellar-mass Kerr black hole (spin 0.9375) in SANE and MAD states, and shows that the outflow power, normalized to an advective accretion rate, lands in the ULX luminosity range. That is a genuine GRMHD check of the MM19 steady-state model, and it comes with a cross-code robustness check using HARMPI, which is more than most proceedings papers bother with. The second genuinely new piece is the ZAMO-frame power decomposition: the profiles show a peak near the horizon that weakens as spin decreases (the BZ signature) and a second peak around 30 rg attributed to the BP mechanism. That's a nice diagnostic and worth having.\n\nThe paper is also transparent about its main limitation, in the body. Section 3.2 states explicitly that there is no radiative cooling and the computed power is 'only the precursor to ULX luminosities or the upper bound.' The problem is that the abstract and conclusion don't carry that caveat. They say the systems produce 'high luminosities like ULXs' and 'high outflow power, well within the observed ULX luminosity range.' That is an overclaim. The simulated quantity is mechanical energy flux—thermal, kinetic, and Poynting—not radiated X-ray luminosity. A mechanical outflow of 1e40 erg/s can emit far less in X-rays if it is radiatively inefficient. The connection to observed ULX luminosity is not yet established.\n\nThe physical normalization is also a hand choice: Mdot_phy = 0.05 Mdot_Edd and MBH = 20 Msun. The power scales linearly with both, so lower values drop the profiles below the ULX band. The paper does not test sensitivity to these choices, and the time-averaged power definitions carry no error bars. That said, the magnetic field strength of ~1e7 G is an emergent outcome of the simulation, not an input, which supports the MM19 field estimate.\n\nOn the citation pattern, the main dependencies are MM19 for the model and an in-preparation paper for the req values; the req values are also displayed in Fig. 3, so the missing reference is not load-bearing. Self-citation here is not a red flag.\n\nVerdict: this is a modest but honest contribution. The core simulation work is solid, the ZAMO diagnostic is useful, and the cross-code check helps. The central ULX claim is defensible only as an upper bound on mechanical power. That is not fatal, but the abstract needs to match the body. I'd take it for a reading group if anyone is working on ULX jet models, and I'd cite the ZAMO decomposition if I needed a clean way to separate BZ and BP in simulation data. A serious editor should send it to review; the authors should be asked to fix the abstract and add a sentence about radiative efficiency.","headline":"GRMHD check of MM19 is real and the ZAMO BZ/BP decomposition is useful, but the abstract overstates mechanical outflow power as ULX luminosity; the body admits it is an upper bound.","tokens_in":6949,"tokens_out":2980,"would_cite":true,"duration_ms":25488,"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":"Highly magnetized advective accretion flows around stellar-mass black holes can produce outflow power in the observed ULX luminosity range, according to GRMHD simulations.","keywords":["ultraluminous X-ray sources","hard state","GRMHD simulations","magnetized advective accretion flows","Blandford-Znajek mechanism","Blandford-Payne mechanism","magnetically arrested disks","outflow power"],"falsifier":"Measure the accretion rate of a hard-state ULX with an independently determined black-hole mass; if the rate is substantially below 0.05 Eddington for a ~20 solar-mass black hole, the normalized outflow power falls below the observed ULX band, contradicting the model.","tokens_in":5877,"feed_emoji":"🌌","tokens_out":8321,"duration_ms":63784,"temperature":0.7,"pith_summary":"This paper argues that hard-state ultraluminous X-ray sources (ULXs) can be explained as highly magnetized, advective accretion flows around stellar-mass black holes, without invoking intermediate-mass black holes or modified Eddington limits. Using general relativistic magnetohydrodynamic (GRMHD) simulations of such flows, the authors show that the outflow power they produce falls inside the observed ULX luminosity range, and that the magnetic fields near the black hole reach about $10^7$ gauss. They also decompose the outflow power in the frame of a zero angular momentum observer and identify two components, one tied to black-hole spin (Blandford-Znajek) and one tied to magnetocentrifugal launching from the disk (Blandford-Payne). Because the simulations omit radiative cooling, the computed power is presented as a precursor or upper bound to the X-ray luminosity rather than the luminosity itself.","feed_headline":"Magnetized black-hole disks can power ultraluminous X-ray outflows","feed_subtitle":"Simulations show a 20-solar-mass black hole with a magnetized disk reaches the observed ULX luminosity band.","key_machinery":"The load-bearing object is the outflow power $P(r)=\\dot{M}(r)-\\dot{E}(r)$, computed from the stress-energy tensor of the GRMHD flow, where $\\dot{M}$ is the mass accretion rate and $\\dot{E}$ the inward energy flux. It is normalized to physical units by multiplying the dimensionless ratio by $\\dot{M}_{\\rm phy} c^2$, with $\\dot{M}_{\\rm phy}=0.05\\,\\dot{M}_{\\rm Edd}$ and a 20-solar-mass black hole. A second object is the ZAMO-frame radial energy flux, which separates the power into a spin-dependent inner peak and an outer disk-launched peak. Together these definitions carry the argument: the first places the simulated flows in the ULX band, and the second attributes the two components to the Blandford-Znajek and Blandford-Payne mechanisms.","core_discovery":"The central discovery claimed is that a 20-solar-mass black hole accreting at $0.05$ of the Eddington rate, surrounded by a highly magnetized advective disk, produces time-averaged outflow power in the ULX band. The power is defined as the difference between the mass accretion rate and the inward energy flux, $P(r)=\\dot{M}(r)-\\dot{E}(r)$, normalized by $\\dot{M}_{\\rm phy} c^2$. In both the standard (SANE) and magnetically arrested (MAD) initial magnetic configurations, the authors find powers in the ULX range, with MAD flows giving higher power; the required magnetic field at the black hole is about $10^7$ gauss, consistent with earlier steady-state calculations. In the ZAMO frame, the radial power profile shows a peak that grows with black-hole spin and ends near the ergosphere, interpreted as Blandford-Znajek spin extraction, and a second peak around $30$ gravitational radii, interpreted as Blandford-Payne magnetocentrifugal outflow. The authors state that the absence of radiative cooling means the simulated power is only the precursor to, or upper bound on, ULX luminosities.","pith_inferences":["If the mechanical outflow power is converted to radiation with even moderate efficiency, these flows could also account for the kinetic power of ULX jets and winds, not only the X-ray band.","The spin dependence of the inner ZAMO peak predicts that higher-spin sources should show stronger BZ-dominated jet power; retrograde-spin simulations could cleanly isolate the BP component.","Because the power scales linearly with the adopted accretion rate and black-hole mass, the model makes a quantitative prediction: a ULX with an accretion rate an order of magnitude lower at the same mass would fall out of the ULX band.","The 2.5-dimensional axisymmetric setup likely exaggerates the magnetic-barrier oscillations in MAD runs; full 3D simulations may reduce the MAD-versus-SANE power gap."],"forward_implications":["Hard-state ULXs can be modeled without intermediate-mass black holes or a modified Eddington limit.","Magnetic field strengths near $10^7$ gauss suffice for ULX-level outflow power, without invoking super-Eddington fields far from the black hole.","Magnetically arrested disks give higher outflow power than standard (SANE) disks, making MAD flows a promising engine for the most luminous hard-state ULXs.","ZAMO-frame power profiles separate spin-powered (Blandford-Znajek) from disk-launched (Blandford-Payne) outflows, providing a diagnostic for jet-launching mechanisms.","Because the simulations exclude radiative cooling, the computed power is an upper bound to the eventual X-ray luminosity, not the luminosity itself."],"supporting_citations":[{"why":"Supplies the steady-state model of hard-state ULXs as magnetized advective flows that this work verifies with GRMHD simulations.","marker":"[7]"},{"why":"Provides the BHAC GRMHD code used for all simulations.","marker":"[8]"},{"why":"Provides the analytic torus initial condition used to set up the accretion disk.","marker":"[9]"},{"why":"Defines the outflow power $P(r)=\\dot{M}-\\dot{E}$ and the SANE/MAD initial magnetic vector potentials used here.","marker":"[10]"},{"why":"Independent GRMHD simulation of the same setup used to check robustness across codes and computational domains.","marker":"[11]"},{"why":"Defines the steady-flow radius used to choose the equilibrium region for power extraction.","marker":"[12]"},{"why":"Documents the high efficiency of MAD flows, used to explain why MAD simulations have higher outflow power.","marker":"[14]"},{"why":"Names the Blandford-Znajek spin-extraction mechanism used to interpret the inner ZAMO power peak.","marker":"[15]"},{"why":"Names the Blandford-Payne magnetocentrifugal mechanism used to interpret the outer ZAMO power peak.","marker":"[16]"}],"fun_headline_variants":["Magnetized black hole disks hit ULX luminosity","Simulations show magnetized disks reach ULX levels","Black hole spin and magnetized disk produce ULX","Magnetized flows around black holes can shine as ULXs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a hard-state ULX accretes at about 5% of the Eddington rate onto a 20-solar-mass black hole, so the mechanical outflow power, computed with radiative cooling omitted, can stand in for the X-ray luminosity.","fun_headline_variants_meta":{"raw":{"variants":["Magnetized black hole disks hit ULX luminosity","Simulations show magnetized disks reach ULX levels","Black hole spin and magnetized disk produce ULX","Magnetized flows around black holes can shine as ULXs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00145,"raw_usage":{"total_tokens":5860,"prompt_tokens":986,"completion_tokens":4874,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":4809}},"tokens_in":602,"tokens_out":4874,"duration_ms":36416,"temperature":1.0,"reasoning_tokens":4809,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:05:36.528697+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the accretion rate of a hard-state ULX with an independently determined black-hole mass; if the rate is substantially below 0.05 Eddington for a ~20 solar-mass black hole, the normalized outflow power falls below the observed ULX band, contradicting the model.","supporting_citations":[{"cited_title":"Mondal, B","cited_arxiv_id":null,"evidence_quote":"Supplies the steady-state model of hard-state ULXs as magnetized advective flows that this work verifies with GRMHD simulations."},{"cited_title":"Porth, H","cited_arxiv_id":null,"evidence_quote":"Provides the BHAC GRMHD code used for all simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the analytic torus initial condition used to set up the accretion disk."},{"cited_title":"Chatterjee, R","cited_arxiv_id":null,"evidence_quote":"Defines the outflow power $P(r)=\\dot{M}-\\dot{E}$ and the SANE/MAD initial magnetic vector potentials used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Independent GRMHD simulation of the same setup used to check robustness across codes and computational domains."},{"cited_title":"Narayan, A","cited_arxiv_id":null,"evidence_quote":"Defines the steady-flow radius used to choose the equilibrium region for power extraction."},{"cited_title":"Tchekhovskoy, R","cited_arxiv_id":null,"evidence_quote":"Documents the high efficiency of MAD flows, used to explain why MAD simulations have higher outflow power."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Names the Blandford-Znajek spin-extraction mechanism used to interpret the inner ZAMO power peak."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Names the Blandford-Payne magnetocentrifugal mechanism used to interpret the outer ZAMO power peak."}],"review_version":1}