{"id":"54107fa0-39a3-4fb2-802f-c8595bb72bba","arxiv_id":"2506.06445","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"TOI-3884's host star rotates every 11.04 days, and its persistent polar spot is offset roughly 13 degrees from the pole, leading to refined stellar inclination and obliquity measurements.","lead":"New multi-band, multi-epoch transit observations of the star TOI-3884 yield a rotation period of about 11.04 days and show that the starspot which crosses every transit sits about 13 degrees from the stellar pole. The refined geometry resolves earlier discrepancies between published studies and sharpens the target for upcoming atmospheric spectroscopy of the transiting super-Neptune.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported spot latitude and geometry rely on a single circular, non-evolving spot; a second spot in the fit causes non-convergence, and the 1.74-day phase offset in Sec 5.5 suggests the assumed spot configuration drifts, so the quoted parameter shifts may be model-induced.","rationale":"The reader's weakest_assumption identifies the same load-bearing assumption: the single circular, non-evolving spot parametrization in Section 4. I agree with the CONDITIONAL verdict. My stress-test adds two specific pieces of evidence: (1) the paper explicitly reports that fitting multiple spots caused MCMC non-convergence, so alternative spot configurations that could break the geometry degeneracy were not adequately explored; and (2) the forward model to the monitoring light curve requires an ad hoc 1.74-day phase shift, which is a substantial longitude offset (~57 deg) and indicates that the spot configuration was not stable across the combined dataset, contra the modeling assumption. The rotation period itself (11.04 d) is well supported independently by GLS, GP, ZTF, and the transit spot-crossing evolution, and the broad picture of a polar spot and misaligned orbit is robust. However, the precise values of i_star, lambda, and theta_spot that constitute the central claim are model-dependent. A concrete two-spot or evolving-spot refit would settle whether the quoted geometry is unique. Therefore the verdict remains CONDITIONAL, as the reader recommended, because the data support the broad picture but not yet the precise parameter values as definitive.","tokens_in":22226,"tokens_out":1490,"duration_ms":14699,"concrete_test":"Refit the three MuSCAT transits with a two-spot model (e.g., one polar spot plus one smaller, lower-latitude spot) or with a spot that is allowed to evolve slightly in latitude and/or radius between transit epochs, using the same priors and noise treatment. If the resulting i_star and lambda shift by more than the quoted 1-sigma uncertainties (~1-9 deg), or if the posterior becomes multimodal across qualitatively different geometries, then the spot-crossing geometry is not uniquely determined by the data and the quoted values should be treated as model-dependent.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central geometry i_star=139.9 and lambda=41.0 is constrained almost entirely by the spot-crossing bumps and their evolution across the three MuSCAT transits (Section 4). The model fixes one circular spot with constant radius, temperature, and latitude over ~40 days, with only longitude evolving at Prot. The paper acknowledges both that multiple-spot MCMC chains failed to converge and that the model's simplicity may introduce systematic biases (Sec 5.6, Sec 6). This is an internal admission that the uniqueness of the quoted geometry has not been demonstrated. The 1.74-day time offset needed to align the forward-modeled rotational modulation with the monitoring light curve (Sec 5.5) is also larger than the Prot uncertainty integrated over 254 days, implying either a spot longitude/configuration change or inconsistent spot properties between epochs. A 1.74-day offset at 11.04 d corresponds to ~57 deg of spot longitude, not a negligible phase shift. If the spot evolved significantly between epochs, the single-spot assumption for the transit epochs becomes unsupported, and the formal uncertainties on i_star, lambda, and theta_spot understate the true model uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents new multi-band, multi-epoch transit photometry of TOI-3884b from MuSCAT3/4 and ground-based r-band monitoring from LCO/Sinistro. The authors report a stellar rotation period of 11.043 +0.054/-0.053 days, detected independently with GLS and Gaussian-process modeling, and use it to interpret the time-variable spot-crossing signals seen in three transits. Modeling the spot as a single circular, time-invariant feature with only its longitude evolving, they derive a spot latitude of -76.8 deg, a spot radius of 0.425 R_star, a spot-photosphere temperature difference of 200 K, and a refined system geometry with i_star = 139.9 deg and lambda = 41.0 deg. They also reanalyze TESS and ZTF data to argue that the spot configuration changed on multi-year timescales and that the spot lay closer to the pole during the TESS epoch. The paper reconciles several previously discrepant parameters, including the projected rotation velocity and the impact parameter, and discusses implications for transmission spectroscopy and M-dwarf magnetic activity.","tokens_in":22524,"tokens_out":7267,"duration_ms":76393,"significance":"If the derived geometry is correct, TOI-3884 is an unusually valuable system: a nearly polar-orbit super-Neptune transiting a mid-M dwarf with a large spot offset from the stellar pole. The paper provides the first robust rotation-period measurement for this star, a multi-band spot temperature constraint, and an updated obliquity that reconciles the contradictory A22 and L23 results. The work is also relevant to upcoming JWST transmission spectroscopy by quantifying the transit light source effect for this target. Strengths include the use of two independent period-search methods, simultaneous four-band transit photometry, and transparent public release of the custom spocon code and data products. However, the central geometric claims rest on a deliberately simplified single-spot model, and the manuscript itself acknowledges that this model may introduce systematic biases and that more complex spot configurations cannot be uniquely recovered. The significance is therefore real but conditional on the robustness of the spot model to plausible violations of its assumptions.","major_comments":[{"comment":"The quoted spot latitude, stellar inclination, and projected obliquity are derived entirely from the evolution of a single circular spot whose radius, temperature, and latitude are held fixed over the ~40-day span between transits. This is the load-bearing assumption of the paper. The authors note that MCMC chains failed to converge when multiple spots were introduced, but non-convergence of a multi-spot model does not by itself validate the single-spot solution. I request a robustness demonstration: for example, fit each of the three transits independently with the spot longitude (and ideally radius and latitude) free, check whether the recovered phases are consistent with a single rotating spot, and run injection-recovery tests with a second spot to quantify how easily the single-spot model could absorb a two-spot configuration. Without such a test, the formal 1-sigma intervals on i_star, lambda, and theta_spot are conditional on the single-spot model and likely underestimate the true model uncertainty.","section":"Sec. 4 and Table 3"},{"comment":"The forward-modeled rotational modulation requires a 1.74-day phase offset to align with the Sinistro light curve. At P_rot = 11.04 d this is roughly 57 degrees of spot longitude, which is comparable to or larger than the ~1.2-day (1-sigma) uncertainty propagated from the period uncertainty over the 254-day gap. The same figure also shows that the model fails to reproduce the lowest observed flux levels, which the authors attribute to additional spots or more complex spot geometries. This indicates that the spot configuration during the monitoring epoch is not identical to the configuration assumed for the transit epochs, and it weakens the assumption of a time-invariant single spot over the full interval. I ask the authors to quantify how much the derived geometry shifts if spot parameters are allowed to evolve between the three transits, or if the forward-model comparison is restricted to the transit epochs rather than extrapolated across the 254-day gap.","section":"Sec. 5.5 and Fig. 7"},{"comment":"The TESS reanalysis is presented as supporting the idea that the spot was closer to the pole during the TESS epoch, but this analysis is not independent of the MuSCAT-derived geometry. The priors in Table 5 are taken directly from the MuSCAT posterior, including i_star ~ N(139.9, 2.0), lambda ~ N(41.0, 9.0), rho_star, T_phot, and P_rot. The derived TESS spot latitude of -88.1 deg and lambda of 26.2 deg are therefore conditioned on the MuSCAT solution. I recommend rerunning the TESS fit with broad, physically motivated priors, or with priors taken only from A22 and L23, and reporting how the recovered spot latitude changes. The conclusion about a more polar spot during the TESS epoch should be reframed as conditional on the MuSCAT geometry unless the fit is shown to be robust to prior choice.","section":"Appendix C and Table 5"},{"comment":"The posterior from the transit fit is multimodal, and the authors impose U(90, 180) on i_star and U(10.8, 11.3) on P_rot after inspecting the MCMC output. This post-hoc selection may be physically justified, because the -40 deg and 140 deg stellar-inclination solutions correspond to the same axis with the opposite pole visible, and because the lower-likelihood period solutions are clearly disfavored. However, the manuscript should state this equivalence explicitly and should show the marginalized posteriors before and after the prior restriction, so the reader can assess how much of the quoted central values and uncertainties comes from the prior choice. At present, Table 3 reports only the restricted solution, which makes the influence of the post-hoc prior difficult to evaluate.","section":"Sec. 4 and Sec. 5.6"},{"comment":"The revised transit geometry (b = 0.402 ± 0.019, rho_star = 11.85 ± 0.30) differs substantially from both A22 and L23, which reported near-zero impact parameters and rho_star of 14.3 ± 1.1 and 15.26 ± 2.04. The authors attribute this to the geometric leverage of the spot-crossing features, while also noting that the spot-model simplicity could introduce systematic bias. Because b and a/R_star are strongly degenerate in ordinary transit fitting, and because the spot model is simplified, I request a control analysis: for example, fit the same light curves with the spot-crossing regions masked, or with a different spot model (e.g., two spots or a non-circular spot), and report whether b and rho_star remain consistent. This would clarify whether the geometric revision is a robust consequence of the data or a consequence of the assumed spot parameterization.","section":"Sec. 5.6 and Table 4"}],"minor_comments":[{"comment":"The abstract contains the typo 'poler orbit' in place of 'polar orbit'; please correct it.","section":"Abstract"},{"comment":"The phrase 'photsphere-spot temperature difference' is missing a letter; it should read 'photosphere-spot temperature difference.'","section":"Sec. 5.3"},{"comment":"The table would be clearer if the 'conditioning priors' (rho_star and T_14) were separated explicitly from the sampled parameters; as written, a reader may mistake them for directly fitted quantities.","section":"Table 3"},{"comment":"The caption notes that the apparent overlap in the spot-latitude contours is due to smoothing and not actual mixing; this is an important clarification, but the figure would be more informative if the two groups were shown without over-smoothing or with the raw density contours.","section":"Appendix A and Fig. 9"},{"comment":"The sentence explaining the 1.74-day offset says it 'can also be explained by error propagation'; given the calculation in the text, it would be clearer to state the propagated uncertainty explicitly (roughly 1.2 days at 1 sigma) so the reader can see that the offset is approximately 1.4 sigma rather than being negligible.","section":"Sec. 5.5"},{"comment":"The manuscript should check for consistency between the statement in Sec. 5.6 that the absolute stellar inclination is 'approximately 40 deg' and the tabulated i_star of 139.9 deg; the relationship between these two conventions is explained, but it is easy for a reader to misread as a contradiction.","section":"Table 4 and Sec. 6"}],"recommendation":"major_revision","confidential_remarks":"This is a valuable observational paper with a clear dataset and a transparent analysis, and the central claims are defensible if the requested robustness tests are added. My main concern is that the paper's headline geometric parameters are extracted from a deliberately simplified spot model whose limitations are acknowledged in the text but not quantitatively evaluated. The requested tests — per-transit spot fits, injection-recovery with multiple spots, a prior-free TESS reanalysis, and control fits with masked spot crossings — are all within the scope of the existing data and should be feasible. The manuscript also contains an unusual acknowledgment about correcting the inclination-vector direction after arXiv posting; the current version appears internally consistent, but I would recommend a final pass to confirm that all sign conventions and the text in Sec. 4 and Appendix A agree."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: TOI-3884 finally has a measured rotation period, and the paper makes a good case that its large spot sits about 13 degrees from the pole rather than exactly on it. The new MuSCAT3/4 transits and LCO monitoring are well reduced, the GLS detection at 11.04 d is secure (FAP < 1e-12, and ZTF shows the same period in several seasons), and the multi-band spot temperature (Tspot ~2950 K, deltaT ~200 K) is consistent with A22 and the empirical spot-temperature relation. The paper is transparent about its degeneracies and about where the model is simplistic. That transparency is a real strength.\n\nThe main soft spot is the single circular, non-evolving spot used for the transit fits. The chains didn't converge with two spots, and the posterior shows a bimodality that the data can't resolve; the quoted i_star = 139.9 and lambda = 41.0 come from a model that is acknowledged to be simple. The post-hoc priors U(90,180) on i_star and the P_rot window are used to pick a branch, though the 11.04-day branch was already favored by log-probability, so that is not circular. The stellar density comes out at 11.85, about 2 sigma lower than L23's 15.26, which is a bit uncomfortable given the prior comes from L23. The TESS reanalysis gives a spot at -88 deg, which suggests the spot configuration does evolve over years.\n\nThe stress-test's 1.74-day phase offset in Sec 5.5 is real, and the paper's claim that P_rot uncertainty alone can explain it is not convincing: over 254 days the accumulated phase error is about 0.1 days, not 1.74. That offset implies the spot changed between the transit epoch and the monitoring epoch, or that the model phase is off. But the stress-test goes too far in saying this undermines the single-spot assumption for the three transits, which span only 40 days. It is evidence of long-term spot evolution, not necessarily of intra-baseline changes. Still, it means the formal uncertainties on i_star, lambda, and spot latitude are probably underestimated.\n\nWho is this for? Anyone working on M-dwarf spot properties, obliquity measurements, or preparing transmission spectroscopy of TOI-3884b. It deserves a serious referee and publication after the spot-model caveats are made even more explicit. I would not use the quoted geometry as a tie-breaker for dynamical theory without independent confirmation from more flexible spot models and contemporaneous monitoring.","headline":"Solid new rotation period and off-pole spot evidence for TOI-3884, with geometry that is plausible but model-limited.","tokens_in":23123,"tokens_out":3699,"would_cite":true,"duration_ms":36928,"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":"Multi-epoch, multi-band transit photometry of TOI-3884 measures the host star's 11.043-day rotation and shows its persistent spot sits about 13 degrees from the stellar pole, yielding stellar inclination $i_\\star = 139.9$ deg and…","keywords":["Stellar rotation","Exoplanet systems","M dwarf stars","Planet hosting stars","Starspots","Multi-color photometry","Transit photometry","Polar orbit"],"falsifier":"Observe the next several transits continuously in at least two bands while simultaneously monitoring the star in r-band, and test whether the 11.043-day rotation model with the published spot latitude ($-76.8$ deg) and longitude ($262$ deg) predicts the exact phase and shape of each spot-crossing bump; a drift in the bump's crossing time beyond the rotation-period uncertainty, or a sharp low-flux minimum like the one seen in 2024–2025 that the single-spot model cannot reproduce, would refute the single fixed-circular-spot picture and the geometry derived from it.","tokens_in":21991,"feed_emoji":"🪐","tokens_out":15744,"duration_ms":140719,"temperature":0.7,"pith_summary":"TOI-3884 is a rare exoplanet system where every transit of a super-Neptune shows a spot-crossing bump, because the planet crosses a large spot near the pole of an M dwarf viewed almost pole-on. This paper uses three epochs of simultaneous four-band transit photometry (MuSCAT3/4) plus ground-based photometric monitoring to show that the spot is not exactly on the pole and that the star rotates with a period of 11.043 days. With the rotation period as a prior, the paper fits all twelve transit light curves with a single circular spot that rotates between transits, deriving a stellar inclination of $i_\\star = 139.9^{+1.2}_{-2.0}$ deg, a projected obliquity of $\\lambda = 41.0^{+3.7}_{-9.0}$ deg, and a spot latitude of about $-77$ deg. These values resolve the earlier disagreement between two previous studies, and they place TOI-3884b among the few known planets on near-polar orbits around low-mass stars. The paper argues that this well-characterized spot geometry matters for interpreting the planet's transmission spectrum, because the unocculted spot contaminates transit depths in a wavelength-dependent way.","feed_headline":"11.04-day spin resolves TOI-3884's disputed geometry","feed_subtitle":"Multi-band transit bumps place the host star's spot off the pole and its planet on a near-polar orbit.","key_machinery":"The analysis is carried by the `fleck` spot-crossing transit model, which places a circular, disk-like spot on the stellar surface with a contrast (temperature), radius, latitude, and longitude, and recomputes the transit light curve as the spot rotates with the stellar period. Its key partner is the rotation period $P_\\mathrm{rot} = 11.043^{+0.054}_{-0.053}$ days, measured two independent ways (a generalized Lomb-Scargle periodogram and a Gaussian-process quasi-periodic fit) from the Sinistro r-band monitoring light curve. The rotation period is the thread that ties the three epochs together: the spot's phase between transits is set by $P_\\mathrm{rot}$ and $i_\\star$, so the time-variable spot-crossing bumps constrain $i_\\star$ and the projected obliquity $\\lambda$, breaking the degeneracy between competing geometric solutions that could not be distinguished before. The spot contrast in each of the four bands is computed from BT-Settl model spectra, so the simultaneous $griz$ photometry also pins down the spot temperature ($\\Delta T = 200^{+11}_{-9}$ K).","core_discovery":"On the paper's own terms, the central discovery is that the persistent spot-crossing signal of TOI-3884b is produced by one large circular spot located at latitude $\\theta_\\mathrm{spot} = -76.8^{+2.5}_{-4.6}$ deg (about 13 degrees from the south pole), longitude $\\phi_\\mathrm{spot}=262^{+6}_{-14}$ deg, radius $0.425^{+0.018}_{-0.011}\\,R_\\star$, and temperature about 200 K cooler than the photosphere, whose longitude advances with the 11.043-day stellar rotation period. The time-variable shape of the spot-crossing bumps, combined with the rotation period measured from the ~5% r-band modulation, breaks the degeneracies that made earlier analyses disagree: it selects the solution with $i_\\star = 139.9^{+1.2}_{-2.0}$ deg and $\\lambda = 41.0^{+3.7}_{-9.0}$ deg, corresponding to a true obliquity of $\\Psi = 61.9^{+2.3}_{-5.6}$ deg. The fit also revises the impact parameter from near zero to $b=0.402^{+0.018}_{-0.019}$, using the spot geometry to break the usual degeneracy between $a/R_\\star$ and $b$. The same model explains why TESS light curves from an earlier epoch showed no rotational modulation: re-analysis of the TESS transits yields a spot at latitude about $-88$ deg, nearly exactly on the pole, implying the spot configuration has evolved over several years.","pith_inferences":["Beyond the paper: if the MuSCAT transit light curves are re-fit with two spots or a spot that evolves between epochs, the recovered $i_\\star$ and $\\lambda$ could shift; such a test would reveal whether the single-spot simplification is the source of the claimed reconciliation.","Beyond the paper: the model's failure to reproduce the sharp minimum in the monitoring light curve hints at an additional lower-latitude spot; a season of simultaneous transits and monitoring would reveal whether such a spot also alters the spot-crossing bumps and biases the obliquity.","Beyond the paper: if polar spots are common on slowly rotating field M dwarfs, as this paper hints, then transmission spectra of aligned-orbit planets around these stars are systematically contaminated by unocculted near-polar spots, and only contemporaneous multi-band monitoring can calibrate the effect.","Beyond the paper: the 1.74-day phase offset between the transit-epoch spot model and the later monitoring light curve could be spot migration rather than period uncertainty; tracking the offset over several seasons would test for longitude drift from differential rotation."],"forward_implications":["TOI-3884b is placed among the small population of Neptune-sized planets with close-to-polar orbits around M dwarfs, with true obliquity $\\Psi = 61.9^{+2.3}_{-5.6}$ deg, giving dynamical models a test case.","The spot temperature of $2952\\pm36$ K (with $\\Delta T = 200^{+11}_{-9}$ K) provides a precise point on the empirical spot-temperature relation for M dwarfs.","The measured chromaticity, with the transit depth about 7% deeper in $g$ than in $z$, is consistent with unocculted-spot contamination at a spot covering fraction of about 20%, so transmission spectra of TOI-3884b must be modeled with the spot configuration taken into account.","The derived $v\\sin i_\\star \\approx 0.89$ km/s indicates that previously reported values of 1.1 and 3.6 km/s are near the measurement limit and should be treated as upper limits.","The TESS-epoch spot appears nearly exactly polar, whereas the 2024–2025 spot is offset by about 13 degrees, indicating spot configuration changes on timescales of several years."],"supporting_citations":[{"why":"It provides the discovery characterization of TOI-3884b and the previous spot-crossing model whose parameter values (for example $i_\\star=47\\pm8.5$ deg and $\\lambda=151\\pm11$ deg) this paper reconciles.","marker":"A22"},{"why":"It supplies the stellar parameters used as Gaussian priors in the transit fits, including effective temperature, mass, radius, and orbital period, and the independent characterization that this paper compares with.","marker":"L23"},{"why":"It provides the fleck package, the spot-crossing transit model that places a circular spot on the stellar surface and is the core of the light-curve fitting.","marker":"B. Morris 2020"},{"why":"It provides the generalized Lomb-Scargle periodogram that first detects the 11.03-day rotation signal in the Sinistro monitoring light curve.","marker":"M. Zechmeister & M. Kürster 2009"},{"why":"It provides the celerite quasi-periodic Gaussian-process kernel used to measure the rotation period and its uncertainty independently.","marker":"D. Foreman-Mackey et al. 2017"},{"why":"It provides the BT-Settl model spectra used by the custom spocon module to compute wavelength-dependent spot contrast for the spot temperature.","marker":"F. Allard 2014"},{"why":"It supplies the geometric relation used to convert the fitted stellar inclination, orbital inclination, and projected obliquity into the true obliquity.","marker":"D. C. Fabrycky & J. N. Winn 2009"},{"why":"It provides the transit light source effect formula used to interpret the wavelength-dependent transit depth as unocculted-spot contamination.","marker":"B. V. Rackham et al. 2018"}],"fun_headline_variants":["11.04-day spin settles TOI-3884's disputed geometry","Off-pole spot plus rotation pin TOI-3884's true obliquity","Spot sits 13° off pole, unlocking TOI-3884's tilt","Multi-band bumps fix TOI-3884's spin and spot size","Persistent polar spot reveals TOI-3884's true spin-orbit geometry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that, over the roughly 40 days spanned by the three transits, the stellar surface is described by one circular spot whose radius, temperature, and latitude/longitude do not change and only its rotation phase advances; if the spot is non-circular, evolves, or coexists with other spots, the derived $i_\\star$, $\\lambda$, and spot properties could shift, a simplification the paper itself cautions can introduce systematic bias.","fun_headline_variants_meta":{"raw":{"variants":["11.04-day spin settles TOI-3884's disputed geometry","Off-pole spot plus rotation pin TOI-3884's true obliquity","Spot sits 13° off pole, unlocking TOI-3884's tilt","Multi-band bumps fix TOI-3884's spin and spot size","Persistent polar spot reveals TOI-3884's true spin-orbit geometry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001294,"raw_usage":{"total_tokens":5426,"prompt_tokens":1233,"completion_tokens":4193,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":849,"completion_tokens_details":{"reasoning_tokens":4091}},"tokens_in":849,"tokens_out":4193,"duration_ms":30727,"temperature":1.0,"reasoning_tokens":4091,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:56:59.005803+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the next several transits continuously in at least two bands while simultaneously monitoring the star in r-band, and test whether the 11.043-day rotation model with the published spot latitude ($-76.8$ deg) and longitude ($262$ deg) predicts the exact phase and shape of each spot-crossing bump; a drift in the bump's crossing time beyond the rotation-period uncertainty, or a sharp low-flux minimum like the one seen in 2024–2025 that the single-spot model cannot reproduce, would refute the single fixed-circular-spot picture and the geometry derived from it.","supporting_citations":[],"review_version":1}