{"id":"f6fc0aaf-67c6-40f0-84ad-283bfa3443c5","arxiv_id":"2508.06488","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using 49 young stars, about one-third show a significant projected tilt between stellar spin and outer protoplanetary disk, indicating primordial misalignment is common.","lead":"Astronomers measured the tilt between the spin axes of 49 young Sun-like stars and their planet-forming disks. About one-third show significant misalignment, suggesting some exoplanet and Solar System obliquities may be set at birth rather than by later planet interactions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"16/49 misalignment count lacks a null-hypothesis calibration; the acknowledged equator-on bias is uncorrected, so 33% may overstate the true primordial misalignment fraction.","rationale":"The paper's central claim reduces to the statement that 16 of 49 systems, after Bayesian inference, have star–disk minimum obliquity posteriors whose MAP exceeds twice their lower 1σ width. The central assertion would be secure if this count were a calibrated estimator of the population misalignment fraction. The robustness grid (Figs. 2–3) is a valuable and honest check of R* and Prot systematics, and it shows the qualitative conclusion survives plausible ±15% perturbations. The authors also transparently discuss the equator-on excess in SD-4, including p = 0.054. These are real strengths. But neither suffices to calibrate the classification step: the Δi posterior is the absolute difference of two noisy inclinations, so under perfect alignment the mode of this posterior will not be exactly zero whenever the individual i* and idisk modes differ, and the 2σ_lower rule has no published false-positive rate. Likewise, the equator-on excess is acknowledged but not propagated into the headline fraction. I therefore see the most load-bearing gap as the absence of an end-to-end null test. The reader's chosen weakest assumption (selection bias) is related but not identical, which is why I mark partial agreement. A simple injection-recovery simulation would settle whether the 16/49 count needs to be revised downward. If the null false-positive expectation is only ~2–3 stars and the selection correction is modest, the 'about one-third' conclusion survives; if the null expectation is higher, the central number weakens. The appropriate verdict remains conditional pending this test, so no change from the reader's verdict is needed.","tokens_in":51687,"tokens_out":8075,"duration_ms":100378,"concrete_test":"Generate a null catalog of 49 systems with true obliquity zero: for each adopted disk inclination, set the true stellar inclination equal to the disk inclination, draw observed vsini*, Prot, and R* from the sample's adopted values and uncertainties, run the exact Bayesian pipeline (Eqs. 2–4), and classify each system with the same MAP > 2σ_lower rule. Repeat at least 1000 times and report the median and 95% range of the number classified as misaligned. If the null expectation is significantly above ~2.5 (i.e., 5% of 49), the 16/49 count and the 33% rate need to be corrected for false positives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline '16/49' is a count of stars whose individual Δi posterior satisfies MAP > 2σ_lower. This is not a null-hypothesis-tested statistic. Because Δi is defined as |i* − idisk|, the posterior is nonnegative and its mode is pulled away from 0 whenever the i* and idisk posterior peaks differ, even if the true star–disk obliquity is zero. With typical Δi uncertainties near 10°, a modest random difference in the two measured inclination modes (e.g., 60° vs 50°) can produce a MAP Δi of ~10° with a lower 1σ width ~5°, satisfying the 2σ criterion. The paper never simulates the all-aligned null ensemble to measure how many of the 49 stars would be misclassified at this threshold. A 10% false-positive rate alone would turn 16/49 into ~11/49 (22%). The acknowledged equator-on excess (14 vs ~8 expected, p = 0.054) is a separate, compounding bias in the same direction: a sample enriched in i* ≈ 90° will systematically increase Δi relative to a random-inclination sample, and no correction/reweighting is applied to the headline 33%. Either effect alone could move the central number below the claimed one-third.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper measures minimum star–disk obliquities (Δi = |i* − idisk|) for 49 isolated, young, Sun-like (F6–M) stars with resolved ALMA protoplanetary disks, using stellar inclinations inferred from vsini*, Prot, and R* through the Bayesian expression in Eq. (2). The central result is that 16/49 systems are classified as misaligned at the MAP > 2σ_lower level, corresponding to a claimed primordial misalignment rate of 33%+7/−6; a 25-cell robustness grid over ±30% perturbations of Prot and R* yields Pmis values between roughly 29% and 71%, with the nominal value at 33%. The paper argues this primordial misalignment distribution is broadly consistent with the obliquity distributions of hot/warm Jupiters and that the Sun's ~6° obliquity could be primordial. The authors explicitly note that Δi is a lower limit on the true obliquity and that the sample contains an overabundance of equator-on stars (14/49 with i* > 80°), but they do not correct the headline rate for either effect.","tokens_in":51991,"tokens_out":5525,"duration_ms":71296,"significance":"If the central claim holds, the paper provides the first statistically meaningful direct measurement of primordial star–disk obliquity in isolated low-mass stars, with direct implications for the origin of exoplanet spin–orbit misalignments and for the Solar System's obliquity. The strengths of the manuscript are its transparent Bayesian formulation, the explicit propagation of uncertainties from vsini*, Prot, and R*, the 25-cell systematic robustness grid, the full posterior and light-curve galleries in the supplementary material, and the use of published/open-source tools. The measurement is not circular: Δi is computed from independent stellar and disk inclination measurements. However, the headline 'one-third' depends on a classification threshold whose false-positive rate is not calibrated, and on an acknowledged selection bias that is not corrected; both issues are load-bearing for the main claim.","major_comments":[{"comment":"The 16/49 count is the number of systems for which the MAP of the non-negative, folded Δi posterior exceeds twice the lower HDI width. This is stated as a 'departure from 0° at 95.4% confidence,' but it is not a calibrated false-positive rate. Under the null hypothesis of true alignment, independent measurement noise in i* and idisk (typical Δi uncertainty ~10°) will produce Δi posteriors whose modes are pulled away from zero for some fraction of the sample. No injection/recovery test or fully-aligned synthetic ensemble is reported, so the expected number of false 'misaligned' classifications is unknown. A modest false-positive rate, e.g. 10%, would change the headline from 16/49 to ~11/49 (22%). Please add a null-hypothesis calibration and report a decontaminated or corrected misalignment fraction.","section":"Main text, first Results paragraph; Methods 'Star-Disk Minimum Obliquity, Δi'"},{"comment":"The paper finds 14/49 stars with i* > 80° versus about 8 expected for an isotropic distribution (p = 0.054) and identifies plausible observational selection toward equator-on stars from both Prot detection and vsini broadening. Because Δi is defined as |i* − idisk|, an excess of equator-on stars mechanically increases the mean and tail of the Δi distribution. This selection effect is acknowledged but not propagated into the 16/49 headline or into the Pmis values of Fig. 3. The manuscript should either reweight by the inclination selection function, report results for a random-inclination subsample, or provide an explicit upper bound on the bias. As written, the central '33%' claim is not robust to this acknowledged bias.","section":"Methods 'Characterization of Systematics and Biases'; Supplementary SD-4"},{"comment":"The 25-cell robustness grid is a strength, but the handling of non-physical cells is unclear. The text states that points with veq < vsini by more than 2σ are excluded from further analysis, yet Fig. 3 prints Pmis for every cell, including cells with many faded/non-physical points (e.g., +30% Prot, −30% R* shows Pmis = 71±10%). If non-physical systems are excluded from the KDE, the displayed Pmis values should be recomputed on the reduced sample and reported as such; otherwise the high-Pmis cells cannot be used to support the statement that 'the misalignment fraction remains substantial.' At minimum, please mark cells for which the non-physical exclusion changes Pmis by more than the quoted uncertainty.","section":"Figs. 2–3 and Methods 'Characterization of Systematics and Biases'"}],"minor_comments":[{"comment":"The abstract says 'about one third of isolated young systems exhibit primordial misalignment' without a confidence interval; the value is later quoted as 33+7/−6%. Please include the binomial/credible interval in the abstract and main-text statement, and phrase the claim as a lower limit given that Δi is a minimum obliquity.","section":"Abstract and main text"},{"comment":"The expression '33^!\"#$%' appears garbled; the intended value is 33+7/−6%. Please correct the typesetting throughout.","section":"Main text, Results paragraph"},{"comment":"The title and abstract use 'Sun-like stars' while the sample is dominated by K and M dwarfs, with a few G stars. Please either justify the phrase by the mass range (0.5–1.5 M☉) or soften it to 'low-mass stars' in the title or abstract.","section":"Title and sample description"},{"comment":"Please define Pmis explicitly. It is used interchangeably as a 'misalignment probability' but it is unclear whether it is the fraction of the KDE above 0°, the fraction of posterior draws satisfying the 2σ_lower criterion, or an integral over the modeled Δi distribution.","section":"Fig. 3 caption"},{"comment":"The notation P)./ and σ2*+ is difficult to read; please ensure all symbols are defined in the main text and that the relation to veq = 2πR*/Prot is stated consistently.","section":"Methods, Eq. (2)"},{"comment":"The equator-on bias discussion is candid and useful, but it should be referenced explicitly in the main-text paragraph that introduces the 16/49 fraction, since it directly qualifies that headline number.","section":"Supplementary SD-4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and addresses an important question, and the observational effort is substantial. My main concern is that the headline 'one-third' is presented as a robust measurement while the two dominant systematics—false positives under the null hypothesis and the acknowledged equator-on selection bias—are not quantitatively corrected. Both are fixable with simulations or reweighting, and the paper would be suitable for publication after such an analysis is added and the claims are softened to 'at least one-third' or a corrected value is reported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What's actually new: this is the first robust statistical measurement of star-disk misalignment for young Sun-like stars, with 49 systems compared to 7 in the previous best effort. The 43 new TESS/K2 rotation periods and the careful Bayesian treatment of i* from vsini, Prot, and R* are solid, reproducible work. The 25-cell robustness grid over ±30% in Prot and R* is a genuine strength: it shows the qualitative conclusion—misalignment is not rare—survives large systematic shifts. The compilation of disk inclinations and stellar properties is transparent and thorough. Credit where it's due: this is a useful dataset and the lower-limit Δi distribution is a real measurement.\n\nThe soft spots are real but not disqualifying. The 16/49 misalignment count uses MAP > 2× lower uncertainty, but there's no simulation of the all-aligned null hypothesis. Since Δi is nonnegative and both i* and idisk have noise, this flag can fire even when the true obliquity is zero. A 10% false-positive rate would drop the count to roughly 11/49, or 22%. That's a load-bearing gap, and it's fixable: draw aligned pairs from the same uncertainty distributions and count how many get flagged. Second, the sample is biased toward equator-on stars (14/49 with i*>80° vs ~8 expected, p=0.054). The authors acknowledge this in SD-4 but never apply a correction or reweighting. Since disk inclinations cluster around 30–75°, pushing i* toward 90° inflates Δi in the same direction as the headline. These two effects could move the true rate below one-third.\n\nThere's also the opposite caveat, which the paper handles honestly: Δi is a lower limit on the true obliquity, so the true misaligned fraction could be higher. That cuts both ways, but it means the abstract's \"about one-third\" overstates the precision of what is really a lower-limit measurement affected by selection effects.\n\nOn balance, the paper deserves a serious referee. The sample is valuable, the analysis is mostly careful, and the authors are transparent about their main systematic concern. A good referee should push for a null-hypothesis calibration, a selection-bias-corrected estimate, and a more careful framing of 33% as a lower limit rather than a point estimate. I'd cite this for the sample and the measured Δi distribution even while being cautious about the exact fraction. It's a good reading-group paper precisely because the biases are illustrative.","headline":"A genuine population-level first for primordial star-disk obliquity, but the 33% headline needs a null-hypothesis test and a selection-bias correction before I'd trust the exact number.","tokens_in":52548,"tokens_out":2823,"would_cite":true,"duration_ms":38298,"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":"About one in three young, isolated Sun-like stars forms with a planet-forming disk tilted relative to its spin axis — evidence that spin-orbit misalignment, possibly including the Sun's 6° tilt, can be set at birth.","keywords":["primordial stellar obliquity","star-disk misalignment","protoplanetary disks","T Tauri stars","stellar rotation","disk inclination","planet formation","spin-orbit alignment"],"falsifier":"Re-measure the stellar inclinations of the same 49 disks with a technique that does not require a measured rotation period or vsini* — for example Zeeman-Doppler imaging of spot patterns or asteroseismology — and recount systems whose minimum obliquity is at least 2σ above zero. If the sample's excess of edge-on stars (14 of 49 versus ~8 expected) is a selection artifact, the fraction should drop below 33%; if the excess is real geometry, the fraction should stay at or above 33%.","tokens_in":51591,"feed_emoji":"🪐","tokens_out":15669,"duration_ms":165213,"temperature":0.7,"pith_summary":"This paper asks whether the tilted orbits seen in many exoplanet systems are fixed at birth or imposed later by gravitational dynamics. It answers with a measurement: among 49 young, isolated, Sun-like stars whose planet-forming disks have been imaged out to tens or hundreds of au, 16 (about one in three) have a disk measurably tilted relative to the star's spin axis. Because the quantity measured — the absolute difference between the sky-projected stellar and disk inclinations — is only a lower limit on each system's true obliquity, the authors frame the 33% figure as a minimum occurrence rate of primordial misalignment. The finding matters because it would mean a large share of the spin-orbit misalignment seen in mature planetary systems, possibly including the Sun's own 6° tilt, is written in during star formation rather than carved out later.","feed_headline":"One in three Sun-like stars is born with a tilted planet disk","feed_subtitle":"Survey of 49 young stars finds 16 tilted disks, hinting the Sun's 6-degree tilt may be set at birth.","key_machinery":"The central quantity is the minimum star–disk obliquity, Δi = |i* − idisk| — the absolute difference between the sky-projected stellar inclination and the outer-disk inclination traced by ALMA at tens to hundreds of au. The stellar inclination i* is recovered with a Bayesian posterior that combines the projected rotation velocity vsini*, a rotation period from TESS/K2 photometry, and the stellar radius, then propagates the full probability distributions into a per-star Δi distribution. A system counts as misaligned when the mode of its Δi distribution sits at least two lower-side uncertainties (95.4% confidence) above 0°. The population is further characterized by a kernel density estimate a","core_discovery":"The paper establishes a 33% primordial misalignment rate: 16 of 49 isolated young Sun-like stars are misaligned at 95.4% confidence. The measured quantity, Δi, is the absolute difference between the sky-projected stellar and disk inclinations — a lower limit on the true obliquity. Most systems cluster near alignment (sample average about 17°; population peak 10°–25°), but the tail extends past 60°. Misalignment shows no significant correlation with stellar mass, temperature, radius, rotation period, or disk inclination. Because the stars are gravitationally isolated, the misalignment is argued to be primordial — plausibly turbulent cloud-core collapse or late infall — and the Sun's 6° obliqu","pith_inferences":["Editorial inference: two opposing biases bracket the true rate — the lower-limit nature of Δi pushes the intrinsic misalignment fraction above 33%, while the sample's apparent excess of edge-on stars (14 of 49 beyond 80° versus ~8 expected, p ≈ 0.05) could push it below; an inclination-unbiased sample is the cleanest way to decide.","Editorial inference: a discriminating test of the late-infall channel would be to split the sample by disk radius or mass — the paper tests only stellar properties — since streamer-fed misalignment predicts systematically larger Δi for more extended outer disks.","Editorial inference: the missing piece is the disk's position angle on the sky and the star's spin direction; resolved CO kinematics plus Zeeman-Doppler imaging of just a handful of the 16 misaligned systems would convert this lower-limit distribution into true obliquities.","Editorial inference: applying the same machinery to younger (Class I) and older (debris-disk) populations would turn this single snapshot into an evolutionary sequence, testing whether primordial misalignment persists or decays across the disk lifetime."],"forward_implications":["A formation-theory benchmark: the misalignment fraction stays substantial (roughly 29–71%) across the paper's 25 systematic tests that shift rotation periods and radii by ±30%, so the one-in-three rate is not an artifact of the adopted stellar parameters.","The Sun's 6° obliquity needs no exotic post-formation mechanism: the pre-main-sequence Sun (about K9 at 5 Myr) falls inside the aligned majority of the sample, making a slightly tilted birth disk a natural explanation.","The star–outer-disk Δi distribution resembles the measured obliquity distributions of hot and warm Jupiters, so a meaningful fraction of mature giant-planet misalignments could be primordial rather than dynamical in origin.","Because Δi is a lower limit on the true obliquity, the high-Δi tail implies some systems could harbor extreme misalignments capable of producing polar or retrograde hot Jupiters, though primordial retrograde disks should remain rare.","Misalignment shows no significant correlation with stellar mass, temperature, radius, rotation, or disk inclination, so it behaves like a random initial condition of the star-disk system rather than a tracer of stellar properties."],"supporting_citations":[{"why":"Supplies the Bayesian framework for turning vsini*, rotation period, and radius into a stellar-inclination posterior.","marker":"10"},{"why":"Provides the analytical posterior expression for stellar inclination used here, plus the differential-rotation uncertainty treatment and the prior Δi-style obliquity methodology.","marker":"11"},{"why":"The largest earlier attempt at primordial star-disk obliquities (seven single stars with broad constraints) that this study extends to 49 systems.","marker":"6"},{"why":"Supplies the hot- and warm-Jupiter obliquity comparison sample and the parametric hierarchical-Bayes modeling approach used to characterize the Δi distribution.","marker":"15"},{"why":"Hydrodynamic simulations of turbulent cloud-core collapse predicting obliquities up to about 80°, the theoretical scale invoked to interpret the observed spread.","marker":"18"},{"why":"Mechanism by which massive planets break inner-outer disk coupling and torque the host star, producing large obliquities (roughly 20°–150°) on timescales matching the sample's ages.","marker":"26"},{"why":"Model constraints on primordial star-disk misalignment produced by late gas infall via accretion streamers, one of the two leading scenarios for the observed misalignment.","marker":"27"},{"why":"Helioseismic determination of the solar rotation axis, anchoring the Sun's ~6° obliquity used in the comparison.","marker":"33"},{"why":"Measurement of the Solar System's invariable plane, anchoring the Sun's ~6° obliquity used in the comparison.","marker":"34"}],"fun_headline_variants":["1 in 3 sunlike stars born with tilted planet disks","Sun's 6-degree tilt may date to birth, not later chaos","Most star disks align, but a third start misaligned","One-third of sunlike stars emerge with wonky disks","Tilted disks at birth: 1 in 3 sunlike stars"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The headline one-in-three rate assumes the 49-star sample is not biased toward edge-on stars in a way that inflates the misalignment count: the paper itself finds 14 of 49 stars nearly edge-on where random orientation predicts about 8 (p ≈ 0.05) and, in the Methods and Supplementary SD-4, discusses but does not correct this selection effect.","fun_headline_variants_meta":{"raw":{"variants":["1 in 3 sunlike stars born with tilted planet disks","Sun's 6-degree tilt may date to birth, not later chaos","Most star disks align, but a third start misaligned","One-third of sunlike stars emerge with wonky disks","Tilted disks at birth: 1 in 3 sunlike stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000192,"raw_usage":{"total_tokens":1199,"prompt_tokens":777,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":521,"completion_tokens_details":{"reasoning_tokens":332}},"tokens_in":521,"tokens_out":422,"duration_ms":4585,"temperature":1.0,"reasoning_tokens":332,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T22:40:28.738051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the stellar inclinations of the same 49 disks with a technique that does not require a measured rotation period or vsini* — for example Zeeman-Doppler imaging of spot patterns or asteroseismology — and recount systems whose minimum obliquity is at least 2σ above zero. If the sample's excess of edge-on stars (14 of 49 versus ~8 expected) is a selection artifact, the fraction should drop below 33%; if the excess is real geometry, the fraction should stay at or above 33%.","supporting_citations":[],"review_version":1}