The radii of small exoplanets form two populations, super-Earths and sub-Neptunes, separated by a gap known as the radius valley. This feature could be produced by the removal of atmospheres by stellar or internal heating or by the lack of an initial envelope. We used transit photometry and radial velocity measurements to detect and characterize four exoplanets orbiting LHS 1903, a red dwarf star in the Milky Way's thick disk. These four planets have orbital periods ranging from 2.2 to 29.3 days and span the radius valley within a single planetary system. The derived densities indicate that LHS 1903 b is rocky, whereas LHS 1903 c and LHS 1903 d have extended atmospheres. The most distant planet from the host star, LHS 1903 e, has no gaseous envelope, indicating that it formed from gas-depleted material.

We report the discovery and characterization of the multiplanetary system around TOI-4311, a K dwarf kinematically between the Galactic thick disc and Hercules stream. TOI-4311 hosts an ultra-short-period super-Earth (P similar to 0.99 d, 1 . 376(-0.080)(+0 . 077) R-circle plus) and a longer period sub-Neptune (P similar to 15 d, 2.47(-0 . 11)(+0.12) R-circle plus) that was first detected in the Transiting Exoplanet Survey Satellite photometry. Using follow-up observations with CHaracterising ExOPlanet Satellite and High Accuracy Radial Velocity Planet Searcher (HARPS), we refine the planetary radius of both planets, derive the mass of planet b ( 4 . 5(-1.4)(+1 . 5) M-circle plus), and confirm the planetary nature of planet c. Intriguingly, a third periodic signal is clearly detected in our HARPS Radial Velocities (RVs) that we cannot link to stellar activity. This signal could be attributed to a third planet (P similar to 38 d, Msin(i) = 26.4(-6.8)(+6 . 3) M-circle plus) in the system; however, with the current photometric data set we do not find a transit. Our dynamical analysis highlights that this potential outer planet would remain stable. Using the precise radius and mass for TOI-4311 b, we model its interior structure and find that it is very dense given the host star's galactic kinematics and chemistry. Hence, this system could challenge current formation theories and provide insights into planet formation across the galaxy.

Changes in Earth’s orbit and axial tilt affect the distribution of solar radiation across the planet’s surface, influencing climate and leaving recognizable imprints in sedimentary layers over time. These recurring patterns, known as Milankovitch cycles, arise from gravitational interactions with other bodies in the Solar System. Their preservation in stratigraphic records allows the construction of continuous and precise age–depth models based on astronomical orbital cycles, helping to compensate for the discontinuities and uncertainties associated with other absolute dating methods.In this study, we examine several tens of stratigraphic records spanning 40–100 Ma to establish and refine astronomical time scales by extracting signals related to Earth’s orbital eccentricity. We apply our newly released tool, AstroGeoFit (https://www.astrogeo.eu/astrogeofit/), which accommodates flexible time-varying sedimentation rates within stratigraphic sequences and quantifies associated uncertainties by combining genetic algorithm optimization with Bayesian approaches.Beyond improving age–depth models at the scale of individual sections or cores, our results contribute towards a more systematic reconstruction of Earth’s eccentricity evolution from 40 to 100 Ma, derived primarily from stratigraphic records, providing intrinsic geological insights into the behavior of planetary orbits in deep time.

Context. One of the main bottlenecks in assessing the accuracy of the Mars orbit is the unknown value of the asteroids in the main asteroid belt. A modeling with 343 asteroids as point masses is used today, with the relative masses fit to observational data. Aims. We propose an innovative method for reducing the number of asteroids implemented as point masses, which in turn reduces the number of parameters to be fit without a significant degradation of the postfit residuals. Methods. We used boosting decision trees to obtain a ranking by relative importance of the 343 asteroids of the current main belt modeling. Results. We were able to remove more than 100 of these asteroids without significantly degrading the postfit residuals and with a significant improvement of the uncertainties. Furthermore, we verified that the postfit masses found with the new modeling (INPOP25c) are consistent with respect to the mass estimation from an independent approach by means of the albedo properties of the asteroids.

We explore a hybrid expansion of the disturbing function in planetary dynamics that combines elements of the classical Laplace and Legendre developments. This formulation retains the structure of the Laplace expansion, but expresses the inverse of the mutual distance as a series whose terms keep an exact dependence on both the eccentricity and the semi-major axis ratio. We use it to construct the first-order secular Hamiltonian of the planar 3-body problem, relevant for modeling the long-term evolution of planetary systems. We assess the convergence of the new expansion numerically and compare it with that of the Laplace and Legendre series across a range of orbital configurations. The results show that the new expansion provides consistent performance across diverse dynamical regimes, bridging the domains of applicability of the two classical approaches.

Understanding the physics of planetary magma oceans has been the subject of growing efforts, in light of the increasing abundance of solar system samples and extrasolar surveys. A rocky planet harboring such an ocean is likely to interact tidally with its host star, planetary companions, or satellites. To date, however, models of the tidal response and heat generation of magma oceans have been restricted to the framework of weakly viscous solids, ignoring the dynamical fluid behavior of the ocean beyond a critical melt fraction. Here we provide a handy analytical model that accommodates this phase transition, allowing for a physical estimation of the tidal response of lava worlds. We apply the model in two settings: the tidal history of the early Earth–Moon system in the aftermath of the giant impact, and the tidal interplay between short-period exoplanets and their host stars. For the former, we show that the fluid behavior of the Earth's molten surface drives efficient early lunar recession to ~25 Earth radii within 10 4 –10 5 yr, in contrast with earlier predictions. For close-in exoplanets, we report on how their molten surfaces significantly change their spin–orbit dynamics, allowing them to evade spin–orbit resonances and accelerating their track toward tidal synchronization from a gigayear to megayear timescale. Moreover, we reevaluate the energy budgets of detected close-in exoplanets, highlighting how the surface thermodynamics of these planets are likely controlled by enhanced, fluid-driven tidal heating, rather than vigorous insolation, and how this regime change substantially alters predictions for their surface temperatures.

Atmospheric thermal tides arise from the diurnal contrast in stellar irradiation. They exert a significant influence on the long-term rotational evolution of rocky planets because they can accelerate the planetary spin, thereby counteracting the decelerating effect of classical gravitational tides. Consequently, equilibrium tide-locked states may emerge, as exemplified by Venus and hypothesised for Precambrian Earth. Quantifying the atmospheric thermal torque and elucidating its dependence on tidal frequency – both in the low- and high-frequency regimes – is therefore essential. In particular, we focus here on the resonance that affected early Earth, which is associated with a forced Lamb wave. Within the framework of linear theory, we develop a new analytical model of the atmospheric response to both gravitational an thermal tidal forcings for two representative vertical temperature profiles that bracket the atmospheres of rocky planets: (i) an isothermal profile (uniform temperature) and (ii) an isentropic profile (uniform potential temperature). Dissipative processes are incorporated via Newtonian cooling. We demonstrate that the isothermal and isentropic cases are governed by the same general closed-form solution, and we derive explicit expressions for the three-dimensional tidal fields (pressure, temperature, density and wind velocities) throughout the spherical atmospheric shell. These results constitute the foundation for two forthcoming papers, in which analytical formulae for the thermotidal torque will be presented and compared with numerical solutions obtained from General Circulation Models (GCMs).

A crucial chemical link between stars and their orbiting exoplanets is thought to exist. If universal, this connection could affect the formation and evolution of all planets. Therefore, this potential vital link needs testing by characterizing exoplanets around chemically-diverse stars. We present the discovery of two planets orbiting the metal-poor, kinematic thick-disc K-dwarf TOI-2345. TOI-2345 b is a super-Earth with a period of 1.05 d and TOI-2345 c is a sub-Neptune with a period of 21 d. In addition to the target being observed in four TESS sectors, we obtained five CHEOPS visits and 26 radial velocities from HARPS. By conducting a joint analysis of all the data, we find TOI-2345 b to have a radius of 1 . 504+ 0 . 047-0.044 R (R) and a mass of 3 . 49 +/- 0 . 85 M-(R); and TOI-2345 c to have a radius of 2.451+0.045-0.046R (R) and a mass of 7.27+2.27-2.45M((R)). To explore chemical links between these planets and their host star, we model their interior structures newly accounting for devolatized stellar abundances. TOI-2345 adds to the limited sample of well-characterized planetary systems around thick disc stars. This system challenges theories of formation and populations of planets around thick disc stars with its Ultra-Short Period super-Earth and the wide period distribution of these two planets spanning the radius valley.

Aims: We aim to observe the transits and occultations of WASP-33b, which orbits a rapidly-rotating δ Scuti pulsator, with the goal of measuring the orbital obliquity via the gravity-darkening effect, and constraining the geometric albedo via the occultation depth. Methods: We observed four transits and four occultations with CHEOPS, and employ a variety of techniques to remove the effects of the stellar pulsations from the light curves, as well as the usual CHEOPS systematic effects. We also performed a comprehensive analysis of low-resolution spectral and Gaia data to re-determine the stellar properties of WASP-33. Results: We measure an orbital obliquity 111.3 +0.2 -0.7 degrees, which is consistent with previous measurements made via Doppler tomography. We also measure the planetary impact parameter, and confirm that this parameter is undergoing rapid secular evolution as a result of nodal precession of the planetary orbit. This precession allows us to determine the second-order fluid Love number of the star, which we find agrees well with the predictions of theoretical stellar models. We are unable to robustly measure a unique value of the occultation depth, and emphasise the need for long-baseline observations to better measure the pulsation periods.

Extrasolar planetary systems commonly exhibit planets on eccentric orbits, with many systems located near or within mean-motion resonances, showcasing a wide diversity of orbital architectures. Such complex systems challenge traditional secular theories, which are limited to first-order approximations in planetary masses or rely on expansions in orbital elements--eccentricities, inclinations, and semi-major axis ratios--that are subject to convergence issues, especially in highly eccentric, inclined, or tightly-packed systems. To overcome these limitations, we develop a numerical approach to second-order perturbation theory based on the Lie transform formalism. Our method avoids the need for expansions in orbital elements, ensuring broader applicability and more robust convergence. We first outline the Hamiltonian framework for the 3-body planetary problem, and apply a canonical transformation to eliminate fast angle dependencies, deriving the secular Hamiltonian up to second order in the mass ratio. We then use the fast Fourier transform algorithm to numerically simulate, in an accurate way, the long-term evolution of planetary systems near or away from mean-motion resonances. Finally, we validate our methods against well-known planetary configurations, such as the Sun-Jupiter-Saturn system, as well as to exoplanetary systems like WASP-148, TIC 279401253 and GJ 876, demonstrating the applicability of our models across a wide range of planetary configurations.

Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the different proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and one sub-Neptune. We observed HIP 29442 with CHEOPS for a total of 9.6 days, which we modelled jointly with two sectors of TESS data to derive planetary radii of 3.410 +/- 0.046, 1.551 +/- 0.045, and 1.538 +/- 0.049 R-circle plus for planets b, c, and d, which orbit HIP 29442 with periods of 13.6, 3.5, and 6.4 days, respectively. For planet d this value deviates by more than 3 sigma from the median value reported in the discovery paper, leading us to conclude that caution is required when using TESS photometry to determine the radii of small planets with low per-transit signal-to-noise ratios and large gaps between observations. Given the high precision of these new radii, combining them with published RVs from ESPRESSO and HIRES provides us with ideal conditions to investigate the internal structure and formation pathways of the planets in the system. We introduced the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then applied hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identified planetary system analogues in a synthetic population generated with the Bern model for planet formation and evolution. Based on this analysis, we find that the planets likely formed on opposing sides of the water iceline from a protoplanetary disk with an intermediate solid mass. We finally report that the observed parameters of the HIP 29442 system are compatible with a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope and with a scenario with water-rich sub-Neptunes.

Rotational dynamics of the Earth, over geological timescales, have profoundly affected local and global climatic evolution, probably contributing to the evolution of life. To better retrieve the Earth’s rotational history, and motivated by the published hypothesis of a stabilized length of day during the Precambrian, we examined the effect of thermal tides on the evolution of planetary rotational motion. The hypothesized scenario is contingent upon encountering a resonance in atmospheric Lamb waves, whereby an amplified thermotidal torque cancels out the opposing torque generated by the oceans and solid interior, driving the Earth into rotational equilibrium. With this scenario in mind, we constructed an ab initio model of thermal tides on rocky planets describing a neutrally stratified atmosphere. The model takes into account dissipative processes with Newtonian cooling and diffusive processes in the planetary boundary layer. We retrieved, from this model, a closed-form solution for the frequency-dependent tidal torque, which captures the main spectral features previously computed using 3D general circulation models. In particular, under longwave heating, diffusive processes near the surface and the delayed thermal response of the ground prove to be responsible for attenuating, and possibly annihilating, the accelerating effect of the thermotidal torque at the resonance. When applied to the Earth, our model prediction suggests the occurrence of the Lamb resonance in the Phanerozoic, but with an amplitude that is insufficient for the rotational equilibrium. Interestingly, though our study was motivated by the Earth’s history, the generic tidal solution can be straightforwardly and efficiently applied in exoplanetary settings.

Context . The TOI-178 system consists of a nearby, late-K-dwarf with six transiting planets in the super-Earth to mini-Neptune regime, with radii ranging from to 2.9 R ⊕ and orbital periods between 1.9 and 20.7 days. All the planets, but the innermost one, form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, thereby providing important anchors for planet formation models. Aims . We aim to improve the characterisation of the architecture of this key system and, in particular, the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we propose to use it as a test bench for the robustness of the planetary mass determination with each technique. Methods . We performed a global analysis of all the available photometry from CHEOPS, TESS and NGTS, and radial velocity from ESPRESSO, using a photo-dynamical modelling of the light curve. We also tried different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by transit-timing variations (TTVs) and radial velocities (RVs). Results . We demonstrate how stellar activity prevents a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis has resulted in a robust mass determination for planets c to 𝑔, with precision of ~ 12% for the mass of planet c, and better than 10% for planets d to 𝑔. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable for the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain.

Context. Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as 'lonely' and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed and well characterised. Aims. Within the contest of a multi-planet system hosting a gas giant on short orbits, we characterise the TOI-1130 system by measuring masses and orbital parameters. This is a two-transiting planet system with a Jupiter-like planet (c) on a 8.35 days orbit and a Neptune-like planet (b) on an inner (4.07 days) orbit. Both planets show strong anti-correlated transit timing variations (TTVs). Furthermore, radial velocity (RV) analysis showed an additional linear trend, a possible hint of a non-transiting candidate planet on a far outer orbit. Methods. Since 2019, extensive transit and radial velocity observations of the TOI-1130 have been acquired using TESS and various ground-based facilities. We present a new photo-dynamical analysis of all available transit and RV data, with the addition of new CHEOPS and ASTEP+ data, which achieve the best precision to date on the planetary radii and masses and on the timings of each transit. Results. We were able to model interior structure of planet b constraining the presence of a gaseous envelope of H/He, while it was not possible to assess the possible water content. Furthermore, we analysed the resonant state of the two transiting planets, and we found that they lie just outside the resonant region. This could be the result of the tidal evolution that the system underwent. We obtained both masses of the planets with a precision of less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively.

HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but the observations show no transit. We find that large ($>22.4$ hours) transit timing variations (TTVs) could explain this non-detection during the CHEOPS observation window. We also investigated the possibility of an incorrect orbital solution, which would have major implications for our knowledge of this system. If $P_\mathrm{d} \neq 278.36$ d, the periods that minimize the eccentricity would be $101.22$ d and $371.14$ d. The shortest orbital period will be tested by TESS, which will observe HIP 41378 in Sector 88 starting in January 2025. Our study shows the importance of a mission like CHEOPS, which today is the only mission able to make long observations (i.e., from space) to track the ephemeris of long-period planets possibly affected by large TTVs.

Earth-like planets in the habitable zone of low-mass stars undergo strong tidal effects that modify their spin states. These planets are expected to host dense atmospheres that can also play an important role in the spin evolution. On one hand, gravitational tides tend to synchronise the rotation with the orbital mean motion, but on the other hand, thermal atmospheric tides push the rotation away and may lead to asynchronous equilibria. Here, we investigate the complete tidal evolution of Earth-like planets by taking into account the effect of obliquity and eccentric orbits. We adopted an Andrade rheology for the gravitational tides and benchmarked the unknown parameters with the present rotation of Venus. We then applied our model to Earth-like planets, and we show that asynchronous rotation can be expected for planets orbiting stars with masses between 0.4 and 0.9 M-circle dot and semi-major axes between 0.2 and 0.7 au. Interestingly, we find that Earth-like planets in the habitable zone of stars with masses similar to 0.8 M-circle dot may end up with an equilibrium rotation of 24 h. We additionally find that these planets can also develop high obliquities, which may help sustain temperate environments.

The AU Microscopii planetary system is only 24 Myr old, and its geometry may provide clues about the early dynamical history of planetary systems. Here, we present the first measurement of the Rossiter-McLaughlin effect for the warm sub-Neptune AU Mic c, using two transits observed simultaneously with the European Southern Observatory's (ESO's) Very Large Telescope (VLT)/Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), CHaracterising ExOPlanet Satellite (CHEOPS), and Next-Generation Transit Survey (NGTS). After correcting for flares and for the magnetic activity of the host star, and accounting for transit-timing variations, we find the sky-projected spin-orbit angle of planet c to be in the range lambda(c) = 67.8(-49.0)(+31.7) degrees (1 sigma). We examine the possibility that planet c is misaligned with respect to the orbit of the inner planet b (lambda(b) = -2.96(-10.30)(+10.44)), and the equatorial plane of the host star, and discuss scenarios that could explain both this and the planet's high density, including secular interactions with other bodies in the system or a giant impact. We note that a significantly misaligned orbit for planet c is in some degree of tension with the dynamical stability of the system, and with the fact that we see both planets in transit, though these arguments alone do not preclude such an orbit. Further observations would be highly desirable to constrain the spin-orbit angle of planet c more precisely.

The CHaracterising ExOPlanet Satellite (CHEOPS) is a partnership between the European Space Agency and Switzerland with important contributions by 10 additional ESA member States. It is the first S-class mission in the ESA Science Programme. CHEOPS has been flying on a Sun-synchronous low Earth orbit since December 2019, collecting millions of short-exposure images in the visible domain to study exoplanet properties. A small yet increasing fraction of CHEOPS images show linear trails caused by resident space objects crossing the instrument field of view. To characterize the population of satellites and orbital debris observed by CHEOPS, all and every science images acquired over the past 3 years have been scanned with a Hough transform algorithm to identify the characteristic linear features that these objects cause on the images. Thousands of trails have been detected. This statistically significant sample shows interesting trends and features such as an increased occurrence rate over the past years as well as the fingerprint of the Starlink constellation. The cross-matching of individual trails with catalogued objects is underway as we aim to measure their distance at the time of observation and deduce the apparent magnitude of the detected objects. As space agencies and private companies are developing new space-based surveillance and tracking activities to catalogue and characterize the distribution of small debris, the CHEOPS experience is timely and relevant. With the first CHEOPS mission extension currently running until the end of 2026, and a possible second extension until the end of 2029, the longer time coverage will make our dataset even more valuable to the community, especially for characterizing objects with recurrent crossings.

The Earth’s rotation has been decelerating throughout its history due to tidal dissipation, but the variation of the rate of this deceleration through time has not been established. We present a detailed analysis of eight geological datasets to constrain the Earth’s rotational history from 650 to 240 Mya. The results allow us to test physical tidal models and point to a staircase pattern in the Earth’s deceleration from 650 to 280 Mya. During this time interval, the Earth–Moon distance increased by approximately 20,000 km and the length of day increased by approximately 2.2 h. Specifically, there are two intervals with high Earth rotation deceleration, 650 to 500 Mya and 350 to 280 Mya, separated by an interval of stalled deceleration from 500 to 350 Mya. The interval with stalled deceleration is attributed mainly to reduced tidal dissipation due to the continent-ocean configuration at the time, not to changes in Earth’s dynamical ellipticity from continental assembly or glaciation. Modeling indicates that, except for the very recent time, tidal dissipation is the main driver for decelerating Earth rotation. One potential implication of our findings is that the Earth’s tidal dissipation, along with Earth’s rotation deceleration, may play a role in the evolving Earth.