{"id":"660ff11e-9c0b-4c51-97dd-c1518f0344c9","arxiv_id":"1908.06135","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The four components of the young binaries S CrA and VV CrA are K7 to M1 stars with variable infrared veiling and high accretion luminosities, placing the systems between the Class I and Class II evolutionary stages.","lead":"This paper used five years of high-resolution near-infrared Keck spectra to measure the spectral types, disk veiling, and hydrogen emission of the four stars in two young binary systems in Corona Australis. A generalist might read it because the results make these heavily obscured stars clear examples of actively accreting, transitional young stars, which matters for how binaries influence planet formation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Epoch-mixed SED photometry leaves the Class I/II classification underdetermined; the only Class I component sits right at the boundary.","rationale":"The paper has strong, partly independent evidence for the qualitative central claim: multi-epoch H-band spectra show variable veiling and Br16 emission in all four components, new spectral types are K7-M1, and the derived accretion luminosities are high. Those results support the broad conclusion that S CrA and VV CrA are young, actively accreting binaries. The more specific 'bridging Class I and Class II' statement, however, is pinned to SED spectral indices assembled from photometry spanning roughly two decades, combined with flux ratios from a single 2015 epoch. Because the paper itself documents large flux-ratio variability, this is the least secure link in the argument. I agree with the reader's weakest_assumption and do not see a reason to change the conditional verdict: the classification should be checked against contemporaneous photometry and given realistic uncertainties. The internal 'Class II' wording for S CrA A in Section 3.7 reinforces that the classification step needs tightening, but it does not overturn the core youth/accretion result.","tokens_in":19255,"tokens_out":10376,"duration_ms":100506,"concrete_test":"Recompute the four spectral indices in Table 4 using, for each star, K- and N-band photometry from a single contemporaneous source: for VV CrA use only the Scicluna et al. (2016) K and 10.35 micron fluxes; for S CrA use only the McCabe et al. (2006) K and 10.8 micron fluxes. Also repeat the calculation using the Prato et al. (2003) K-band rows in Table 8 as an alternative epoch. If VV CrA B's alpha drops below 0.3, or if S CrA A's classification changes, then the Class I/bridging conclusion is an artifact of mixing epochs and should be reworded; if alpha remains above 0.3 under both alternatives, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 'bridging Class I and Class II' claim rests on the four spectral indices in Table 4, computed in Section 3.7 by combining 2015 NIRC2 flux ratios (Table 3) with 2MASS/WISE, McCabe et al. (2006), and Scicluna et al. (2016) photometry from different epochs (Table 8). The paper itself shows that component flux ratios are strongly time-variable: VV CrA's J ratio changed from 0.25 in 1996 to 0.015 in 2015 (Section 2.2). No uncertainties are quoted for the alpha values. VV CrA B has alpha = 0.37, just above the Class I/flat-spectrum boundary of 0.3, and the text in Section 3.7 even notes that its SED 'appears relatively flat.' Recomputing alpha with the contemporaneous Scicluna et al. (2016) K- and N-band fluxes instead of the mixed-epoch combination could plausibly move it below 0.3, removing the only Class I object and weakening the specific 'bridging' conclusion. The section also calls S CrA A a 'Class II source' on the basis of silicate emission even though its alpha = -0.17 is flat-spectrum, an internal inconsistency that should be resolved before the classification is used as evidence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents multi-epoch Keck/NIRSPEC H-band spectroscopy of all four components of the young binary systems S CrA and VV CrA, supplemented by archival J- and K-band spectra and 2015 NIRC2 adaptive-optics photometry. The authors determine spectral types for the four stars, measure H-band veiling at multiple epochs, measure Br16 equivalent widths and correct them for veiling, derive Paβ and Brγ accretion luminosities and mass accretion rates, and construct component-resolved SEDs from new and archival photometry. The central conclusion is that S CrA and VV CrA are young binary systems bridging the Class I and Class II evolutionary stages, with high accretion luminosities and variable emission lines.","tokens_in":19523,"tokens_out":5639,"duration_ms":51049,"significance":"If the results hold, the paper supplies the first secure spectral classification of all four components, a valuable multi-epoch veiling record, and one of the few component-resolved SEDs for such close binaries. The careful template-based spectral typing, the explicit treatment of the J- and K-band equivalent widths as lower limits, and the use of archival plus new data are strengths. However, the evolutionary-stage claim rests on epoch-mixed SED construction, and the veiling-versus-equivalent-width correlation is partly algebraic, so the central conclusions are not yet established at the level of rigor claimed in the abstract.","major_comments":[{"comment":"The four spectral indices that drive the \"bridging Class I/II\" conclusion are computed by combining 2015 AO flux ratios (Table 3) with 2MASS/WISE, Prato et al. (2003), McCabe et al. (2006), and Scicluna et al. (2016) photometry from different epochs, and no uncertainty is quoted for any α. The paper itself shows that secondary-to-primary flux ratios change by large factors (VV CrA J ratio: 0.25 in 1996 vs 0.015 in 2015, Section 2.2), so the assumption that the 2015 ratios describe the earlier photometric epochs is not justified. Because VV CrA B's α = 0.37 sits only 0.07 above the Class I/flat-spectrum boundary, a plausible epoch-related shift could move it into the flat-spectrum or Class II category and remove the only Class I object. Please recompute α from single-epoch ratios or quantify how α varies when the ratios are allowed to vary across their observed range.","section":"Section 3.7, Tables 4 and 8"},{"comment":"The text states that silicate emission in S CrA A \"agrees with its spectral index to verify its classification as a Class II source,\" but α = -0.17 lies in the flat-spectrum range (0.3 to -0.3) by the paper's own definition. This is an internal inconsistency in the classification scheme. Please either classify S CrA A as flat-spectrum and discuss the consequence for the evolutionary-stage claim, or provide a definition that makes α = -0.17 consistent with Class II.","section":"Section 3.7, Table 4"},{"comment":"The plotted \"veiling-corrected\" Br16 equivalent width is computed as EW_obs × (1 + rH), and the observed Br16 EWs in Table 5 are nearly constant for each star. The plotted quantity is therefore essentially a linear function of rH by construction, so the reported small p-values for VV CrA A and B do not by themselves demonstrate a physical correlation between accretion emission and NIR veiling. Please add a regression of the observed (uncorrected) EW against rH, or a permutation test that randomizes the observed EWs, to show that the correlation is not introduced by the correction factor.","section":"Section 3.6, Table 7, Figure 4"},{"comment":"The accretion luminosities and mass accretion rates are labeled lower limits because the Paβ and Brγ EWs are not veiling-corrected and because the line EWs, photometry, and flux ratios come from epochs separated by as much as ~13 years, yet Section 5 states the \"high accretion luminosity\" result as a firm conclusion. The text acknowledges a possible ~1 mag uncertainty from the epoch mismatch (citing Gahm et al. 2008); please propagate this into Lacc and Ṁ and state explicitly whether the \"high accretion\" conclusion survives when the input fluxes are varied by that amount.","section":"Section 3.5, Table 6"}],"minor_comments":[{"comment":"The text refers to a \"single epoch (2003/05/23 UT)\" for the Paβ and Brγ measurements, but Table 2 lists no observation on that date; the closest entries are 2003/05/14 and 2003/09/07. Please correct the date or the table.","section":"Section 3.4 and Table 2"},{"comment":"The row label \"VV CrA\" should read \"VV CrA A\" for consistency with the other rows.","section":"Table 6"},{"comment":"The spectral-index equation is written as α = d log(λFλ)/d log(λ) = αλ2 − αλ1; the second equality omits the denominator and should be written as (log(λ2F2) − log(λ1F1))/(log λ2 − log λ1), with the flux convention specified explicitly.","section":"Section 3.7"},{"comment":"The note says α is measured \"between 10.5 and 2.3 μm,\" which reverses the order of the wavelengths given in the text (λ1 ≈ 10.5 μm, λ2 ≈ 2.3 μm); please make the convention consistent.","section":"Table 4 note"},{"comment":"The facility is listed as \"Keck:II (ESI)\" in the acknowledgments, but the observations were made with NIRSPEC; please correct this.","section":"Facility line"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is worth knowing for the measurements, not for the classification. The authors get first spectral types for all four components of S CrA and VV CrA (K7, M1, M1, M0) and a nice multi-epoch H-band veiling series. Both are solid and will be used by anyone working on these systems. The accretion luminosities are lower limits but standard, and they support youth.\n\nThe soft spots are in the SED analysis. The alpha indices in Table 4 combine 2015 flux ratios with photometry from 2MASS/WISE, McCabe 2006, and Scicluna 2016, which are not contemporaneous. The paper itself shows the secondary/primary J ratio for VV CrA changed from 0.25 to 0.015 over two decades. Without error bars, VV CrA B sitting at alpha = 0.37, just above the Class I boundary, is not convincing. Also, the text calls S CrA A a Class II source after showing its alpha = -0.17 is flat-spectrum; that is a direct inconsistency. I could not reproduce the quoted alpha values from Table 8 with the stated definition, so a careful recalculation is needed.\n\nSection 3.6's veiling-EW correlation is another soft spot: the corrected EW is obtained by multiplying observed EW by (1 + r), and with observed EW roughly constant the plotted relation is mostly algebraic. The p-values are not telling you what the text says they tell you. This is secondary, not fatal.\n\nThe central qualitative claim—young, actively accreting T Tauri systems with variable veiling and emission, transitional between Class I and Class II—probably survives. There is enough independent evidence in the spectral types, line emission, and SED shapes. But the specific 'bridging' classification depends too much on the epoch-mixed alpha values.\n\nI would send it to review. A good referee will ask for contemporaneous photometry for the SEDs, error propagation on alpha, and a rewrite of the classification text. The measurements themselves deserve to be published.\n\nRegards.","headline":"First spectral types and multi-epoch veiling for S/VV CrA are genuinely useful, but the SED classification is too shaky to support the 'Class I/II bridging' claim without revisions.","tokens_in":20067,"tokens_out":10976,"would_cite":true,"duration_ms":93528,"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":"All four stars in the S CrA and VV CrA binaries are actively accreting and bridge the Class I–Class II evolutionary stages","keywords":["T Tauri stars","young binary stars","near-infrared spectroscopy","veiling","accretion luminosity","spectral energy distributions","Corona Australis","Class I/Class II transition"],"falsifier":"Take new adaptive-optics JHK and mid-infrared imaging of both systems in a single epoch, measure the component flux ratios and fluxes simultaneously, and recompute the SED spectral indices; if all four $\\alpha$ values fall below $-0.3$, the flat-spectrum and Class I classifications, and with them the bridging-stage conclusion, would be overturned.","tokens_in":19030,"feed_emoji":"⭐","tokens_out":13102,"duration_ms":111423,"temperature":0.7,"pith_summary":"This paper uses five years of high-resolution near-infrared spectra to characterize all four stars in the young binary systems S CrA and VV CrA in the Corona Australis star-forming region. It establishes the first spectral types for the components (K7 for S CrA A, M1 for S CrA B and VV CrA A, M0 for VV CrA B) and shows that their near-infrared veiling and hydrogen emission lines vary strongly from epoch to epoch. Combining these spectra with archival photometry, the authors place the four stars at roughly the same young evolutionary stage, with accretion luminosities higher than those typical of classical T Tauri stars, and conclude that the systems are bridging the Class I and Class II phases. If correct, these closely matched stars offer a controlled comparison for how individual stellar properties and shared environments shape disk accretion and veiling.","feed_headline":"All four stars in two young binaries are still actively accreting","feed_subtitle":"New near-infrared spectra give the first spectral types for the four components and show their disk accretion varies with time.","key_machinery":"The argument rests on simultaneous spectral typing and veiling measurement using the depth ratio of two Fe I lines at $1.562$ and $1.563\\,\\mu$m flanking an unresolved OH doublet near $1.5625\\,\\mu$m in the observed H-band order. Because the OH/Fe ratio is temperature-sensitive in low-mass stars, matching a veiled template to the observed line ratios yields both the spectral type and the H-band veiling at each epoch. Around this core sit the SED spectral index $\\alpha = d\\log(\\lambda F_\\lambda)/d\\log \\lambda$, measured between 2.3 and 10.5 $\\mu$m, which assigns each component to Classes I, flat-spectrum, or II, and the $\\mathrm{Pa}\\beta$ and $\\mathrm{Br}\\gamma$ line-luminosity relations that convert measured equivalent widths into accretion luminosities and mass accretion rates. The same H-band order also contains the $\\mathrm{Br}_{16}$ emission line, whose veiling-corrected equivalent widths are compared with the measured veiling across epochs.","core_discovery":"The central discovery is that this paper establishes the first spectral types for all four components of S CrA and VV CrA and finds that the systems are not classical Class II T Tauri stars but objects in transition between the Class I and Class II embedded stages. All four components are actively accreting from optically thick circumstellar disks, with high, time-variable near-infrared veiling and hydrogen emission lines. The SED spectral indices place S CrA A, S CrA B, and VV CrA A in the flat-spectrum category ($\\alpha \\approx -0.17$, $-0.17$, and $0.13$) and VV CrA B in Class I ($\\alpha = 0.37$), with silicate absorption in VV CrA consistent with an earlier embedded stage. The paper also finds that the infrared companion VV CrA B has an earlier spectral type (M0) than VV CrA A (M1), making it the more massive component despite being fainter at shorter wavelengths because of its high extinction.","pith_inferences":["If the slope of the veiling versus $\\mathrm{Br}_{16}$ equivalent-width relation traces the inner dust truncation radius relative to the corotation radius, then the systematic difference between S CrA and VV CrA could be tested directly by measuring inner-disk radii with infrared or submillimeter interferometry in both systems.","Because component flux ratios varied by factors of several to ten between 1996 and 2015, SED-based evolutionary classifications of close young binaries are only as reliable as the epoch matching between the resolved flux ratios and the broadband photometry; future surveys should prioritize contemporaneous measurements.","The Fe I/OH line-ratio method demonstrated here could be applied to other embedded young binaries to map the Class I–Class II transition across a larger population, provided similar high-resolution H-band spectra are available."],"forward_implications":["If the classification is right, S CrA A, S CrA B, and VV CrA A are flat-spectrum sources rather than Class II, and VV CrA B is Class I, so the two binaries occupy the evolutionary window in which infalling envelopes are giving way to circumstellar disks.","The four components share similar masses, ages, and temperatures, so differences in veiling, accretion, and emission-line strength between them point to systemic or environmental factors rather than gross stellar parameters.","Because VV CrA B is M0 and hence more massive than its M1 primary, the nominal 'secondary' in this infrared companion system is actually the more massive star, hidden by extinction.","The mass accretion rates derived here, roughly $10^{-8}\\,M_\\odot\\,\\mathrm{yr}^{-1}$ and exceeding typical T Tauri values, imply these are relatively young, vigorously accreting systems.","The epoch-to-epoch changes in veiling and $\\mathrm{Br}_{16}$ emission mean that single-epoch snapshots misrepresent the accretion state, so multi-epoch monitoring is required to characterize such systems."],"supporting_citations":[{"why":"Supplies the earlier low-resolution NIR spectroscopy, K-band veiling, and 1996 flux ratios for S CrA that this work recalibrates and extends.","marker":"Prato et al. (2003)"},{"why":"Describes the Fe I/OH line-ratio method used here for simultaneous spectral typing and veiling measurement.","marker":"Prato (2007)"},{"why":"Provides one set of H-band spectral templates used for matching the observed line ratios.","marker":"Bender et al. (2005)"},{"why":"Provides the H-band spectral templates, including the K7 star GJ 281, used to classify the components.","marker":"Prato et al. (2002)"},{"why":"Supplies the mid-infrared photometry for the S CrA component SEDs.","marker":"McCabe et al. (2006)"},{"why":"Supplies the mid-infrared photometry and previous disk modeling for VV CrA, including the comparison accretion rate for VV CrA A.","marker":"Scicluna et al. (2016)"},{"why":"Provides the empirical relation between Paβ/Brγ line luminosity and accretion luminosity used to estimate accretion rates.","marker":"Alcalá et al. (2017)"},{"why":"Provides the pre-main-sequence isochrones used to derive ages, masses, radii, and stellar luminosities.","marker":"Baraffe et al. (2015)"},{"why":"Provides the median SEDs of young stars of matching spectral type used for comparison in the SED plots.","marker":"Fang et al. (2013)"},{"why":"Identified VV CrA B as an infrared companion, the classification this paper refines and ages.","marker":"Koresko et al. (1997)"}],"fun_headline_variants":["Young binaries reveal first spectral types: all four stars still feeding","Four stars, two systems: bridging the gap between embedded and mature","S CrA and VV CrA: first spectra show four stars bridging Class I-II","Two baby binaries caught mid-makeover: all four still accreting"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 2015 secondary-to-primary flux ratios, which the paper itself shows changed dramatically over earlier epochs, are representative of the epochs when the older JHK and mid-infrared photometry used to build the SEDs was taken; if not, the derived spectral indices could cross the class boundaries.","fun_headline_variants_meta":{"raw":{"variants":["Young binaries reveal first spectral types: all four stars still feeding","Four stars, two systems: bridging the gap between embedded and mature","S CrA and VV CrA: first spectra show four stars bridging Class I-II","Two baby binaries caught mid-makeover: all four still accreting"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00093,"raw_usage":{"total_tokens":3977,"prompt_tokens":932,"completion_tokens":3045,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":2964}},"tokens_in":548,"tokens_out":3045,"duration_ms":20116,"temperature":1.0,"reasoning_tokens":2964,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:55:18.590360+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take new adaptive-optics JHK and mid-infrared imaging of both systems in a single epoch, measure the component flux ratios and fluxes simultaneously, and recompute the SED spectral indices; if all four $\\alpha$ values fall below $-0.3$, the flat-spectrum and Class I classifications, and with them the bridging-stage conclusion, would be overturned.","supporting_citations":[{"cited_title":"2005, AJ, 129, 402, 10.1086/426331","cited_arxiv_id":null,"evidence_quote":"Provides one set of H-band spectral templates used for matching the observed line ratios."},{"cited_title":"2002, APJ, 569, 863, 10.1086/339397","cited_arxiv_id":null,"evidence_quote":"Provides the H-band spectral templates, including the K7 star GJ 281, used to classify the components."},{"cited_title":"M., Prato , L., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the mid-infrared photometry for the S CrA component SEDs."},{"cited_title":"S., van Boekel, R., et al","cited_arxiv_id":null,"evidence_quote":"Provides the median SEDs of young stars of matching spectral type used for comparison in the SED plots."}],"review_version":1}