{"id":"f64e0993-0223-4b9a-b997-e555ebbb39a2","arxiv_id":"2505.02227","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Gaia23bab is a prototypical EXor, with an accretion rate around 2e-7 solar masses per year and hydrogen emission consistent with temperatures of 5000-12500 K and densities of 1e8-1e11 cm^-3.","lead":"Astronomers caught a young star, Gaia23bab, during its 2023 outburst and measured how rapidly it is feeding from the disk around it. The measured feeding rate and the temperatures and densities in its hot gas make it a textbook example of an EXor, a young star that grows in repeated bursts.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The IRTF-epoch extinction is fitted to force tracer agreement, so the multi-tracer accretion rate and the subsequent excitation analysis are not fully independent; the concern is real but does not overturn the central estimate.","rationale":"The reader's conditional verdict matches the actual vulnerability. I read the central claim as the quantitative accretion rate plus the EXor classification. The quantitative claim is externally supported by the LBT value of 2.5 ± 0.6e-7, and the NOT epoch independently gives values bracketing 2e-7, but the IRTF multi-tracer agreement is achieved by optimizing A_V rather than by an independent extinction measurement. The systematic A_V uncertainty is therefore absent from the quoted 2.0 ± 0.5e-7, even though a plausible range of A_V near 2.7 to 4.1 changes the near-IR line luminosities by only about 10 to 30 percent. This makes the rounded central value robust, but the error bar and the description of independence are overstated. The same circularly chosen A_V is then used for the excitation analysis, which is an additional weakness; independent of that, the Balmer and Paschen decrements imply densities that differ by three orders of magnitude, so the excitation-based confirmation of the EXor classification is not as secure as the abstract implies. The classification itself still rests on the accretion rate, the recurrence of outbursts, and the spectral appearance, so the paper does not need rejection. The conditional verdict should stand, with a requirement to marginalize over A_V and to quote systematic errors in the accretion rate and excitation parameters. This matches the reader's weakest_assumption on the fitted extinction.","tokens_in":1940,"tokens_out":881,"duration_ms":145373,"concrete_test":"Recompute the IRTF-epoch accretion rate from the Table 2 line fluxes with A_V fixed to the independent quiescent value 3.2 ± 0.5 (and, as a sensitivity bound, 4.0), instead of the best-fit 3.6. If the weighted mean Macc remains within about 30% of 2.0e-7 and all tracers are consistent at the fixed A_V, then the fitted A_V is not load-bearing; if the mean shifts by more than about 30% or the tracer-to-tracer scatter exceeds 2σ at every fixed A_V, the claimed multi-tracer rate is partly an artifact of the A_V choice. The same fixed-A_V recomputation should be carried through to the Sect. 3.5 excitation fits to see whether the reported temperature and density ranges change.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weak point is the determination of A_V for the IRTF epoch in Sect. 3.4. The paper states that different A_V values were considered because the accretion rates derived from the different lines are supposed to be the same for the right A_V. Thus A_V = 3.6 is obtained by requiring the Alcalá et al. (2017) tracers to agree, and the same dereddened fluxes are then used in Sect. 3.5 to derive temperatures and densities. This makes the description of the accretion rate as measured independently with several line tracers too strong: at the IRTF epoch the inter-tracer agreement is imposed by the fitting procedure, and at the NOT epoch only Paβ and Brγ are available, giving rates of 1.9e-7 and 3.2e-7 that differ by a factor of 1.7. The central value is nevertheless not fragile: 3.6 lies within 1σ of the independently fitted quiescent A_V = 3.2 ± 0.5 (Sect. 3.3), the near-IR tracers have small A_λ/A_V, and Giannini et al. (2024) independently report 2.5 ± 0.6e-7. The real issue is that the quoted 0.5e-7 error omits the systematic A_V contribution and that the excitation analysis inherits a circularly chosen A_V.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a multi-wavelength photometric and spectroscopic study of the young eruptive star Gaia23bab during its 2023 outburst, triggered by a Gaia alert. The authors combine optical and near-infrared photometry, NIR spectra from NOT and IRTF, and archival data to characterize the outburst and previous activity. They derive revised stellar parameters (M1-type, Teff = 3630 K, L* = 0.72 L_sun, R* = 2.3 R_sun, M* = 0.40 M_sun), measure accretion rates at two epochs using several line tracers (Ca II, He I, Pa, and Br lines), and compare hydrogen line ratios to Case B theory and Kwan-Fischer models to infer excitation temperatures and densities. The main quantitative claim is that the accretion rate during the 2023 outburst was approximately 2 x 10^-7 M_sun/yr, consistent across epochs and with the independent LBT estimate of Giannini et al. (2024), and that the derived excitation conditions support the classification of Gaia23bab as a prototypical EXor.","tokens_in":24866,"tokens_out":2998,"duration_ms":37589,"significance":"If the central claim holds, the paper provides a well-observed example of an EXor with a stable accretion rate across three epochs of an outburst and a detailed comparison of multiple hydrogen-series diagnostics. The analysis combines original data with previously published LBT observations, and the revised stellar parameters (from a quiescent LBT spectrum) improve on earlier estimates. The authors are transparent about their assumptions, including the fitted extinction and the acknowledged limitations of Case B theory in T Tauri environments. The paper also documents a previous 2017 outburst and compares the two events, contributing to the growing sample of Gaia-alerted eruptive stars. However, the quantitative accretion rate relies on a fitted A_V at the IRTF epoch, and the quoted uncertainty omits several systematic contributions, so the precision of the central value is likely overstated.","major_comments":[{"comment":"The IRTF-epoch extinction A_V = 3.6 mag is determined by requiring that accretion rates from different lines agree, and the same dereddened fluxes are then used in Sect. 3.5 for the hydrogen-line excitation analysis. This makes the excitation analysis and the 'multiple independent tracers' description partially circular. The quoted error of 0.5 x 10^-7 M_sun/yr is only the standard deviation among tracers at a fixed A_V; it does not include the systematic contribution from the assumed A_V uncertainty, the distance uncertainty (900 ± 45 pc), or the uncertainties in R* and M*. I recommend propagating these and reporting a total systematic uncertainty, or at least explicitly stating that the quoted error is statistical only.","section":"Sect. 3.4"},{"comment":"The statement in the abstract and Sect. 4 that the accretion rate was 'measured independently with several line tracers' is stronger than the data support. At the NOT epoch, only Paβ and Brγ are available, and they give accretion rates of 1.9 x 10^-7 and 3.2 x 10^-7 M_sun/yr, which differ by a factor of about 1.7. Even at the IRTF epoch, the inter-tracer agreement is imposed by the fitted A_V. I suggest reporting the NOT-epoch rates separately and rephrasing the claim to 'consistent with about 2 x 10^-7 M_sun/yr across epochs' rather than implying fully independent multi-tracer agreement at each epoch.","section":"Sect. 3.4"},{"comment":"The Balmer decrement is best fitted with n_H = 10^8 cm^-3, while the Paschen decrement at the same epochs requires n_H = 10^11 cm^-3, a discrepancy of three orders of magnitude. The paper presents this without discussing whether it indicates two physically distinct emitting regions, model limitations, or an error in the extinction correction. Since this discrepancy bears on the claim that the derived temperatures and densities confirm the EXor classification, the authors should discuss its origin and its implications for the robustness of the excitation conditions.","section":"Sect. 3.5"}],"minor_comments":[{"comment":"The sentence 'A temperature of 3750 K is unlikely since it is below the stellar temperature (see Sect. 3.3)' is incorrect: the stellar Teff quoted in Sect. 3.3 is 3630 K, so 3750 K is above it. The intended argument (that 3750 K is close to or below the photospheric temperature) still holds, but the wording should be corrected.","section":"Sect. 3.5.2"},{"comment":"The notation '1.9 ∼ 10−7 M⊙ yr−1 and 3.2 ∼ 10−7 M⊙ yr−1' should use multiplication signs (e.g., 1.9 × 10^-7) for consistency with the rest of the paper.","section":"Sect. 3.4 and Table 2"},{"comment":"The right-hand panel of Figure 3 uses 'JD−2450000' as the x-axis label, while the left-hand panel uses 'MJD−57000'; please use a consistent time axis for both panels.","section":"Fig. 3"},{"comment":"In the text, '109 cm−3' and '1010 cm−3' should be typeset as '10^9 cm^-3' and '10^10 cm^-3' to avoid confusion.","section":"Sect. 3.5.1"},{"comment":"The service BHTOM is first mentioned in the text as 'BHTOM service' and in the facilities line as 'BHTOM.space'; please unify the spelling.","section":"Sect. 2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript overlaps with Giannini et al. (2024), which already confirmed the EXor classification and measured an accretion rate of (2.5 ± 0.6) x 10^-7 M_sun/yr at the LBT epoch. The present paper adds two more epochs and revised stellar parameters, which is a useful extension, but the novelty relative to the earlier work should be sharpened in the introduction. The central numerical claim is robust, but the omitted systematic error and the partially circular extinction determination require attention before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take on arXiv:2505.02227. This is a follow-up to Giannini et al. 2024 on the EXor Gaia23bab. What is actually new: revised stellar parameters (M1, L=0.72 Lsun, M=0.4 Msun), accretion rates at two additional epochs during the 2023 outburst, and the first excitation temperature/density estimates from HI line ratios. The accretion rate at the NOT and IRTF epochs is consistent with ~2e-7 Msun/yr, matching the LBT value from Giannini et al. That consistency across three epochs and several tracers is the paper's real strength. The light-curve comparison of the 2017 and 2023 outbursts is a nice bonus.\n\nThe soft spots are where you'd expect in a single-object study like this. The IRTF-epoch A_V = 3.6 is not measured independently; it is chosen so the accretion rates from different lines agree. That same dereddened flux set is then fed into the excitation analysis. So the multi-tracer agreement at that epoch is partly imposed by construction, and the quoted ±0.5e-7 error does not include the A_V systematic. The excitation results are fragile in other ways too: the authors themselves note Case B has been questioned for T Tauri stars, and the Balmer decrement prefers nH=1e8 while Paschen prefers 1e11—three orders of magnitude apart. I'd treat the temperatures and densities as indicative, not measured.\n\nThat said, the central claim holds up. A_V=3.6 is within 1σ of the independently fitted quiescent A_V=3.2±0.5, near-IR tracers have low extinction sensitivity, and the independent LBT rate agrees. The circularity concern is real but does not overturn the accretion rate. The paper is honest about its assumptions, and the data tables are thorough.\n\nWho gets value: anyone working on EXor statistics, accretion diagnostics, or Gaia-alert follow-up. It is exactly the kind of well-documented single-object study that the episodic accretion community needs to build samples. It does not introduce new physics or methodology, and it does not resolve a long-open question, but it is a solid benchmark.\n\nRecommendation: send it to peer review. A referee should ask for error bars that include A_V systematics and a discussion of why the Balmer and Paschen densities disagree, but neither issue warrants rejection.","headline":"Careful follow-up of a single EXor gives a solid ~2e-7 Msun/yr accretion rate, but the excitation temperatures and densities are shakier than the abstract suggests.","tokens_in":25753,"tokens_out":2318,"would_cite":true,"duration_ms":27138,"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":"This paper shows that during the 2023 outburst of the young star Gaia23bab, the mass accretion rate was roughly $2\\times10^{-7}$ solar masses per year at two independent epochs, and that hydrogen line excitation conditions mark the system…","keywords":["EXors","episodic accretion","young stellar objects","eruptive variables","accretion rate","hydrogen recombination lines","Gaia photometric alerts","pre-main-sequence stars"],"falsifier":"A direct, method-independent measurement of the extinction toward Gaia23bab—for instance, fitting the quiescent optical and near-infrared spectrum with a stellar atmosphere plus a dust screen, or comparing an X-ray hydrogen column with the optical reddening—would settle whether $A_V = 3.6$ mag is right; a value outside $3.2 \\pm 0.5$ mag would rescale all extinction-corrected line luminosities and could move the inferred accretion rate and excitation temperatures outside the EXor ranges claimed here.","tokens_in":24403,"feed_emoji":"⭐","tokens_out":9903,"duration_ms":115007,"temperature":0.7,"pith_summary":"This paper studies a young star that the Gaia satellite flagged as brightening by two magnitudes in March 2023, and asks what the outburst reveals about how such stars build their mass. By measuring the strengths of several hydrogen and metallic emission lines in spectra taken near maximum light and during the fade, it derives the rate at which disk material falls onto the star: about $2\\times10^{-7}$ solar masses per year, a value that stays the same at two epochs and matches a third measurement from earlier work. The paper then compares the hydrogen line ratios to models of ionized gas and accretion flows, finding temperatures of roughly 5000--12500 K and densities of $10^8$--$10^{11}$ particles per cubic centimeter, values typical of EX Lupi-type eruptive stars. In short, the paper makes the case that Gaia23bab is a prototypical EXor and that its 2023 outburst is a well-characterized instance of episodic accretion.","feed_headline":"Star pulled a steady 0.0000002 solar masses per year in 2023 outburst","feed_subtitle":"Hydrogen-line ratios across three epochs mark Gaia23bab as a prototypical EXor young star.","key_machinery":"Multi-epoch near-infrared spectroscopy is the workhorse: Paschen and Brackett hydrogen lines, the Ca II triplet, He I, and O I are detected in spectra taken near maximum and during the fading phase, and the ratios of their fluxes are the observables that carry the argument. The accretion rate follows from converting extinction-corrected line luminosities ($L_{\\rm line}=4\\pi d^2 f_{\\rm line}$) to accretion luminosities using the empirical relations of Alcalá et al. (2017) and then applying the magnetospheric accretion formula $\\dot{M}_{\\rm acc}=1.25\\,L_{\\rm acc}R_\\star/(GM_\\star)$ with an assumed inner-disk radius of $5R_\\star$. The excitation analysis compares observed Brackett decrements to Case B recombination predictions (Hummer & Storey 1987), in which the plasma is opaque to Lyman-$\\alpha$ and optically thin in higher lines so that line ratios depend mainly on temperature and density, and compares Balmer and Paschen decrements to the T Tauri wind and accretion models of Kwan & Fischer (2011) and Edwards et al. (2013), with $\\chi^2$ minimization selecting best-fit temperature and density values.","core_discovery":"The paper establishes that during the 2023 outburst of Gaia23bab the mass accretion rate was $\\sim 2.0\\times10^{-7}\\,M_\\odot\\,\\mathrm{yr}^{-1}$, measured at two independent epochs with several tracers, and consistent with a third epoch from earlier work; the accretion rate did not change significantly while the star faded. The hydrogen Brackett, Paschen, and Balmer line ratios, compared with Case B recombination and accretion-flow models, yield excitation temperatures of roughly 5000--12500 K and hydrogen densities of $10^8$--$10^{11}\\,\\mathrm{cm}^{-3}$, with the Ca II line ratios pointing to a denser, cooler component ($n_H\\sim10^{12}\\,\\mathrm{cm}^{-3}$, $T\\lesssim7500$ K) likely associated with the disk boundary layer. On this basis the authors conclude that Gaia23bab is a prototypical EXor, matching EXor behavior in both its accretion rate and its emitting-gas conditions.","pith_inferences":["A testable extension follows from the extinction assumption: if future independent measurements fix $A_V$, the same multi-line method can be turned into a calibration of EXor accretion rates that does not depend on a self-consistency choice.","The Paγ/He I ratio above 1, unusual for classical T Tauri stars, may serve as a quick spectroscopic tag for EXor-type eruptive stars; checking it on a larger sample would show whether it separates burst-phase accretors from quiescent ones.","Because the 2017 and 2023 bursts look similar in amplitude and color, Gaia23bab is a candidate for predicting recurrence from a single burst; monitoring through the next decline would test whether the color loop and the disappearance of CO emission repeat.","If the accretion rate stayed constant while the continuum faded, then the optical decline is driven more by changing extinction or disk emission than by a drop in accretion power; this can be checked by correlating the color evolution with the accretion-tracer fluxes."],"forward_implications":["The accretion rate was essentially unchanged across three epochs (near maximum, mid-fade, and later fade), so the 2023 outburst maintained a roughly constant mass-feeding rate rather than spiking and decaying.","Hydrogen line excitation conditions overlap those measured in other EXors, so Gaia23bab joins the class not only by accretion rate but by the physical state of its emitting gas.","The two documented outbursts (2017 and 2023) show similar amplitudes and color behavior, implying the underlying disk or accretion instability repeats on similar timescales.","The Ca II line ratios imply a higher-density, cooler component near the disk boundary layer, suggesting the accretion flow in Gaia23bab is structured rather than a single uniform medium.","The non-detection of the CO bandhead and the weakening of Na I at the later epoch show that the disk's near-infrared emission faded faster than the accretion tracers, revealing which spectral features track the declining phase of an EXor outburst."],"supporting_citations":[{"why":"Gives the 900 ± 45 pc distance and cluster membership used to convert fluxes to luminosities, and first noted the EXor candidacy.","marker":"(Kuhn et al. 2023)"},{"why":"Provided the earlier LBT spectrum, the third-epoch accretion rate (2.5 ± 0.6) × 10−7 M⊙ yr−1, and the initial confirmation of the EXor nature.","marker":"(Giannini et al. 2024)"},{"why":"Supplies the empirical relations that turn extinction-corrected line luminosities into accretion luminosities.","marker":"(Alcalá et al. 2017)"},{"why":"Gives the magnetospheric accretion formula and the assumed 5R★ inner-disk radius used to derive the accretion rate.","marker":"(Hartmann et al. 1998)"},{"why":"Provides the Case B recombination line-ratio grids used for the Brackett series excitation fits.","marker":"(Hummer & Storey 1987)"},{"why":"Provides Balmer and Paschen emission models as well as Ca II, He I, and O I ratio grids used for the temperature and density analysis.","marker":"(Kwan & Fischer 2011)"},{"why":"Gives the Paschen decrement model predictions used in the excitation fits.","marker":"(Edwards et al. 2013)"},{"why":"Provides the evolutionary tracks from which the stellar mass 0.40 M⊙ is read off.","marker":"(Siess et al. 2000)"},{"why":"Supplies the bolometric correction and effective temperature for the M1 spectral type used to set the stellar parameters.","marker":"(Pecaut & Mamajek 2013)"}],"fun_headline_variants":["Steady accretion marks Gaia23bab as textbook EXor","Gaia23bab's 2023 outburst fits EXor mold","Hydrogen lines peg Gaia23bab as classic EXor","Steady 2e-7 solar masses/yr: Gaia23bab's EXor signature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the star is reddened by an extinction of 3.6 magnitudes, a value chosen because it makes the accretion rates from different lines agree, and the same correction is then used to deredden the hydrogen line fluxes that set the temperatures and densities.","fun_headline_variants_meta":{"raw":{"variants":["Steady accretion marks Gaia23bab as textbook EXor","Gaia23bab's 2023 outburst fits EXor mold","Hydrogen lines peg Gaia23bab as classic EXor","Steady 2e-7 solar masses/yr: Gaia23bab's EXor signature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000591,"raw_usage":{"total_tokens":2833,"prompt_tokens":1070,"completion_tokens":1763,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":686,"completion_tokens_details":{"reasoning_tokens":1687}},"tokens_in":686,"tokens_out":1763,"duration_ms":16339,"temperature":1.0,"reasoning_tokens":1687,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:57:41.076815+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct, method-independent measurement of the extinction toward Gaia23bab—for instance, fitting the quiescent optical and near-infrared spectrum with a stellar atmosphere plus a dust screen, or comparing an X-ray hydrogen column with the optical reddening—would settle whether $A_V = 3.6$ mag is right; a value outside $3.2 \\pm 0.5$ mag would rescale all extinction-corrected line luminosities and could move the inferred accretion rate and excitation temperatures outside the EXor ranges claimed here.","supporting_citations":[],"review_version":1}