{"id":"464d7316-e0ed-42d1-a645-e214a632b250","arxiv_id":"2506.11998","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":16,"one_line_summary":"Photometry and transit-shape modeling of TOI-3884 reveal an 11-day stellar rotation period and a polar starspot, confirming a misaligned orbit (true obliquity about 77 degrees) for the hot Neptune TOI-3884 b.","lead":"New ground-based photometry of the M dwarf TOI-3884 shows its star rotates every 11.02 days and that its large starspot sits near the stellar pole. This confirms that the transiting planet TOI-3884 b travels on a strongly misaligned orbit, crossing over the polar spot, and makes the system a benchmark for studying polar starspots on M dwarfs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline obliquity ψ⋆=77.4° is prior-dominated: with the alternative stellar density of Mori et al. (2025), the same spot-crossing data yield ψ≈118°, and the transit duration cannot break the b–ρ⋆ degeneracy.","rationale":"The paper's central claim is twofold: (1) TOI-3884 b is on a misaligned orbit, and (2) the true obliquity is ψ⋆ = 77.4°, with a polar spot at 80.5° latitude. The first claim is robust: the 11.020-day rotation period is clearly detected in the Tierras photometry, the spot-crossing phases shift between transits as expected for a rotating spot, and the independent Mori et al. (2025) analysis also finds a strongly misaligned orbit and a polar spot. The second, quantitative claim is the load-bearing one for the paper's status as a benchmark system, and it is exactly where the analysis is least secure. The reader's weakest-assumption analysis identified the same soft spot: the adopted M⋆/R⋆ priors set ρ⋆, which sets b, and ψ⋆ shifts by ~40° under the alternative stellar parameterization. My proposed test would settle whether the Tierras data themselves, freed from those priors, actually prefer b≈0.03 over b≈0.40; absent that test, the quoted 1σ range on ψ⋆ is not a fair representation of model uncertainty. The single-spot model's failure on Transit 258 (§3.2, Fig. 5) is a secondary concern: it shows the spot surface is more complex than modeled, but it is confined to one of six transits and does not directly change ψ⋆ in the same drastic way; it reinforces the need for a robustness check but is not the primary reason the numerical obliquity is uncertain. The ZTF persistence claim is explicitly qualified by the authors (only r-band FAP < 1%, §4.2) and is not required for the misalignment conclusion. I therefore agree with the reader's conditional verdict: accept the qualitative misalignment, but withhold full confidence in the specific obliquity until the density/impact-parameter degeneracy is resolved.","tokens_in":20327,"tokens_out":8043,"duration_ms":182600,"concrete_test":"Re-run the joint Tierras + FLWO 1.2-m + LCO SSO fit of §3.2 with the stellar mass and radius priors replaced by a wide uniform prior on ρ⋆ (e.g., 5–25 g cm⁻³) and a wider eccentricity prior (e.g., N(0.06, 0.10)), keeping all other priors identical. Compare the posterior on b and ψ⋆: if it contains both (b≈0.03, ψ⋆≈77°) and (b≈0.40, ψ⋆≈118°) with comparable support (e.g., Δln Z < 2, or a bimodal ψ⋆ posterior spanning both values within the 95% credible interval), then the headline obliquity is not robust to the stellar-density assumption and the reported 1σ errors are understated. If the wide-prior fit still selects b≈0.03 and ψ⋆≈77° with >95% posterior mass, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, ψ⋆ = 77.4° ± 2.5° (abstract; Table 2), is inherited from the Gaussian priors placed on stellar mass and radius in §3.2 (M⋆ = 0.298 ± 0.018 M⊙, R⋆ = 0.302 ± 0.012 R⊙, from Libby-Roberts et al. 2023). Those priors fix ρ⋆ ≈ 14.4 g cm⁻³ and hence a/R⋆ = 25.06, which drives the impact parameter to b = 0.029 and ip to 89.93°. The true obliquity is then ψ⋆ = 77.4° from i⋆ = 22.3° and λ⋆ = 305.1°. Mori et al. (2025) independently derive b = 0.40 and ψ⋆ ≈ 118° from a different stellar parameterization. The authors defend their b ≈ 0 solution in §4.5 by noting that b = 0.40 requires ρ⋆ ≈ 11.85 g cm⁻³, lower than the LR23 value of 15.3 ± 2.0 g cm⁻³. But the Tierras transit duration is insensitive to this difference: for fixed orbital period, transit duration scales as ρ⋆⁻¹ᐟ³√(1−b²), and (ρ⋆=14.4, b=0.03) vs (ρ⋆=11.85, b=0.40) give nearly identical durations. With the eccentricity free (e = 0.042 ± 0.044), the duration constraint is even weaker. Therefore the b ≈ 0 solution is not demanded by the light-curve shape; it is imposed by the adopted stellar priors. A lower density shifts the chord outward, requires a different i⋆ to reproduce the same spot-crossing phases, and changes ψ⋆ by ~40°, outside the paper's quoted error bars. The qualitative misalignment and polar spot are independently corroborated, so the title-level claim is safe, but the specific benchmark obliquity is prior-dominated, not data-dominated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents Tierras Observatory photometry of TOI-3884 spanning 2024 November 13 to 2025 May 30, including six transits of TOI-3884 b and sparse out-of-transit monitoring. A Lomb-Scargle analysis finds a stellar rotation period of 11.020 +/- 0.015 days, with supporting but weaker archival ZTF detections. The authors simultaneously fit the rotational modulation and the transit shapes with a starry model containing one starspot, stellar rotation, and a Keplerian orbit, obtaining a near-pole-on star (i_star = 22.3 deg) with a large polar spot (r_spot = 31.2 deg, latitude 80.5 deg) and a strongly misaligned orbit (lambda_star = 305.1 deg, psi_star = 77.4 deg). They discuss a positive/negative pole degeneracy, apply the model to archival TESS transits, provide JWST spot-phase predictions, and compare with the independent analysis of Mori et al. (2025). The qualitative picture of a misaligned orbit crossing a persistent polar spot is well supported by the data, but the numerical obliquity depends on the adopted stellar mass and radius priors.","tokens_in":20718,"tokens_out":7667,"duration_ms":92957,"significance":"If the result holds, TOI-3884 b would be a benchmark system: a hot Neptune around an M dwarf on a strongly misaligned orbit crossing a long-lived polar spot, with JWST observations already scheduled. The paper has several concrete strengths: an independent 11-day rotation detection from Tierras and ZTF, a global simultaneous model with 37 parameters over 2019 points (reduced chi-squared = 0.94), an explicit treatment of the positive/negative pole degeneracy, and falsifiable JWST spot-longitude predictions. The central quantitative claim is not yet robust, however. The quoted psi_star = 77.4 deg arises from Gaussian priors on the stellar mass and radius that fix the impact parameter, and the transit duration alone cannot break the b-rho_star degeneracy because the duration scales as rho_star^-1/3 sqrt(1-b^2). The paper therefore needs a robustness demonstration before the specific obliquity value can be adopted as a benchmark.","major_comments":[{"comment":"The headline value psi_star = 77.4 deg +2.3/-2.5 is prior-dominated rather than data-dominated. The Gaussian priors on M_star and R_star from Libby-Roberts et al. (2023) fix rho_star ~ 14.4 g cm^-3 and hence a/R_star ~ 25.06, which drives b = 0.029 and ip = 89.93 deg. The argument in Section 4.5 that Mori et al.'s b = 0.40 is disfavored because it requires rho_star = 11.85 g cm^-3 below the LR23 value is not decisive: for a fixed orbital period, the transit duration scales as rho_star^-1/3 sqrt(1-b^2), so (rho_star=14.4, b=0.03) and (rho_star=11.85, b=0.40) produce nearly identical durations, and the fitted eccentricity (0.042 +/- 0.044) weakens the duration constraint further. The authors should demonstrate directly that the Tierras light curves select b~0 in the absence of the LR23 rho_star prior, or marginalize psi_star over the full allowed stellar-density range, before quoting the obliquity with +/-2.5 deg errors. As written, the abstract's numerical claim is not robust to the b-rho_star degeneracy.","section":"Section 3.2, Table 2, and Section 4.5"},{"comment":"The single-spot model is explicitly a poor fit to Transit 258 in both the Tierras and FLWO g'-band data, and the two-spot test in Section 3.2 locks all stellar and planetary parameters to the one-spot best-fit values rather than treating the second spot as part of the full global model. Because the spot-crossing morphology is the diagnostic that constrains i_star, lambda_spot, and psi_star, an epoch containing an additional spot crossing implies that the one-spot parameter uncertainties may be underestimated. I recommend either including a second spot in the global fit (or marginalizing over its presence) or explicitly checking that the psi_star posterior is unchanged when Transit 258 is removed from the fit.","section":"Section 3.2 and Figure 5"}],"minor_comments":[{"comment":"The comparison with Mori et al. reports two mutually inconsistent values of the true obliquity in consecutive paragraphs: psi = 118.1 deg +5.6/-2.3 in the first paragraph and psi = 61.9 deg +2.3/-5.6 in the second. Please correct the typo and reconcile the text, because the current wording makes it difficult to assess the claimed 4.6-sigma discrepancy.","section":"Section 4.5"},{"comment":"The caption contains a duplicated phrase, 'with i_star = 22.3 deg, lambda_star = 305.1 deg, ip = 89.934 deg, and and lambda_spot = 80.5 deg'; the extra 'and' should be removed.","section":"Figure 2 caption"},{"comment":"The text refers to 'LSO SSO 1-m' in the description of the fitted data sets; this should read 'LCO SSO 1-m' for consistency with Section 2.2 and Table 1.","section":"Section 3.2"},{"comment":"The archival ZTF support for a seven-year spot lifetime is weak: the i-band periodogram peak has FAP = 19% and the g-band peak has FAP = 4.7%, with only the r-band peak significant at FAP = 0.046%. The abstract's statement that the spot 'has persisted for at least seven years' is stronger than the evidence; the body text's 'suggests the possibility' is appropriately hedged.","section":"Section 4.2 and Abstract"},{"comment":"The MCMC description reports 62 walkers and 100,000 steps but does not provide a convergence diagnostic (e.g., autocorrelation time or Gelman-Rubin statistic). Adding one would strengthen the reported uncertainties.","section":"Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The Tierras dataset is valuable and the qualitative result is likely correct, but the headline obliquity is not yet robust to the stellar-density/impact-parameter degeneracy. I would not accept the current numerical value as final; a robustness test marginalizing over stellar priors or a direct model comparison with the Mori et al. parameterization should be required. I see no citation or novelty concerns, and the paper fits the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nTwo things to know about Tamburo et al. The qualitative result is secure: TOI-3884 b crosses a large spot near the pole of its M dwarf host on a substantially misaligned orbit. The quantitative headline — ψ⋆ = 77.4° ± 2.5° — is not. It is set largely by the Gaussian priors on stellar mass and radius adopted from Libby-Roberts et al. (2023), and it is 4.6σ away from the true obliquity derived by Mori et al. (2025) from a similar spot-crossing analysis.\n\nWhat is genuinely new here is the independent dataset. The Tierras photometry gives a clean 11.020 ± 0.015 day rotation period, the joint model of rotational modulation and six transit spot-crossings is well executed (reduced χ² = 0.94), and the ZTF archival check, while formally significant only in r-band, is a reasonable attempt to constrain spot persistence. The paper also provides practical JWST phase predictions. It is clearly written, transparent about its assumptions, and honest about the discrepancies with Mori et al. Credit is due for that.\n\nThe soft spot, as the stress-test note says, is the b–ρ⋆ degeneracy. Transit duration is insensitive to the difference between the paper's solution (b ≈ 0.03, ρ⋆ ≈ 14.4 g cm⁻³) and Mori et al.'s (b ≈ 0.40, ρ⋆ ≈ 11.85 g cm⁻³) — the scaling gives nearly identical durations. The paper's defense, that ρ = 11.85 is lower than the LR23 value of 15.3, leans on stellar models with their own systematics. If the true density is lower, the same spot-path data would require a different stellar inclination, shifting ψ by roughly 40°, outside the quoted error bars. So the misalignment is robust; the specific obliquity is not.\n\nMinor issues: the single-spot model fits five of six transits well, but Transit 258 needs an ad hoc second spot; the ZTF persistence claim is only significant in r-band; and no data or code appear to be released. None of these change the qualitative conclusion.\n\nThis paper deserves a serious referee. The tension with Mori et al. is scientifically important and referees can require the authors to either run an alternative fit with a lower stellar density prior or present ψ as a function of ρ⋆. I would not cite the headline obliquity as a settled number until that degeneracy is resolved.\n\nMy recommendation: send it to peer review, but ask the authors to address the density degeneracy head-on.","headline":"A well-observed confirmation of a misaligned orbit over a polar spot on TOI-3884, but the headline obliquity is prior-dominated and should be treated as conditional on stellar parameters.","tokens_in":21456,"tokens_out":5904,"would_cite":false,"duration_ms":61761,"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":"TOI-3884 b travels on a misaligned orbit and crosses a large polar starspot on its M-dwarf host.","keywords":["exoplanets","starspots","stellar obliquity","M dwarfs","transit photometry","spot-crossing events","stellar rotation","polar starspots"],"falsifier":"Measure the star's projected rotation velocity $v \\sin i_\\star$ from a high-resolution spectrum: the model predicts $0.57 \\pm 0.04$ km/s, so a value well outside that range would falsify the derived spin geometry, and an independent stellar-density determination from Gaia parallax and stellar models would settle whether the impact parameter is near $0.03$ or near $0.40$, which changes the headline obliquity.","tokens_in":19995,"feed_emoji":"🪐","tokens_out":9292,"duration_ms":97934,"temperature":0.7,"pith_summary":"This paper establishes that the transiting planet TOI-3884 b is not orbiting in its host star's equatorial plane. The authors combine long-term photometric monitoring from the Tierras Observatory with six resolved transit light curves to measure a stellar rotation period of $P_\\mathrm{rot} = 11.020 \\pm 0.015$ days and to show that the persistent bump inside every transit is the planet passing over a single large spot sitting at $80.5^\\circ$ latitude, very close to the star's pole. The best-fit geometry has the star's spin axis tilted so that one pole points nearly at the observer ($i_\\star = 22.3^\\circ$) and the planet's orbit misaligned from the stellar spin by $\\psi_\\star = 77.4^\\circ$. If correct, TOI-3884 becomes the first M-dwarf benchmark in which spot-crossing variations confirm a misaligned orbit and expose a long-lived polar spot to repeated observation.","feed_headline":"TOI-3884 b crosses a polar starspot on a 77-degree tilted orbit","feed_subtitle":"Transit-shape changes over an 11-day rotation prove the orbit is misaligned and map a long-lived polar spot.","key_machinery":"The load-bearing element is a simultaneous model of a rotating spotted stellar surface and the planet's transits, computed with the starry package's spherical-harmonic expansion of surface brightness. A single top-hat spot with radius, contrast, latitude, longitude, and smoothing is placed on the star; as the star rotates with period $P_\\mathrm{rot}$, the spot moves in and out of the transit chord, changing the shape and timing of the spot-crossing bump between epochs, while the same spot produces the out-of-transit sinusoidal modulation. Fitting both data sets together breaks the degeneracy that a static spot model cannot resolve and directly constrains the angle between the stellar spin axis and the planetary orbit normal, $\\psi_\\star$.","core_discovery":"The central claim is that TOI-3884 b has a strongly misaligned orbit and that the spot-crossing events seen in every transit occur because the planet passes over a large spot located very close to the visible rotational pole of its M4 host. The authors fit the star's sinusoidal 1% rotational modulation and the epoch-to-epoch changes in transit shape simultaneously, yielding $P_\\mathrm{rot} = 11.020 \\pm 0.015$ days, a spot radius of $31.2^\\circ$, a spot latitude of $80.5^\\circ \\pm 1.2^\\circ$, and a true stellar obliquity of $\\psi_\\star = 77.4^{+2.3}_{-2.5}$ degrees with a stellar inclination of $i_\\star = 22.3^{+1.8}_{-1.6}$ degrees. This rules out the alternative that the planet's orbital period is synchronized with stellar rotation so the same spot is always under the transit chord. The model also explains why some earlier TESS transits showed no obvious spot crossing: at certain rotation phases the spot rotates out of the transit chord. Archival Zwicky Transient Facility photometry shows an approximately 11-day signal across roughly seven years, suggesting the polar spot is long-lived.","pith_inferences":["A test the authors do not run: compare the spot longitudes predicted for the two JWST programs against the actual transit shapes; a systematic offset would reveal spot migration or differential rotation, while agreement would extend the spot's stability baseline.","If the polar spot really survives for years on a star rotating as slowly as 11 days, it would weaken the usual link between polar spots and rapid rotation and motivate dynamo models for slowly rotating convective stars.","The same spot-crossing technique could be applied to other M-dwarf planets with persistent transit bumps, giving an obliquity measurement path for cool stars where Rossiter–McLaughlin spectroscopy is difficult."],"forward_implications":["The aligned-scenario alternative is ruled out because the stellar rotation period is not an integer fraction of the orbital period: $P_\\mathrm{rot}/P_\\mathrm{orb} = 2.4249 \\pm 0.0033$.","TOI-3884 becomes a benchmark for studying polar starspot evolution on an M dwarf, with photometric evidence that the spot has persisted for at least seven years.","The model's predicted spot longitudes (Table 3) give observers a direct handle on how starspot contamination will affect the JWST Cycle 3 transmission spectra of TOI-3884 b.","The planet joins the small population of misaligned hot Neptunes around cool stars that favor nearly polar orbits, which may point toward disk torquing or secular perturbations by an unseen companion.","If the spot remains stable, future transits will cross it at different angles, allowing a map of the stellar pole and direct constraints on spot latitude drift and spot lifetime."],"supporting_citations":[{"why":"Supplies the stellar mass and radius priors, the transit ephemeris, and the polar-spot scenario that this paper confirms.","marker":"J. E. Libby-Roberts et al. 2023"},{"why":"Provides the discovery transit parameters and the limb-darkening central values used in the modeling.","marker":"J. M. Almenara et al. 2022"},{"why":"Provides the starry spherical-harmonic surface model used to fit the rotating spotted star and transits.","marker":"R. Luger et al. 2019"},{"why":"Provides the emcee MCMC sampler used to explore the 37-parameter model.","marker":"D. Foreman-Mackey et al. 2013"},{"why":"Provides the ZTF archival photometry that suggests the polar spot has persisted for about seven years.","marker":"E. C. Bellm et al. 2018"},{"why":"Presents the independent analysis whose higher impact parameter and different obliquity value are compared and reconciled.","marker":"M. Mori et al. 2025"},{"why":"Establishes the spot-crossing method for identifying misaligned planetary orbits, which this paper applies to an M dwarf.","marker":"R. Sanchis-Ojeda & J. N. Winn 2011"},{"why":"Supplies the SPHINX M-dwarf spectral models used to compute spot contrast from sampled stellar and spot temperatures.","marker":"A. R. Iyer et al. 2023"}],"fun_headline_variants":["Spot-crossing pattern reveals planet's 77-degree tilt","TOI-3884 b's tilted orbit confirmed by polar starspot","Misaligned exoplanet crosses starspot on M dwarf","Polar starspot tracks prove TOI-3884 b's orbit tilt","Persistent polar spot confirms misaligned orbit of TOI-3884 b"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline obliquity value assumes the star's mass and radius (and therefore its density) follow the Gaussian priors from the earlier Libby-Roberts analysis; if the true density is lower, the transit chord would sit farther from the star's center and the derived tilt would shift from about $77^\\circ$ to about $118^\\circ$, so the precise number is only as trustworthy as those priors.","fun_headline_variants_meta":{"raw":{"variants":["Spot-crossing pattern reveals planet's 77-degree tilt","TOI-3884 b's tilted orbit confirmed by polar starspot","Misaligned exoplanet crosses starspot on M dwarf","Polar starspot tracks prove TOI-3884 b's orbit tilt","Persistent polar spot confirms misaligned orbit of TOI-3884 b"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000791,"raw_usage":{"total_tokens":3575,"prompt_tokens":1125,"completion_tokens":2450,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":741,"completion_tokens_details":{"reasoning_tokens":2356}},"tokens_in":741,"tokens_out":2450,"duration_ms":19596,"temperature":1.0,"reasoning_tokens":2356,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T01:00:45.520422+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the star's projected rotation velocity $v \\sin i_\\star$ from a high-resolution spectrum: the model predicts $0.57 \\pm 0.04$ km/s, so a value well outside that range would falsify the derived spin geometry, and an independent stellar-density determination from Gaia parallax and stellar models would settle whether the impact parameter is near $0.03$ or near $0.40$, which changes the headline obliquity.","supporting_citations":[{"cited_title":"Multi-band, Multi-epoch Photometry of the Spot-crossing System TOI-3884: Refined System Geometry and Spot Properties","cited_arxiv_id":"2506.06445","evidence_quote":"Presents the independent analysis whose higher impact parameter and different obliquity value are compared and reconciled."}],"review_version":1}