{"id":"809496dc-1da8-4123-ae8b-688ec5a0122c","arxiv_id":"2608.11010","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"TESS photometry of three late O-type supergiants reveals evenly spaced frequencies interpreted as rotationally split gravity modes, with a few residual frequencies plausibly matching Rossby wave predictions.","lead":"This paper uses TESS brightness data to look for periodic pulsations and rotation signals in three O-type supergiant stars. It finds evenly spaced frequency patterns it links to rotation, and it tentatively identifies some leftover frequencies as Rossby waves, which may explain the flickering of massive stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on interpreting the fitted equal spacing Δf as rotational splitting, but the paper's own frequency tables show all detected frequencies as integer/half-integer multiples of Δf, which a co-rotating spot or harmonic series explains without invoking g-mode splitting or Rossby…","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the fitted frequency spacing Δf is assumed to be rotational splitting, and every derived rotation period and Rossby-mode identification depends on that assumption. I agree with that diagnosis and add two concrete stress points. First, the degeneracy with a corotating spot or harmonic series is not merely a theoretical worry; the paper's own Tables D.1 and D.2 show that essentially all frequencies in HD 192639 and HD 195592 are integer or half-integer multiples of Δf, and Sect. 6 explicitly lists corotating features and the superperiod effect as possible explanations. Second, the derived rotation period for HD 192639 is internally inconsistent with the adopted vsini=82 km/s: P_rot=16.1 d with R=19.8 R☉ gives v_rot≈62 km/s, which cannot satisfy v_rot sin i = 82 km/s for any inclination. This does not disprove the Rossby interpretation for the other two stars, but it removes one of three supporting cases and indicates that the conversion from Δf to f_rot is not robust. The paper is honest and tentative, and a conditional verdict is appropriate. The proposed model-comparison test would settle whether the equal spacings and r-mode matches survive when the spot/harmonic alternative is fitted with the same data and window function.","tokens_in":24999,"tokens_out":4613,"duration_ms":54247,"concrete_test":"For each star, perform a model comparison on the full set of detected frequencies: (A) the g-mode rotational-splitting plus Rossby-mode model with f_rot derived from Eq. (5); (B) a harmonic/co-rotating spot model with frequencies f_k = k·f_rot, including TESS orbital sidebands at ±0.074 d^-1. Use the actual TESS window function to generate synthetic periodograms under both models and compute the Bayesian evidence or a bootstrap false-alarm rate for the claimed equal spacings. For HD 192639, impose the physical constraint v_rot ≥ vsini (equivalently P_rot ≤ 2πR/vsini) in both fits; if the best-fit f_rot violates this bound, the rotational-splitting identification for that star is ruled out.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the interpretation of the fitted equal spacing Δf as rotational splitting via Eq. (5), which fixes f_rot and then enters every Rossby-mode identification through Eq. (4). This step is not secured against the alternatives the authors themselves list in Sect. 6: a long-lived co-rotating spot or a harmonic series produces integer and half-integer multiples of a fundamental frequency, and Tables D.1–D.2 indeed show that essentially all detected frequencies in HD 192639 and HD 195592 are integer or half-integer multiples of Δf. For high-ℓ r-modes, Eq. (4) reduces to σ ≈ mΩ = m·Δf once f_rot is derived from Δf, so the graphical matches in Figs. 4, 7, and 10 partly echo the same underlying pattern rather than provide independent confirmation. No false-alarm probability is given for the r-mode matches, and for HD 188001 the spacing Δf=0.0576 d^-1 is only comparable to 1.5/ΔT=0.055, so the 'constant spacing' is barely resolved. Moreover, the derived rotation periods are not all internally consistent: for HD 192639, R=19.8 R☉ and P_rot=16.1 d give v_rot≈62 km/s, which is below the adopted vsini=82 km/s; since sin i≤1 this is impossible at any inclination, yet Table 3 lists i≈90°. The r-mode interpretation for that star therefore rests on an internally inconsistent rotation estimate. If Δf is not rotational splitting, all three Rossby identifications lose their basis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyses TESS light curves of three late O-type supergiants (HD 188001, HD 192639, and HD 195592) using generalized Lomb-Scargle periodograms and weighted wavelet Z-transforms, with iterative pre-whitening and Monte Carlo frequency uncertainties. The authors identify a set of low-frequency peaks in each star, note an approximately constant frequency spacing Δf, interpret that spacing as rotational splitting via Eq. (5), derive rotation periods, classify most frequencies as l=1 or l=2 g-modes with the help of LPCODE/LP-PUL models, and then identify residual frequencies as Rossby-wave candidates using the dispersion relation Eq. (4). The central claim is that the detected evenly spaced frequency sets are rotationally split g-modes and that the residual pattern is consistent with Rossby waves, which would make the derived rotation periods and inclination angles physically meaningful.","tokens_in":25434,"tokens_out":4067,"duration_ms":36645,"significance":"If the central claim were robust, the paper would provide the first systematic evidence for Rossby waves in late O-type supergiants and would connect the low-frequency photometric variability of these stars to coherent rotation-driven oscillations rather than purely stochastic red noise. The authors use two complementary period-search methods, explicitly quote Monte Carlo frequency uncertainties, exploit multi-sector TESS observations spanning about five years, and compare against non-rotating LPCODE/LP-PUL pulsation models; they also repeatedly use cautious phrasing such as 'potential evidence' and 'tentative'. These are genuine strengths. However, the identification chain from Δf to f_rot to the r-mode assignments is not yet statistically or physically secured, and the derived parameters contain an internal inconsistency for HD 192639. The result is therefore plausible but not established at the level needed for its strong conclusions.","major_comments":[{"comment":"The rotation solution for HD 192639 is internally inconsistent: with R=19.8 R_sun and P_rot=16.1 d, the equatorial velocity is 2πR/P ≈62 km/s, which is less than the adopted vsin i=82 km/s. Since sin i ≤1, this is impossible at any inclination, yet Table 3 lists i≈90°. The r-mode identifications for this star, built on f_rot=0.0622 d^-1 from Eq. (5), therefore cannot be considered supported by the current analysis.","section":"§5.2, Table 3"},{"comment":"The load-bearing assumption that the fitted spacing Δf is rotational splitting is not tested against the alternatives listed in §6. Tables D.1 and D.2 show that essentially all detected frequencies are integer or half-integer multiples of Δf; a long-lived co-rotating spot or a harmonic series produces the same pattern without invoking g-mode splitting. For HD 188001, Δf=0.0576 d^-1 is only comparable to 1.5/ΔT=0.055 d^-1, so the 'constant spacing' is barely resolved. Because f_rot is derived from Δf via Eq. (5) and enters every r-mode identification through Eq. (4), this ambiguity is load-bearing for the central claim.","section":"§5.1, Eq. (5)"},{"comment":"The r-mode matches are not an independent confirmation: for high ℓ, Eq. (4) reduces to σ≈mΩ = mΔf once f_rot is set by Δf, so the graphical agreement in Figs. 4, 7, and 10 partly reproduces the same fitted spacing. No false-alarm probability is given for the matches, and several candidate frequencies are detected in only one sector (e.g., f2 for HD 188001 in sector 14, Table C.1). A quantitative test, such as a Monte Carlo false-alarm estimate over the searched (ℓ,m) grid, is needed before the 'potential evidence' can be distinguished from chance alignment.","section":"§6, Figs. 4, 7, 10"},{"comment":"The detection threshold was lowered to S/N≥4 because the recommended TESS threshold of S/N≥5 could exclude low-frequency RWs. While this is transparently stated, no estimate is given of how many spurious frequencies this admits, and the residual noise window is a broad 5 d^-1. Given that several r-mode candidates are single-sector detections and the rotation periods have large uncertainties, the significance of the r-mode identifications should be quantified under the adopted threshold.","section":"§3.1"}],"minor_comments":[{"comment":"The abstract names the third target as HD 159952, while the body and tables consistently use HD 195592.","section":"Abstract"},{"comment":"The term 'thermodes' is introduced without a definition; its relationship to toroidal modes should be clarified at first use.","section":"§4.1"},{"comment":"The note says 'The obtained period is in column 4', but in Table D.2 the period appears in column 5; the column numbering in the note should be corrected.","section":"Table D.2"},{"comment":"The text refers to 'f15' as identified in Table D.1, but Table D.1 does not use the f_i labels from Table C.2; aligning the row labels between the two tables would improve readability.","section":"§5.2 and Table D.1"}],"recommendation":"major_revision","confidential_remarks":"The analysis is potentially publishable, but the HD 192639 rotation-period inconsistency and the lack of a false-alarm test are serious enough that I would not recommend acceptance in the current form. If the authors can resolve the rotation-period contradiction and provide a statistical test distinguishing Δf as rotational splitting from harmonic or spot models, the paper could become suitable for A&A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. This is an honest, careful observational study of three late O supergiants with TESS light curves. What's new: first systematic frequency decomposition of these stars in TESS, with period lists from two complementary methods, phase-folded diagrams, and g-mode identification against LPCODE models. The authors are appropriately cautious, calling the Rossby evidence 'potential' and 'tentative' throughout, and they do list alternative explanations in the discussion. That intellectual honesty is real.\n\nThe soft spot is the load-bearing interpretation. The entire Rossby identification depends on treating the fitted equal spacing Δf as rotational splitting via Eq. (5). But the authors' own tables show that essentially all detected frequencies are integer or half-integer multiples of Δf. That is exactly the pattern produced by a long-lived co-rotating spot or a harmonic series, and the paper does not quantitatively exclude those alternatives. For HD 188001, the spacing is barely resolved, just above the 1.5/ΔT criterion. For HD 192639, there is an internal inconsistency: with R=19.8 R☉ and P_rot=16.1 d, the equatorial velocity is about 62 km/s, but the adopted vsini is 82 km/s, so the star cannot appear equator-on (i≈90°) as listed. That error undercuts the r-mode identification for that star. Throughout, there is no false-alarm probability for the r-mode matches, and the S/N threshold was lowered to 4 to keep candidates alive.\n\nNone of this makes the paper worthless. The frequency lists and g-mode matches are useful, and the red-noise connection is plausible. But the Rossby result is not independently established. If the spacing is not rotational splitting, all three Rossby identifications lose their basis.\n\nI would send this to peer review: the observational work deserves referee scrutiny, and the authors need to be pushed on the statistical front and on the spot/harmonic alternative. The paper is a solid candidate after major revision, not as is. For my own work, I would not cite the Rossby detection yet; I might cite the frequency lists. I would bring it to reading group as a case study in how a tentative detection can be vulnerable to interpretation.","headline":"Careful TESS analysis of three O supergiants, but the Rossby-wave interpretation rests on a frequency spacing that could equally be a spot or harmonic series.","tokens_in":25963,"tokens_out":2588,"would_cite":false,"duration_ms":23933,"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":"TESS light curves of three late O-type supergiants reveal rotationally split g-modes and candidate Rossby waves.","keywords":["Rossby waves","O-type supergiants","stellar pulsations","g-modes","rotational splitting","TESS photometry","asteroseismology","red noise"],"falsifier":"A direct measurement of the rotation period of HD 188001, HD 192639, or HD 195592—from time-resolved spectroscopy of line-profile variations, from a longer TESS baseline that resolves the frequency spacing well beyond $1.5/\\Delta T$, or from an independent asteroseismic analysis—that disagrees with the 14.5, 16.1, or 6.1 day values would falsify the rotational-splitting interpretation and, with it, the Rossby-mode identifications that depend on the derived rotation frequency.","tokens_in":24834,"feed_emoji":"🌊","tokens_out":11245,"duration_ms":91001,"temperature":0.7,"pith_summary":"This paper reports that the low-frequency variability seen in TESS light curves of three late O-type supergiants—HD 188001, HD 192639, and HD 195592—is not just stochastic red noise but contains coherent, rotation-related oscillations. The periodograms show a nearly constant frequency spacing in each star, which the authors interpret as rotational splitting of g-mode pulsations (mostly $\\ell=1$ or $\\ell=2$). Using the splitting formula $\\Delta f = m(1-C_\\ell) f_{\\rm rot}$, they derive rotation periods of roughly 14.5, 16.1, and 6.1 days, and they find that the residual independent frequencies are consistent with the dispersion relation of global Rossby waves, $\\sigma = -2m\\Omega/[\\ell(\\ell+1)] + m\\Omega$. If this interpretation survives longer baselines, it would mean that Rossby waves are excited in evolved massive stars and that part of the broadband 'red noise' seen in such stars is actually made of discrete rotationally driven modes.","feed_headline":"Rossby waves may explain low-frequency flicker in O-type supergiants","feed_subtitle":"TESS light curves reveal evenly spaced g-mode quintuplets and residual frequencies matching Rossby waves.","key_machinery":"The load-bearing identities are the rotational-splitting formula $\\Delta f = m(1-C_\\ell) f_{\\rm rot}$ (with $C_\\ell=0.166$ for two stars and $0.5$ for the third) and the inertial-frame Rossby dispersion relation $\\sigma = -2m\\Omega/[\\ell(\\ell+1)] + m\\Omega$, where Rossby waves are large-scale waves in a rotating fluid whose restoring force is the Coriolis force. The splitting formula converts a measured constant frequency spacing into a rotation frequency, and the dispersion relation is used to match residual observed frequencies to specific $(\\ell,m)$ Rossby modes through a graphical comparison. The period search combines a generalised Lomb-Scargle periodogram with a weighted wavelet Z-transform, followed by iterative pre-whitening to isolate significant frequencies, and the g-mode identifications are anchored by stellar evolution models and adiabatic pulsation calculations.","core_discovery":"The authors' central claim is that each of the three stars shows an approximately uniform frequency spacing in its TESS periodogram, which they identify as the rotational splitting of g-mode oscillations (mostly $\\ell=1$ or $\\ell=2$), and that after accounting for these modes plus harmonics and combination frequencies, the remaining independent frequencies match the inertial-frame Rossby-wave dispersion relation $\\sigma = -2m\\Omega/[\\ell(\\ell+1)] + m\\Omega$ for modes such as $(\\ell,m)=(3,2)$, $(3,3)$, and $(5,4)$. From the splitting they derive rotation periods of $14.5\\pm4.2$ days for HD 188001, $16.1\\pm2.6$ days for HD 192639, and $6.1\\pm1.0$ days for HD 195592, with inclination angles near $90^\\circ$ for the first two and about $20^\\circ$ for the third. They further argue that the co-existence of g-modes and Rossby candidates, plus a pattern of frequencies at integer and half-integer multiples of the splitting, points to rotational modulation and suggests that Rossby waves may contribute to the red-noise component observed in massive supergiants.","pith_inferences":["If Rossby waves are truly excited in these stars, they could act as a mechanism for angular momentum transport in the near-surface layers of massive stars, a process the paper mentions but does not model.","Applying the same analysis to a larger TESS sample of late O-type and B-type supergiants would test whether the constant-spacing pattern and Rossby candidates are a general property of evolved massive stars rather than a coincidence in three objects.","Spectroscopic follow-up targeting line-profile variability at the predicted Rossby frequencies could independently confirm the modes, because Rossby waves perturb surface temperature and pressure.","The HD 188001 result is the least secure because its spacing ($0.0576$ d$^{-1}$) is only marginally above the restrictive Rayleigh resolution ($0.055$ d$^{-1}$); an independent measurement of its rotation period would settle whether the quoted $14.5$ day period is real."],"forward_implications":["If the frequency spacing is rotational splitting, the derived rotation periods ($14.5$, $16.1$, and $6.1$ days) and inclination angles give direct photometric rotation constraints for these stars.","The residual frequencies matching the Rossby dispersion relation imply that the low-frequency variability of O-type supergiants is not purely stochastic; part of it is a set of discrete, rotationally driven modes, which could contribute to the observed red noise.","The g-modes that fit the adopted stellar models are predicted by late main-sequence evolutionary tracks, so these 'supergiants' may still be core hydrogen-burning stars whose supergiant appearance is set by their winds.","A longer TESS baseline, with frequency resolution better than $1.5/\\Delta T$, would test whether the equally spaced frequency sets and the Rossby candidates persist and sharpen."],"supporting_citations":[{"why":"Provides the Rossby-wave dispersion relation $\\sigma = -2m\\Omega/[\\ell(\\ell+1)]$ used to identify candidate modes.","marker":"Zaqarashvili et al. 2021"},{"why":"Supplies the rotational-splitting relation $\\Delta f = m(1-C_\\ell)f_{\\rm rot}$ that converts the observed spacing into rotation periods.","marker":"Dziembowski & Goode 1992"},{"why":"Establishes the photometric detectability of Rossby modes and the mode-visibility arguments used for candidate identifications.","marker":"Saio et al. 2018"},{"why":"Provides the adopted effective temperature, gravity, radius, and $v\\sin i$ for all three targets.","marker":"Gormaz-Matamala et al. 2022"},{"why":"Provides the projected rotation velocity for HD 192639 used in the inclination-angle estimate.","marker":"Holgado et al. 2022"},{"why":"Derives the correction to the Rossby dispersion relation and notes that high-order modes approach $\\sigma\\simeq m\\Omega$, supporting the observed $\\sigma=m\\Omega$ pattern.","marker":"Papaloizou & Pringle 1978"}],"fun_headline_variants":["Rossby waves may explain flicker in O-type supergiants","TESS finds g-modes and Rossby wave hints in O supergiants","Gravity modes and possible Rossby waves in late O stars","Low-frequency O-star variability: Rossby waves at play?","O-type supergiants show g-modes and Rossby wave candidates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the nearly constant frequency spacing seen in each star's periodogram is rotational splitting of g-modes, rather than a corotating spot, binary motion, or an instrumental artifact; if that premise fails, the rotation periods and every Rossby-mode identification built on them collapse.","fun_headline_variants_meta":{"raw":{"variants":["Rossby waves may explain flicker in O-type supergiants","TESS finds g-modes and Rossby wave hints in O supergiants","Gravity modes and possible Rossby waves in late O stars","Low-frequency O-star variability: Rossby waves at play?","O-type supergiants show g-modes and Rossby wave candidates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000465,"raw_usage":{"total_tokens":2411,"prompt_tokens":1127,"completion_tokens":1284,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":743,"completion_tokens_details":{"reasoning_tokens":1191}},"tokens_in":743,"tokens_out":1284,"duration_ms":12043,"temperature":1.0,"reasoning_tokens":1191,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:06:58.301148+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the rotation period of HD 188001, HD 192639, or HD 195592—from time-resolved spectroscopy of line-profile variations, from a longer TESS baseline that resolves the frequency spacing well beyond $1.5/\\Delta T$, or from an independent asteroseismic analysis—that disagrees with the 14.5, 16.1, or 6.1 day values would falsify the rotational-splitting interpretation and, with it, the Rossby-mode identifications that depend on the derived rotation frequency.","supporting_citations":[],"review_version":1}