{"id":"adb08e1b-0648-4fbf-83d8-10683aa31bea","arxiv_id":"2505.10460","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"A non-conservative radiation field around a luminous inner corona generates magnetic fields up to roughly 10^5 Gauss, about five percent of equipartition, over the viscous infall time.","lead":"This paper calculates how eddies in the radiation field around a black hole accretion disk can generate magnetic fields above the disk. The authors find that with a bright, compact inner corona the fields could reach about 100,000 Gauss and become dynamically important.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 10^5 G claim rests on linear, v=0 induction equations integrated to the viscous time, a regime the paper itself says those equations do not cover.","rationale":"The reader identified the missing nonlinear saturation, and I agree. The contradiction is internal, not merely a question of consensus. Equation 9 is derived by dropping the Hall term as higher order and setting v = 0, and the paragraph after Eq. 13 explicitly limits Eqs. 11–13 to short timescales. Yet the abstract and Section 4.2 use t = t_v, the viscous infall time, to reach 10^5 G. A few percent of equipartition is by definition the regime where Lorentz and Hall terms can no longer be neglected, so the computed endpoint lies outside the validity domain of the equations actually solved. The motionless assumption is also in tension with adopting v_r from Eq. 24 as the relevant time: a radially moving fluid element should have a ∇×(v×B) term and a changing source along its trajectory. These are not small quantitative corrections; they determine whether growth saturates far below the claimed values or at all. The underlying mechanism, generation of seed fields by the curl of the radiation force, is physically plausible, and the Keplerian-disc result may survive, but the central quantitative claim requires a self-consistent nonlinear calculation. Since the headline result is unsupported by the paper's own equations, the REJECT verdict is appropriate and no adjustment to the reader's verdict is needed.","tokens_in":11660,"tokens_out":6484,"duration_ms":68411,"concrete_test":"Integrate Eq. 7 numerically for the fiducial case of Section 4.2 (10 Msun, m = 0.1, lc = 1, rs = 3rg) with the Hall term -(σ_c/(c e n_e))∇×((∇×B)×B) and the advection term ∇×(v×B) included, using v_r from Eq. 24 and a range of densities n_e = 10^10–10^16 cm^-3, and record Bmax(t). If the peak field saturates below ~10^5 G, or if reaching ~10^5 G requires t much longer than t_v(r = 3rg) with α = 0.1, the linear-extrapolation claim fails. An independent analytical check is to evaluate the Hall saturation field at r ~ 3rg and compare it with the claimed 10^5 G for the same n_e.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The 10^5 G result is obtained by inserting t = t_v (Eq. 24) into the linear solutions Eqs. 11–13. But those solutions follow from Eq. 9, which is derived under two explicit restrictions stated in Section 2: (a) the plasma is essentially at rest (v → 0), and (b) terms nonlinear in B, notably the Hall term in Eq. 7, can be dropped because B is small. The paper itself warns after Eq. 13 that Eqs. 11–13 are valid only for short timescales, and footnote 1 states that the Hall-term assumption fails once B becomes significant. The claimed endpoint is precisely that regime: a few percent of equipartition is dynamically significant. Using the viscous infall time also requires the plasma to move radially (Eq. 24), contradicting v → 0; in a moving plasma the induction equation should include ∇×(v×B) and an advected, time-varying source. Neither the saturation level of the Hall term nor the advective modification is calculated at any point. Thus the headline field strength is an extrapolation of a linear, stationary-plasma solution into a regime the equations of the paper do not cover, rather than a solution of the model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that magnetic fields around black hole accretion discs are generated by the non-conservative nature of the radiation force field, i.e., by ∇×F_rad ≠ 0. The authors compute the radiation flux components from a geometrically thin, optically thick Keplerian disc plus an optional inner corona, and then evaluate the linear induction equation dB/dt = -(c/e)∇×f_rad (Eq. 9) to obtain the time-dependent field components (Eqs. 11-13). For a disc without a corona they recover fields of a few Gauss, consistent with earlier work. When a compact, Eddington-luminosity inner corona (lc = 1, rs = 3rg) is added, they report maximum fields up to ~10^5 G, a few percent of equipartition, developing over the viscous infall time t_v (Eq. 24). The paper argues that such fields can be dynamically significant and may affect disc evolution, jet launching, and polarization signatures.","tokens_in":11918,"tokens_out":3607,"duration_ms":34855,"significance":"If the 10^5 G result were correct, the mechanism would provide a new, radiation-driven route to magnetizing the inner regions of black hole accretion flows, with potential implications for jet launching and X-ray polarization. The paper contains a relatively detailed numerical calculation of the radiation flux components and their derivatives, and it clearly states the model assumptions and the parameter dependences. However, the central quantitative claim rests on an extrapolation of a linear, stationary-plasma solution into a regime that the paper itself identifies as requiring additional terms. The significance is therefore conditional on whether that extrapolation can be justified, and the present manuscript does not provide such a justification.","major_comments":[{"comment":"The headline field strength of 10^5 G is obtained by inserting t = t_v (Eq. 24) into the linear solutions Eqs. (11)–(13), which are derived from Eq. (9) under the assumptions v → 0 and small B. However, the text immediately after Eq. (13) states that Eqs. (11)–(13) “are valid for short time scales,” and footnote 1 states that the neglect of the Hall term “fails at later times, when the magnitude of the magnetic field becomes significant.” The endpoint of the integration, a field at a few percent of equipartition, is precisely that regime. The paper provides no calculation of the Hall-term saturation or of the advective term ∇×(v×B), so the 10^5 G value is an unsupported extrapolation rather than a solution of the model.","section":"Sec. 2, Eqs. (11)–(13) and Sec. 4.2"},{"comment":"Equation (9) assumes a plasma that is essentially at rest (v → 0), so that the induction term ∇×(v×B) is negligible. In Sec. 3.1 the available growth time is instead taken as the viscous infall time t_v = R/v_r, with v_r given by Eq. (24). A plasma that advects radially on the viscous timescale is not at rest, and the induction equation for such a flow should contain ∇×(v×B) and an advection-modified source term. The paper does not solve this full equation; it linearly superposes a static-plasma source onto a moving-plasma timescale. This inconsistency directly affects the claimed growth time and final field strength.","section":"Sec. 2, Eq. (9), and Sec. 3.1, Eq. (24)"},{"comment":"The paper argues that the generated field reaches ~5% of equipartition and is therefore dynamically significant, affecting disc and jet evolution. At such field strengths the Lorentz force modifies the momentum balance of Eq. (1) that underlies the radiation source term, and the Hall term in Eq. (7) becomes comparable to the radiation source. The authors acknowledge this only in footnote 1 and do not estimate the saturation value or the timescale at which the linear approximation breaks down. Without a nonlinear or saturation calculation, the claim that such magnitudes develop within the viscous timescale is not established.","section":"Sec. 4.2, Fig. 3"}],"minor_comments":[{"comment":"The phrase “few percentage of equipartition” should be “a few percent of equipartition”; it appears in both the abstract and Sec. 4.2.","section":"Abstract and Sec. 4.2"},{"comment":"The first author's surname is typeset as “Vy as” instead of “Vyas.”","section":"Author list"},{"comment":"The monotonic scalings Bmax ∝ 1/r_s^2 and Bmax ∝ l_c shown in Fig. 5 are direct consequences of B ∝ t ∇×F with F ∝ L_c/A_c; the text should state explicitly that these are consistency checks rather than new predictions.","section":"Fig. 5 and Sec. 4.2"},{"comment":"The assumption of constant angular velocity for the inner corona (Eq. 20) is introduced with only a citation to McKinney & Narayan (2007); given that the corona here is a radiative, weakly magnetized region, a more quantitative justification is needed for this velocity profile.","section":"Eq. (20) and Sec. 3"},{"comment":"The conclusion claims that the vertical component Bz dominates, but the paper does not compare the relative magnitudes of Bz, Br, and Bφ quantitatively; a quantitative statement would strengthen the claim about jet launching.","section":"Sec. 5, Fig. 3"}],"recommendation":"reject","confidential_remarks":"The manuscript is a straightforward application of a known linear induction mechanism to a new source geometry (an inner corona with a specific velocity profile), and the numerical evaluation of the radiation flux derivatives is careful. The fatal issue is that the headline 10^5 G result is obtained by integrating the linear, v=0 equations over the viscous timescale, a regime the authors themselves state the equations do not cover. A revision that computes the nonlinear saturation or restricts the claims to the short-time, weak-field regime could be reconsidered, but the present version does not, in my view, meet the standard for publication. I also note that the parameter scalings in Fig. 5 are essentially predetermined by linearity, so the predictive content of the paper beyond the geometry is limited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper says a compact, luminous, rotating corona can generate poloidal magnetic fields near black hole accretion discs that reach ~10^5 G, a few percent of equipartition, on the viscous infall time. The underlying mechanism is not new—non-conservative radiation fields generating fields via curl f_rad is the cosmic battery idea from Bisnovatyi-Kogan & Blinnikov, Ando et al., Shiromoto et al. What is new is the explicit computation of the full radiation flux and its curl for a two-component disc: a standard Keplerian disc plus a rigidly rotating inner corona. For the Keplerian-only case they reproduce the old weak-fields result (a few Gauss). For the corona case they find a large amplification. That is a legitimate and physically motivated extension, and the plots and calculation seem careful.\n\nThe problem is the central quantitative claim. The field evolution is obtained from dB/dt = -(c/e) ∇×f_rad, which is Eq. 9. As the paper itself states, that equation assumes the plasma is essentially at rest and that the Hall term (∝ j×B) is negligible because B is small. The authors then integrate this linear equation up to the viscous timescale t_v, using a radial velocity that implies the plasma is moving. They write after Eq. 13 that Eqs. 11–13 are valid only for short timescales, and footnote 1 says the Hall-term assumption fails once B becomes significant. The claimed endpoint—~5% of equipartition—is exactly that regime. So the headline field strengths are a linear extrapolation into a regime the model's own equations do not cover. No saturation calculation, no advection term, no Hall-term estimate is given. This is not a minor detail; it is the load-bearing step for the dynamical-significance claim.\n\nSome smaller things: the corona rotation law v_ϕ ∝ r_d is an assumption (constant angular velocity) that affects the curl; the viscous timescale depends on the unspecified α; and the \"motionless plasma\" assumption is in tension with using t_v at all. But those are secondary. The main issue stands.\n\nTo their credit, the authors are explicit about the caveats—footnote 1 and the note after Eq. 13 are honest. The mechanism is plausible and the setup is a genuine extension. But as it stands, the paper does not justify the 10^5 G fields. A referee could reasonably ask for a saturation-level estimate by keeping the Hall term and the advective terms, or a time-dependent numerical integration of the full induction equation. That is a substantial revision, not a copy-edit.\n\nMy recommendation: send it to peer review. The idea is worth refereeing; the flaw is specific and fixable, and the authors have shown they know where the boundary lies. It should not be accepted in its current form, but a desk rejection would waste a useful extension.","headline":"The corona extension is a real idea, but the 10^5 G claim comes from applying a short-timescale, motionless-plasma equation for a full viscous time.","tokens_in":12437,"tokens_out":3776,"would_cite":false,"duration_ms":31911,"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":"The curl of radiation from a rotating inner corona can magnetize a black hole accretion disc to about 100,000 Gauss within the viscous timescale.","keywords":["magnetic field generation","radiation force","accretion disc corona","black hole accretion","charge separation","battery effect","equipartition","jet launching"],"falsifier":"Simulate the same configuration—a 10-solar-mass black hole, a thin disc at 0.1 Eddington accretion, and a corona at $r_s=3r_g$ radiating at the Eddington luminosity—using the full induction equation including Hall, advection, and back-reaction terms, and read off the maximum field at the viscous time; if it is not within an order of magnitude of $10^5$ G, the central claim fails.","tokens_in":11435,"feed_emoji":"🧲","tokens_out":11348,"duration_ms":100604,"temperature":0.7,"pith_summary":"This paper claims that radiation itself can be the engine that magnetizes black hole accretion discs. Because the radiation force on the plasma is not a pure gradient, its curl drives charge separation, and the induction equation grows a magnetic field that is linear in time. For a standard thin Keplerian disc the effect is weak, producing only a few Gauss. With a compact inner corona radiating near the Eddington limit, however, the same mechanism produces fields of order $10^5$ Gauss—a few per cent of equipartition with the gas pressure—within the viscous infall time. If correct, this radiation-driven battery is a significant magnetization channel alongside the magneto-rotational instability, with direct consequences for jet launching and for X-ray polarization.","feed_headline":"Inner corona drives black-hole disc fields to 100,000 Gauss","feed_subtitle":"Radiation's curl separates charge above the disc and yields near-equipartition fields within the viscous timescale.","key_machinery":"The load-bearing object is the curl of the radiation force, $\\nabla\\times\\mathbf{f}_{\\rm rad}$, which acts as a battery source in the induction equation. For Thomson scattering, $\\mathbf{f}_{\\rm rad}=\\sigma_T\\mathbf{F}/c$, so the source is proportional to $\\nabla\\times\\mathbf{F}$, the rotation of the radiation flux. The calculation evaluates the three flux components and all six spatial derivatives by integrating the Lorentz-boosted specific intensity of the disc and corona over a thin disc plane, using a pseudo-Newtonian Keplerian profile for the outer disc and a constant-angular-velocity, solid-body profile for the inner corona. The resulting field components grow linearly with time as $(\\sigma_T/e)\\,t$ times combinations of flux derivatives. The key geometric fact is that the rotating corona at small radius produces a much larger flux curl than the Keplerian disc alone, which is what amplifies the field by orders of magnitude.","core_discovery":"The paper's central claim is that the non-conservative component of the radiation field above an accretion disc is a viable magnetic-field source, and that an inner corona makes it dynamically important. Starting from the electron momentum equation and Maxwell's equations, the authors derive the source term $d\\mathbf{B}/dt = -(c/e)\\nabla\\times\\mathbf{f}_{\\rm rad}$, with $\\mathbf{f}_{\\rm rad} = \\sigma_T\\mathbf{F}/c$. They then compute all components of the radiation flux and their coordinate derivatives for a geometrically thin disc consisting of an outer Keplerian part and an inner corona with solid-body rotation. A bare Keplerian disc yields fields of a few Gauss, while adding an Eddington-luminous corona inside $r_s\\simeq 3\\,r_g$ raises the peak field to about $10^5$ G near the inner disc, roughly 5% of the equipartition value, on the viscous timescale. The paper presents this as a demonstration that radiation-driven fields from a corona can reach dynamically significant strengths in realistic times.","pith_inferences":["A direct test of the paper's rotation assumption would be to rerun the flux calculation with a Keplerian inner corona ($v_\\phi\\propto r_d^{-1/2}$) instead of the solid-body profile; if the $10^5$ G result is materially reduced, the mechanism hinges on the coronal velocity law.","The paper's linear-growth caveat implies that the physical end state may be a field pinned near the few-per-cent equipartition level rather than continuing to grow; even at that level it could act as a coherent seed for magneto-rotational turbulence in the inner disc.","Applying the same battery to an active galactic nucleus, where coronae are routinely observed and viscous times are much longer, would suggest larger integrated fields—but the saturation question becomes more important because there is more time for back-reaction to act.","Observationally, a steady, axisymmetric vertical field near the disc–corona interface should produce stable X-ray polarization aligned with the jet axis, distinguishing this mechanism from turbulent fields; this prediction is testable with current and upcoming X-ray polarimeters."],"forward_implications":["A luminous inner corona at $r_s\\simeq 3$–$10\\,r_g$ makes radiation-driven fields dynamically significant, reaching about $10^5$ G and a few per cent of equipartition within the viscous timescale.","The corona-driven field is dominated by the vertical component $B_z$ near the disc–corona interface, the geometry that is a prerequisite for magnetically launched outflows and jets.","The produced field scales linearly with corona luminosity and inversely with the square of corona radius, so compact and Eddington-luminous coronae are the configurations in which this mechanism matters.","Because these fields decay faster with height than magneto-rotational-instability fields, their observational imprint would be inner-disc spectral breaks or radial variations in X-ray polarization rather than a volume-filling disc field."],"supporting_citations":[{"why":"Establishes the Poynting–Robertson radiation battery and the few-Gauss baseline for a Keplerian disc that this paper extends, and supplies the viscous infall velocity used for the growth timescale.","marker":"Bisnovatyi-Kogan & Blinnikov 1977"},{"why":"Generalizes the radiation battery to finite conductivity and finite disc extent; the standard-disc result this work compares against.","marker":"Bisnovatyi-Kogan et al. 2002"},{"why":"Introduces the cosmic battery and defines the equipartition field strength that the corona case approaches within a few per cent.","marker":"Contopoulos & Kazanas 1998"},{"why":"Provides the short-timescale induction equation with the curl-of-radiation-force source term that is the paper's starting point.","marker":"Ando et al. 2010"},{"why":"Applies the same radiation-curl battery in another astrophysical context, supporting the general validity of the source term.","marker":"Shiromoto et al. 2014"},{"why":"Supplies the standard thin-disc intensity profile used to compute the Keplerian disc radiation flux.","marker":"Shakura & Sunyaev 1973"},{"why":"Gives the pseudo-Newtonian potential from which the Keplerian and coronal velocity profiles are derived.","marker":"Paczyński & Wiita 1980"},{"why":"Motivates the constant-angular-velocity inner corona, the assumption that produces the large flux curl.","marker":"McKinney & Narayan 2007"},{"why":"Radiation-MHD simulations motivating the homogeneous specific intensity assumed for the corona.","marker":"Jiang et al. 2019"}],"fun_headline_variants":["Corona's radiation curl yields 100,000 Gauss fields","Radiation curl magnetizes black hole accretion discs","Non-conservative radiation drives strong black hole disc fields","Radiation field curl seeds near-equipartition fields in discs","Magnetic fields from radiation curl reach 100,000 Gauss with corona"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes linear growth continues for the full viscous infall time, even though the paper states its growth equation is only valid while the field is weak; if back-reaction slows growth once the field becomes dynamically significant, the predicted peak strength is an overestimate.","fun_headline_variants_meta":{"raw":{"variants":["Corona's radiation curl yields 100,000 Gauss fields","Radiation curl magnetizes black hole accretion discs","Non-conservative radiation drives strong black hole disc fields","Radiation field curl seeds near-equipartition fields in discs","Magnetic fields from radiation curl reach 100,000 Gauss with corona"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000927,"raw_usage":{"total_tokens":4017,"prompt_tokens":1035,"completion_tokens":2982,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":651,"completion_tokens_details":{"reasoning_tokens":2898}},"tokens_in":651,"tokens_out":2982,"duration_ms":22321,"temperature":1.0,"reasoning_tokens":2898,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:08:35.211989+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the same configuration—a 10-solar-mass black hole, a thin disc at 0.1 Eddington accretion, and a corona at $r_s=3r_g$ radiating at the Eddington luminosity—using the full induction equation including Hall, advection, and back-reaction terms, and read off the maximum field at the viscous time; if it is not within an order of magnitude of $10^5$ G, the central claim fails.","supporting_citations":[{"cited_title":"S., Lovelace, R","cited_arxiv_id":null,"evidence_quote":"Generalizes the radiation battery to finite conductivity and finite disc extent; the standard-disc result this work compares against."},{"cited_title":"2010, ApJ, 716, 1566, doi: 10.1088/0004-637X/716/2/1566","cited_arxiv_id":null,"evidence_quote":"Provides the short-timescale induction equation with the curl-of-radiation-force source term that is the paper's starting point."},{"cited_title":"2014, ApJ, 782, 108, doi: 10.1088/0004-637X/782/2/108","cited_arxiv_id":null,"evidence_quote":"Applies the same radiation-curl battery in another astrophysical context, supporting the general validity of the source term."}],"review_version":1}