{"id":"6187f34b-bddb-478f-87c3-c3a4bd87ec53","arxiv_id":"2501.02272","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"TOI-6038 A b is a newly discovered dense sub-Saturn planet at the boundary between the Neptunian ridge and savanna, with mass 78.5 Earth masses and radius 6.41 Earth radii.","lead":"We report the discovery and characterization of TOI-6038 A b, a dense sub-Saturn planet with a mass of 78.5 Earth masses and a radius of 6.41 Earth radii, orbiting a bright F-type star in a wide binary system. It sits at the boundary between the Neptunian ridge and savanna, making it a useful probe of planetary migration and atmospheric evolution.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass measurement rests on 29 of 43 RVs, with 14 seeing-excluded points never shown or tested; if those shift K, the headline mass is not robust.","rationale":"The paper is a standard discovery; the measured radius rests on TESS photometry with speckle and RV confirmation, and the density is robust to the low-mass stellar solution. The weakest link in the central mass claim is the undocumented RV cut. The reader's weakest assumption (internal structure) is model-dependent but does not threaten the headline mass/radius. The paper also lacks a fully reproducible data/code package, but that is a transparency issue, not a correctness threat. A single concrete check, refitting with the excluded RVs, would settle whether the mass measurement is biased. Until that is shown, conditional acceptance is appropriate.","tokens_in":27575,"tokens_out":8995,"duration_ms":93366,"concrete_test":"Obtain the 14 excluded RVs with uncertainties from the authors; refit the circular single-planet model with EXOFASTv2 on all 43 points (or on the 29 plus a jitter model) and compare K and MP. If K shifts by more than ~2.8 m/s (1 sigma) or the best-fit period changes, the mass claim is not robust; if K and MP remain within 1 sigma, the exclusion is benign and the central claim stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.3 states that fourteen additional PARAS-2 spectra were excluded from the RV analysis 'due to poor seeing conditions', but the excluded points are not listed, their phase coverage is not shown, and no fit including them is reported. The reported semi-amplitude K=23.5±2.8 m/s is the sole basis for the headline mass MP=78.5±9.9 M_Earth. A data-quality cut is normally benign if it is random in orbital phase, but the paper does not demonstrate this. If the excluded points preferentially sample certain phases (e.g., the pre-recoating epoch with larger errors, or poor-seeing nights that avoid or track the RV maximum), K could be biased by more than its quoted uncertainty. This is a selection-bias risk in the central mass claim. The internal-structure model dependence raised by the reader affects the secondary core-mass/envelope-fraction interpretation, not the measured mass and radius, so the RV cut is the more load-bearing gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"TOI-6038 A b is presented as a newly discovered dense sub-Saturn orbiting a bright, metal-rich late F star in a wide binary system. The analysis combines nine TESS transits from Sectors 18 and 58 with 29 PARAS-2 radial velocities, and uses EXOFASTv2 to jointly model the photometry, RVs, SED, and MIST isochrones. The adopted circular-orbit, high-stellar-mass solution yields P = 5.8267311 days, R_P = 6.41 R_Earth, M_P = 78.5 M_Earth, and a bulk density of 1.62 g/cm3, placing the planet at the boundary between the Neptunian ridge and savanna. Internal-structure modeling gives a core mass of about 58 M_Earth and a 26% H/He envelope mass fraction. The paper discusses tidal circularization, a possible von Zeipel-Lidov-Kozai migration pathway driven by the binary companion, and future atmospheric and orbital follow-up opportunities.","tokens_in":27781,"tokens_out":19637,"duration_ms":186322,"significance":"If the derived parameters hold, the paper delivers a precise mass and radius for a sub-Saturn in a sparsely populated region of parameter space, with a density consistent with the high-density ridge population and a wide binary companion that is relevant for testing high-eccentricity migration scenarios. The analysis has several strengths: contamination is quantified with TESS-cont; the planetary origin of the RV signal is supported by in-phase TESS photometry and null activity-indicator periodograms; speckle imaging and Gaia RUWE rule out close stellar companions; and both circular/eccentric and high/low stellar-mass solutions are compared. The two main weaknesses are the undocumented exclusion of 14 radial velocities from the fit and a numerical inconsistency in the quoted Kozai timescale; both are fixable and do not, by themselves, invalidate the central mass and radius measurement.","major_comments":[{"comment":"The manuscript states that fourteen additional PARAS-2 spectra were excluded from the RV analysis 'due to poor seeing conditions', but these points are never listed, their orbital phases are not shown, and no fit including them is reported. Since the reported semi-amplitude K = 23.5 ± 2.8 m/s is the sole basis for the headline mass M_P = 78.5 ± 9.9 M_Earth, a phase-dependent selection could bias K by more than its formal uncertainty. Please provide the excluded epochs and RVs (or at minimum their phases, uncertainties, and a seeing/SNR criterion), and either include them in a fit with appropriate weights or demonstrate explicitly that their exclusion is phase-independent and does not shift K.","section":"§2.3, Table 4"},{"comment":"The quoted von Zeipel-Lidov-Kozai timescale is numerically inconsistent with the stated parameters. Using the paper's formula t_vZLK = (a_comp/a_1,0)^3 (m_primary/m_comp) t_P,0 with a_comp = 3217 AU, a_1,0 = 5 AU, m_primary/m_comp ≈ 1.5, and t_P = 5.8267 days gives t_vZLK ≈ 6 × 10^6 yr, not 4 × 10^9 yr; the quoted 4 Gyr would correspond to a_1,0 ≈ 0.5 AU rather than 5 AU. Please correct the numerical value or the assumed initial semi-major axis and re-evaluate the sentence claiming that the timescale is 'comparable to the age of the system'. The qualitative conclusion based on epsilon_GR ≈ 3 may be unaffected, but the stated timescale is wrong as written.","section":"§4.2"}],"minor_comments":[{"comment":"The low-mass stellar solution is reported only in the appendix and is not discussed in the main text. Since it gives M_P ≈ 72.5 M_Earth versus 78.5 M_Earth for the adopted high-mass solution, a sentence in §3.6 stating that the low-mass solution is within the quoted 68% interval and leaves the density and ridge/savanna classification unchanged would help readers judge the robustness of the result.","section":"§3.6, Tables 2 and 5"},{"comment":"The internal-structure results (M_core ≈ 58 M_Earth, f_env ≈ 0.26) are computed with one core mass-radius relation (Otegi et al. 2020) and one envelope model (Chen & Rogers 2016). Please state explicitly that these values are model-dependent, and if possible quote the spread from an alternative envelope model (for example Lopez & Fortney 2014) rather than only the formal fit uncertainty.","section":"§4.1, Table 3"},{"comment":"The TESS aperture contains a 10.4% flux contribution from TOI-6038 B (Star 1). Although the PDCSAP light curves are corrected for crowding and the paper reasonably does not include an extra dilution factor, a fit with a free dilution parameter constrained by the TESS-cont estimate would directly test the sensitivity of R_P and hence of the density to residual contamination errors.","section":"§2.1"},{"comment":"There is a minor wording error: 'not additional relevant peaks pop up' should read 'no additional relevant peaks pop up.'","section":"§3.5"},{"comment":"The sentence 'TOI-6038 is only the fifth known multiple star system...' should refer to TOI-6038 A, since TOI-6038 names the whole binary system.","section":"§3.2"},{"comment":"The quantity labeled 'envelope radius' appears to be the envelope thickness R_P - R_core rather than a radius in the usual sense; consider renaming it to 'envelope thickness' for clarity.","section":"§4.1, Table 3"}],"recommendation":"major_revision","confidential_remarks":"The central characterization appears sound and the paper is a useful addition to the sparse sample of well-measured sub-Saturns. The main risk to the headline mass is the undocumented RV cut; the Kozai timescale issue is a clear numerical error that should be corrected. Both concerns are addressable in a revision, so I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid discovery paper in the normal TESS + radial-velocity mold. It reports a new sub-Saturn, TOI-6038 A b, with a mass of 78.5 ± 9.9 M_Earth and radius 6.41 ± 0.20 R_Earth from 29 PARAS-2 RVs and nine TESS transits. The planet sits at the boundary between the Neptunian ridge and savanna, with a density (1.62 g/cm3) that fits the dense ridge population. There's also a wide binary companion at ~3200 AU, which makes it a reasonable target for future obliquity and atmospheric follow-up.\n\nWhat's new is the object itself; the methods are entirely standard. That's not a criticism. The paper does the usual checks well: TESS contamination analysis with TESS-cont, speckle imaging, CCF activity indicators, and a BIC/AIC comparison that favors a circular orbit. The authors are honest about the migration story: their own vZLK estimate suggests the wide companion is probably insufficient to drive high-eccentricity migration, and they leave disk migration and undetected companions open. That's appropriate framing, not overclaiming.\n\nThe main soft spot is the excluded radial velocities. The paper says fourteen additional spectra were discarded for poor seeing, but they are never listed, their phase coverage is never shown, and no fit including them is reported. The semi-amplitude K of 23.5 ± 2.8 m/s is the sole basis for the headline mass. A poor-seeing cut is usually benign if it's random in phase, but the paper doesn't demonstrate that. If those points preferentially sample part of the orbit, K could shift by more than the quoted uncertainty. This is a legitimate referee question, not a demonstrated flaw. The authors need to show the excluded points and run a sensitivity fit.\n\nThe internal structure interpretation — a ~58 M_Earth core and 26% envelope — depends on the chosen two-layer models and should be treated as illustrative, not robust. The measured mass and radius are the secure result. I also wished the paper shipped code, but the RVs are in the table and the light curves are public, so that's minor.\n\nThis paper deserves peer review. It's a useful data point for the ridge/savanna population and a good target for follow-up. I'd send it to review with a request to address the RV exclusion. I'd bring it to a reading group focused on exoplanet demographics; otherwise it's not a must-read.","headline":"A solid, standard sub-Saturn characterization; the unexplained exclusion of 14 RVs is the one gap that needs a referee's attention.","tokens_in":28396,"tokens_out":3616,"would_cite":true,"duration_ms":32973,"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":"The paper establishes TOI-6038 A b as a dense, core-dominated sub-Saturn at the boundary between the Neptunian ridge and savanna, and shows the wide binary companion alone probably cannot explain its orbit.","keywords":["sub-Saturn","Neptunian ridge","Neptunian savanna","radial velocity","TESS photometry","internal structure","high-eccentricity tidal migration","wide binary"],"falsifier":"Measure the Rossiter-McLaughlin effect in the stellar lines during transit: a large spin-orbit misalignment would contradict the paper's conclusion that the wide companion cannot drive high-eccentricity migration and would reopen that channel, while a small misalignment would favor disk-driven, coplanar migration.","tokens_in":27414,"feed_emoji":"🪐","tokens_out":8106,"duration_ms":69769,"temperature":0.7,"pith_summary":"The paper reports the discovery and characterization of TOI-6038 A b, a sub-Saturn exoplanet orbiting a bright, metal-rich late F-type star that is itself the primary of a wide binary system. Combining TESS transit photometry with PARAS-2 radial velocities, the authors establish a mass of $78.5^{+9.5}_{-9.9}\\,M_\\oplus$, a radius of $6.41^{+0.20}_{-0.16}\\,R_\\oplus$, and a bulk density of $1.62^{+0.23}_{-0.24}\\,\\mathrm{g\\,cm^{-3}}$, placing the planet at the border between the Neptunian ridge and the Neptunian savanna. The high density matches the ridge population, which has been proposed to arrive at short periods through high-eccentricity tidal migration, yet a first-order dynamical analysis says the wide stellar companion cannot by itself drive such migration. This makes the system a test bed for distinguishing dynamically driven migration, disk-driven migration, and in-situ formation for planets in this transition regime.","feed_headline":"Dense sub-Saturn found on the edge of the Neptunian ridge","feed_subtitle":"At 5.83 days and 1.62 g/cm3 it sits where two migration histories compete.","key_machinery":"The argument is carried by joint modeling of TESS photometry and PARAS-2 radial velocities within EXOFASTv2, which simultaneously fits the transit, radial-velocity, spectral energy distribution, and MIST stellar models to produce stellar and planetary parameters; a circular orbit is adopted because AIC/BIC strongly favor it. The population placement depends on the period-radius and density-period diagrams defined by Castro-González et al. (2024a,b), which separate the Neptunian ridge from the desert and savanna. Internal-structure claims rest on the photoevolver code, using the Otegi et al. (2020) rocky-core mass-radius relation and the Chen & Rogers (2016) envelope model to solve for core mass and envelope fraction. The migration assessment uses the vZLK timescale and the short-range force parameters $\\epsilon_\\mathrm{GR}$, $\\epsilon_\\mathrm{rot}$, and $\\epsilon_\\mathrm{td}$ from Liu et al. (2015) to evaluate whether the wide stellar binary can excite the eccentricity needed for high-eccentricity tidal migration.","core_discovery":"TOI-6038 A b is a dense sub-Saturn with $M_P = 78.5^{+9.5}_{-9.9}\\,M_\\oplus$, $R_P = 6.41^{+0.20}_{-0.16}\\,R_\\oplus$, and a circular orbit with period $P = 5.8267311^{+0.0000074}_{-0.0000068}$ days around a $1.291^{+0.066}_{-0.060}\\,M_\\odot$ late F-type star. The paper's central claim is that this object occupies the outer boundary of the Neptunian ridge, the orbital-period over-density of super-Neptunes and sub-Saturns at $3.2$--$5.7$ days, and that its density of $1.62^{+0.23}_{-0.24}\\,\\mathrm{g\\,cm^{-3}}$ places it with the dense ridge population ($\\rho_P \\simeq 1.5$--$2.0\\,\\mathrm{g\\,cm^{-3}}$) rather than the low-density savanna. Internal-structure modeling yields a rocky/iron core of about $58\\,M_\\oplus$ and an H/He envelope mass fraction $f_\\mathrm{env} \\approx 0.26$, so the planet is core-dominated. The authors further argue that vZLK oscillations induced by the wide companion at $\\approx 3217$ AU are probably suppressed by relativistic apsidal precession, leaving high-eccentricity migration by an undetected companion or early disk-driven migration as open formation routes.","pith_inferences":["If a future Rossiter-McLaughlin measurement finds a small spin-orbit angle, the dynamically hot migration channel would be largely ruled out for this planet and disk-driven migration would become the leading explanation, implying that the ridge population may contain many dynamically cold members.","The position of this planet at the edge of the ridge makes the boundary an empirical probe: discovering more planets just inside and outside 5.7 days would reveal whether the ridge-savanna transition in density is sharp or continuous.","A metastable-helium observation that sets an upper limit on atmospheric escape would independently test the thin-envelope structure; escape rates much larger than expected for a core-dominated planet would call the internal-structure model into question.","The authors' preliminary result that GR precession suppresses vZLK-driven migration for a formation radius near 5 AU suggests that similar transition planets in wide binaries could serve as clean tests of whether distant stellar companions ever drive close-in migration."],"forward_implications":["TOI-6038 A b becomes a benchmark for the ridge-savanna transition: any successful theory of Neptunian-planet migration and evolution should reproduce a dense, core-dominated planet on a 5.83-day circular orbit in a wide binary.","The density dichotomy between ridge and savanna is reinforced by this object, supporting the view that the two regimes correspond to different formation or migration sub-populations rather than a single mass-loss sequence.","With a bright host star, the planet is observable for spin-orbit (Rossiter-McLaughlin) and atmospheric-escape (e.g., metastable helium) studies, which can discriminate between the proposed migration histories.","The derived core mass of about $58\\,M_\\oplus$ and envelope fraction of $0.26$ provide a constraint on the timing of gas accretion and disk dissipation for sub-Saturns.","The system adds to the small sample of hot sub-Saturns in multiple star systems, enabling statistical tests of whether wide companions trigger high-eccentricity migration."],"supporting_citations":[{"why":"Supplies EXOFASTv2, the joint-fitting suite used to model TESS photometry and PARAS-2 radial velocities together.","marker":"Eastman et al. (2019)"},{"why":"Defines the Neptunian ridge, desert, and savanna boundaries used to place the planet in period-radius space.","marker":"Castro-González et al. (2024a)"},{"why":"Establishes the density dichotomy between ridge and savanna and provides the TESS-cont contamination tool used to evaluate nearby sources.","marker":"Castro-González et al. (2024b)"},{"why":"Provides the rocky core mass-radius relation used by photoevolver to derive the core mass.","marker":"Otegi et al. (2020)"},{"why":"Supplies the envelope structure model used to compute envelope radius and mass fraction as a function of irradiation and age.","marker":"Chen & Rogers (2016)"},{"why":"Documents the PARAS-2 radial velocity pipeline and speckle reduction steps, and reports the similarly dense ridge planet TOI-6651 b used for comparison.","marker":"Baliwal et al. (2024)"},{"why":"Gives the short-range force parameters used to assess vZLK suppression in the binary system.","marker":"Liu et al. (2015)"},{"why":"Provides the tidal circularization timescale equation used to argue the planet's orbit may not yet be circularized.","marker":"Adams & Laughlin (2006)"}],"fun_headline_variants":["Dense sub-Saturn straddles the Neptunian ridge edge","TOI-6038 A b: dense world at a migration crossroads","Edge-of-ridge planet challenges migration scenarios","Dense sub-Saturn in the Neptunian ridge transition zone"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The measured mass and radius are secure, but the inferred core mass of $58\\,M_\\oplus$ and envelope fraction of $0.26$ assume a two-layer planet whose core obeys the Otegi et al. (2020) rocky mass-radius relation and whose envelope follows Chen & Rogers (2016); if the true interior composition deviates from these models, the compositional claims would change.","fun_headline_variants_meta":{"raw":{"variants":["Dense sub-Saturn straddles the Neptunian ridge edge","TOI-6038 A b: dense world at a migration crossroads","Edge-of-ridge planet challenges migration scenarios","Dense sub-Saturn in the Neptunian ridge transition zone"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000937,"raw_usage":{"total_tokens":4249,"prompt_tokens":1426,"completion_tokens":2823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1042,"completion_tokens_details":{"reasoning_tokens":2747}},"tokens_in":1042,"tokens_out":2823,"duration_ms":21143,"temperature":1.0,"reasoning_tokens":2747,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:13:50.018734+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Rossiter-McLaughlin effect in the stellar lines during transit: a large spin-orbit misalignment would contradict the paper's conclusion that the wide companion cannot drive high-eccentricity migration and would reopen that channel, while a small misalignment would favor disk-driven, coplanar migration.","supporting_citations":[],"review_version":1}