{"id":"3d065a7d-98bd-4c64-a47f-6401dd960868","arxiv_id":"2506.04486","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Adding viscous shear heating to Ae disk models makes Hα emission much easier to produce for spectral types A2 and cooler, and this dependence is proposed as a new way to measure the disk viscosity parameter α, though the current data give inconsistent values.","lead":"This paper adds frictional heating from the disk's own rotation to models of gas disks around A-type stars, and shows that for stars cooler than spectral type A2 this heating strongly boosts the predicted hydrogen emission. It then tries to read off the disk's viscosity parameter α from how the number of detected Ae stars falls with temperature, but the two tests give contradictory answers.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Survey incompleteness and the mapping of noiseless model Hα detectability to LAMOST S/N≈90 spectra are unmodeled, so the count-based α≤0.1 constraint in Eq. (10) is not established.","rationale":"The paper has two independent pieces: the Bedisk/Beray modeling with shear heating, and the comparison with the LAMOST CAe sample. The modeling result—that shear heating raises disk temperatures and increases the Hα emission fraction for A2 and later spectral types—is a legitimate, parameter-swept prediction; I do not see an internal inconsistency in Eqs. (3)-(6) or in the radiative-equilibrium treatment that would undermine this. The load-bearing weakness is in the second piece, the calibration to observed CAe numbers and EW distributions. The reader's weakest assumption identifies exactly this: the equivalence between noiseless model classifications and LAMOST catalog membership, plus the treatment of raw LAMOST A-star counts as unbiased population fractions. If weak A3/A4 emission is missed at S/N≈90, the observed N_e_A34/N_e_A01 = 15/109 is too low, so the derived E_A34 ≈ 12% and the resulting α ≤ 0.1 are biased. This is a one-way bias that would also weaken the EW-CDF method for the same spectral types. The paper honestly notes the internal inconsistency between the count-based and EW-based α values and attributes it to small-number statistics; the selection-function issue is a more concrete and testable explanation. My proposed test—injecting noise into the model profiles and recomputing the completeness—would settle whether the count-based constraint survives. Because the modeling contribution stands on its own and the observational calibration is clearly flagged by the authors as preliminary, the appropriate verdict remains CONDITIONAL, unchanged from the reader's assessment.","tokens_in":20384,"tokens_out":7350,"duration_ms":71884,"concrete_test":"Add Gaussian noise to the computed noiseless Hα profiles for A2-A4 models at R=2000 to match the LAMOST S/N≈90 and continuum normalization of Section 4; rerun the Section 2.3 peak-finding classification and measure the completeness of weak-emission models (e.g., EW < 1 Å) as a function of spectral type and inclination. If a substantial fraction of models classified as emitting in the noiseless case are lost at S/N=90, then N_e_A34 is systematically incomplete and the count-based α ≤ 0.1 constraint from Eq. (10) is invalid; one should then either add a detection-efficiency term to Eq. (9) or restrict the comparison to EWs above a reliable detection limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The proposed α diagnostic is only as strong as the mapping between model Hα classifications and the LAMOST CAe catalog. Section 2.3 defines 'detectable Hα emission' by peak-finding on noiseless, R=2000 profiles, while Section 4.1 uses raw LAMOST A-star counts (f_A34/f_A01 = 0.86) and assumes D_c2/D_c1 = 1 in Eqs. (8)-(10). LAMOST is a targeted survey with magnitude/color selection toward the anticenter, and the CAe sample of Anusha et al. (2021) has S/N ≈ 90 at Hα; no completeness or selection model is included. Because A3/A4 CAe stars have systematically weaker, narrower Hα emission (Figure 5), a noiseless peak-finder will flag many weak profiles that real LAMOST spectra at S/N≈90 would not. This biases N_e_A34 downward; the inferred E_A34 ≈ 12% and the count-based α ≤ 0.1 then reflect survey incompleteness rather than disk viscosity. The same selection issue affects the EW CDF fits in Section 4.2, where the A3/A4 samples are the very objects most likely to be incomplete. Until this selection function is modeled, neither the α ≤ 0.1 nor the α ≈ 1.0 conclusion is secure, so the central claim that this dependence can be used to constrain α is not yet demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper models the thermal structure of gaseous circumstellar disks around main-sequence A-type stars (classical Ae stars) using the Bedisk/Beray code suite, adding a volumetric viscous shear-heating term d(R,Z) = (3/2) α P Ω (Eq. 6) for a Keplerian disk with α-viscosity. On a grid of 174,240 star+disk models spanning spectral types A0–A5, disk density parameters, inclination, and α = 0.01–1.0, the authors compute Hα profiles and classify them as emission or absorption. They find that for A2 and later spectral types, shear heating significantly raises disk temperatures and increases the fraction of models with detectable Hα emission (Figures 3 and 4). The paper then attempts to constrain α using the LAMOST DR5 CAe sample of Anusha et al. (2021): a count-based comparison (Section 4.1) favors α ≤ 0.1, while fitting Hα equivalent-width CDFs (Section 4.2) favors α ≈ 1.0; the paper concludes that these results are inconsistent and attributes this to small-number statistics. The central proposal is that the spectral-type dependence of Hα emission can be used as a radiative diagnostic of disk viscosity.","tokens_in":20685,"tokens_out":4358,"duration_ms":40857,"significance":"If the shear-heating mechanism is correct, the paper offers a novel, observationally accessible probe of the viscosity parameter α in cool Ae disks, complementing traditional viscous-timescale estimates used for classical Be stars. The theoretical modeling is thorough: the grid is large, the radiative-equilibrium calculations are described in detail, and the paper is commendably explicit about the internal inconsistency between the two proposed constraints. The comparison to observations, however, rests on assumptions about survey completeness and statistical significance that are not yet demonstrated; as it stands, the paper establishes the plausibility of the diagnostic rather than a secure measurement of α. The honest reporting of the disagreement between the count-based and EW-based constraints is a strength, but it also means the central claim is not yet fully supported.","major_comments":[{"comment":"The count-based constraint α ≤ 0.1 equates the noiseless, R = 2000 model Hα emission classification of Section 2.3 with actual detection in the LAMOST DR5 sample of Anusha et al. (2021), which has S/N ≈ 90 at Hα, without modeling LAMOST's target selection or the detection threshold for weak, narrow emission lines. Because A3/A4 CAe stars have systematically weaker Hα emission (Figure 5), a noiseless peak-finder will classify many model profiles as 'emitting' that would not be detected in real LAMOST spectra, biasing N_e_A34/B01 downward. The inferred E_A34 ≈ 12% and the resulting α ≤ 0.1 may therefore reflect survey incompleteness rather than disk viscosity. The authors should either model the LAMOST selection function and S/N-dependent detection probability or clearly state the range of completeness corrections that would change the inferred α.","section":"Section 4.1, Eqs. (8)–(10)"},{"comment":"The claim that only α = 1.0 yields a consistent set of fits across spectral types A0–A4 is based on a grid search over 1512 (μ_ρ, σ_ρ, μ_n, σ_n) combinations per spectral type, with the 'common region' in Figure 8 identified visually rather than by a formal statistical test. No multiple-testing correction is applied, so the appearance of an overlap region at α = 1.0 could be a chance coincidence given the large number of trials. The authors should quantify the significance of the overlap, for example by performing a permutation test that scrambles spectral-type labels or by computing the expected number of chance overlaps under a null model.","section":"Section 4.2, Figure 8"},{"comment":"The two constraints derived in Sections 4.1 and 4.2 are mutually inconsistent (α ≤ 0.1 from counts versus α ≈ 1.0 from EW CDFs), and the α = 1.0 model predicts A5 CAe stars that are not observed. The paper attributes this to small-number statistics, but the inconsistency directly undermines the stated conclusion that the dependence of Hα on shear heating 'can be used to constrain' α. The authors should provide a quantitative joint assessment, such as a likelihood-based combination of both constraints, or explicitly conclude that current data do not yet provide a reliable α measurement and specify what sample sizes or observations would resolve the tension.","section":"Section 5, Conclusions"},{"comment":"The assumption D_c2/D_c1 ≈ 1 (that the fraction of A-type stars hosting disks is the same for A0–A1 and A3–A4) is not justified beyond a plausibility argument. If the disk-formation efficiency varies with spectral type—for example, because stellar rotation rates or wind strengths change across this range—the inferred E_A34 and hence the α constraint would shift. The authors should discuss the sensitivity of Eq. (10) to plausible ranges of D_c2/D_c1.","section":"Section 4.1, Eq. (9)"}],"minor_comments":[{"comment":"There is a typo in 'its' versus 'it's' in the paragraph discussing spectral-type re-estimation: 'to estimate it's equivalent width' should be 'its equivalent width'.","section":"Section 4"},{"comment":"The legend entry 'α=0.0' appears inconsistent with the text, which lists α = 0.01 as the lowest value; the figure caption or the text should be corrected to match.","section":"Figure 3"},{"comment":"The sentence 'this comparison is be repeated over the 1512 combinations' contains a grammatical error ('is be') and should read 'this comparison is repeated'.","section":"Section 4.2"},{"comment":"The table entries are log10 of KS probabilities; the caption should state this explicitly and explain why negative values are used, as the current description is terse.","section":"Table 3"},{"comment":"The exclusion of parameter combinations that cannot produce samples of 20 or more emission-line stars should be discussed in the main text, because it could bias the reported numbers of fitting models, particularly for late spectral types with small α.","section":"Section 4.2, footnote 12"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and addresses an interesting question. The main concern is that both observational constraints are presented without a rigorous treatment of survey completeness and multiple testing, and the two constraints disagree. The authors are transparent about this disagreement, which speaks well of their scientific honesty, but the central claim is not yet established. I would encourage a revision that adds a selection-function model, formal significance estimates for the EW-CDF overlap, and a joint assessment of the two constraints. No issues with the theoretical derivation or the computational grid were found."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — read this one for the modeling, not for the α measurement. The genuinely new thing is that they add the volumetric shear heating term d = (3/2)αPΩ to Bedisk's radiative equilibrium and show it changes predicted disk temperatures and Hα detectability for A2 and later. That is a real, clearly explained advance, and the grid is enormous and fully specified (174,240 profiles). The shear heating derivation is standard and consistent, and the paper is honest about the limitations.\n\nThe soft spots are all in the comparison to LAMOST. First, 'detectable Hα emission' is defined by peak-finding on noiseless, R=2000 model profiles, but the LAMOST CAe sample has S/N≈90 and unknown target-selection completeness. Because A3/A4 emission is systematically weaker and narrower, the count-based α≤0.1 (Eq. 10) probably reflects survey incompleteness as much as disk viscosity. That is a load-bearing problem for the count argument, and the stress-test note is right about it. Second, Section 4.2's statistics are sloppy: logP ≥ -0.60 is called the 5% level, but that is P≈25%, and they search 1512 parameter combinations per panel without any multiple-comparison correction. That inflates the number of 'fits' and weakens the α=1.0 preference. Third, the two methods yield inconsistent α values (≤0.1 from counts, ≈1.0 from EW CDFs), and α=1.0 models predict A5 emitters that are not seen. The paper acknowledges this, but it means the proposed α diagnostic is plausible, not demonstrated.\n\nWhat holds up: the shear-heating modeling is a solid contribution that should stand regardless of the calibration outcome. The paper is also refreshingly direct about its own inconsistencies. With a corrected threshold, a proper multiple-testing treatment, and a model of LAMOST's selection function, the diagnostic could become useful. As is, it deserves a serious referee but needs revision.","headline":"Solid shear-heating models for cool Ae disks; the α constraint is plausible but not yet demonstrated, so referee it but expect revision.","tokens_in":21285,"tokens_out":2268,"would_cite":true,"duration_ms":23306,"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":"Viscous shear heating—not just starlight—can power Hα emission in the disks of the coolest Classical Ae stars, and the strength of that emission can be used to estimate the disk viscosity parameter α.","keywords":["Classical Ae stars","circumstellar disks","Hα emission","shear heating","disk viscosity","Shakura-Sunyaev alpha parameter","decretion disks","LAMOST survey"],"falsifier":"A high signal-to-noise Hα survey of a large, spectroscopically complete sample of A3–A5 stars in the same survey fields, counting emission-line stars down to equivalent widths well below the current detection limit, would decide between the count-based α≤0.1 conclusion and the EW-based α≈1 conclusion; in particular, finding numerous A3/A4 CAe stars or any A5 CAe star would falsify the low-α count constraint, while finding none at high signal-to-noise would challenge the high-α EW fit.","tokens_in":20104,"feed_emoji":"⭐","tokens_out":7731,"duration_ms":65709,"temperature":0.7,"pith_summary":"Classical Ae (CAe) stars are A-type main-sequence stars whose Hα emission betrays a dust-free circumstellar disk, but for the coolest members the central star's ultraviolet radiation is too weak to explain the observed emission. This paper asks whether viscous shear heating—the energy dissipated as the disk's Keplerian rotation shreds itself through viscosity—supplies the missing warmth. Using a radiative-equilibrium disk code with a volumetric shear heating term $d(R,Z)=\\tfrac{3}{2}\\alpha P\\Omega$, the authors find that shear heating becomes important at spectral type A2 and later, and that the presence and strength of Hα emission grow with the viscosity parameter α. They propose that the steep decline in CAe numbers toward A3/A4, together with the distribution of Hα equivalent widths, can be used as an ensemble constraint on α, and they attempt such constraints with a 159-star survey sample. A sympathetic reader would care because this is a new, purely radiative route to a quantity—disk viscosity—that previously could only be estimated from the timing of disk variability.","feed_headline":"Cool Ae star disks may need viscous shear to shine in Hα","feed_subtitle":"Hα emission after spectral type A2 could put a number on disk viscosity.","key_machinery":"The load-bearing object is the volumetric shear-heating rate $d(R,Z)=\\frac{3}{2}\\alpha P\\Omega$, inserted into the radiative-equilibrium solver Bedisk (Sigut & Jones 2007) as an extra source in the energy balance. It derives from the standard viscous dissipation rate $D(R)=\\frac{1}{2}\\nu\\Sigma(R\\,d\\Omega/dR)^2$ with $\\nu=\\alpha c_s H$ and Keplerian rotation, and its magnitude is set entirely by the α parameter; because shear heating and stellar photoionization scale differently with gas density, the ratio of heating to cooling changes with spectral type, which is why the effect switches on near A2. Around this rate the paper builds a grid of 174,240 Hα profiles over spectral type, disk density parameters, α, and inclination, automatically classifying each profile as emission or absorption at R=2000 to define what 'would be observed as a CAe star'.","core_discovery":"On the paper's own terms, the central claim is that viscous shear heating is a necessary thermal ingredient in the coolest Classical Ae star disks, and that its strength is readable from Hα emission. The authors compute two-dimensional temperature structures of gaseous, axisymmetric disks around main-sequence A0–A5 stars with the Bedisk code, adding to the usual radiative heating by the stellar photoionizing field a viscous dissipation rate $d(R,Z)=(3/2)\\alpha P\\Omega$ (with $\\alpha$ the Shakura–Sunyaev viscosity parameter, $P$ the gas pressure, and $\\Omega$ the Keplerian angular velocity). For A0 and A1 disks, even $\\alpha=1$ changes little; for A2 and later, shear heating substantially raises disk temperatures and increases both the fraction of model disks that produce detectable Hα emission and the strength of that emission. The paper therefore proposes that the observed decline in CAe incidence from A0 to A4 and the Hα equivalent-width distributions can jointly constrain α. Applying the proposal to the survey sample of 159 CAe stars, they find a tension: the count-based decline favors $\\alpha \\lesssim 0.1$, while matching the equivalent-width distributions across spectral types favors $\\alpha\\approx 1.0$ (with no consistent solution for $\\alpha\\le 0.3$). The paper reads this as a promising but not yet settled diagnostic, limited by small numbers of A3/A4 stars and by the absence of A5 members.","pith_inferences":["If the EW-based preference for α≈1 survives larger samples, it would imply that CAe disks are more viscous at late spectral types than the α≈0.1–0.3 values typically inferred for classical Be disks, suggesting a spectral-type trend in the viscosity mechanism.","The same radiative method could be extended to shell stars, whose Hα central absorption is also sensitive to disk temperature and which are observed to linger to A5–A7; matching shell-star statistics would give an independent cross-check.","The count-versus-EW discrepancy might dissolve if disk incidence is not equal across spectral types, e.g., if slower rotation among later A stars makes disks rarer; measuring CAe disk incidence independently would separate this from the viscosity constraint.","With larger samples from later survey data releases, the α≈1 fit's 5%-level KS overlap over many trial distributions could be re-tested with proper multiple-testing control, since 1512 parameter combinations were searched."],"forward_implications":["If shear heating is required, purely radiatively heated disk models systematically underpredict Hα emission for A2–A4 stars, so any thermal model of cool Ae disks must include a viscous heating term.","The observed CAe fraction versus spectral type becomes a direct readout of α: larger α keeps more late-A disks hot enough to emit Hα.","For A0/A1 stars Hα is insensitive to α, so those stars cannot constrain viscosity; constraints come almost entirely from A2 and later.","An α near 1 predicts detectable Hα emission from A5 disks, which the current survey sample does not contain; the absence of A5 CAe stars is therefore a test of the high-α end.","Because the count-based and EW-based constraints currently disagree (α≤0.1 versus α≈1), the paper's own conclusion is that more cool CAe stars, especially at A3 and later, are needed to resolve the ensemble α."],"supporting_citations":[{"why":"Supplies the 159-star CAe sample whose spectral-type counts and Hα equivalent widths are the observational constraints.","marker":"Anusha et al. (2021)"},{"why":"Defines the α-prescription $\\nu=\\alpha c_s H$ that the paper adopts for the viscosity and whose parameter α the paper aims to constrain.","marker":"Shakura & Sunyaev (1973)"},{"why":"Provides the viscous decretion disk model and the Keplerian rotation assumption, and sets the numerical factor in the shear heating rate.","marker":"Lee et al. (1991)"},{"why":"Gives the standard viscous dissipation expression $D(R)=\\frac{1}{2}\\nu\\Sigma(R\\,d\\Omega/dR)^2$ from which the volumetric heating rate is derived.","marker":"Pringle (1981)"},{"why":"Describes the Bedisk radiative-equilibrium code into which the shear heating term is inserted.","marker":"Sigut & Jones (2007)"},{"why":"Documents the Beray code used to compute the 174,240 Hα profiles.","marker":"Sigut (2018)"},{"why":"Supplies the disk density parameter ranges and the Gaussian distribution scheme used to generate synthetic Hα samples.","marker":"Sigut & Ghafourian (2023)"},{"why":"Provides the main-sequence mass and effective temperature calibration and spectral-type assignments used for the central star models.","marker":"Gray (2022)"}],"fun_headline_variants":["Hα emission in cool Ae stars may pin down disk viscosity","Viscous shear heats cool Ae disks, Hα reveals its strength","Cool Ae star Hα output could measure disk viscosity","Ae disks shine in Hα once shear heating kicks in","Hα from cool Ae disks offers a viscosity diagnostic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The premise that would sink the central claim is the equivalence between 'detectable Hα emission' as defined by peak-finding on noise-free model profiles at R=2000 and actual membership in the CAe catalog, with raw survey A-star counts and equal disk incidence across spectral types treated as unbiased.","fun_headline_variants_meta":{"raw":{"variants":["Hα emission in cool Ae stars may pin down disk viscosity","Viscous shear heats cool Ae disks, Hα reveals its strength","Cool Ae star Hα output could measure disk viscosity","Ae disks shine in Hα once shear heating kicks in","Hα from cool Ae disks offers a viscosity diagnostic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000266,"raw_usage":{"total_tokens":1651,"prompt_tokens":1027,"completion_tokens":624,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":643,"completion_tokens_details":{"reasoning_tokens":540}},"tokens_in":643,"tokens_out":624,"duration_ms":6034,"temperature":1.0,"reasoning_tokens":540,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:44:08.391667+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high signal-to-noise Hα survey of a large, spectroscopically complete sample of A3–A5 stars in the same survey fields, counting emission-line stars down to equivalent widths well below the current detection limit, would decide between the count-based α≤0.1 conclusion and the EW-based α≈1 conclusion; in particular, finding numerous A3/A4 CAe stars or any A5 CAe star would falsify the low-α count constraint, while finding none at high signal-to-noise would challenge the high-α EW fit.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the Bedisk radiative-equilibrium code into which the shear heating term is inserted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the Beray code used to compute the 174,240 Hα profiles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the main-sequence mass and effective temperature calibration and spectral-type assignments used for the central star models."}],"review_version":1}