{"id":"cf6400fb-0e86-4d4a-8dac-2d0789fe6a98","arxiv_id":"2412.08557","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"CHEOPS transit photometry of WASP-33b confirms the previously detected nodal precession and yields a Love number consistent with stellar models, but fails to measure the occultation depth.","lead":"Using four CHEOPS transits, the authors show that the transit path of hot Jupiter WASP-33b is shifting over time, consistent with the known nodal precession of its orbit. They also derive a stellar Love number that agrees with theoretical models, but they cannot extract a reliable occultation depth or albedo.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The photometric 'confirmation' of precession is prior-driven: CHEOPS cannot distinguish +b from -b, so the plotted agreement with the Watanabe curve assumes the sign that the model predicts.","rationale":"The paper's central claim, as stated in the abstract and Section 8, is that CHEOPS photometry confirms the nodal precession of WASP-33b and that this demonstrates photometry-alone detectability. What has to be true for that claim is that the CHEOPS b measurements discriminate between the precession model and alternatives such as constant b or the opposite sign of b. The manuscript itself identifies the sign degeneracy in note 7, and the reader correctly isolates this as the weakest assumption. My reading of Table 5 and Fig. 8 supports that: the reported b is an unsigned magnitude, the plotted sign is taken from the Watanabe et al. (2022) model, and the model curve is not refit to the new points. If the sign were free, the CHEOPS points would be equally consistent with positive b, placing them off the model curve, so the photometry alone does not confirm the model. I credit the authors for disclosing this limitation, for the honest injection-recovery failure on the occultations, and for the careful re-derivation of stellar parameters and the independent PIPE reduction. Those strengths make the paper a solid consistency study, but the 'confirmation' and 'photometry alone' language should be softened or the sign-dependence explicitly carried into the conclusions. The reader's CONDITIONAL verdict is therefore appropriate and does not need to change.","tokens_in":30868,"tokens_out":9101,"duration_ms":100396,"concrete_test":"Refit the four CHEOPS transits (individually or as the two pairs used in Section 7.3) with the TLCM gravity-darkened model in two runs, one constraining b > 0 and one constraining b < 0, keeping all other settings identical to case (iii) of Table 5. Compare the two branches by Delta BIC and by the location of the resulting b(t) values relative to the Watanabe et al. (2022) curve; if the sign branches fit equally well (Delta BIC < 3) and only the negative branch lies on the model, the photometric confirmation is prior-driven rather than established by the CHEOPS data alone.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 7.3 (note 7) states that transit light-curve fits cannot distinguish +b from -b, and Fig. 8 plots the CHEOPS b values as negative because that sign is expected from the Watanabe et al. (2022) Doppler-tomography model. The fiducial CHEOPS fit (case iii of Table 5) reports b = 0.064+0.075-0.058 as a positive value, i.e. an unsigned magnitude from the photometric fit. If the true sign at the CHEOPS epochs were positive, the same photometric fits would place the points on the opposite side of b=0, where the Watanabe curve is not located. The agreement of the CHEOPS points with the precession curve is therefore not an independent photometric test of the model; it is a consistency check between |b| and a curve whose sign was supplied by the model. Section 7.3 also states that the authors do not refit the precession model to the CHEOPS data, only overplot the Watanabe curve. The abstract's 'confirm' and Section 8's claim that the work demonstrates that precession 'can be detected from photometry alone' are stronger than what sign-free photometry can establish. The authors disclose the degeneracy clearly, which is to their credit, but the headline claim should be conditioned on the DT sign prior.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes four CHEOPS transits and four occultations of WASP-33b, a hot Jupiter orbiting a rapidly rotating delta Scuti star. The transit light curves are modeled with a gravity-darkening prescription, yielding measurements of the transit impact parameter and stellar inclination; the main fiducial fit fixes the sky-projected obliquity lambda to the Doppler-tomography value of Watanabe et al. (2022). The authors report that the CHEOPS impact parameters are consistent with the nodal precession model of Watanabe et al. (2022), combine that model's J2 with a newly derived stellar rotation period to obtain the fluid Love number k2 (0.0099 +/- 0.0012 or 0.0074 +/- 0.0011), and find agreement with theoretical stellar models. They also re-derive stellar parameters from spectra and Gaia data, characterize pulsations from TESS, and show through extensive tests that the occultation depth cannot be reliably extracted from the CHEOPS data.","tokens_in":31076,"tokens_out":4531,"duration_ms":47620,"significance":"If the nodal precession confirmation were robust, the paper would demonstrate a purely photometric detection channel for spin-orbit precession, with direct relevance to PLATO. The study is careful in several respects: the pulsation treatment is systematic and cross-checked with multiple frequency lists; the occultation analysis includes injection-recovery tests and honestly concludes that no reliable depth can be measured; the stellar parameters are re-determined with modern data; and the sign degeneracy of the impact parameter is explicitly acknowledged in Section 7.3. However, the central claim of photometrically confirming nodal precession is weakened by the fact that the CHEOPS b values are sign-degenerate, so the plotted agreement with the Watanabe et al. (2022) curve assumes the sign supplied by that model. The paper is a useful consistency check and a source of improved stellar parameters, but the headline claim is stronger than the data support.","major_comments":[{"comment":"The transit light-curve fits cannot distinguish +b from -b, as stated in note 7 of Section 7.3, and the CHEOPS points are plotted with negative b because that sign is expected from the Doppler-tomography model of Watanabe et al. (2022). The agreement of the CHEOPS impact parameters with the precession curve is therefore not an independent photometric confirmation of nodal precession; it is a consistency check between |b| and a curve whose sign was supplied by the model. The abstract's 'confirm' and Section 8's claim that precession 'can be detected from photometry alone' should be tempered to reflect this prior dependence, or the analysis should be extended to show that the photometry can discriminate between the two signs when combined with the model.","section":"Section 7.3 and Fig. 8"},{"comment":"The fiducial CHEOPS measurement b = 0.064+0.075-0.058 is consistent with zero and with a broad range of precession phases, and the authors state that they do not refit the precession model to the CHEOPS data. As presented, the photometric measurement adds little quantitative constraint on the precession parameters, and the word 'confirm' overstates the evidential weight. I recommend adding a quantitative comparison, for example a likelihood ratio between the Watanabe et al. (2022) precession model and a constant-b model evaluated with the CHEOPS points, or an explicit statement of the posterior probability that the CHEOPS b values lie on the model curve given the photometric uncertainties.","section":"Table 5, case (iii), and Section 7.3"},{"comment":"The headline obliquity psi = 111.3+0.2-0.7 deg is quoted from the fit in which lambda is fixed to the Watanabe et al. (2022) Doppler-tomography value, so the quoted uncertainty excludes both the uncertainty in that fixed lambda and the systematic choice among the four degenerate photometric solutions shown in Fig. 6. The free gravity-darkening fit (case ii of Table 5) yields lambda = -148+22-17 deg and psi = 139+12-15 deg with substantially larger errors, and the 'correct' cluster in Fig. 6 is evidently identified using the DT value. The precision of the headline obliquity is therefore largely inherited from the DT prior rather than from the CHEOPS photometry; this should be stated explicitly and the error propagation should include the lambda uncertainty.","section":"Section 7.1 and Eq. 3"}],"minor_comments":[{"comment":"The sentence beginning 'Tthe best ephemeris for planning future observations' contains a typo; 'Tthe' should be 'The'.","section":"Table 5, note"},{"comment":"The word 'temperatute' in the first sentence should be 'temperature'.","section":"Section 3.1"},{"comment":"Please define what is meant by the 'correct' solution in the caption, since the identification of this cluster is used to select the photometric solution in the free gravity-darkening fit.","section":"Fig. 6 caption"},{"comment":"The statement that the Claret (2023) values 'must be multiplied by two' would be easier to verify if the definition of k2 used by Claret were given explicitly, for instance in a short equation or by citing the specific formula.","section":"Section 7.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is transparent about the b sign degeneracy and the occultation difficulties, which is commendable. The main issue is that the abstract and Section 8 claim a photometric confirmation of nodal precession despite the fact that the plotted agreement relies on the sign of b being taken from the Doppler-tomography model. A revised version that softens the central claim and adds a quantitative assessment of the evidence would be appropriate for the journal. The stellar parameter update and the Love number comparison remain valuable even if the precession 'confirmation' is reframed as a consistency check."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, transparent CHEOPS analysis, but the abstract sells it as a stronger confirmation of nodal precession than the data actually support. The four transits show the expected gravity-darkened asymmetry and the fitted |b| is consistent with the Watanabe et al. (2022) precession curve, but photometry cannot distinguish a transit north of the stellar equator from one south of it. The plotted agreement therefore assumes the sign that the model predicts. The authors disclose this in note 7 of Section 7.3, which is to their credit, but the phrase 'can be detected from photometry alone' in Section 8 is not supported without that prior.\n\nWhat is genuinely new and useful: the CHEOPS light curves, the redetermined stellar parameters (Teff, R*, M*, age) from spectra plus Gaia, a cleaned 108-frequency TESS pulsation list, an updated ephemeris from four new transit timings, and a k2 estimate that combines the Watanabe J2 with a new rotation period. The occultation section is a model of honesty: over 200 fits, depth values that scatter wildly, and injection–recovery tests that fail to recover a 500 ppm signal. That is exactly how to report a null result.\n\nSoft spots, in order of importance. First, the confirmation claim: the sign degeneracy means the photometry alone would look equally consistent with the mirror image of the precession curve, so the test is really a consistency check on |b| once the sign is taken from the DT prior. Second, the headline obliquity psi = 111.3+0.2−0.7 comes from a fit with lambda fixed to the Watanabe value; the free fit gives psi = 139+12−15, which is 2–3 sigma off. The paper flags this, but the abstract’s precision is attached to a prior-anchored fit. Third, the Love number is not a fresh measurement of J2; it reuses the precession model’s J2, so the agreement with theory is partly by construction. Fourth, beta is fixed at 0.23 with a known uncertainty, which adds a modest systematic.\n\nNone of these are fatal. The central consistency claim holds up and the limitations are mostly in the text. The paper will be useful to anyone working on gravity-darkened transits, hot-Jupiter obliquities, or PLATO forecasts. It deserves peer review and probably acceptance after the authors temper the abstract and the closing claim. I’d cite it for the pulsation list and the honest occultation failure analysis.","headline":"A transparent CHEOPS analysis that makes a useful confirmation but overstates it: photometry alone cannot fix the sign of b, so the 'photometric detection' of precession leans on the Doppler-tomography prior.","tokens_in":32230,"tokens_out":3495,"would_cite":true,"duration_ms":36639,"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":"CHEOPS photometry confirms WASP-33b's orbit is precessing, and the derived stellar Love number matches theory.","keywords":["nodal precession","WASP-33b","transit impact parameter","gravity darkening","stellar Love number","CHEOPS photometry","hot Jupiter","stellar pulsations"],"falsifier":"A decisive test would be a Doppler-tomography observation of a WASP-33 transit taken at an epoch where the precession model predicts a significantly negative $b$; Doppler tomography measures both the magnitude and the sign of $b$, so if the sign turned out to be positive, the CHEOPS points would not lie on the published precession curve and the photometric confirmation would fail.","tokens_in":30600,"feed_emoji":"🪐","tokens_out":10839,"duration_ms":105275,"temperature":0.7,"pith_summary":"The paper aims to show that CHEOPS transit photometry can independently confirm the nodal precession of the hot Jupiter WASP-33b, and that the precession can be turned into a measurement of the host star's second-order fluid Love number. The four new transit light curves give an impact parameter that falls on the secular precession curve established by earlier Doppler-tomography work, and combining the precession model's quadrupole moment with a new stellar rotation period yields $k_{2,\\star}$ values in agreement with theoretical stellar models. This matters because a Love number measured from precession is a direct probe of how mass is distributed inside the star, and because detecting precession photometrically opens the way for future wide-field space photometry to find more such systems. The same fits measure a true obliquity of $111.3^{+0.2}_{-0.7}$ degrees, consistent with previous measurements, and the paper redetermines stellar parameters and updates the ephemeris. The occultation observations, by contrast, do not yield a robust eclipse depth because the star's pulsations dominate the signal.","feed_headline":"CHEOPS confirms precession in WASP-33 system","feed_subtitle":"New transit fits peg WASP-33's fluid Love number at 0.0099, matching stellar models.","key_machinery":"The load-bearing mechanism is the nodal precession of WASP-33b: the star's rapid rotation makes it oblate, and the resulting gravitational quadrupole moment $J_2$ torques the tilted planetary orbit, producing a roughly 700-year sinusoidal drift of the transit impact parameter $b$. The paper fits the CHEOPS transits with a gravity-darkened transit model, which is necessary because the oblate, gravity-darkened stellar disc makes the light curve asymmetric; from those fits it extracts $b$ and the stellar inclination $i_\\star$. It then converts the precession model's $J_2$ into the second-order fluid Love number $k_{2,\\star}$, a dimensionless measure of how centrally concentrated the stellar mass is, through $k_{2,\\star} = 3J_2(\\Omega_{\\rm crit}/\\Omega_{\\rm rot})^2$, using the rotation period implied by $i_\\star$.","core_discovery":"The paper's central claim is that space-based transit photometry alone confirms the nodal precession of the hot Jupiter WASP-33b. Fitting four CHEOPS transits with a gravity-darkened model gives a transit impact parameter $b$ that follows the secular curve predicted by the published nodal-precession model, whose roughly 700-year cycle carries $b$ from positive to negative values over the observational window. The precession model supplies the stellar gravitational quadrupole moment $J_2 = (1.36^{+0.15}_{-0.12})\\times 10^{-4}$; combining this with the rotation period derived from the fitted stellar inclination gives a second-order fluid Love number $k_{2,\\star}=0.0099\\pm0.0012$ (or $0.0074\\pm0.0011$ using an alternative inclination), in agreement with theoretical stellar models. The same fits measure the true stellar obliquity $\\psi = 111.3^{+0.2}_{-0.7}$ degrees, consistent with previous Doppler-tomography values. The paper also redetermines the stellar radius and mass, updates the ephemeris, and finds that the occultation depth cannot be robustly extracted because the stellar pulsations dominate the eclipse signal.","pith_inferences":["Beyond the paper: a simultaneous Doppler-tomography and photometry campaign would turn the sign ambiguity of $b$ into a decisive test, since transit photometry alone cannot distinguish northern from southern transits.","Beyond the paper: the failed occultation recovery suggests that eclipse photometry of pulsating hot Jupiters needs long, same-band out-of-eclipse baselines or simultaneous spectroscopy; pulsation frequencies borrowed from another instrument at another epoch do not solve the problem.","Beyond the paper: if precession-based Love numbers like this one can be measured for several hot Jupiters around rapid rotators, the ensemble would map how stellar internal structure depends on mass and age, a test stellar models have not yet faced at this precision."],"forward_implications":["Nodal precession in a hot-Jupiter system can be confirmed from photometry alone, without the spectroscopic Doppler-tomography measurements that previously detected it.","Long-baseline space photometry of bright stars should be able to detect and track such precession in real time, extending the sample beyond the few systems known today.","The measured Love number $k_{2,\\star}=0.0099\\pm0.0012$ (or $0.0074\\pm0.0011$) agrees with theoretical stellar models, supporting the models' internal-structure predictions for this type of star.","The updated ephemeris gives transit-time predictions with 1-sigma uncertainty below 90 seconds for the rest of this decade and below 150 seconds through the 2030s.","The transit fits with the obliquity fixed to the Doppler-tomography value constrain the stellar inclination to $i_\\star = 100.5\\pm8.4$ degrees."],"supporting_citations":[{"why":"Supplies the nodal-precession model, the precession rate, and the $J_2$ value that the CHEOPS $b$ measurements are compared against.","marker":"Watanabe et al. (2022)"},{"why":"Discovery paper establishing WASP-33b, the rapid rotation, and the original Doppler-tomography measurements of obliquity and impact parameter.","marker":"Collier Cameron et al. (2010)"},{"why":"Provides the formula connecting $J_2$, rotation, and the second-order fluid Love number used to derive $k_2$.","marker":"Ragozzine & Wolf (2009)"},{"why":"Gives the theoretical stellar-model Love numbers against which the measured $k_2$ is judged.","marker":"Claret (2023)"},{"why":"Establishes the gravity-darkened transit light-curve asymmetry that lets the fits extract obliquity and inclination.","marker":"Barnes (2009)"},{"why":"Provides the light-curve modeller used for the gravity-darkened transit fits and wavelet noise treatment.","marker":"Csizmadia (2020)"},{"why":"Supplies one of the TESS pulsation frequency lists used for pulsation removal and the occultation injection-recovery tests.","marker":"von Essen et al. (2020)"},{"why":"A prior gravity-darkened TESS fit whose treatment of the gravity-darkening coefficient is followed, and a photometric $b$ point at the same epoch.","marker":"Dholakia et al. (2022)"},{"why":"Provides TESS pulsation frequencies and a gravity-darkened fit giving another photometric $b$ measurement for comparison.","marker":"Kálmán et al. (2022)"},{"why":"Supplies the reference ephemeris against which the new transit timings are compared.","marker":"Ivshina & Winn (2022)"}],"fun_headline_variants":["CHEOPS nails WASP-33b's precession and stellar Love number","WASP-33b's precession pinned down by CHEOPS transit fits","CHEOPS confirms WASP-33b precession, Love number matches theory","WASP-33b's nodal precession confirmed: Love number from CHEOPS","CHEOPS pinpoints WASP-33b precession, stellar Love number"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The confirmation assumes the sign of the transit impact parameter is negative, taken from the earlier precession model, because transit photometry cannot distinguish a transit across the northern half of the stellar disc from one across the southern half.","fun_headline_variants_meta":{"raw":{"variants":["CHEOPS nails WASP-33b's precession and stellar Love number","WASP-33b's precession pinned down by CHEOPS transit fits","CHEOPS confirms WASP-33b precession, Love number matches theory","WASP-33b's nodal precession confirmed: Love number from CHEOPS","CHEOPS pinpoints WASP-33b precession, stellar Love number"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000543,"raw_usage":{"total_tokens":2646,"prompt_tokens":1037,"completion_tokens":1609,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":1502}},"tokens_in":653,"tokens_out":1609,"duration_ms":13297,"temperature":1.0,"reasoning_tokens":1502,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:45:58.388720+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a Doppler-tomography observation of a WASP-33 transit taken at an epoch where the precession model predicts a significantly negative $b$; Doppler tomography measures both the magnitude and the sign of $b$, so if the sign turned out to be positive, the CHEOPS points would not lie on the published precession curve and the photometric confirmation would fail.","supporting_citations":[{"cited_title":"2022, , 512, 4404","cited_arxiv_id":null,"evidence_quote":"Supplies the nodal-precession model, the precession rate, and the $J_2$ value that the CHEOPS $b$ measurements are compared against."},{"cited_title":"& Wolf , A","cited_arxiv_id":null,"evidence_quote":"Provides the formula connecting $J_2$, rotation, and the second-order fluid Love number used to derive $k_2$."},{"cited_title":"2023, , 674, A67","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical stellar-model Love numbers against which the measured $k_2$ is judged."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the gravity-darkened transit light-curve asymmetry that lets the fits extract obliquity and inclination."},{"cited_title":"2020, , 496, 4442","cited_arxiv_id":null,"evidence_quote":"Provides the light-curve modeller used for the gravity-darkened transit fits and wavelet noise treatment."},{"cited_title":"C., et al","cited_arxiv_id":null,"evidence_quote":"Supplies one of the TESS pulsation frequency lists used for pulsation removal and the occultation injection-recovery tests."},{"cited_title":"2022, , 925, 185","cited_arxiv_id":null,"evidence_quote":"A prior gravity-darkened TESS fit whose treatment of the gravity-darkening coefficient is followed, and a photometric $b$ point at the same epoch."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the reference ephemeris against which the new transit timings are compared."}],"review_version":1}