{"id":"40a89212-ec99-4bda-ac5b-d74f99d27490","arxiv_id":"2607.26803","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"RY Lup's 3.75-day RV and LSD variations are compatible with a ~0.04 au binary (q~0.6, e~0.23, i~13°) or with multi-spot activity, but both interpretations have serious physical tensions.","lead":"High-resolution ESPRESSO spectra of the T Tauri star RY Lup show 3.75-day line-profile variations that can be fit either as a close-in binary or as cool spots. Neither model is ruled out; the binary solution conflicts with tidal circularization and the measured v sin i, so the paper calls for denser multi-technique follow-up.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Dual-Gaussian centroids after selective epoch exclusion may not faithfully trace two Keplerian components; the binary solution is under-constrained by the retained data.","rationale":"The Reader correctly isolates the dual-Gaussian centroid extraction and post-hoc epoch selection as the weakest link. The paper itself is admirably cautious (abstract and conclusions leave both hypotheses open and flag the geometric/tidal tensions), so the CONDITIONAL verdict and the framing “compatible with” rather than “evidence for” already match the data quality. My concrete test simply operationalises that concern: if the orbital elements are stable under controlled re-inclusion of the discarded points, the binary solution gains weight; if they shift, the quantitative parameters should be down-weighted further. No stronger internal inconsistency or methodological flaw is required to keep the verdict at CONDITIONAL.","tokens_in":26465,"tokens_out":569,"duration_ms":12525,"concrete_test":"Re-run the MCMC orbital fit of Sect. 3.2 on three controlled subsets: (i) all epochs with dual-Gaussian solutions (no ω-pinning exclusions), (ii) only Season 1+4 (lowest activity), (iii) a bootstrap that randomly restores 50 % of the previously excluded complex-morphology epochs. If the recovered (e, i, K1/K2) change by more than their quoted 1-σ uncertainties, or if the low-e mode disappears, the published binary parameters are selection-dependent and the claim must remain strictly “compatible with”.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim (a~0.04 au SB2 with q~0.65, e~0.23, i~13°) rests on dual-Gaussian centroids extracted from LSD profiles (Sect. 3.2). Several epochs with complex morphology (esp. Season 2/3) are excluded from the orbital fit, and two further points (φ=0.2, 0.7) are dropped specifically “to constrain ω”. The MCMC posterior is bimodal (Fig. C.1); the low-e solution is preferred partly because it is “more physically motivated”. The retained RVs still produce an orbit whose inclination is incompatible with the measured v sin i~32 km s⁻¹ under any plausible stellar radius and spin-orbit state (Sect. 4.1), and whose eccentricity should have been erased by tidal circularisation on ~1–5 Myr timescales. Thus the assumption that the dual-Gaussian centroids faithfully sample Keplerian motion of two stars, rather than residual activity or blending, is the single most load-bearing and least secure step.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper presents a 327-day VLT/ESPRESSO campaign on the classical T Tauri star RY Lup, using LSD line profiles to test whether the known ~3.75 d photometric period arises from a close spectroscopic binary or from cool surface spots. Dual-Gaussian decomposition of the LSD profiles yields RV curves that are fitted with an MCMC Keplerian SB2 model (P≈3.75 d, q≈0.65, e≈0.23, a sini≈0.01 au), which, combined with ALMA dynamical masses, implies a nearly face-on binary (i≈13°) strongly misaligned with the inner and outer discs. In parallel, SpotCCF multi-spot modelling reproduces the profile deformations with 1–3 cool spots, but recovers a ~15% epoch-to-epoch scatter in v sin i and no significant rotational period. The authors conclude that neither hypothesis is ruled out and call for denser RV plus interferometric monitoring.","tokens_in":26730,"tokens_out":1820,"duration_ms":44616,"significance":"RY Lup is a well-studied transition disc with a large cavity and known inner–outer disc misalignment; a secure close companion would be an important dynamical driver and a useful benchmark for cavity-carving models. The work brings high-resolution ESPRESSO time series, standard LSD and ROTFIT analyses, a full MCMC orbital solution, and a Bayesian multi-spot comparison to the same dataset, and it is appropriately cautious in not claiming a definitive detection. The explicit confrontation of the binary geometry with v sin i, tidal circularisation timescales, and GRAVITY inner-rim size is a strength. If the binary solution can be placed on firmer footing (or cleanly rejected), the result would matter for both disc architecture and the interpretation of RV jitter in active PMS stars.","major_comments":[{"comment":"Sect. 3.2, Table 1 notes, and Fig. 4: Several epochs (2023 Jan 15, Mar 24, Mar 27, and the complex Season 2/3 profiles) are excluded from the orbital fit because of complex LSD morphology, and two further points (φ≈0.2 and 0.7; Apr 14 and Mar 2) are dropped specifically “to constrain the argument of periastron.” These cuts are load-bearing for the adopted (P, K1, K2, e, ω) solution in Table 2. The manuscript needs objective, pre-stated rejection criteria (e.g., a quantitative goodness-of-fit or peak-separation threshold) and a sensitivity test showing how P, e, q, and a sini change when the dropped points are restored or when only morphology-based exclusions are applied. Without that, it is unclear whether the dual-Gaussian centroids faithfully sample two Keplerian components.","section":"Sect. 3.2, Table 2, Fig. 4"},{"comment":"Sect. 3.2 and Appendix C (Figs. C.1–C.2): The MCMC posterior is explicitly bimodal, with a low-e and a higher-e (e~0.5) family. The low-e solution is adopted partly because it is “more physically motivated” given the short period. Preferring one mode on physical grounds while the retained solution still has e≈0.23—which the same section later notes should circularise on ~1–5 Myr timescales (Zahn & Bouchet 1989)—is internally in tension. Please report both modes with relative evidence (Δχ² or Bayesian odds), and either justify why e≈0.23 remains acceptable for a ~few-Myr system or treat eccentricity as poorly constrained.","section":"Sect. 3.2, Appendix C"},{"comment":"Sect. 4.1: Under the adopted SB2 masses and i≈13°, spin–orbit synchronisation at P=3.75 d predicts v sin i≈5–6 km s⁻¹, an order of magnitude below the measured ≈32 km s⁻¹ (ROTFIT and SpotCCF). The text correctly flags this as strongly disfavouring a low-inclination binary, yet the abstract and conclusions still present the variations as “compatible with” a close-in binary at ~0.04 au. The abstract and final bullet list should be rebalanced so that the geometric and tidal inconsistencies are stated at the same prominence as the formal orbital fit; “compatible” overstates the current evidence.","section":"Sect. 4.1, Abstract, Sect. 5"},{"comment":"Sect. 3.3 and 4.2: The claim that “spots alone cannot explain the observed variability” rests on (i) ~15% dispersion in retrieved v sin i and (ii) non-detection of a period above 10% FAP in the SpotCCF time series. SpotCCF treats spots as black (T=0 K) and fits each epoch independently with free v sin i; both choices can inflate apparent v sin i scatter and destroy phase coherence. A controlled test with v sin i fixed (or tightly prior-constrained) to the ROTFIT value, and/or with a realistic spot–photosphere contrast, is needed before concluding that activity is insufficient. If fixed-v sin i models still fail systematically, that would substantially strengthen the argument.","section":"Sect. 3.3, Sect. 4.2"}],"minor_comments":[{"comment":"Table 1: Several RV entries are quoted as round numbers with large or identical errors (e.g., 5.0±2.0, −10.0±2.0). Clarify whether these are fit failures, manual assignments, or heavily constrained solutions, and ensure the orbital fit does not silently treat them as high-weight measurements.","section":"Table 1"},{"comment":"Fig. 3: State explicitly which RV series (RV1, RV2, or both) enters the GLS, the frequency grid and oversampling, and whether the 1% FAP is analytic or bootstrap-based. The bootstrapping test (50 draws of 80%) is useful; consider showing the period distribution.","section":"Sect. 3.2, Fig. 3"},{"comment":"Sect. 3.1 / Table A.1: Mean v sin i from ROTFIT is 32.3±1.1 km s⁻¹, while literature values are 16–25 km s⁻¹. A short discussion of why ESPRESSO+ROTFIT yields a higher value (resolution, veiling treatment, template grid) would help non-specialist readers.","section":"Sect. 3.1"},{"comment":"Fig. 9 is schematic only; if retained, label the adopted i_binary, i_inner, i_outer and a values so the misalignment claim is visually quantitative.","section":"Fig. 9"},{"comment":"Typos/notation: “a a Markov Chain” (Sect. 3.2); “GRA VITY” spacing is inconsistent; mass-function lines in Table 2 give M1 sin³i = 0.010 and M2 sin³i = 0.070, which appears swapped relative to q=M2/M1=0.65 (primary should be the more massive component). Please check.","section":"Sect. 3.2, Table 2"},{"comment":"Data availability points to a Zenodo DOI; ensure the deposited material includes the per-epoch LSD profiles, dual-Gaussian centroids, and SpotCCF posteriors needed to reproduce Figs. 4–6.","section":"Sect. 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a careful, well-executed exploratory study and is appropriate in scope for A&A. My major_revision recommendation is driven by the need to harden the data-selection and abstract wording around the SB2 solution, not by any sense that the work is unsound. If the authors supply the requested sensitivity tests and rebalance the abstract, this could move to accept after one round. The M1/M2 sin³i swap in Table 2 looks like a simple labelling error but should be caught before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new thing here is a real ESPRESSO LSD time series on RY Lup plus a side-by-side binary vs SpotCCF treatment. They recover the known 3.75 d period in the RVs, publish dual-Gaussian elements, multi-spot solutions, and R'_HK, and they do not over-claim: abstract and conclusions leave both hypotheses open and flag the geometric and tidal problems themselves.\n\nWhat they do well is the observational craft. ROTFIT parameters, LSD construction, GLS, MCMC, and Bayesian spot comparison are standard and cleanly reported. The Zenodo deposit is there. The discussion of the face-on binary vs v sin i ~32 km/s, the circularisation timescale, and the 15% v sin i scatter under spots is frank. That honesty is useful; a lot of T Tauri RV papers bury the tension.\n\nThe soft spot is real but already mostly owned by the authors. The SB2 solution rests on dual-Gaussian centroids after dropping complex-morphology epochs and two more points to pin ω. The MCMC is bimodal; they prefer the low-e mode partly on physical grounds. The retained orbit (a~0.04 au, q~0.65, e~0.23, i~13°) fights both the measured broadening and tidal expectations. So the quantitative binary parameters are under-constrained by the retained data and should stay framed as “compatible with,” not “evidence for.” Spot modelling fits the bumps better but produces unstable v sin i and no coherent period, so it does not cleanly win either. Sparse sampling is the practical limit.\n\nThis is for people who work on T Tauri RVs, transition-disc companions, or activity mitigation. It is a careful case study, not a definitive detection. I would send it to referees; the data products and the dual analysis deserve the airtime, with the usual request to keep the binary numbers provisional. Worth a look if you care about this system or the method; not urgent if you do not.","headline":"Solid ESPRESSO time series and honest dual analysis of RY Lup; the binary numbers are under-constrained and physically tense, but the paper already says so.","tokens_in":27395,"tokens_out":510,"would_cite":true,"duration_ms":10430,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"ESPRESSO line profiles of RY Lup show a 3.75-day signal that fits either a close-in binary or cool spots, and neither explanation is ruled out.","keywords":["T Tauri stars","spectroscopic binaries","stellar activity","protoplanetary discs","radial velocities","least-squares deconvolution","transition discs","RY Lup"],"falsifier":"A dense, multi-epoch campaign that simultaneously measures radial velocities and interferometric inner-disc geometry (or multi-wavelength line-profile contrasts) would show either a stable Keplerian orbit independent of wavelength or wavelength-dependent spot contrasts that erase the binary solution.","tokens_in":27345,"feed_emoji":"⭐","tokens_out":983,"duration_ms":19181,"temperature":0.7,"pith_summary":"RY Lup is a young star with a large cleared cavity in its disc and a long history of photometric variability at 3.75 days. High-resolution ESPRESSO spectra, averaged into high-signal line profiles, show the same period in radial-velocity shifts. Those shifts can be fit as a double-lined spectroscopic binary with separation about 0.04 au and mass ratio near 0.6, but the implied nearly face-on orbit sits strongly misaligned with both the inner and outer discs and retains an eccentricity that is hard to keep at such a short period. The same line deformations can instead be reproduced by one-to-three cool spots on a single rapidly rotating star, yet the retrieved rotational broadening still wanders by about 15 percent and no clean rotational period appears in the spot parameters. The paper therefore leaves both the binary and the pure-spot interpretations open and calls for denser, multi-wavelength monitoring that can separate orbital Doppler motion from temperature-driven spot contrast.","feed_headline":"RY Lup's 3.75-day signal fits binary or spots—neither wins","feed_subtitle":"ESPRESSO line profiles leave a close-in companion and cool-spot models both viable","key_machinery":"Least-squares deconvolution (LSD) mean line profiles, fitted once as dual Gaussians for a double-lined binary orbit and once with the SpotCCF multi-spot code that freezes no v sin i, so that stability of the recovered broadening becomes a diagnostic of whether spots alone can explain the data.","core_discovery":"The ESPRESSO least-squares-deconvolution profiles of RY Lup contain a coherent ~3.75-day radial-velocity signal that is formally compatible with a close-in spectroscopic binary (a ~ 0.04 au, q ~ 0.65, e ~ 0.23, i ~ 13°), yet the identical profile deformations are also reproduced by multi-spot models; a 15 percent scatter in retrieved v sin i and the absence of a significant period in the spot parameters leave neither hypothesis excluded.","pith_inferences":["The tension between a short-period eccentric orbit and tidal theory already hints that any real binary is either younger than expected or still being torqued by residual disc material.","A joint binary-plus-spots model is the most natural next step; the two effects are unlikely to be mutually exclusive in an accreting T Tauri star.","Wavelength-dependent LSD or simultaneous optical–infrared radial velocities would break the degeneracy because spot contrast falls toward the red while orbital Doppler shifts do not."],"forward_implications":["If the binary solution is real, the system is a nearly face-on close binary strongly misaligned with both its inner (~50°) and outer (~70°) discs, implying a highly warped inner disc.","The 0.04 au separation lies inside the 0.12 au dust rim, so the companion could sculpt the innermost disc even while the 60 au millimetre cavity still needs additional clearing agents.","The retained eccentricity e ~ 0.23 at P = 3.75 d conflicts with standard tidal circularisation timescales for fully convective pre-main-sequence stars.","If spots dominate, the required filling factors and the persistent v sin i scatter imply large-scale, rapidly evolving surface features that pure rotational modulation cannot track with the present cadence."],"fun_headline_variants":["RY Lup 3.75-day RV signal fits binary or spots—neither excluded","ESPRESSO profiles leave RY Lup companion and spot models both open","Close-in binary or cool spots? RY Lup line data can't decide","RY Lup's coherent 3.75-day signal still compatible with both hypotheses","Misaligned binary at 0.04 au or spots: RY Lup RV case unresolved"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That the two Gaussian centroids extracted from the LSD profiles (after dropping several complex epochs) truly track the Keplerian motion of two stars rather than residual activity or blending.","fun_headline_variants_meta":{"raw":{"variants":["RY Lup 3.75-day RV signal fits binary or spots—neither excluded","ESPRESSO profiles leave RY Lup companion and spot models both open","Close-in binary or cool spots? RY Lup line data can't decide","RY Lup's coherent 3.75-day signal still compatible with both hypotheses","Misaligned binary at 0.04 au or spots: RY Lup RV case unresolved"]},"model":"grok-4.5","effort":"low","cost_usd":0.003495,"raw_usage":{"total_tokens":1293,"prompt_tokens":957,"num_sources_used":0,"completion_tokens":105,"cost_in_usd_ticks":34948000,"prompt_tokens_details":{"text_tokens":957,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":231,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":957,"tokens_out":105,"duration_ms":5564,"temperature":1.0,"reasoning_tokens":231,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T20:35:33.856918+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A dense, multi-epoch campaign that simultaneously measures radial velocities and interferometric inner-disc geometry (or multi-wavelength line-profile contrasts) would show either a stable Keplerian orbit independent of wavelength or wavelength-dependent spot contrasts that erase the binary solution.","supporting_citations":[],"review_version":1}