Search

In this note we study an iterative belief propagation (IBP) algorithm and demonstrate it's ability to solve sparse combinatorial optimization problems. Similar to simulated annealing (SA), our IBP algorithm attempts to sample from the Boltzmann distribution of the objective function but also uses belief propagation (BP) to improve convergence.

We propose a general framework for a hybrid continuous-discrete algorithm that integrates continuous-time deterministic dynamics with Metropolis-Hastings steps to combine search dynamics with and without detailed balance. Our purpose is to study the non-equilibrium dynamics that leads to the ground state of rugged energy landscapes in this general setting. Our results show that MH-driven dynamics reach ``easy'' ground states faster, indicating a stronger bias in the non-equilibrium dynamics of the algorithm with reversible transition probabilities. To validate this, we construct a set of Ising problem instances with a controllable bias in the energy landscape that makes one degenerate solution more accessible than another. The constructed hybrid algorithm demonstrates significant improvements in convergence and ground-state sampling accuracy, achieving a 100x speedup on GPUs compared to simulated annealing, making it well-suited for large-scale applications.

K2-24 is a planetary system composed of two transiting low-density Neptunians locked in an almost perfect 2:1 resonance and showing large TTVs, i.e., an excellent laboratory to search for signatures of planetary migration. Previous studies performed with K2, Spitzer and RV data tentatively claimed a significant non-zero eccentricity for one or both planets, possibly high enough to challenge the scenario of pure disk migration through resonant capture. With 13 new CHEOPS light curves (seven of planet -b, six of planet -c), we carried out a global photometric and dynamical re-analysis by including all the available literature data as well. We got the most accurate set of planetary parameters to date for the K2-24 system, including radii and masses at 1% and 5% precision (now essentially limited by the uncertainty on stellar parameters) and non-zero eccentricities $e_b=0.0498_{-0.0018}^{+0.0011}$, $e_c=0.0282_{-0.0007}^{+0.0003}$ detected at very high significance for both planets. Such relatively large values imply the need for an additional physical mechanism of eccentricity excitation during or after the migration stage. Also, while the accuracy of the previous TTV model had drifted by up to 0.5 days at the current time, we constrained the orbital solution firmly enough to predict the forthcoming transits for the next ~15 years, thus enabling an efficient follow-up with top-level facilities such as JWST or ESPRESSO., Comment: 19 pages, 11 figures, 8 tables. Accepted for publication in A&A on September 4, 2024. Typos corrected

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. Within the contest of multi-planet system hosting gas-giant on short orbits, we characterise TOI-1130 system by measuring masses and orbital parameters. This is a 2-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. 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 that achieve the best precision to date on the planetary radii and masses and on the timings of each transit. 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 less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively., Comment: 22 pages, 10 figures, Accepted for publication by Astronomy and Astrophysics

Co-orbital objects, also known as trojans, are frequently found in simulations of planetary system formation. In these configurations, a planet shares its orbit with other massive bodies. It is still unclear why there have not been any co-orbitals discovered thus far in exoplanetary systems or even pairs of planets found in such a 1:1 mean motion resonance. Reconciling observations and theory is an open subject in the field. The main objective of the TROY project is to conduct an exhaustive search for exotrojans using diverse observational techniques. In this work, we analyze the radial velocity time series informed by transits, focusing the search around low-mass stars. We employed the alpha-test method on confirmed planets searching for shifts between spectral and photometric mid-transit times. This technique is sensitive to mass imbalances within the planetary orbit, allowing us to identify non-negligible co-orbital masses. Among the 95 transiting planets examined, we find one robust exotrojan candidate with a significant 3-sigma detection. Additionally, 25 exoplanets show compatibility with the presence of exotrojan companions at a 1-sigma level, requiring further observations to better constrain their presence. For two of those weak candidates, we find dimmings in their light curves within the predicted Lagrangian region. We established upper limits on the co-orbital masses for either the candidates and null detections. Our analysis reveals that current high-resolution spectrographs effectively rule out co-orbitals more massive than Saturn around low-mass stars. This work points out to dozens of targets that have the potential to better constraint their exotrojan upper mass limit with dedicated radial velocity observations. We also explored the potential of observing the secondary eclipses of the confirmed exoplanets to enhance the exotrojan search., Comment: 41 pages, 9 figures, 6 tables. Accepted in Astronomy & Astrophysics

A systematic, population-level discrepancy exists between the densities of exoplanets whose masses have been measured with transit timing variations (TTVs) versus those measured with radial velocities (RVs). Since the TTV planets are predominantly nearly resonant, it is still unclear whether the discrepancy is attributed to detection biases or to astrophysical differences between the nearly resonant and non resonant planet populations. We defined a controlled, unbiased sample of 36 sub-Neptunes characterised by Kepler, TESS, HARPS, and ESPRESSO. We found that their density depends mostly on the resonant state of the system, with a low probability (of $0.002_{-0.001}^{+0.010}$) that the mass of (nearly) resonant planets is drawn from the same underlying population as the bulk of sub-Neptunes. Increasing the sample to 133 sub-Neptunes reveals finer details: the densities of resonant planets are similar and lower than non-resonant planets, and both the mean and spread in density increase for planets that are away from resonance. This trend is also present in RV-characterised planets alone. In addition, TTVs and RVs have consistent density distributions for a given distance to resonance. We also show that systems closer to resonances tend to be more co-planar than their spread-out counterparts. These observational trends are also found in synthetic populations, where planets that survived in their original resonant configuration retain a lower density; whereas less compact systems have undergone post-disc giant collisions that increased the planet's density, while expanding their orbits. Our findings reinforce the claim that resonant systems are archetypes of planetary systems at their birth.

Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the 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 a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TESS data to derive planetary radii of $3.410\pm0.046$, $1.551\pm0.045$ and $1.538\pm0.049$ R$_\oplus$ for planets b, c and d, which orbit HIP 29442 with periods of 13.6, 3.5 and 6.4 days. 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 S/N 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 introduce the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then apply hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identify 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 both a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope as well as a scenario with water-rich sub-Neptunes., Comment: 30 pages, 17 figures, accepted for publication in A&A

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All 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, hence providing important anchors for planet formation models. 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 aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try 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 each method. We show how stellar activity is preventing us from obtaining 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 resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. 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 during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain.

We propose a novel algorithm that extends the methods of ball smoothing and Gaussian smoothing for noisy derivative-free optimization by accounting for the heterogeneous curvature of the objective function. The algorithm dynamically adapts the shape of the smoothing kernel to approximate the Hessian of the objective function around a local optimum. This approach significantly reduces the error in estimating the gradient from noisy evaluations through sampling. We demonstrate the efficacy of our method through numerical experiments on artificial problems. Additionally, we show improved performance when tuning NP-hard combinatorial optimization solvers compared to existing state-of-the-art heuristic derivative-free and Bayesian optimization methods., Comment: Accepted to ICML2024

The TOI-421 planetary system contains two sub-Neptune-type planets and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. We jointly analysed photometric data of three TESS sectors and six CHEOPS visits as well as 156 radial velocity data points to retrieve improved planetary parameters. We also searched for TTVs and modelled the interior structure of the planets. Finally, we simulated the evolution of the primordial H-He atmospheres of the planets using two different modelling frameworks. We determine the planetary radii and masses of TOI-421 b and c to be $R_{\rm b} = 2.64 \pm 0.08 \, R_{\oplus}$, $M_{\rm b} = 6.7 \pm 0.6 \, M_{\oplus}$, $R_{\rm c} = 5.09 \pm 0.07 \, R_{\oplus}$, and $M_{\rm c} = 14.1 \pm 1.4 \, M_{\oplus}$. We do not detect any statistically significant TTV signals. Assuming the presence of a hydrogen-dominated atmosphere, the interior structure modelling results in both planets having extensive envelopes. While the modelling of the atmospheric evolution predicts for TOI-421 b to have lost any primordial atmosphere that it could have accreted at its current orbital position, TOI-421 c could have started out with an initial atmospheric mass fraction somewhere between 10 and 35%. We conclude that the low observed mean density of TOI-421 b can only be explained by either a bias in the measured planetary parameters (e.g. driven by high-altitude clouds) and/or in the context of orbital migration. We also find that the results of atmospheric evolution models are strongly dependent on the employed planetary structure model., Comment: 21 pages, 12 figures, accepted for publication in A&A

Co-orbital exoplanets are a by-product of the models of formation of planetary systems. However, none have been detected in nature thus far. Although challenging, the observation of co-orbital exoplanets would provide valuable information on the formation of planetary systems as well as on the interactions between planets and their host star. After a brief review of the stability and formation issues of co-orbital systems, some observational methods dedicated to their detection are presented., Comment: Hans Deeg \S Juan Antonio Belmonte. Handbook of Exoplanets, 2nd Edition, Springer International Publishing AG, part of Springer Nature, In press

TOI-1227 b is an 11 Myr old validated transiting planet in the middle of its contraction phase, with a current radius of 0.85 R$_J$. It orbits a low-mass pre-main sequence star (0.170 M$_\odot$, 0.56 R$_\odot$) every 27.4 days. The magnetic activity of its young host star induces radial velocity jitter and prevents good measurements of the planetary mass. We gathered additional transit observations of TOI-1227 b with space- and ground-based telescopes, and we detected highly significant transit-timing variations (TTVs). Their amplitude is about 40 minutes and their dominant timescale is longer than 3.7 years. Their most probable origin is dynamical interactions with additional planets in the system. We modeled the TTVs with inner and outer perturbers near first and second order resonances; several orbital configurations provide an acceptable fit. More data are needed to determine the actual orbital configuration and eventually measure the planetary masses. These TTVs and an updated transit chromaticity analysis reinforce the evidence that TOI-1227 b is a planet., Comment: 12 pages, 9 figures, accepted for publication in A&A

The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. Mass estimates derived from a preliminary radial velocity (RV) dataset suggest that the planetary densities do not decrease in a monotonic way with the orbital distance to the star, contrary to what one would expect based on simple formation and evolution models. To improve the characterisation of this key system and prepare for future studies (in particular with JWST), we perform a detailed photometric study based on 40 new CHEOPS visits, one new TESS sector, as well as previously published CHEOPS, TESS, and NGTS data. First we perform a global analysis of the 100 transits contained in our data to refine the transit parameters of the six planets and study their transit timing variations (TTVs). We then use our extensive dataset to place constraints on the radii and orbital periods of potential additional transiting planets in the system. Our analysis significantly refines the transit parameters of the six planets, most notably their radii, for which we now obtain relative precisions $\lesssim$3%, with the exception of the smallest planet $b$ for which the precision is 5.1%. Combined with the RV mass estimates, the measured TTVs allow us to constrain the eccentricities of planets $c$ to $g$, which are found to be all below 0.02, as expected from stability requirements. Taken alone, the TTVs also suggest a higher mass for planet $d$ than the one estimated from the RVs, which had been found to yield a surprisingly low density for this planet. However, the masses derived from the current TTV dataset are very prior-dependent and further observations, over a longer temporal baseline, are needed to deepen our understanding of this iconic planetary system., Comment: 20 pages, 5 figures, 9 tables. Accepted for publication in A&A

The distribution of close-in exoplanets is shaped by a complex interplay between atmospheric and dynamical processes. The Desert-Rim Exoplanets Atmosphere and Migration (DREAM) program aims at disentangling those processes through the study of the hot Neptune desert, whose rim hosts planets that are undergoing, or survived, atmospheric evaporation and orbital migration. In this first paper, we use the Rossiter-McLaughlin Revolutions (RMR) technique to investigate the orbital architecture of 14 close-in planets ranging from mini-Neptune to Jupiter-size and covering a broad range of orbital distances. While no signal is detected for the two smallest planets, we were able to constrain the sky-projected spin--orbit angle of six planets for the first time, to revise its value for six others, and, thanks to constraints on the stellar inclination, to derive the 3D orbital architecture in seven systems. These results reveal a striking three-quarters of polar orbits in our sample, all being systems with a single close-in planet but of various stellar and planetary types. High-eccentricity migration is favored to explain such orbits for several evaporating warm Neptunes, supporting the role of late migration in shaping the desert and populating its rim. Putting our measurements in the wider context of the close-in planet population will be useful to investigate the various processes shaping their architectures., Comment: 25 pages and 12 figures (plus Appendix)

The HD108236 system was first announced with the detection of four small planets based on TESS data. Shortly after, the transit of an additional planet with a period of 29.54d was serendipitously detected by CHEOPS. In this way, HD108236 (V=9.2) became one of the brightest stars known to host five small transiting planets (R$_p$<3R$_{\oplus}$). We characterize the planetary system by using all the data available from CHEOPS and TESS space missions. We use the flexible pointing capabilities of CHEOPS to follow up the transits of all the planets in the system, including the fifth transiting body. After updating the host star parameters by using the results from Gaia eDR3, we analyzed 16 and 43 transits observed by CHEOPS and TESS, respectively, to derive the planets physical and orbital parameters. We carried out a timing analysis of the transits of each of the planets of HD108236 to search for the presence of transit timing variations. We derived improved values for the radius and mass of the host star (R$_{\star}$=0.876$\pm$0.007 R$_{\odot}$ and M$_{\star}$=0.867$_{-0.046}^{+0.047}$ M$_{\odot}$). We confirm the presence of the fifth transiting planet f in a 29.54d orbit. Thus, the system consists of five planets of R$_b$=1.587$\pm$0.028, R$_c$=2.122$\pm$0.025, R$_d$=2.629$\pm$0.031, R$_e$=3.008$\pm$0.032, and R$_f$=1.89$\pm$0.04 [R$_{\oplus}$]. We refine the transit ephemeris for each planet and find no significant transit timing variations for planets c, d, and e. For planets b and f, instead, we measure significant deviations on their transit times (up to 22 and 28 min, respectively) with a non-negligible dispersion of 9.6 and 12.6 min in their time residuals. We confirm the presence of planet f and find no significant evidence for a potential transiting planet in a 10.9d orbital period, as previously suggested. Full abstract in the PDF file., Comment: 18 Figures and 25 pages. Accepted for publication in A&A

We report the characterisation of a multi-planetary system composed of five exoplanets orbiting the K-dwarf HD~23472 (TOI-174). In addition to the two super-Earths that were previously confirmed, we confirm and characterise three Earth-size planets in the system using ESPRESSO radial velocity observations. The planets of this compact system have periods of $P_d \sim 3.98\,$, $P_e \sim 7.90\,$, $P_f \sim 12.16\,$, $P_b \sim 17.67,\,$ and $P_c \sim 29.80\,$days and radii of $R_d \sim 0.75\,$ , $R_e \sim 0.82,$, $R_f \sim 1.13\,$, $R_b \sim 2.01,\,$ and, $R_c \sim 1.85\,$ \REarth. Because of its small size, its proximity to planet d's transit, and close resonance with planet d, planet e was only recently found. The planetary masses were estimated to be $M_d =0.54\pm0.22$, $M_e =0.76\pm0.30$, $M_f =0.64_{-0.39}^{+0.46}$, $M_b = 8.42_{-0.84}^{+0.83}$, and $M_c = 3.37_{-0.87}^{+0.92}$ \MEarth. These planets are among the lightest planets, with masses measured using the radial velocity method, demonstrating the very high precision of the ESPRESSO spectrograph. We estimated the composition of the system's five planets and found that their gas and water mass fractions increase with stellar distance, suggesting that the system was shaped by irradiation. The high density of the two inner planets ($\rho_d = 7.5_{-3.1}^{+3.9}$ and $\rho_e = 7.5_{-3.0}^{+3.9}\, \mathrm{g.cm^{-3}}$) indicates that they are likely to be super-Mercuries. This is supported by the modelling of the internal structures of the planets, which also suggests that the three outermost planets have significant water or gas content If the existence of two super-Mercuries in the system is confirmed, this system will be the only one known to feature two super-Mercuries, making it an excellent testing bed for theories of super-Mercuries formation. (Abridged), Comment: 24 pages, 13 figures, Published in A&A

We experimentally demonstrate that networks of identical photonic spiking neurons based on coupled degenerate parametric oscillators can show various chimera states, in which, depending on their local synchronization and desynchronization, different kinds of spiking dynamics can develop in a self-organized manner. Even when only a static interaction is implemented, through synchronized inputs from connected neurons, the spiking mode of photonic neurons can be spontaneously and adaptively changed between the Class-I and Class-II modes classified by A. L. Hodgkin. This spontaneous spiking-mode shift induces a significant change in the spiking frequency despite the all neurons having the same natural spiking frequency, which encourages the generation of chimera states. Controllability and self-organized flexibility of the spiking modes in the present system allow us to create an experimental platform to explore the nature of chimera states in neuromorphic spiking neural networks., Comment: 26 pages, 12 figures

S-type planets, which orbit one component of multiple-star systems, place strong constraints on the planet formation and evolution models. A notable case study is Kepler-444, a triple-star system whose primary is orbited by five planets smaller than Venus in a compact configuration, and for which the stellar binary companion revolves around the primary on a highly eccentric orbit. Having access to the most precise up-to-date masses and orbital parameters is highly valuable to understand formation and evolution processes. We provide the first full dynamical exploration of this system, with the goal to refine those parameters. The planetary system does not appear in any of low-order two or three-planet mean-motion resonances (MMR). We provide the most precise up-to-date dynamical parameters for the planets and the stellar binary companion, using an approach that makes use of the Numerical Analysis of Fundamental Frequencies (NAFF) fast chaos indicator. The orbit of the latter is constrained by new observations from HIRES and Gaia, and also by the stability analysis. This update further challenges the planets formation processes. We also test the dynamical plausibility of a sixth planet in the system, following hints observed in the Hubble Space Telescope (HST) data. We find that this putative planet could exist over a broad range of masses, and with an orbital period roughly comprised between 12 and 20 days. We note an overall good agreement of the system with short-term orbital stability. This suggests that a diverse range of planetary system architectures could be found in multiple-star systems, potentially further challenging the planet formation models., Comment: 13 pages, 9 figures. Accepted for publication in A&A

We present the discovery of two exoplanets transiting TOI-836 (TIC 440887364) using data from TESS Sector 11 and Sector 38. TOI-836 is a bright ($T = 8.5$ mag), high proper motion ($\sim\,200$ mas yr$^{-1}$), low metallicity ([Fe/H]$\approx\,-0.28$) K-dwarf with a mass of $0.68\pm0.05$ M$_{\odot}$ and a radius of $0.67\pm0.01$ R$_{\odot}$. We obtain photometric follow-up observations with a variety of facilities, and we use these data-sets to determine that the inner planet, TOI-836 b, is a $1.70\pm0.07$ R$_{\oplus}$ super-Earth in a 3.82 day orbit, placing it directly within the so-called 'radius valley'. The outer planet, TOI-836 c, is a $2.59\pm0.09$ R$_{\oplus}$ mini-Neptune in an 8.60 day orbit. Radial velocity measurements reveal that TOI-836 b has a mass of $4.5\pm0.9$ M$_{\oplus}$ , while TOI-836 c has a mass of $9.6\pm2.6$ M$_{\oplus}$. Photometric observations show Transit Timing Variations (TTVs) on the order of 20 minutes for TOI-836 c, although there are no detectable TTVs for TOI-836 b. The TTVs of planet TOI-836 c may be caused by an undetected exterior planet.

Transit Timing Variations (TTVs) can provide useful information on compact multi-planetary systems observed by transits, by putting constraints on the masses and eccentricities of the observed planets. This is especially helpful when the host star is not bright enough for radial velocity follow-up. However, in the past decades, numerous works have shown that TTV-characterised planets tend to have a lower densities than RV-characterised planets. Re-analysing 34 Kepler planets in the super-Earth to sub-Neptunes range using the RIVERS approach, we show that at least part of these discrepancies was due to the way transit timings were extracted from the light curve, which had a tendency to under-estimate the TTV amplitudes. We recover robust mass estimates (i.e. low prior dependency) for 23 of the planets. We compare these planets the RV-characterised population. A large fraction of these previously had a surprisingly low density now occupy a place of the mass-radius diagram much closer to the bulk of the known planets, although a slight shift toward lower densities remains, which could indicate that the compact multi-planetary systems characterised by TTVs are indeed composed of planets which are different from the bulk of the RV-characterised population. These results are especially important for obtaining an unbiased view of the compact multi-planetary systems detected by Kepler, TESS, and the upcoming PLATO mission.

Over the past years, the amount of detected multi-planet systems significantly grew, an important sub-class of which being the compact configurations. A precise knowledge of them is crucial to understand the conditions with which planetary systems form and evolve. However, observations often leave these systems with large uncertainties, notably on the orbital eccentricities. This is especially prominent for systems with low-mass planets detected with Radial Velocities (RV), the amount of which is more and more important in the exoplanet population. It is becoming a common approach to refine these parameters with the help of orbital stability arguments. Such dynamical techniques can be computationally expensive. In this work we use an alternative procedure faster by orders of magnitude than classical N-body integration approaches. We couple a reliable exploration of the parameter space with the precision of the Numerical Analysis of Fundamental Frequencies (NAFF, Laskar 1990) fast chaos indicator. We also propose a general procedure to calibrate the NAFF indicator on any multi-planet system without additional computational cost. This calibration strategy is illustrated on HD 45364, in addition to yet-unpublished measurements obtained with the HARPS and CORALIE high-resolution spectrographs. We validate the calibration approach on HD 202696. We test the performances of this stability-driven approach on two systems with different architectures. First we study HD 37124, a 3-planet system composed of planets in the Jovian regime. Then, we analyse HD 215152, a compact system of four low-mass planets. We demonstrate the potential of the NAFF stability-driven approach to refine the orbital parameters and planetary masses. We stress the importance of undertaking systematic global dynamical analyses on every new multi-planet system discovered., Comment: 14 pages, 14 figures. Accepted for publication in A&A

Context: TOI-2076 is a transiting three-planet system of sub-Neptunes orbiting a bright (G = 8.9 mag), young ($340\pm80$ Myr) K-type star. Although a validated planetary system, the orbits of the two outer planets were unconstrained as only two non-consecutive transits were seen in TESS photometry. This left 11 and 7 possible period aliases for each. Aims: To reveal the true orbits of these two long-period planets, precise photometry targeted on the highest-probability period aliases is required. Long-term monitoring of transits in multi-planet systems can also help constrain planetary masses through TTV measurements. Methods: We used the MonoTools package to determine which aliases to follow, and then performed space-based and ground-based photometric follow-up of TOI-2076 c and d with CHEOPS, SAINT-EX, and LCO telescopes. Results: CHEOPS observations revealed a clear detection for TOI-2076 c at $P=21.01538^{+0.00084}_{-0.00074}$ d, and allowed us to rule out three of the most likely period aliases for TOI-2076 d. Ground-based photometry further enabled us to rule out remaining aliases and confirm the $P=35.12537\pm0.00067$ d alias. These observations also improved the radius precision of all three sub-Neptunes to $2.518\pm0.036$, $3.497\pm0.043$, and $3.232\pm0.063$ $R_\oplus$. Our observations also revealed a clear anti-correlated TTV signal between planets b and c likely caused by their proximity to the 2:1 resonance, while planets c and d appear close to a 5:3 period commensurability, although model degeneracy meant we were unable to retrieve robust TTV masses. Their inflated radii, likely due to extended H-He atmospheres, combined with low insolation makes all three planets excellent candidates for future comparative transmission spectroscopy with JWST., Comment: 18 pages, 7 figure, 8 tables. Accepted for publication in Astronomy & Astrophysics. Photometry available on CDS/Vizier. Python package presented for modelling duotransiting planet candidates available at https://github.com/hposborn/MonoTools. Modelling & figure-creation code available at https://github.com/hposborn/TOI2076

Transit Timing Variations (TTVs) can provide useful information for systems observed by transit, by putting constraints on the masses and eccentricities of the observed planets, or even constrain the existence of non-transiting companions. However, TTVs can also prevent the detection of small planets in transit surveys, or bias the recovered planetary and transit parameters. Here we show that Kepler-1972 c, initially the "not transit-like" false positive KOI-3184.02, is an Earth-sized planet whose orbit is perturbed by Kepler-1972 b (initially KOI-3184.01). The pair is locked in a 3:2 Mean-motion resonance, each planet displaying TTVs of more than 6h hours of amplitude over the duration of the Kepler mission. The two planets have similar masses $m_b/m_c =0.956_{-0.051}^{+0.056}$ and radii $R_b=0.802_{-0.041}^{+0.042}R_{Earth}$, $R_c=0.868_{-0.050}^{+0.051}R_{Earth}$, and the whole system, including the inner candidate KOI-3184.03, appear to be coplanar. Despite the faintness of the signals (SNR of 1.35 for each transit of Kepler-1972 b and 1.10 for Kepler-1972 c), we recovered the transits of the planets using the RIVERS method, based on the recognition of the tracks of planets in river diagrams using machine learning, and a photo-dynamic fit of the lightcurve. Recovering the correct ephemerides of the planets is essential to have a complete picture of the observed planetary systems. In particular, we show that in Kepler-1972, not taking into account planet-planet interactions yields an error of $\sim 30\%$ on the radii of planets b and c, in addition to generating in-transit scatter, which leads to mistake KOI3184.02 for a false positive. Alleviating this bias is essential for an unbiased view of Kepler systems, some of the TESS stars, and the upcoming PLATO mission., Comment: arXiv admin note: text overlap with arXiv:2111.06825

We report the discovery and characterisation of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in TESS photometry. To characterise the system, we performed and retrieved CHEOPS, TESS, and ground-based photometry, HARPS high-resolution spectroscopy, and Gemini speckle imaging. We characterise the host star and determine $T_{\rm eff, \star}=4734\pm67$ K, $R_{\star}=0.726\pm0.007$ $R_{\odot}$, and $M_{\star}=0.748\pm0.032$ $M_{\odot}$. We present a novel detrending method based on PSF shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of $P_{\rm b}=6.44387\pm0.00003$ d, a radius of $R_{\rm b}=2.59\pm0.04$ $R_{\oplus}$, and a mass of $M_{\rm b}=13.5_{-1.8}^{+1.7}$ $M_{\oplus}$, whilst TOI-1064 c has an orbital period of $P_{\rm c}=12.22657^{+0.00005}_{-0.00004}$ d, a radius of $R_{\rm c}=2.65\pm0.04$ $R_{\oplus}$, and a 3$\sigma$ upper mass limit of 8.5 ${\rm M_{\oplus}}$. From the high-precision photometry we obtain radius uncertainties of $\sim$1.6%, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterised sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further RVs are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass-radius space, and it allows us to identify a trend in bulk density-stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets., Comment: 30 pages, 24 figures, 6 tables including the Appendix; accepted for publication in MNRAS

Transit timing variations (TTVs) can provide useful information for systems observed by transit, as they allow us to put constraints on the masses and eccentricities of the observed planets, or even to constrain the existence of non-transiting companions. However, TTVs can also act as a detection bias that can prevent the detection of small planets in transit surveys that would otherwise be detected by standard algorithms such as the Boxed Least Square algorithm (BLS) if their orbit was not perturbed. This bias is especially present for surveys with a long baseline, such as Kepler, some of the TESS sectors, and the upcoming PLATO mission. Here we introduce a detection method that is robust to large TTVs, and illustrate its use by recovering and confirming a pair of resonant super-Earths with ten-hour TTVs around Kepler-1705. The method is based on a neural network trained to recover the tracks of low-signal-to-noise-ratio(S/N) perturbed planets in river diagrams. We recover the transit parameters of these candidates by fitting the light curve. The individual transit S/N of Kepler-1705b and c are about three times lower than all the previously known planets with TTVs of 3 hours or more, pushing the boundaries in the recovery of these small, dynamically active planets. Recovering this type of object is essential for obtaining a complete picture of the observed planetary systems, and solving for a bias not often taken into account in statistical studies of exoplanet populations. In addition, TTVs are a means of obtaining mass estimates which can be essential for studying the internal structure of planets discovered by transit surveys. Finally, we show that due to the strong orbital perturbations, it is possible that the spin of the outer resonant planet of Kepler-1705 is trapped in a sub- or super-synchronous spin-orbit resonance.