{"id":"fb564e54-c5fd-4a9a-876e-13ff77feef42","arxiv_id":"2412.18772","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In the FRADO model, a photon-flux-dependent cloud emissivity makes broad-line profiles insensitive to the radial scaling of the disk mass-loss rate, while under uniform emissivity a steeper r^-3/2 scaling better matches observed single-peaked profiles.","lead":"This paper tests how the radial decline of mass loss from an accretion disk changes the shapes of broad emission lines in active galactic nuclei, using the FRADO dusty-outflow model. It finds that the assumed emissivity of the emitting clouds influences the line profiles more than the radial outflow scaling does.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim hinges on a binary photon-flux cut whose position and shape are untested; shifting or smoothing the threshold could restore or reverse the s=-0.5 vs s=-1.5 profile difference.","rationale":"The reader identified the binary photon-flux threshold as the weakest assumption, and I agree. This is the single point on which the central claim pivots: without the step-function cut, the two outflow scalings produce clearly different line profiles (double-peaked vs single-peaked, Figure 2), and the paper's stated preference for the steeper scaling is based on that difference. With the cut, the profiles become nearly identical, leading to the conclusion that emissivity conditions matter more. Since the cut is borrowed from an external photoionization study and is applied as a hard threshold with no sensitivity analysis, the central claim is not yet secure. The paper's own discussion acknowledges the need for a more sophisticated emissivity treatment, which further supports this concern. The concrete test I propose would directly probe whether the conclusion survives reasonable variations in the threshold or its smoothness. Because the paper is explicitly preliminary and the reader's CONDITIONAL verdict already reflects this uncertainty, my concern does not change the verdict; it reinforces the conditionality.","tokens_in":11822,"tokens_out":1687,"duration_ms":16835,"concrete_test":"Re-run the 2.5D FRADO line-profile calculation for s = -0.5 and s = -1.5 with log(phi) thresholds of 17, 18, and 19 cm^-2 s^-1, and also with a smooth emissivity weight w(phi) = phi / (phi + phi_0) for phi_0 spanning 10^17 to 10^19 cm^-2 s^-1. For each variant, compute the RMS difference and FWHM between the two radial scalings. If the profiles differ by more than about 10% in any reasonable threshold placement or smoothing, then the conclusion that emissivity treatment dominates over outflow radial scaling is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central conclusion - that the emissivity treatment dominates over the radial outflow scaling - rests on the step-function photon-flux threshold introduced in Section 4 and applied in Figure 4. Clouds contribute only if log(phi) ≳ 17-18 cm^-2 s^-1, a hard cut adopted from [37]. This binary threshold is what makes the s = -0.5 and s = -1.5 profiles nearly identical: it selects clouds from a narrow radial range, washing out the dynamical differences visible in Figure 2. However, no sensitivity analysis is provided. The threshold is stated as a range (17-18), not a fixed value, and the real emissivity of H-beta or MgII is a smooth function of ionization parameter, column density, and spectral shape, not a step in photon flux alone. If the effective threshold sits lower or higher, or if emissivity rises gradually with flux, the inner-cloud selection changes and the two outflow scalings may again produce distinguishable profiles. The paper itself concedes in Section 5 that 'a more sophisticated treatment of the cloud emissivity' is needed. Because the headline claim is precisely that cloud physical conditions outweigh the radial outflow behavior, the unvalidated binary cut is the most load-bearing assumption in the argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the 2.5D FRADO model to compute BLR cloud distributions and Doppler-broadened emission-line profiles for two radial scalings of the disk outflow rate, s = -0.5 and s = -1.5. With uniform emissivity the two scalings give different profiles (double-peaked for s = -0.5 versus single-peaked for s = -1.5); with a photon-flux threshold (log φ ≈ 17–18) the profiles become similar. The paper concludes that the emission-line shape depends more critically on cloud physical conditions and the adopted emissivity treatment than on the radial outflow scaling, while still tentatively favoring the steeper s = -1.5 scaling for the mean SDSS quasar benchmark.","tokens_in":12101,"tokens_out":4823,"duration_ms":39035,"significance":"If the central conclusion survives scrutiny, it is a useful caution against overinterpreting line-profile shape as a direct diagnostic of the radial mass-loss profile in the BLR, and it motivates a more detailed treatment of cloud emissivity. The paper's strengths include a self-consistent calculation of the launching radii (no arbitrary inner/outer radius inputs), the use of an external photoionization constraint for the photon-flux threshold, and an explicit comparison of two outflow scalings. However, the quantitative support for the secondary preference for s = -1.5 is not yet established, and the central emissivity conclusion rests on a single binary threshold with no sensitivity analysis. The paper therefore presents a plausible and interesting exploratory result, but the evidence is currently provisional.","major_comments":[{"comment":"The derivation in Eqs. (4)–(5) gives ˙M_z(r) ∝ r^{α-1}; with the quoted α ≈ -0.5, the predicted index is s = -1.5, not s = -0.5. The shallower case s = -0.5 is never derived from the stated optically thin or dust-opacity scalings (Eqs. 2, 3, and 5). Since the paper presents both s values as physically motivated outflow scenarios, the s = -0.5 branch needs either an explicit derivation or a clear label as an exploratory ad hoc case; otherwise the comparison in Figures 1, 2, and 4 is not between two model predictions.","section":"Section 3, Eq. (5)"},{"comment":"The preference for s = -1.5 is asserted from a visual comparison with the 'overall shape of the low-ionization broad emission profile' of the SDSS mean quasar, but the observed composite profile is not shown and no quantitative goodness-of-fit or profile metric is provided. The statement that s = -1.5 'better reproduces observed profiles' is therefore not supported by the evidence presented. An overlay with the Vanden Berk et al. composite, or a quantitative measure such as peak-to-wing ratio or FWHM comparison, is needed to justify this secondary claim.","section":"Section 4, Fig. 2; Section 5, Discussion"},{"comment":"The central conclusion that the emissivity treatment dominates over the outflow scaling rests on a binary photon-flux threshold: Section 4 quotes log φ ≳ 17–18 cm⁻² s⁻¹, but Figure 3 and Section 5 use φ ≳ 10^18 cm⁻² s⁻¹, so the exact threshold used in the profile calculation is not specified reproducibly. No sensitivity analysis is performed: varying the threshold within the stated range, smoothing the step function, or adopting a gradually rising emissivity could restore differences between the s = -0.5 and s = -1.5 profiles. Given that Section 5 itself concedes that 'a more sophisticated treatment of the cloud emissivity' is needed, the headline claim is currently provisional.","section":"Section 4, Fig. 4; Section 5, Discussion"}],"minor_comments":[{"comment":"Several grammatical issues should be corrected, e.g., 'how the mass ejection rate contribute' and 'Overally' in Section 2.","section":"Abstract and throughout"},{"comment":"'Thompson scattering' should be 'Thomson scattering'.","section":"Section 3"},{"comment":"The model is referred to inconsistently as '2.5 FRADO' and '2.5D FRADO'; use the latter throughout.","section":"Section 4"},{"comment":"Some reference entries contain malformed arXiv identifiers (e.g., 'arXiv:astro-ph/astro-ph/0306389') and one reference (ref. 19) is cited only as an arXiv e-print without journal details; these should be cleaned up.","section":"References"},{"comment":"The manuscript uses placeholder journal metadata ('Universe 2024, 1, 0' and a placeholder DOI); these need to be completed before formal submission.","section":"Front matter"}],"recommendation":"major_revision","confidential_remarks":"The paper is built largely on the author's own FRADO model series, which is not disqualifying, but the argument would be considerably stronger with an observed-profile overlay and a sensitivity test of the photon-flux threshold. The inconsistency between the derived s = -1.5 and the adopted s = -0.5 should be resolved, either by deriving both scalings or by clearly labeling the shallower case as exploratory. The manuscript also carries placeholder MDPI metadata that should be completed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a preliminary parameter study, not a breakthrough, but it does one thing that is genuinely new and useful. Naddaf takes his 2.5D FRADO model and compares two radial scalings for the disk outflow rate, s=-0.5 and s=-1.5, and shows that when you impose a photon-flux threshold on cloud emission (log phi ~ 17-18 cm^-2 s^-1), the line profile becomes nearly independent of which scaling you pick. That is an honest result: it tells BLR modelers that emissivity treatment can be a bigger lever than the radial outflow law. The paper is also transparent about its own limitations, explicitly flagging that a more sophisticated emissivity treatment is needed and that radiative transfer is ignored.\n\nThe soft spots are real but fixable. First, the s=-0.5 case is never derived; the paper's own Eq. 5 with alpha ~ -0.5 gives s=-1.5, so the shallower index appears as an ad hoc alternative. Second, the preference for s=-1.5 in the uniform-emissivity case rests on a qualitative comparison to the SDSS composite, but the observed profile is not plotted, so the reader cannot check the claim. Third, and most important: the central conclusion that emissivity dominates over outflow scaling is entirely dependent on the binary photon-flux cut. The threshold is quoted as a range (17-18), not a fixed value, and there is no sensitivity analysis. If the real emissivity rises smoothly with flux, or the effective threshold sits lower, the two scalings may again give distinguishable profiles. The paper's own Section 5 concession is effectively an admission of this.\n\nThe citation pattern is self-heavy, but that is reasonable since the model is the author's own; the external threshold is from Pandey et al. (2023) and is used appropriately. No data or code are shipped, but the model is described well enough that a competent person could reproduce the calculation.\n\nWho should read this: anyone building BLR models, especially those working with radiation-pressure or failed-wind pictures. It is a small step, but it identifies a real degeneracy.\n\nI would send it to referees, because the central question is clear and the weakness is specific and addressable. But I would tell the author to add a derivation or justification for s=-0.5, plot the observed comparison profile, and test the sensitivity to the threshold before I would trust the headline claim.","headline":"A modest but honest FRADO parameter study; the new flux-threshold result is real but hinges on a step-function cut that the paper itself admits is too simple.","tokens_in":12613,"tokens_out":2406,"would_cite":false,"duration_ms":21089,"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":"Low-ionization line shapes trace cloud illumination more than the disk-outflow radial law.","keywords":["broad line region","outflow rate","active galactic nuclei","radiation pressure","emission line profiles","FRADO model","dusty winds","photon flux threshold"],"falsifier":"Compute the profiles with a photoionization-based continuous emissivity, assigning each cloud a line emissivity that rises gradually with incident flux and depends on density, and compare the $s=-0.5$ and $s=-1.5$ cases; if the profiles diverge again, the paper's central claim fails. Observationally, a sample of quasars with independent kinematic measures of outflow steepness, such as blueshifted absorption features, that shows line profiles tracking outflow steepness would contradict the conclusion.","tokens_in":11623,"feed_emoji":"🔭","tokens_out":8568,"duration_ms":73677,"temperature":0.7,"pith_summary":"The paper asks whether the radial steepness of mass loss from an accretion disk controls the shapes of low-ionization broad emission lines in active galactic nuclei. Using the FRADO model, a failed radiatively accelerated dusty outflow that forms the broad line region (BLR), it compares two radial outflow laws, $\\dot{M}_z(r)\\propto r^{-0.5}$ and $\\propto r^{-1.5}$, for a benchmark quasar. When every cloud contributes equally to the line, the two laws give different profiles: double-peaked for the shallower scaling and single-peaked for the steeper one. But when only clouds with a photon flux above about $10^{18}\\ \\mathrm{cm^{-2}\\,s^{-1}}$ (the illumination needed to form H$\\beta$/MgII) are allowed to emit, both scalings produce the same single-peaked profile. The paper concludes that the observed line shape depends more critically on the clouds' physical conditions and the adopted emissivity than on the radial behavior of the disk outflow, while still favoring the $r^{-3/2}$ scaling as the better match to the mean quasar.","feed_headline":"Cloud illumination, not outflow law, sets quasar line shapes","feed_subtitle":"A photon-flux cutoff erases differences between two radial disk-outflow models, so cloud conditions dominate the profile.","key_machinery":"The carrying object is FRADO (Failed Radiatively Accelerated Dusty Outflow), a model upgraded to 2.5D in which dusty clumps are lifted from the surface of a Shakura-Sunyaev accretion disk by radiation pressure, lose their dust as they rise, and fall back ballistically; the model self-consistently outputs the BLR cloud distribution with no free launching parameters. The paper varies the radial outflow efficiency through the power-law index $s$ in $\\dot M_z(r)\\propto r^s$: $s=-0.5$ comes from an energy-wise approach with dust opacity scaling as $r^{-0.5}$, and $s=-1.5$ comes from single-scattering momentum balance. Emission-line profiles are then computed by Doppler-shifting cloud velocities at the mean viewing angle of $39^\\circ$ under two emissivity assumptions: uniform emission from all clouds, or emission only from clouds meeting the photon-flux threshold $\\log\\phi\\gtrsim17$-$18\\ \\mathrm{cm^{-2}\\,s^{-1}}$. The threshold is the mechanism that makes the two outflow scalings converge, and it is what carries the argument that cloud conditions, not the outflow law, dominate the line shape.","core_discovery":"The central claim is that the radial power-law index of the disk mass-loss rate is a secondary factor in shaping low-ionization broad emission lines; the primary factor is how cloud emissivity is treated. In the 2.5D FRADO setup, the shallower scaling $s=-0.5$ spreads clouds over a wider, slower BLR and, under uniform emissivity, yields a double-peaked line, while $s=-1.5$ concentrates clouds near the black hole and yields a single-peaked line. Imposing the physically motivated photon-flux condition $\\log\\phi\\gtrsim17$-$18\\ \\mathrm{cm^{-2}\\,s^{-1}}$ for H$\\beta$/MgII formation erases this distinction: only clouds launched from a narrow radial range near the black hole are bright enough to contribute, so both scalings give nearly identical single-peaked profiles. The paper therefore concludes that the shapes of low-ionization lines are dictated more by the illumination and dynamics of the inner clouds than by the specific radial form of the outflow, with the steeper $r^{-3/2}$ scaling still preferred for the mean quasar benchmark.","pith_inferences":["I read the flux-threshold result as implying that the observable low-ionization BLR is effectively a narrow, well-illuminated annulus; if true, reverberation-mapping size measurements trace the illumination-limited annulus rather than the full radial extent of the outflow.","A direct extension is to model H$\\beta$ and MgII together: because the two lines form under different threshold fluxes, the profile difference between them should be a diagnostic of the radial outflow law.","The coincidence that $r^{-3/2}$ matches the Keplerian orbital-frequency scaling hints at a feedback loop in which returning failed-wind clouds disturb the disk surface and trigger new cloud launches; a hydrodynamical simulation of cloud re-impact could test whether such feedback self-regulates the outflow rate."],"forward_implications":["Observed single-peaked low-ionization profiles cannot uniquely determine the radial outflow law, because the photon-flux threshold makes both $s=-0.5$ and $s=-1.5$ produce the same shape.","The steeper $r^{-3/2}$ scaling predicts a more compact, higher-velocity BLR, so BLR size and velocity dispersion measurements can still constrain the outflow efficiency even when profile shape cannot.","Across the transition to the dust-free inner region, the outflow launching efficiency should drop much more steeply, with power-law indices below $-2$ and down to roughly $-3.5$, concentrating high-ionization emission near the black hole.","For the mean SDSS quasar benchmark, the model yields broad lines with FWHM near 5000$-$6000 km/s and a slight blueshift, placing it in Population B2 (the moderate-accretion tile of the Eigenvector 1 classification) with negligible gravitational redshift."],"supporting_citations":[{"why":"Supplies the 2.5D FRADO cloud distributions, velocities, and self-consistent launching radii used to compute the line profiles.","marker":"[15]"},{"why":"Defines the analytical FRADO framework that the paper extends to compare outflow scalings.","marker":"[22]"},{"why":"Provides the photon-flux threshold for H$\\beta$/MgII formation that the second emissivity scenario adopts.","marker":"[37]"},{"why":"Gives the mean composite quasar spectrum that serves as the observational benchmark for the line shapes.","marker":"[34]"},{"why":"Establishes that the benchmark parameters represent typical SDSS quasars.","marker":"[35]"},{"why":"Supplies the Shakura-Sunyaev thin-disk flux $F(r)\\propto r^{-3}$ from which the outflow scalings are derived.","marker":"[26]"},{"why":"Provides the steeper radial outflow decline from disk-wind models that the paper compares against the FRADO scalings.","marker":"[20]"}],"fun_headline_variants":["Disk outflow law plays second fiddle to cloud illumination","Photon cutoff erases outflow differences in quasar lines","Cloud emissivity dominates over disk outflow profile in BLR","For low-ionization lines, cloud conditions beat outflow radial law","Radial outflow efficiency is not the key to quasar line profiles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The key load-bearing premise is the step-function emissivity threshold: a cloud either contributes to the line once its photon flux exceeds about $10^{18}\\ \\mathrm{cm^{-2}\\,s^{-1}}$ or not at all; if real emissivity varies smoothly with flux, or if the threshold sits at a different value, the radial outflow scaling would regain a visible role.","fun_headline_variants_meta":{"raw":{"variants":["Disk outflow law plays second fiddle to cloud illumination","Photon cutoff erases outflow differences in quasar lines","Cloud emissivity dominates over disk outflow profile in BLR","For low-ionization lines, cloud conditions beat outflow radial law","Radial outflow efficiency is not the key to quasar line profiles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000482,"raw_usage":{"total_tokens":2373,"prompt_tokens":926,"completion_tokens":1447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":1365}},"tokens_in":542,"tokens_out":1447,"duration_ms":10523,"temperature":1.0,"reasoning_tokens":1365,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:29:42.830171+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the profiles with a photoionization-based continuous emissivity, assigning each cloud a line emissivity that rises gradually with incident flux and depends on density, and compare the $s=-0.5$ and $s=-1.5$ cases; if the profiles diverge again, the paper's central claim fails. Observationally, a sample of quasars with independent kinematic measures of outflow steepness, such as blueshifted absorption features, that shows line profiles tracking outflow steepness would contradict the conclusion.","supporting_citations":[{"cited_title":"The picture of BLR in 2.5-D FRADO: Dynamics & Geometry","cited_arxiv_id":"2102.00336","evidence_quote":"Supplies the 2.5D FRADO cloud distributions, velocities, and self-consistent launching radii used to compute the line profiles."},{"cited_title":"Broad-line region in active galactic nuclei: Dusty or dustless?","cited_arxiv_id":"2310.05089","evidence_quote":"Provides the photon-flux threshold for H$\\beta$/MgII formation that the second emissivity scenario adopts."},{"cited_title":"Black holes in binary systems","cited_arxiv_id":null,"evidence_quote":"Supplies the Shakura-Sunyaev thin-disk flux $F(r)\\propto r^{-3}$ from which the outflow scalings are derived."}],"review_version":1}