Planetary systems

arXiv:2607.07182v1 Announce Type: new Abstract: The LTT3780 system hosts two known exoplanets-LTT3780b, a rocky super-Earth, and LTT3780c, a temperate sub-Neptune-orbiting a nearby M dwarf on opposite sides of the radius valley. LTT3780c has been proposed as a candidate Hycean world, making the system an important target for astrobiological investigation, particularly in light of recent JWST atmospheric observations. Although biosignature and technosignature searches both seek evidence of life beyond Earth, these approaches have historically been pursued independently. Well-characterized exoplanet systems provide an opportunity to combine these complementary search strategies. In this work, we conducted radio technosignature observations of the LTT3780 system using both the Allen Telescope Array (ATA) and the Karl G. Jansky Very Large Array (VLA). The two facilities provide complementary observational capabilities, with the ATA optimized for wide-band multi-beam post-processing analyses and the VLA enabling high-sensitivity real-time interferometric searches. Across approx 30 hr of total observing time, we searched for narrowband Doppler-drifting signals in the frequency range approx 1--10 GHz. After applying comprehensive radio-frequency interference mitigation and multi-beam consistency tests, no candidate signals consistent with astrophysical or technosignature origins were identified. We place minimum detectable effective isotropic radiated power limits of 4.7 X 10^12--3.6 X 10^13W across the observed bands and facilities. Although no technosignatures were detected, this work demonstrates how complementary observation and analysis strategies can be applied to exoplanets of astrobiological interest and serves as a pathfinder for future combined biosignature and technosignature investigations.

arXiv:2607.02660v1 Announce Type: new Abstract: Protoplanetary disks often host substructures such as rings and gaps, which trace key processes in planet formation and dust evolution. However, narrow rings may remain unresolved due to limited observational resolution, hiding critical information about early planetesimal formation and the amount of dust present. We apply the azimuthal brightness modulation method, based on the modulation produced by multiple unresolved optically thick rings embedded in an optically thin background Scardoni+2024, to multi-wavelength ALMA observations of CI Tau to identify unresolved rings and constrain their geometry and optical depth. We analysed CI Tau archival ALMA continuum observations in bands 3, 6, and 7, extracting azimuthal brightness profiles along narrow annuli and comparing them with forward modeled synthetic observations of inclined disks containing unresolved rings. We detect the azimuthal signature at $\sim22$ au in all three bands, consistent with unresolved, optically thick rings embedded in a lower optical depth background. Multi wavelength modelling constrains the rings' geometry and optical depth, consistent with conditions expected for streaming instability and early planetesimal formation. Our results demonstrate the applicability of this azimuthal signature technique to real disks, reveal fine scale dust substructures in CI Tau, and illustrate a new method to study the early stages of planet formation below the nominal resolution limit.

arXiv:2607.00192v1 Announce Type: new Abstract: The disruption and accretion of planetary material onto white dwarfs is expected to be inherently dynamic and stochastic, potentially driving variability in the accretion rate and therefore the shape and depth of the photospheric metal absorption lines. This paper presents an 18-year optical spectroscopic monitoring campaign of five warm (11,000-23,000K) polluted white dwarfs with sinking timescales of days-months, observed using Magellan/MIKE and SALT/HRS to directly test this prediction. At four of the five systems, no statistically significant variability is detected over baselines of 15-18 years corresponding to hundreds to thousands of diffusion timescales, with inferred accretion rates stable to within 15-30% (1$\sigma$) showing remarkably stable accretion on decadal timescales. This implies that either the processes maintaining the accretion of the disrupted planetary material are stable on the same timescales, or that currently uncharacterized photospheric processes act to smooth observable abundance variations on these timescales. The one exception, WD 0106$-$328, shows statistically significant variability in the 4481A Mg II doublet from the ground-based data. Yet no significant equivalent width or abundance changes are seen between two Hubble Space Telescope ultraviolet spectra taken in 2016 and 2025, despite probing a larger set of transitions. This may imply that the ground-based observations witnessed a stochastic excursion from a stable baseline accretion rate, rather than a sustained change in the bulk accretion rate.

arXiv:2606.24757v1 Announce Type: new Abstract: Earth and Venus represent two evolutionary outcomes arising from initially molten 'magma ocean' periods, followed by lifetimes of chemical and geophysical divergence. Their physics is common to all rocky planets and is accessible to simulations that adopt coupled interior-atmosphere modelling approaches. Our understanding of planet histories and interpretation of current states is dependent on this modelling, yet existing codes vary in their approximations. Here, we present the first results from the Coupled atmospHere Interior modeL Intercomparison (CHILI) project; benchmarking planetary evolution codes in the context of Earth and Venus to identify key model sensitivities. Our 'nominal' Earth models predict magma ocean solidification timescales within 4 Myr of thermal evolution, and are consistent with empirical constraints on Earth's early history. Venus scenarios exhibit more diverse behaviours where prolonged magma ocean stages can be conditionally sustained for 50 Myr. Cooling timescales correlate with initial hydrogen and carbon budgets, but model-specific treatments of volatile partitioning and vertical energy transport introduce substantial inter-model variance. Different parametrisations of mantle geodynamics, convection, melting curves, rheological properties, and radiative transfer give rise to divergent evolutionary behaviours. Discrepancies in atmospheres generated by magma ocean outgassing underscore these differences, although C-H-O compositions with surface pressures exceeding 100 bar are favoured. This intercomparison identifies critical sensitivities in volatile partitioning, escape processes, mantle viscosity, and melting. Validating these treatments is essential for enabling deep insight into the early histories of the Solar System's terrestrial planets, and for drawing meaningful interpretations from ongoing observational exoplanet campaigns.

arXiv:2606.30644v1 Announce Type: new Abstract: Hydrogen-silicate miscibility can significantly alter the interior structure and thermal evolution of sub-Neptunes. We consider the interplay between this miscibility and stellar-driven atmospheric escape. We find that, for the first $\sim 100$ Myrs, sub-Neptunes store most of their hydrogen content within their miscible interiors, protecting it from escape. As hydrogen is removed from the top of the atmosphere, more hydrogen is exsolved from the miscible interior, resupplying the envelope mass and delaying the planet's contraction when compared with models that do not account for miscibility. Regardless of miscibility, atmospheric escape reproduces the young planet observations from TESS, and we highlight the emergence of the primordial Neptune desert at short orbital periods. We construct a population-level test for the prevalence of miscible sub-Neptunes which exploits their slower radial contraction. We find that $\sim 70-100$ observed young sub-Neptunes with ages $\lesssim 100$ Myrs are required to answer this question.

arXiv:2606.26381v1 Announce Type: new Abstract: The PDS 70 cavity hosts two confirmed directly imaged protoplanets and a third inner planet candidate (~13au). Despite its Keplerian motion, its unusually blue spectrum challenges a planetary interpretation. We further investigate the presence and nature of a third inner planet using new SPHERE and GRAVITY+ observations. Using the star-hopping strategy, we obtained coronagraphic IRDIS polarimetric observations in the H-band, and non-coronagraphic observations with IRDIFS in the YJHK-bands. We also searched for a planetary signal with GRAVITY in the 4UT configuration. We consistently detect two elongated inner emissions with SPHERE: the previously proposed planet candidate and another feature that appears to share the same orbit while leading it by ~120$^\circ$. Both features show dust-scattered-light spectra but different colors, possibly indicating different grain sizes. Such configuration is consistent with co-orbital dust accumulated at the stable Lagrangian regions of a distinct and yet undetected planet. GRAVITY yields a marginal (3$\sigma$) detection at the predicted location along the same orbit ($\rho=76.2\pm0.29$mas, PA=$226.50\pm0.21^\circ$), and consistent with a ~3$M_{\rm Jup}$ planet. This new planet candidate is aligned with a narrow shadow that we detect in the outer disk. We also detect polarized emission very close to the star likely arising from the inner disk ($i\sim50^\circ$, PA $\sim135^\circ$). The apparent embedding of the two dust clumps within it motivates an inner-disk origin as an alternative scenario. We conclude that the previously reported third planet candidate traces a dust clump either trailing an unseen planet on the same orbit or a rotating substructure within the inner disk. Further observations are needed to test these scenarios. Confirming the new GRAVITY planet candidate would support co-orbital substructures as indirect tracers of protoplanets.

arXiv:2606.26026v1 Announce Type: new Abstract: Giant impacts, collisions between planet-sized bodies, play an important role in planet and moon formation. As we enter a new era of large-scale surveys, such as \textit{Gaia} and LSST at the Vera C. Rubin observatory, there is potential to directly observe the remnants produced in such events and gain insights into the process of planet formation. Here, by modelling the emission and cooling of a series of giant impact remnants, we show that giant impacts are detectable as a sudden brightening followed by gradual dimming in the optical and near-infrared. Giant impacts between Earth-composition planets were simulated using the smoothed particle hydrodynamics code {\small SWIFT}, with total colliding masses ranging from 0.2 to 4 $M_{\oplus}$. By constraining the location of the photic surface of the post-impact bodies produced by the simulations, the initial luminosities of post-impact bodies were found to be between $5\times10^{-5}$ and $10^{-1}$ solar luminosities, with luminosity falling roughly exponentially on a timescale between 1 and 2000 days. Based on our results, along with estimates for planet and giant impact occurrence rates, we anticipate that between 0 and 14 terrestrial giant impacts will be observed in the full release of \textit{Gaia} epoch photometry in DR4, with at least a comparable number found by LSST. Identifying the remnants of multiple giant impacts will offer a powerful constraint on the frequency of giant impacts in the galaxy and hence the role of such collisions in planet formation.

arXiv:2606.18338v1 Announce Type: cross Abstract: The search for life beyond Earth will depend on detecting faint signatures in the atmospheres of potentially habitable exoplanets. Interpreting those signatures requires understanding the host planet's climate: the same molecule may signal life on one planet and abiotic chemistry on another. Global climate models (GCMs) provide this understanding, but individual runs can require up to millions of core-hours and substantial domain expert time. Machine-learning emulators could remove this bottleneck, but progress has been limited by the absence of a curated, multi-model exoclimate dataset. We introduce ThousandWorlds, an ML-ready benchmark for exoclimate emulation and for the broader regime of low-data, multi-simulator, parameter-to-field regression. The dataset contains approximately 1800 simulations from five GCMs, mapping eight planet parameters to 3D atmospheric fields including temperature, humidity, winds, clouds, and radiation. Three nested subsets define progressively harder challenges: single-simulator regression, multi-simulator regression with complete observations, and multi-simulator regression with structured missingness. We propose two evaluation protocols: one for ranking methods, and one that measures performance relative to the disagreement between GCMs themselves. We evaluate seven baselines spanning simple methods, deep learning, and Gaussian processes. GP-based methods perform best, suggesting that ThousandWorlds exposes a regime where off-the-shelf deep learning does not yet succeed. Data: https://doi.org/10.57967/hf/8695. Code: https://github.com/edstevenson/ThousandWorlds.

arXiv:2606.18338v1 Announce Type: cross Abstract: The search for life beyond Earth will depend on detecting faint signatures in the atmospheres of potentially habitable exoplanets. Interpreting those signatures requires understanding the host planet's climate: the same molecule may signal life on one planet and abiotic chemistry on another. Global climate models (GCMs) provide this understanding, but individual runs can require up to millions of core-hours and substantial domain expert time. Machine-learning emulators could remove this bottleneck, but progress has been limited by the absence of a curated, multi-model exoclimate dataset. We introduce ThousandWorlds, an ML-ready benchmark for exoclimate emulation and for the broader regime of low-data, multi-simulator, parameter-to-field regression. The dataset contains approximately 1800 simulations from five GCMs, mapping eight planet parameters to 3D atmospheric fields including temperature, humidity, winds, clouds, and radiation. Three nested subsets define progressively harder challenges: single-simulator regression, multi-simulator regression with complete observations, and multi-simulator regression with structured missingness. We propose two evaluation protocols: one for ranking methods, and one that measures performance relative to the disagreement between GCMs themselves. We evaluate seven baselines spanning simple methods, deep learning, and Gaussian processes. GP-based methods perform best, suggesting that ThousandWorlds exposes a regime where off-the-shelf deep learning does not yet succeed. Data: https://doi.org/10.57967/hf/8695. Code: https://github.com/edstevenson/ThousandWorlds.

arXiv:2605.26142v1 Announce Type: new Abstract: White dwarf stars, the endpoint of stellar evolution for 97% of stars in our Milky Way, offer a unique and powerful window into the bulk elemental composition of rocky exoplanetary bodies. Up to 50% of single white dwarfs are observed with photospheric metal lines from accreted exoplanetary bodies (called 'polluted' white dwarfs), and spectroscopic observations reveal the bulk composition of this material. High-resolution (R>15,000) UV spectra are essential for detecting many elements present in the material, such as the volatile elements imperative for habitability studies (C, N, O, P, S) and key rock-forming elements required to constrain interior structure (e.g. Fe, Si, Mg, Al, Ni). HST, through its COS and STIS spectrographs, remains the only facility capable of performing this science in the near future. Looking to the next decade, the scientific case for continued HST UV observations of polluted white dwarfs is compelling on three fronts (i) as a standalone to enable the bulk composition of exoplanetary material to be measured in a statistically significant sample, (ii) as essential groundwork for the Habitable Worlds Observatory (HWO), and (iii) in a powerful synergy with JWST, to enable characterization of the bulk mineralogy and bulk elemental composition of exoplanetary material. This white paper argues that continued UV spectroscopic capabilities with HST is a high-return investment for white dwarf and exoplanet science, and preserving and prioritizing HST's UV capabilities through at least 2035 is crucial to maximize the scientific return from HST, JWST, and HWO.

arXiv:2605.26152v1 Announce Type: new Abstract: Almost every known planet host will evolve into a white dwarf, and the surviving planetary material will continue to orbit this stellar remnant. Asteroids perturbed onto star-grazing orbits will become disrupted, forming an accretion disk which causes "enrichment" of the otherwise pure hydrogen or helium atmosphere. Measurements of these photospheric abundances give detailed insights into the interior compositions of exo-planetesimals with an accuracy not possible for intact exoplanets around main sequence stars. This method has revealed the diversity of rocky material in our solar neighborhood, including primitive, chondritic planetesimals, fragments of planetary cores, and even analogues of Kuiper belt objects. The planetesimal abundances can be used as an input to interior structure models. The far-ultraviolet is a key wavelength range for this field because it contains strong transitions for almost every element of interest, many of which are undetectable using ground-based optical spectroscopy. Without the FUV, we will no longer have access to the C, N, P, S content of exoplanetary bodies and thus will no longer be able to probe how volatiles interact with refractories, which is crucial to understanding planet formation-and even the origin of life. The medium resolution and high sensitivity of COS on HST has been indispensable in determining the compositions of dozens of exo-planetesimals. However, the only two medium resolution FUV-capable spectrographs are currently onboard HST, with no plans for replacements until the 2040s. An extension to the HST mission is critical for the field of white dwarf planetary systems, because the loss of FUV capability would leave us blind to volatiles. Boosting the orbit of HST would allow us to measure volatile abundances, determine the rocky planetary occurrence rate, investigate differentiation, and probe for photospheric abundance variability.

arXiv:2606.04283v1 Announce Type: new Abstract: The close-in small planet population may be formed either with hydrogen/helium dominated envelopes or with water-rich interiors. Both scenarios reproduce the present day planet population in mass, radius, and periods, and are difficult to differentiate with the mature planet demographic. Hydrogen/Helium `gas-dwarfs' have low mean molecular weight atmospheres, while `water-worlds' have envelopes that are significantly heavier, and as such these two scenarios have different evolution tracks that diverge in the first ~50 Myr of their evolution. We show that a low Earth orbit multi-band photometric survey mission, within the scope of the NASA Small Explorers Program (SMEX), can determine the frequency of young close-in planets at the 5% level and definitively differentiate between the competing `gas-dwarf' and `water-world' hypotheses. We simulate a 2.5 year mission capable of simultaneous multi-band near-ultraviolet (NUV), optical, and near infrared (NIR) wide field photometry. Such a mission would perform a photometric survey of 30 different stare-fields selected to probe the young star population. The mission will yield ~100 transiting planets in young star clusters and associations with ages <50 Myr. In comparison, only 20 such planets are known from K2 and TESS today.

arXiv:2605.04840v1 Announce Type: new Abstract: Magma ocean crystallisation sets up the early structure and long-term evolution of terrestrial planets. Recent seismic evidence signals the presence of a silicate layer at the base of Mars' mantle. Magma-ocean crystallisation and subsequent overturn has been invoked as a hypothesis for this layer's origin. However, while a magma ocean existed in both Earth and Mars, there is no seismic evidence for a basal layer in present-day Earth. In this study, we apply a parameterized-convection model to study whether the effect of partial melting in the growing mantle on overlying magma ocean composition can explain this discrepancy. Melts from the mantle buffer the crystallising magma ocean, limiting progressive differentiation, iron enrichment and the density anomaly of the overturned layer. This buffering is more efficient for larger planets with more vigorous mantle convection and for planets that are originally less enriched in iron. Consequently, a shallow magma ocean is more iron enriched and denser on Mars than on Earth, providing an explanation for the Mars-Earth difference in present-day structure of the mantle. We also predict a dichotomy in terrestrial-exoplanet interior structures, with a population with small, stratified mantles and another with large, mostly-homogeneous mantles.

arXiv:2605.03009v1 Announce Type: new Abstract: The dust observed in debris disks is the result of a collisional cascade initiated from $\sim$ km-sized parent bodies. Using near-infrared to sub-millimeter observations, we can probe particle sizes spanning 2-3 orders of magnitude, and with sufficient angular resolution we can follow the dynamics of these dust particles. Observations taken as part of the ALMA ARKS program allowed for a detailed comparison with near-infrared scattered light observations, at unprecedented resolution. The comparison between the two wavelength regimes reveals that for most gas-bearing debris disks, the distribution of small dust grains peaks outward of the distribution of large dust grains. In this paper we investigate whether gas-dust interactions can explain such radial offsets. We perform numerical simulations and compute surface brightness profiles at several wavelengths to assess which parameters drive these radial offsets. We find that while larger gas masses lead to more efficient outward radial drift, the resulting radial offset strongly depends on the optical depth of the disk, as the drift efficiency directly competes with the particles' collisional lifetime. We also find that increasing the relative number of $\mu$m-sized dust grains usually yields a larger radial offset between scattered light and millimeter observations. Finally, we show that mid-infrared observations can complement near-infrared and sub-millimeter images, and we discuss the formation of secondary rings at near-infrared wavelengths. The angular resolution achieved by the ARKS program has opened a new avenue to study the dynamics of dust particles in debris disks, revealing unexpected differences between the appearance of the disks scattered light and thermal emission. We showed that gas-dust interactions can explain the observed radial offsets and provide pointers as to which parameters have the most significant impact.

arXiv:2512.05816v2 Announce Type: replace Abstract: Sub-Neptunes are the most common type of detected exoplanet, yet their observed masses and radii are degenerate with several interior structures. One possibility is that sub-Neptunes have silicate/iron interiors and H$_2$-dominated atmospheres ($\mu$<3.8 g mol$^{-1}$), i.e., they are 'gas dwarfs'. If gas dwarfs have molten interiors, interactions between their magma oceans and atmospheres will produce distinct observational signatures. These signatures may break the degeneracy in interior structure, while providing insight into their interior processes, history, and population trends. We expect all such planets are born molten, but under what conditions do they remain molten today? We use the coupled interior-climate evolution model, PROTEUS, to estimate the 'solidification shoreline': the instellation flux boundary (as a function of stellar $T_{\rm eff}$) that separates molten gas dwarfs from solidified ones. Our results show that 98% of detected sub-Neptunes occupy a region of parameter space consistent with their having permanent magma oceans, if they are gas dwarfs. While mantle $f{\rm O}_2$ and bulk volatile C/H ratio both influence magma ocean cooling, planets with oxidising mantles and carbon-rich atmospheres are likely to have high mean-molecular weight atmospheres ($\mu$>3.8 g mol$^{-1}$) and are thus outside the scope of this study. Therefore, most detected sub-Neptunes, if they are gas dwarfs, have permanent magma oceans. This result motivates further research into the interactions between molten interiors and overlying atmospheres, and campaigns to identify unambiguous signatures of these interactions.

arXiv:2511.16734v2 Announce Type: replace Abstract: We present measurements of key protoplanetary disk properties inferred from parametric models of ALMA 12CO spectral line visibilities. We derive gas-disk radii, integrated fluxes, optically thick emission layers, and brightness temperature profiles for the disk population of the old (4 - 14 Myr) Upper Scorpius star-forming region. We measure CO emission sizes for 37 disks with bright CO J=3-2 emission (S/N > 10 on the integrated flux; out of the 83 disks with CO detections), finding that the median radius containing 90% of the flux is ~84 au, with radii spanning from 23 up to 243 au. We report a correlation between the 12CO brightness temperatures and stellar luminosities, with a Pearson coefficient of 0.6, and we use it to prove that the 12CO optically thick emission layer primarily emanates from a region below the super-heated dust, which is optically thin to the stellar irradiation. Moreover, we derive 33 CO emission surface height profiles, finding a median aspect ratio ~ 0.16 in a range from ~0.01 up to ~0.45 over the sample. Finally, we comment on the multiple systems in our sample, of which only some were already known. These results re-affirm how it is possible to derive bulk disk properties by modeling moderate angular resolution ALMA visibilities.

arXiv:2604.19743v1 Announce Type: new Abstract: Disks (Keplerian or otherwise, particulate or fluid) are often assumed to have densities that drop off vertically as Gaussians. Recent mm-wave imaging of circumstellar debris disks contradicts this assumption, revealing vertical profiles in dust that resemble Lorentzians. As part of the ARKS ALMA Large Program, we calculate how Lorentzians and Gaussians define an evolutionary sequence for disks of gravitationally scattering (viscously stirring) particles. When orbits are crossing and eccentricities $e \gg$ inclinations $i$, each scattering changes a particle's inclination by $\pm \,\Delta i \propto i$. A random walk with fixed steps in $\Delta i/i = \Delta \ln i$ produces a log normal $i$ distribution, whose thick tail at large $i$ leads to thick Lorentzian tails in density. This result holds independent of the origin of the large eccentricities; what matters is that relative motions parallel to the disk midplane are faster than perpendicular motions. After enough scatterings, $i$ comes into equipartition with $e$, $\Delta i$ stops exponentiating, and the vertical density profile relaxes to a Gaussian. We estimate the numbers and masses of perturbers needed to stir themselves and observable dust grains in Lorentzian and Gaussian debris disks imaged by ARKS. The big bodies may be sufficiently few in number as to be collisionless, in which case their masses range from the Moon to several Earths. But if Pluto-sized or smaller, the big body stirrers may be so numerous and collide so frequently that they can source the collisional cascades that produce observable dust.

arXiv:2604.15035v2 Announce Type: replace Abstract: Following the 2022 alert of a TESS object of interest transiting TOI-5624 (a G7 V star $\sim$100 pc away), a CHEOPS campaign in 2023 detected four planetary signals at $P_b\approx3.4$, $P_c\approx7.9$, $P_d\approx13.7$, and $P_e\approx21.5$ days, later confirmed by additional TESS and CHEOPS photometry in 2024-2025. After analysing the TESS & CHEOPS photometric data, we extracted and modelled the HARPS-N & SOPHIE RV time series using two independent methodologies both within an MCMC framework. We further integrated the N-body equations of motion, while simultaneously fitting the transit times and the detrended RVs, to dynamically characterise the system. We present the discovery of four transiting sub-Neptunes with radii of $R_b=2.314\pm0.035 R_{\oplus}$, $R_c=2.474\pm0.042 R_{\oplus}$, $R_d=3.584_{-0.050}^{+0.051} R_{\oplus}$, and $R_e=3.247_{-0.043}^{+0.042} R_{\oplus}$ and masses of $M_b=9.4\pm1.4 M_{\oplus}$, $M_c=4.8\pm1.9 M_{\oplus}$, $M_d=4.9\pm2.2 M_{\oplus}$, and $M_e=8.9_{-3.0}^{+2.9} M_{\oplus}$. Our photometric analysis reveals that the outermost transiting planet TOI-5624 e shows significant TTVs. We find a robust Keplerian signal in the RV time series close to the 2:1 period commensurability with TOI-5624 e, which explains the TTV pattern exhibited by TOI-5624 e according to our dynamical analysis. We label this non-transiting planet as TOI-5624 f and find its minimum mass to be $M_f\sin{i_f}=13.0\pm3.7 M_{\oplus}$. Among the known systems hosting more than four planets, the remarkable precision with which the radii have been measured (<1.7%) and the firm assessment (>3$\sigma$) of the mass for at least three planets has been previously reached only for TRAPPIST-1. Additional photometric observations will enable a better sample of the TTV modulation and a more robust dynamical determination of the masses.

arXiv:2505.10539v2 Announce Type: replace Abstract: The first transmission spectrum of the habitable-zone sub-Neptune K2-18 b with JWST has opened a new avenue for atmospheric characterisation of temperate low-mass exoplanets. The observations led to inferences of methane and carbon dioxide, as well as of dimethyl sulfide (DMS) and/or dimethyl disulfide (DMDS), both potential biosignatures. In the present work we conduct a broad and agnostic search for other chemical species in the atmosphere of K2-18 b. Our exploration includes 661 molecules, spanning a wide range of trace gases, including biotic, abiotic, and anthropogenic gases on Earth. We investigate possible preference for any of these gases, compared to a model only including the previously-detected CH$_4$ and CO$_2$, using three metrics: (a) preference in the JWST mid-infrared (MIR) spectrum, (b) preference in the JWST near-infrared spectrum, for species preferred in MIR and (c) plausible sources of production. We find that only DMS consistently results in Bayes factors $\ln B \geq 2.0$ across the datasets considered independently, though in the near-infrared this depends on detector offsets, as previously reported. The threshold of $\ln B \geq 2.0$ is motivated by the conventional threshold of $\ln B \geq 2.5$ for moderate preference, allowing for an empirical uncertainty of 0.5. A few other gases also provide comparable fits to a subset of the data or only with some of the retrieval codes used, but with limited known plausible sources. Our study highlights the need for further observations to distinguish between possible trace gases in the atmosphere of K2-18 b and theoretical work to establish their plausible sources.

arXiv:2604.15035v1 Announce Type: new Abstract: Following the 2022 alert of a TESS object of interest transiting TOI-5624 (a G7 V star $\sim$100 pc away), a CHEOPS campaign in 2023 detected four planetary signals at $P_b\approx3.4$, $P_c\approx7.9$, $P_d\approx13.7$, and $P_e\approx21.5$ days, later confirmed by additional TESS and CHEOPS photometry in 2024-2025. After analysing the TESS & CHEOPS photometric data, we extracted and modelled the HARPS-N & SOPHIE RV time series using two independent methodologies both within an MCMC framework. We further integrated the N-body equations of motion, while simultaneously fitting the transit times and the detrended RVs, to dynamically characterise the system. We present the discovery of four transiting sub-Neptunes with radii of $R_b=2.314\pm0.035 R_{\oplus}$, $R_c=2.474\pm0.042 R_{\oplus}$, $R_d=3.584_{-0.050}^{+0.051} R_{\oplus}$, and $R_e=3.247_{-0.043}^{+0.042} R_{\oplus}$ and masses of $M_b=9.4\pm1.4 M_{\oplus}$, $M_c=4.8\pm1.9 M_{\oplus}$, $M_d=4.9\pm2.2 M_{\oplus}$, and $M_e=8.9_{-3.0}^{+2.9} M_{\oplus}$. Our photometric analysis reveals that the outermost transiting planet TOI-5624 e shows significant TTVs. We find a robust Keplerian signal in the RV time series close to the 2:1 period commensurability with TOI-5624 e, which explains the TTV pattern exhibited by TOI-5624 e according to our dynamical analysis. We label this non-transiting planet as TOI-5624 f and find its minimum mass to be $M_f\sin{i_f}=13.0\pm3.7 M_{\oplus}$. Among the known systems hosting more than four planets, the remarkable precision with which the radii have been measured (<1.7%) and the firm assessment (>3$\sigma$) of the mass for at least three planets has been previously reached only for TRAPPIST-1. Additional photometric observations will enable a better sample of the TTV modulation and a more robust dynamical determination of the masses.