{"id":"97d947d6-fe91-408c-910c-a50bf8245e86","arxiv_id":"1908.07256","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The 65.767-day brightness cycle of the Herbig AeBe star R CrA is best explained by a close binary (about 3.0 and 2.3 solar masses) partially hidden by a circumbinary disk, with secular dimming from disk precession.","lead":"R CrA, a bright young star, varies with a 65.8-day cycle whose cause was unknown. This paper argues the cycle comes from a close pair of stars moving behind the rim of a surrounding disk, and that the star's long-term dimming comes from the disk slowly tipping over.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Binary masses and the 61-yr M-dwarf period rest entirely on the assumed W2 luminosity proxy; the internal cross-check shares the same 17% reprocessing assumption.","rationale":"The reader identified the W2 luminosity proxy as the weakest assumption, and this is also the most load-bearing condition for the paper's quantitative claims. The masses, binary separation, and the predicted M-dwarf period all scale from the single-star luminosity derived from the W2 proxy; if that calibration is invalid for R CrA, the central numbers in the abstract and conclusions move. The paper's internal check using extinction-corrected magnitudes is not fully independent because Eq. (1) uses the same fixed 17% reprocessed fraction, so it does not remove the concern. The proposed RV monitoring test is the cleanest way to settle the matter: it directly measures the dynamical masses without relying on the W2 proxy, and it is explicitly suggested by the authors as future work. Since the paper already frames the scenario as an interpretation supported by a consistency check rather than a dynamical confirmation, the appropriate verdict remains CONDITIONAL; the reader's conditional assessment already accounts for this weakness, so no change to the verdict is needed.","tokens_in":21418,"tokens_out":8389,"duration_ms":88675,"concrete_test":"Obtain high-resolution spectra over at least two contiguous 65.767-day cycles and fit a double-lined Keplerian orbit. The model predicts RV semi-amplitudes of roughly ±70 km/s and mass ratio q = 0.77; the measured K1 and K2, combined with the assumed i ≈ 80°, directly yield dynamical masses. If the total mass disagrees with 5.34 ± 0.8 M_sun by more than 1σ, the W2-based luminosity scale is implicated; if it agrees, the central masses survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative core of the central claim — M_A = 3.02 ± 0.43, M_B = 2.32 ± 0.35, total 5.34 ± 0.8 M_sun, separation 0.56 au, and the 61 ± 4 yr companion period — is not dynamical. It is obtained by converting the WISE W2 flux into L_star via M_V = M_W2 + 4.45, which assumes a universal 17% reprocessed fraction and similar inner-disk flaring for all HAeBe stars (Sect. 2). The paper's own cross-check with extinction-corrected M_V partially shares this assumption: the dJ disk-contamination correction (Eq. 1) is computed from the same 17% reprocessing model, so agreement between M_V(W2) and M_V(abs) does not independently validate the fraction for R CrA. If the true reprocessed fraction differs because the central object is a 0.56 au binary or because the disk is seen nearly edge-on, the single-star mass 3.5 M_sun and all derived binary parameters shift, and the precession/triple interpretation loses its quantitative anchor.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes archival and new photometric data of the Herbig AeBe star R CrA, spanning more than a century, together with a single-epoch SINFONI spectrum. It derives stellar parameters using a WISE W2-based luminosity proxy with a fixed 17% disk-reprocessed fraction, identifies a stable 65.767-day period in the optical light curve, and interprets the phased light curve as extinction of a close binary by a circumbinary disk viewed at grazing inclination. From this model the authors derive component masses of 3.02 and 2.32 Msun, a total mass of 5.34 Msun, and a 0.56-au separation, and they propose that the secular dimming is caused by precession of a slightly inclined disk driven by an M-dwarf companion, making R CrA a triple system. The single SINFONI radial velocity of the primary is consistent with the model prediction at the observed phase, but no full orbital RV solution is presented.","tokens_in":21740,"tokens_out":8139,"duration_ms":77816,"significance":"If the interpretation holds, this would be a rare and valuable case of a young intermediate-mass binary embedded in a circumbinary disk, with a well-documented century-long photometric record linking binary/dust geometry to secular disk precession. The paper's strengths include the assembly of the long photometric baseline, the stability of the 65.767-day period across many datasets, the simple physical model that reproduces the phased curves, and the explicit prediction of ±70 km/s RV variations that can be tested with future spectroscopy. However, the quantitative masses and the triple-system scenario rest on a series of model assumptions (universal W2 reprocessing fraction, adopted inclination of 80 degrees, mass-luminosity scaling) and on one radial-velocity epoch, so the central claim is not yet dynamically secured.","major_comments":[{"comment":"The luminosity of R CrA, and hence every mass quoted in the paper, is derived from the assumption that the WISE W2 flux is a universal proxy for total stellar luminosity with a fixed 17% disk-reprocessed fraction. The agreement between MV(W2) and MV(abs) in Fig. 2 is not an independent validation of this fraction for R CrA, because the dJ disk-contamination correction in Eq. (1) is computed from the same assumed 17% reprocessing. R CrA is also included in the Table 1 sample used to set the MV = MW2 + 4.45 zero point, and the inner disk in a 0.56-au binary or a near-edge-on geometry need not have the same reprocessing fraction as the calibration sample. If the 17% assumption fails for R CrA, the single-star mass 3.5 ± 0.5 Msun, the binary masses 3.02/2.32 Msun, the total mass 5.34 Msun, and the 61-yr companion period all shift, weakening the quantitative anchor of the triple-system interpretation.","section":"Sect. 2, Eq. (1) and Fig. 2"},{"comment":"The model has an explicit degeneracy between the inclination of the system and the disk scale height h0, and the adopted inclination of 80 degrees is acknowledged in the text as arbitrary. The absence of eclipses only excludes inclinations larger than about 86 degrees, so the specific geometry claim that the orbit is seen at an angle grazing the circumbinary disk is not uniquely determined by the photometry. The paper should state what constraints, if any, actually set the inclination and should present h0 and the derived flaring angle as conditional on this assumption rather than as measured parameters.","section":"Sect. 5, Modelling the light curve"},{"comment":"The binary masses are not measured dynamically: they are scaled from the single-star mass using the photometric luminosity ratio LB/LA = 0.49 and a pre-main-sequence mass-luminosity relation, and the single SINFONI radial velocity at phase 0.459 is consistent with the model but provides only one point on the orbit. The quoted uncertainties (MA = 3.02 ± 0.43 Msun, MB = 2.32 ± 0.35 Msun) therefore reflect only the propagation of the luminosity error, not the model uncertainty in the mass-luminosity relation or the unresolved degeneracies in the light-curve model. The paper's own call for future RV monitoring underscores that the binary is not yet confirmed by an orbital solution.","section":"Sect. 5, masses; Sect. 6, SINFONI radial velocity"},{"comment":"The claimed period of the M-star companion, 61 ± 4 yr, is computed from Kepler's law using the assumed semimajor axis of 27-28 au and the model-dependent total mass of 5.34 Msun, and the statement that this agrees with the short-term oscillations in phase and h0 is circular because the same total mass entered the period calculation. The ~50-yr spacing of the oscillations in Fig. 8 is not established by a periodogram, so the precession/triple-system conclusion should be presented as a tentative scenario until an RV orbit or a measured astrometric period is available. The concluding sentence 'R CrA is then a triple star' overstates the current evidence.","section":"Sect. 5, final paragraph; Sect. 7"}],"minor_comments":[{"comment":"The text at the end of Sect. 5 refers to 'the short-term oscillation in phase and h0 values shown in Fig. 7', but Fig. 7 is the schematic view of the extinction model; the secular changes are displayed in Fig. 8. Please correct the cross-reference.","section":"Sect. 5, figure cross-reference"},{"comment":"The quantity dJ is not defined; please state whether it is a magnitude offset, a flux ratio, or a color excess, and give its units.","section":"Sect. 2, Eq. (1)"},{"comment":"The statement that the same h0 describes data in all REM bands and hence that the extinction is grey would be more convincing if the per-band fits or a plot of h0 versus wavelength were shown, since h0 is a geometric scale height and its wavelength independence is the key evidence for grey extinction.","section":"Sect. 4.3 / Sect. 5"},{"comment":"The entry 'Temperature log Teff = 3.98 K' has an inappropriate unit: log Teff is dimensionless when Teff is expressed in kelvin. Please remove the 'K' or clarify the notation.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for A&A and presents a plausible interpretation, but the quantitative claims (masses, 61-yr period) go beyond what the data and model assumptions can support. The authors should either obtain dynamical evidence or substantially soften the conclusions and provide a sensitivity analysis of the W2 proxy and inclination assumptions. The single SINFONI RV is a positive test but is not sufficient to confirm the binary dynamically."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the R CrA paper. The core is a simple geometric model: a 65.767-day binary partially eclipsed by a circumbinary disk, with the disk precessing under the influence of an M-dwarf companion. The model reproduces the phased light curves from AAVSO, SuperWASP, ASAS, and REM over a century of data, and the wavelength dependence (strong optical variability, nearly flat K-band) is a nice, clean argument that the modulation is extinction, not accretion or scattering. The authors are also honest about their assumptions — they openly call the 80° inclination arbitrary and note the lack of a dynamical orbit.\n\nWhat is genuinely new: they confirm the period found by Percy et al. (2010), add a physical interpretation, and propose a luminosity calibration for HAeBe stars using WISE W2 as a proxy for bolometric luminosity. That calibration uses a fixed 17% reprocessed fraction, and the cross-check with extinction-corrected magnitudes is not independent — the dJ correction in Eq. (1) is derived from the same 17% model. The stress-test note is right: for an edge-on disk or a binary, the reprocessed fraction could differ, and every derived quantity — the 3.5 M⊙ single-star mass, the 3.02/2.32 M⊙ split, the 0.56 au separation, the 61-yr companion period — shifts accordingly.\n\nAnother soft spot: the precession interpretation is a consistency check, not a prediction. The secular increase in h0 and the phase wobbles are fit after the fact, and the 1500-yr precession period from Zhu (2019) is compared with the observed ~2° change over 100 years. It works, but it is not a unique explanation. The single SINFONI radial velocity agrees with the model, but one RV point does not confirm a binary.\n\nNone of this is fatal. The paper is carefully argued, the model is testable (RV amplitudes ±70 km/s, GAIA photocenter shifts), and the authors flag the weak links. It deserves a serious referee and publication after revision, with the luminosity calibration's circularity made more explicit. A reader working on HAeBe variability, circumbinary disks, or the Coronet region will get real value. I would bring it to reading group and cite it, with a caveat about the masses.","headline":"A plausible and well-tested circumbinary-disk model for R CrA's 65.8-day variability, but the binary masses and the 61-yr triple period are anchored to a luminosity proxy that is not independently validated.","tokens_in":22287,"tokens_out":4501,"would_cite":true,"duration_ms":43427,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"R CrA is a close binary of about 3.0 and 2.3 solar masses whose 65.767-day orbit grazes a circumbinary disk, and its century-long dimming is caused by disk precession driven by an M-dwarf companion, making the system a triple.","keywords":["Herbig AeBe stars","circumbinary disks","stellar multiplicity","disk precession","R Coronae Australis","protoplanetary disks","variable stars","pre-main-sequence stars"],"falsifier":"Measure the system's radial velocity every few days over one 65.767-day orbit. The model predicts periodic variations of order ±70 km/s, with both components contributing near quadrature; finding no such periodic signal, or a velocity curve whose phase does not match the photometric double maximum at the assumed systemic velocity of $-2$ km/s, would rule out the binary interpretation.","tokens_in":21216,"feed_emoji":"⭐","tokens_out":12317,"duration_ms":104905,"temperature":0.7,"pith_summary":"The paper argues that the 65.767-day periodic variability of the Herbig AeBe star R CrA, stable for more than a century, comes from a central binary whose orbit is seen almost edge-on, so that each of the two stars is periodically hidden behind the near rim of a circumbinary disk. It assigns masses of 3.02 and 2.32 solar masses to the two components and attributes the star's long-term fading over the last 120 years to a progressive increase of absorption by that disk. The proposed driver is precession of a slightly inclined disk, induced by an M-dwarf companion found in high-contrast imaging, which would make R CrA a triple system hosting a disk. The case matters because R CrA is one of the nearest Herbig AeBe stars and a bright, observable example of how disk-companion interactions shape what we see from young intermediate-mass stars.","feed_headline":"R CrA's 66-day cycle is a binary grazing a circumbinary disk","feed_subtitle":"If correct, the Herbig AeBe star is a triple system whose precessing disk explains a century of fading.","key_machinery":"The load-bearing mechanism is a photometric model of a two-star system partially hidden by a circumbinary disk: two point sources in circular Keplerian orbit around a common center of mass, viewed at an assumed inclination of 80°, with extinction described by a Gaussian vertical profile $\\exp(-(z/h_0)^2)$, where $h_0$ is the half-thickness of the absorbing slab at the binary's semimajor axis. The two unequal maxima map to the times when each component is farthest from the disk's near side, and the fitted value of $h_0$ controls the amplitude and shape of the curve; its secular increase is the specific evidence for precession. An independent scaling step sets the stellar masses: the absolute magnitude in the 4.6 µm band is taken as a proxy for total stellar luminosity through the calibrated relation $M_V = M_{W2} + 4.45$, equivalent to an inner disk that reprocesses 17% of the stellar light, and this converts photometry into luminosity, mass, and age before the binary split.","core_discovery":"The central claim is that the observed optical light curve of R CrA, with two unequal maxima phased to 65.767±0.007 days, is produced by two stars in a circular Keplerian orbit seen just grazing a co-planar circumbinary disk: as each component moves behind the dusty near side of the disk, its light is attenuated. From the ratio of the two maxima the authors derive a luminosity ratio $L_B/L_A = 0.49\\pm0.07$ and, via a pre-main-sequence mass-luminosity relation, a mass ratio of $0.77\\pm0.05$. Splitting the single-star mass of $3.5\\pm0.5\\,M_\\odot$ in that ratio gives $M_A = 3.02\\pm0.43\\,M_\\odot$ and $M_B = 2.32\\pm0.35\\,M_\\odot$, a total of $5.34\\pm0.8\\,M_\\odot$, and a binary separation of $0.56$ au. Fitted to light-curve segments spanning 120 years, the same model shows the disk's effective vertical height $h_0$ increasing with time, which the paper reads as a slow increase of extinction caused by precession of the outer disk; from a companion mass of $0.25\\,M_\\odot$ at about 30 au, the analytic precession period is about 1500 yr and the implied orbital period of that companion is $61\\pm4$ yr, consistent with the companion's measured radial velocity. R CrA would therefore be a triple system: a close intermediate-mass binary plus a distant low-mass companion, all surrounded by a circumbinary disk.","pith_inferences":["Inference: if the 4.6 µm luminosity proxy holds for other embedded Herbig AeBe stars, the same two-step method (proxy luminosity, then binary light-curve modeling) could recover masses and multiplicities for objects whose optical photometry is heavily extincted, where single-star mass estimates are usually unreliable.","Inference: the binary hypothesis predicts a photocenter wobble of a few milliarcseconds on the 65.767-day orbit; precise epoch astrometry over one cycle could confirm the companion even when spectral lines are too blended to decompose.","Inference: if the precession interpretation is correct, the extinction toward R CrA should oscillate on the roughly 61-yr companion period rather than increase monotonically; the current trend could reverse within a few decades, which would distinguish precession from a true thickening of the disk.","Inference: systems like R CrA may be common - any Herbig AeBe star with a stable period of tens of days, an unequal double-maximum light curve, and near-constant near-infrared brightness is a candidate grazing circumbinary binary, and the model gives a direct way to search for such objects in long-baseline photometric surveys."],"forward_implications":["The central object of R CrA should be listed as a binary with a 65.767-day period, 0.56 au separation, and total mass of 5.34 solar masses, rather than a single Herbig AeBe star, changing its inferred evolutionary state and disk geometry.","A radial-velocity campaign across one orbit should reveal variations of order ±70 km/s, with both stellar components contributing near quadrature; the absence of such a signal would falsify the binary interpretation.","The secular trend in $h_0$ predicts continued slow fading and a measurable change of the disk's effective inclination over the coming decades, consistent with a precession period of about 1500 yr.","Because the model requires grey extinction (the same $h_0$ in $g'r'i'z'$), the circumbinary dust grains must be larger than a few microns, tying the variability to the high $R_V = 4.7$ reddening toward R CrA.","The M-dwarf companion's orbit should have a period near 61 yr, so its position angle and radial velocity should continue to change in a way that can be tracked astrometrically."],"supporting_citations":[{"why":"Provides the M-dwarf companion discovery, the adopted distance of 150 pc, and the infrared spectra whose radial velocity supports the binary and companion orbits.","marker":"Mesa et al. (2019)"},{"why":"Identified the stable 66-day period in the long-term visual light curve that this paper re-determines as 65.767±0.007 days.","marker":"Percy et al. (2010)"},{"why":"Documents the KH 15D analogue of a binary partially hidden by a circumbinary disk, the observational pattern this model extends.","marker":"Aronow et al. (2018)"},{"why":"Supplies the RVb-star variable-extinction scenario used to interpret the R CrA light curve as disk attenuation rather than scattered light.","marker":"Waelkens et al. (1991)"},{"why":"Gives the analytic disk-precession formula used to estimate the ~1500 yr precession period and the few-degree-per-century change in disk inclination.","marker":"Zhu (2019)"},{"why":"Defines the Herbig AeBe sample whose lower envelope sets the MV = MW2 + 4.45 calibration and the 17% reprocessing fraction.","marker":"Hamann & Persson (1992)"},{"why":"Provides the larger comparison sample used to validate the derived ages, masses, and extinctions.","marker":"Vioque et al. (2018)"},{"why":"Supports the inner-disk temperature and sublimation picture that makes near-infrared disk emission roughly proportional to stellar luminosity.","marker":"Lazareff et al. (2017)"}],"fun_headline_variants":["R CrA's 66-day cycle: binary grazing a circumbinary disk","R CrA's century-long fade: a triple system with a precessing disk","R CrA's variability: a close binary plus a precessing disk","R CrA's 66-day pattern reveals binary and disk","Binary and precessing disk solve R CrA's century-long fade"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole luminosity and mass scale - the single-star value near 3.5 solar masses, its split into 3.02 and 2.32 solar masses, the 0.56 au separation, and the 61-year companion period - rests on the assumption that the 4.6 µm brightness of any Herbig AeBe star is a universal proxy for that star's total brightness, with always the same 17% of stellar light reprocessed by the inner disk.","fun_headline_variants_meta":{"raw":{"variants":["R CrA's 66-day cycle: binary grazing a circumbinary disk","R CrA's century-long fade: a triple system with a precessing disk","R CrA's variability: a close binary plus a precessing disk","R CrA's 66-day pattern reveals binary and disk","Binary and precessing disk solve R CrA's century-long fade"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00079,"raw_usage":{"total_tokens":3613,"prompt_tokens":1207,"completion_tokens":2406,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":823,"completion_tokens_details":{"reasoning_tokens":2309}},"tokens_in":823,"tokens_out":2406,"duration_ms":16644,"temperature":1.0,"reasoning_tokens":2309,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:22:08.776753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the system's radial velocity every few days over one 65.767-day orbit. The model predicts periodic variations of order ±70 km/s, with both components contributing near quadrature; finding no such periodic signal, or a velocity curve whose phase does not match the photometric double maximum at the assumed systemic velocity of $-2$ km/s, would rule out the binary interpretation.","supporting_citations":[{"cited_title":"2019, , 624, A4","cited_arxiv_id":null,"evidence_quote":"Provides the M-dwarf companion discovery, the adopted distance of 150 pc, and the infrared spectra whose radial velocity supports the binary and companion orbits."},{"cited_title":"R., Grynko , S., Seneviratne , R., & Herbst , W","cited_arxiv_id":null,"evidence_quote":"Identified the stable 66-day period in the long-term visual light curve that this paper re-determines as 65.767±0.007 days."},{"cited_title":"A., Herbst , W., Hughes , A","cited_arxiv_id":null,"evidence_quote":"Documents the KH 15D analogue of a binary partially hidden by a circumbinary disk, the observational pattern this model extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the RVb-star variable-extinction scenario used to interpret the R CrA light curve as disk attenuation rather than scattered light."},{"cited_title":"2019, , 483, 4221","cited_arxiv_id":null,"evidence_quote":"Gives the analytic disk-precession formula used to estimate the ~1500 yr precession period and the few-degree-per-century change in disk inclination."},{"cited_title":"& Persson , S","cited_arxiv_id":null,"evidence_quote":"Defines the Herbig AeBe sample whose lower envelope sets the MV = MW2 + 4.45 calibration and the 17% reprocessing fraction."},{"cited_title":"2017, , 599, A85","cited_arxiv_id":null,"evidence_quote":"Supports the inner-disk temperature and sublimation picture that makes near-infrared disk emission roughly proportional to stellar luminosity."}],"review_version":1}