Showing 1–50 of 64 results for author: Delisle, J

The ANTARESS workflow I. Optimal extraction of spatially resolved stellar spectra with high-resolution transit spectroscopy

Authors: V. Bourrier, J. -B. Delisle, C. Lovis, H. M. Cegla, M. Cretignier, R. Allart, K. Al Moulla, S. Tavella, O. Attia, D. Mounzer, V. Vaulato, M. Steiner, T. Vrignaud, S. Mercier, X. Dumusque, D. Ehrenreich, J. V. Seidel, A. Wyttenbach, W. Dethier, F. Pepe

Abstract: High-resolution spectrographs open a detailed window onto the atmospheres of stars and planets. As the number of systems observed with different instruments grows, it is crucial to develop a standard in analyzing spectral time series of exoplanet transits and occultations, for the benefit of reproducibility. Here, we introduce the ANTARESS workflow, a set of methods aimed at processing high-resolu… ▽ More High-resolution spectrographs open a detailed window onto the atmospheres of stars and planets. As the number of systems observed with different instruments grows, it is crucial to develop a standard in analyzing spectral time series of exoplanet transits and occultations, for the benefit of reproducibility. Here, we introduce the ANTARESS workflow, a set of methods aimed at processing high-resolution spectroscopy datasets in a robust way and extracting accurate exoplanetary and stellar spectra. While a fast preliminary analysis can be run on order-merged 1D spectra and cross-correlation functions (CCFs), the workflow was optimally designed for extracted 2D echelle spectra to remain close to the original detector counts, limit the spectral resampling, and propagate the correlated noise. Input data from multiple instruments and epochs were corrected for relevant environmental and instrumental effects, processed homogeneously, and analyzed independently or jointly. In this first paper, we show how planet-occulted stellar spectra extracted along the transit chord and cleaned from planetary contamination provide a direct comparison with theoretical stellar models and enable a spectral and spatial mapping of the photosphere. We illustrate this application of the workflow to archival ESPRESSO data, using the Rossiter-McLaughlin effect Revolutions (RMR) technique to confirm the spin-orbit alignment of HD\,209458b and unveil biases in WASP-76b's published orbital architecture. Because the workflow is modular and its concepts are general, it can support new methods and be extended to additional spectrographs to find a range of applications beyond the proposed scope. In a companion paper, we will present how planet-occulted spectra can be processed further to extract and analyze planetary spectra decontaminated from the star, providing clean and direct measurements of atmospheric properties. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 30 pages and 30 figures (plus Appendix). Accepted for publication in A&A

Resonant sub-Neptunes are puffier

Authors: Adrien Leleu, Jean-Baptiste Delisle, Remo Burn, André Izidoro, Stéphane Udry, Xavier Dumusque, Christophe Lovis, Sarah Millholland, Léna Parc, François Bouchy, Vincent Bourrier, Yann Alibert, João Faria, Christoph Mordasini, Damien Ségransan

Abstract: 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 resona… ▽ More 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. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Photo-dynamical characterisation of the TOI-178 resonant chain

Authors: A. Leleu, J. -B. Delisle, L. Delrez, E. M. Bryant, A. Brandeker, H. P. Osborn, N. Hara, T. G. Wilson, N. Billot, M. Lendl, D. Ehrenreich, H. Chakraborty, M. N. Günther, M. J. Hooton, Y. Alibert, R. Alonso, D. R. Alves, D. R. Anderson, I. Apergis, D. Armstrong, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, M. P. Battley, W. Baumjohann , et al. (82 additional authors not shown)

Abstract: 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 ev… ▽ More 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. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

A resonant sextuplet of sub-Neptunes transiting the bright star HD 110067

Authors: R. Luque, H. P. Osborn, A. Leleu, E. Pallé, A. Bonfanti, O. Barragán, T. G. Wilson, C. Broeg, A. Collier Cameron, M. Lendl, P. F. L. Maxted, Y. Alibert, D. Gandolfi, J. -B. Delisle, M. J. Hooton, J. A. Egger, G. Nowak, M. Lafarga, D. Rapetti, J. D. Twicken, J. C. Morales, I. Carleo, J. Orell-Miquel, V. Adibekyan, R. Alonso , et al. (127 additional authors not shown)

Abstract: Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial con… ▽ More Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial conditions and environment. Those in resonance (with their orbital periods related by a ratio of small integers) are particularly valuable because they imply a system architecture practically unchanged since its birth. Here, we present the observations of six transiting planets around the bright nearby star HD 110067. We find that the planets follow a chain of resonant orbits. A dynamical study of the innermost planet triplet allowed the prediction and later confirmation of the orbits of the rest of the planets in the system. The six planets are found to be sub-Neptunes with radii ranging from 1.94 to 2.85 Re. Three of the planets have measured masses, yielding low bulk densities that suggest the presence of large hydrogen-dominated atmospheres. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Published in Nature on November 30, 2023. Supplementary Information can be found in the online version of the paper in the journal

Journal ref: Nature 623, 932-937 (2023)

Exploring the Brown Dwarf Desert with Precision Radial Velocities and Gaia DR3 Astrometric Orbits

Authors: N. Unger, D. Ségransan, D. Barbato, J. -B. Delisle, J. Sahlmann, B. Holl, S. Udry

Abstract: Context. The observed scarcity of brown dwarfs in close orbits (within 10 au) around solar-type stars poses significant questions about the origins of these substellar companions. These questions impact our broader understanding of planetary formation processes. However, to resolve these formation mechanisms, accurate observational constraints are essential. Most of the brown dwarfs have been disc… ▽ More Context. The observed scarcity of brown dwarfs in close orbits (within 10 au) around solar-type stars poses significant questions about the origins of these substellar companions. These questions impact our broader understanding of planetary formation processes. However, to resolve these formation mechanisms, accurate observational constraints are essential. Most of the brown dwarfs have been discovered by radial velocity surveys, but this method introduces uncertainties due to its inability to determine the orbital inclination, leaving the true mass-and thus their true nature-unresolved. This highlights the crucial role of astrometric data, helping us distinguish between genuine brown dwarfs and stars. Aims. We aim to refine the mass estimates of massive companions to solar-type stars, mostly discovered through radial velocity measurements and subsequently validated using Gaia DR3 astrometry, to gain a clearer understanding of their true mass and occurrence rates. Methods. We selected a sample of 31 sources with substellar companion candidates validated by Gaia DR3 and with available radial velocities. Using the Gaia DR3 solutions as prior information, we performed an MCMC fit with the available radial velocity measurements to integrate these two sources of data and thus obtain an estimate of their true mass. Results. Combining radial velocity measurements with Gaia DR3 data led to more precise mass estimations, leading us to reclassify several systems initially labeled as brown dwarfs as low-mass stars. Out of the 32 analyzed companions, 13 are determined to be stars, 17 are sub-stellar, and 2 have inconclusive results with the current data. Importantly, using these updated masses, we reevaluated the occurrence rate of brown dwarf companions (13-80 M$_{jup}$ ) on close orbits (<10 au) in the CORALIE sample, determining a tentative occurrence rate of $0.8^{+0.3}_{-0.2}\%$ △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Accepted by A&A

Journal ref: A&A 680, A16 (2023)

A review of planetary systems around HD 99492, HD 147379 and HD 190007 with HARPS-N

Authors: M. Stalport, M. Cretignier, S. Udry, A. Anna John, T. G. Wilson, J. -B. Delisle, A. S. Bonomo, L. A. Buchhave, D. Charbonneau, S. Dalal, M. Damasso, L. Di Fabrizio, X. Dumusque, A. Fiorenzano, A. Harutyunyan, R. D. Haywood, D. W. Latham, M. López-Morales, V. Lorenzi, C. Lovis, L. Malavolta, E. Molinari, A. Mortier, M. Pedani, F. Pepe , et al. (4 additional authors not shown)

Abstract: The Rocky Planet Search (RPS) program is dedicated to a blind radial velocity (RV) search of planets around bright stars in the Northern hemisphere, using the high-resolution echelle spectrograph HARPS-N installed on the Telescopio Nazionale Galileo (TNG). The goal of this work is to revise and update the properties of three planetary systems by analysing the HARPS-N data with state-of-the-art s… ▽ More The Rocky Planet Search (RPS) program is dedicated to a blind radial velocity (RV) search of planets around bright stars in the Northern hemisphere, using the high-resolution echelle spectrograph HARPS-N installed on the Telescopio Nazionale Galileo (TNG). The goal of this work is to revise and update the properties of three planetary systems by analysing the HARPS-N data with state-of-the-art stellar activity mitigation tools. The stars considered are HD 99492 (83Leo B), HD 147379 (Gl617 A) and HD 190007. We employ a systematic process of data modelling, that we selected from the comparison of different approaches. We use YARARA to remove instrumental systematics from the RV, and then use SPLEAF to further mitigate the stellar noise with a multidimensional correlated noise model. We also search for transit features in the Transiting Exoplanets Survey Satellite (TESS) data of these stars. We report on the discovery of a new planet around HD 99492, namely HD 99492 c, with an orbital period of 95.2 days and a minimum mass of msin i = 17.9 M_Earth, and refine the parameters of HD 99492 b. We also update and refine the Keplerian solutions for the planets around HD 147379 and HD 190007, but do not detect additional planetary signals. We discard the transiting geometry for the planets, but stress that TESS did not exhaustively cover all the orbital phases. The addition of the HARPS-N data, and the use of advanced data analysis tools, has allowed us to present a more precise view of these three planetary systems. It demonstrates once again the importance of long observational efforts such as the RPS program. Added to the RV exoplanet sample, these planets populate two apparently distinct populations revealed by a bimodality in the planets minimum mass distribution. The separation is located between 30 and 50 M_Earth. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 27 pages, 26 figures (13 in Appendix); Accepted for publication in A&A

Journal ref: A&A 678, A90 (2023)

DREAM II. The spin-orbit angle distribution of close-in exoplanets under the lens of tides

Authors: O. Attia, V. Bourrier, J. -B. Delisle, P. Eggenberger

Abstract: The spin-orbit angle, or obliquity, is a powerful observational marker that allows us to access the dynamical history of exoplanetary systems. Here, we have examined the distribution of spin-orbit angles for close-in exoplanets and put it in a statistical context of tidal interactions between planets and their stars. We confirm the observed trends between the obliquity and physical quantities dire… ▽ More The spin-orbit angle, or obliquity, is a powerful observational marker that allows us to access the dynamical history of exoplanetary systems. Here, we have examined the distribution of spin-orbit angles for close-in exoplanets and put it in a statistical context of tidal interactions between planets and their stars. We confirm the observed trends between the obliquity and physical quantities directly connected to tides, namely the stellar effective temperature, the planet-to-star mass ratio, and the scaled orbital distance. We further devised a tidal efficiency factor combining critical parameters that control the strength of tidal effects and used it to corroborate the strong link between the spin-orbit angle distribution and tidal interactions. In particular, we developed a readily usable formula to estimate the probability that a system is misaligned, which will prove useful in global population studies. By building a robust statistical framework, we reconstructed the distribution of the three-dimensional spin-orbit angles, allowing for a sample of nearly 200 true obliquities to be analyzed for the first time. This realistic distribution maintains the sky-projected trends, and additionally hints toward a striking pileup of truly aligned systems. The comparison between the full population and a pristine subsample unaffected by tidal interactions suggests that perpendicular architectures are resilient toward tidal realignment, providing evidence that orbital misalignments are sculpted by disruptive dynamical processes that preferentially lead to polar orbits. On the other hand, star-planet interactions seem to efficiently realign or quench the formation of any tilted configuration other than for polar orbits, and in particular for antialigned orbits. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: Accepted in A&A

Journal ref: A&A 674, A120 (2023)

A statistical model of stellar variability. I. FENRIR: a physics-based model of stellar activity, and its fast Gaussian process approximation

Authors: Nathan C. Hara, Jean-Baptiste Delisle

Abstract: The detection of terrestrial planets by radial velocity and photometry is hindered by the presence of stellar signals. Those are often modeled as stationary Gaussian processes, whose kernels are based on qualitative considerations, which do not fully leverage the existing physical understanding of stars. Our aim is to build a formalism which allows to transfer the knowledge of stellar activity int… ▽ More The detection of terrestrial planets by radial velocity and photometry is hindered by the presence of stellar signals. Those are often modeled as stationary Gaussian processes, whose kernels are based on qualitative considerations, which do not fully leverage the existing physical understanding of stars. Our aim is to build a formalism which allows to transfer the knowledge of stellar activity into practical data analysis methods. In particular, we aim at obtaining kernels with physical parameters. This has two purposes: better modelling signals of stellar origin to find smaller exoplanets, and extracting information about the star from the statistical properties of the data. We consider several observational channels such as photometry, radial velocity, activity indicators, and build a model called FENRIR to represent their stochastic variations due to stellar surface inhomogeneities. We compute analytically the covariance of this multi-channel stochastic process, and implement it in the S+LEAF framework to reduce the cost of likelihood evaluations from $O(N^3)$ to $O(N)$. We also compute analytically higher order cumulants of our FENRIR model, which quantify its non-Gaussianity. We obtain a fast Gaussian process framework with physical parameters, which we apply to the HARPS-N and SORCE observations of the Sun, and constrain a solar inclination compatible with the viewing geometry. We then discuss the application of our formalism to granulation. We exhibit non-Gaussianity in solar HARPS radial velocities, and argue that information is lost when stellar activity signals are assumed to be Gaussian. We finally discuss the origin of phase shifts between RVs and indicators, and how to build relevant activity indicators. We provide an open-source implementation of the FENRIR Gaussian process model with a Python interface. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: Submitted to Astronomy \& Astrophysics

The CORALIE survey for southern extrasolar planets XIX. Brown dwarfs and stellar companions unveiled by radial velocity and astrometry

Authors: D. Barbato, D. Ségransan, S. Udry, N. Unger, F. Bouchy, C. Lovis, M. Mayor, F. Pepe, D. Queloz, N. C. Santos, J. B. Delisle, P. Figueira, M. Marmier, E. C. Matthews, G. Lo Curto, J. Venturini, G. Chaverot, M. Cretignier, J. F. Otegi, M. Stalport

Abstract: A historical planet-search on a sample of 1647 nearby southern main sequence stars has been ongoing since 1998 with the CORALIE spectrograph at La Silla Observatory, with a backup subprogram dedicated to the monitoring of binary stars. We review 25 years of CORALIE measurements and search for Doppler signals consistent with stellar or brown dwarf companions to produce an updated catalog of both kn… ▽ More A historical planet-search on a sample of 1647 nearby southern main sequence stars has been ongoing since 1998 with the CORALIE spectrograph at La Silla Observatory, with a backup subprogram dedicated to the monitoring of binary stars. We review 25 years of CORALIE measurements and search for Doppler signals consistent with stellar or brown dwarf companions to produce an updated catalog of both known and previously unpublished binary stars in the planet-search sample, assessing the binarity fraction of the stellar population and providing perspective for more precise planet-search in the binary sample. We perform new analysis on the CORALIE planet-search sample radial velocity measurements, searching for stellar companions and obtaining orbital solutions for both known and new binary systems. We perform simultaneous radial velocity and proper motion anomaly fits on the subset of these systems for which Hipparcos and Gaia astrometry measurements are available, obtaining accurate estimates of true mass for the companions. We find 218 stars in the CORALIE sample to have at least one stellar companion, 130 of which are not yet published in the literature and for which we present orbital solutions. The use of proper motion anomaly allow us to derive true masses for the stellar companions in 132 systems, which we additionally use to estimate stability regions for possible planetary companions on circumprimary or circumbinary orbits. Finally, we produce detection limit maps for each star in the sample and obtain occurrence rates of $0.43^{+0.23}_{-0.11}\%$ and $12.69^{+0.87}_{-0.77}\%$ for brown dwarf and stellar companions respectively in the CORALIE sample. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 34 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 674, A114 (2023)

Two temperate Earth-mass planets orbiting the nearby star GJ1002

Authors: A. Suárez Mascareño, E. González-Álvarez, M. R. Zapatero Osorio, J. Lillo-Box, J. P. Faria, V. M. Passegger, J. I. González Hernández, P. Figueira, A. Sozzetti, R. Rebolo, F. Pepe, N. C. Santos, S. Cristiani, C. Lovis, A. M. Silva, I. Ribas, P. J. Amado, J. A. Caballero, A. Quirrenbach, A. Reiners, M. Zechmeister, V. Adibekyan, Y. Alibert, V. J. S. Béjar, S. Benatti , et al. (20 additional authors not shown)

Abstract: We report the discovery and characterisation of two Earth-mass planets orbiting in the habitable zone of the nearby M-dwarf GJ~1002 based on the analysis of the radial-velocity (RV) time series from the ESPRESSO and CARMENES spectrographs. The host star is the quiet M5.5~V star GJ~1002 (relatively faint in the optical, $V \sim 13.8$ mag, but brighter in the infrared, $J \sim 8.3$ mag), located at… ▽ More We report the discovery and characterisation of two Earth-mass planets orbiting in the habitable zone of the nearby M-dwarf GJ~1002 based on the analysis of the radial-velocity (RV) time series from the ESPRESSO and CARMENES spectrographs. The host star is the quiet M5.5~V star GJ~1002 (relatively faint in the optical, $V \sim 13.8$ mag, but brighter in the infrared, $J \sim 8.3$ mag), located at 4.84 pc from the Sun. We analyse 139 spectroscopic observations taken between 2017 and 2021. We performed a joint analysis of the time series of the RV and full-width half maximum (FWHM) of the cross-correlation function (CCF) to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with the stellar activity. We detect the signal of two planets orbiting GJ~1002. GJ~1002~b is a planet with a minimum mass $m_p \sin i $ of 1.08 $\pm$ 0.13 M$_{\oplus}$ with an orbital period of 10.3465 $\pm$ 0.0027 days at a distance of 0.0457 $\pm$ 0.0013 au from its parent star, receiving an estimated stellar flux of 0.67 $F_{\oplus}$. GJ~1002 c is a planet with a minimum mass $m_p \sin i $ of 1.36 $\pm$ 0.17 M$_{\oplus}$ with an orbital period of 21.202 $\pm$ 0.013 days at a distance of 0.0738 $\pm$ 0.0021 au from its parent star, receiving an estimated stellar flux of 0.257 $F_{\oplus}$. We also detect the rotation signature of the star, with a period of 126 $\pm$ 15 days. GJ~1002 is one of the few known nearby systems with planets that could potentially host habitable environments. The closeness of the host star to the Sun makes the angular sizes of the orbits of both planets ($\sim$ 9.7 mas and $\sim$ 15.7 mas, respectively) large enough for their atmosphere to be studied via high-contrast high-resolution spectroscopy with instruments such as the future spectrograph ANDES for the ELT or the LIFE mission. △ Less

Submitted 25 January, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 23 pages, 20 figures, 3 tables ; Replacement to fix some typos

L 363-38 b: a planet newly discovered with ESPRESSO orbiting a nearby M dwarf star

Authors: Lia F. Sartori, Christophe Lovis, Jean-Baptiste Delisle, Monika Lendl, Gabriele Cugno, Anna Boehle, Felix Dannert, Andrea Krenn, Jonas L. Gubler, Sascha P. Quanz

Abstract: Context. Planets around stars in the solar neighbourhood will be prime targets for characterisation with upcoming large space- and ground-based facilities. Since large-scale exoplanet searches will not be feasible with such telescopes, it is crucial to use currently available data and instruments to find possible target planets before next generation facilities come online. Aims. We aim at detec… ▽ More Context. Planets around stars in the solar neighbourhood will be prime targets for characterisation with upcoming large space- and ground-based facilities. Since large-scale exoplanet searches will not be feasible with such telescopes, it is crucial to use currently available data and instruments to find possible target planets before next generation facilities come online. Aims. We aim at detecting new extrasolar planets around stars in the solar neighbourhood by blind radial velocity (RV) search with ESPRESSO. Our target sample consist of nearby stars (d < 11 pc) with little (< 10) or no previous RV measurements. Methods. We use 31 radial velocity measurements obtained with ESPRESSO at the VLT between December 2020 and February 2022 of the nearby M dwarf star (M_star = 0.21 M_sun, d = 10.23 pc) L 363-38 to derive the orbital parameters of the newly discovered planet. In addition, we use TESS photometry and archival VLT/NaCo high contrast imaging data to put further constraints on the orbit inclination and the possible planetary system architecture around L 363-38. Results. We present the detection of a new extrasolar planet orbiting the nearby M dwarf star L 363-38. L 363-38 b is a planet with minimum mass mp sin(i) = 4.67+/-0.43 M_Earth orbiting its star with a period P = 8.781+/-0.007 d, corresponding to a semi-major axis a = 0.048+/-0.006 AU, which is well inside the inner edge of the habitable zone. We further estimate a minimum radius rp sin(i) = 1.55 - 2.75 R_Earth and an equilibrium temperature Teq = 330K. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 12 pages, 11 figures

Journal ref: A&A 670, A42 (2023)

Analytical determination of orbital elements using Fourier analysis. II. Gaia astrometry and its combination with radial velocities

Authors: J. -B. Delisle, D. Ségransan

Abstract: The ESA global astrometry space mission Gaia has been monitoring the position of a billion stars since 2014. The analysis of such a massive dataset is challenging in terms of the data processing involved. In particular, the blind detection and characterization of single or multiple companions to stars (planets, brown dwarfs, or stars) using Gaia astrometry requires highly efficient algorithms. In… ▽ More The ESA global astrometry space mission Gaia has been monitoring the position of a billion stars since 2014. The analysis of such a massive dataset is challenging in terms of the data processing involved. In particular, the blind detection and characterization of single or multiple companions to stars (planets, brown dwarfs, or stars) using Gaia astrometry requires highly efficient algorithms. In this article, we present a set of analytical methods to detect and characterize companions in scanning space astrometric time series as well as via a combination of astrometric and radial velocity time series. We propose a general linear periodogram framework and we derive analytical formulas for the false alarm probability (FAP) of periodogram peaks. Once a significant peak has been identified, we provide analytical estimates of all the orbital elements of the companion based on the Fourier decomposition of the signal. The periodogram, FAP, and orbital elements estimates can be computed for the astrometric and radial velocity time series separately or in tandem. These methods are complementary with more accurate and more computationally intensive numerical algorithms (e.g., least-squares minimization, Markov chain Monte Carlo, genetic algorithms). In particular, our analytical approximations can be used as an initial condition to accelerate the convergence of numerical algorithms. Our formalism has been partially implemented in the Gaia exoplanet pipeline for the third Gaia data release. Since the Gaia astrometric time series are not yet publicly available, we illustrate our methods on the basis of Hipparcos data, together with on-ground CORALIE radial velocities, for three targets known to host a companion: HD 223636 (HIP 117622), HD 17289 (HIP 12726), and HD 3277 (HIP 2790). △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: Accepted in A&A

Journal ref: A&A 667, A172 (2022)

A CHEOPS-enhanced view of the HD3167 system

Authors: V. Bourrier, A. Deline, A. Krenn, J. A. Egger, A. C. Petit, L. Malavolta, M. Cretignier, N. Billot, C. Broeg, H. -G. Florén, D. Queloz, Y. Alibert, A. Bonfanti, A. S. Bonomo, J. -B. Delisle, O. D. S. Demangeon, B. -O. Demory, X. Dumusque, D. Ehrenreich, R. D. Haywood, S. B Howell, M. Lendl, A. Mortier, G. Nigro, S. Salmon , et al. (70 additional authors not shown)

Abstract: Much remains to be understood about the nature of exoplanets smaller than Neptune, most of which have been discovered in compact multi-planet systems. With its inner ultra-short period planet b aligned with the star and two larger outer planets d-c on polar orbits, the multi-planet system HD 3167 features a peculiar architecture and offers the possibility to investigate both dynamical and atmosphe… ▽ More Much remains to be understood about the nature of exoplanets smaller than Neptune, most of which have been discovered in compact multi-planet systems. With its inner ultra-short period planet b aligned with the star and two larger outer planets d-c on polar orbits, the multi-planet system HD 3167 features a peculiar architecture and offers the possibility to investigate both dynamical and atmospheric evolution processes. To this purpose we combined multiple datasets of transit photometry and radial velocimetry (RV) to revise the properties of the system and inform models of its planets. This effort was spearheaded by CHEOPS observations of HD 3167b, which appear inconsistent with a purely rocky composition despite its extreme irradiation. Overall the precision on the planetary orbital periods are improved by an order of magnitude, and the uncertainties on the densities of the transiting planets b and c are decreased by a factor of 3. Internal structure and atmospheric simulations draw a contrasting picture between HD 3167d, likely a rocky super-Earth that lost its atmosphere through photo-evaporation, and HD 3167c, a mini-Neptune that kept a substantial primordial gaseous envelope. We detect a fourth, more massive planet on a larger orbit, likely coplanar with HD 3167d-c. Dynamical simulations indeed show that the outer planetary system d-c-e was tilted, as a whole, early in the system history, when HD 3167b was still dominated by the star influence and maintained its aligned orbit. RV data and direct imaging rule out that the companion that could be responsible for the present-day architecture is still bound to the HD\,3167 system. Similar global studies of multi-planet systems will tell how many share the peculiar properties of the HD3167 system, which remains a target of choice for follow-up observations and simulations. △ Less

Submitted 19 September, 2022; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 22 pages, 23 pages, accepted for publication in A&A (18 August 2022). Updated author list in new version

Global dynamics and architecture of the Kepler-444 system

Authors: M. Stalport, E. C. Matthews, V. Bourrier, A. Leleu, J. -B. Delisle, S. Udry

Abstract: 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… ▽ More 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. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 13 pages, 9 figures. Accepted for publication in A&A

Journal ref: A&A 667, A128 (2022)

Gaia Data Release 3: Summary of the content and survey properties

Authors: Gaia Collaboration, A. Vallenari, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran , et al. (431 additional authors not shown)

Abstract: We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photom… ▽ More We present the third data release of the European Space Agency's Gaia mission, GDR3. The GDR3 catalogue is the outcome of the processing of raw data collected with the Gaia instruments during the first 34 months of the mission by the Gaia Data Processing and Analysis Consortium. The GDR3 catalogue contains the same source list, celestial positions, proper motions, parallaxes, and broad band photometry in the G, G$_{BP}$, and G$_{RP}$ pass-bands already present in the Early Third Data Release. GDR3 introduces an impressive wealth of new data products. More than 33 million objects in the ranges $G_{rvs} < 14$ and $3100 <T_{eff} <14500 $, have new determinations of their mean radial velocities based on data collected by Gaia. We provide G$_{rvs}$ magnitudes for most sources with radial velocities, and a line broadening parameter is listed for a subset of these. Mean Gaia spectra are made available to the community. The GDR3 catalogue includes about 1 million mean spectra from the radial velocity spectrometer, and about 220 million low-resolution blue and red prism photometer BPRP mean spectra. The results of the analysis of epoch photometry are provided for some 10 million sources across 24 variability types. GDR3 includes astrophysical parameters and source class probabilities for about 470 million and 1500 million sources, respectively, including stars, galaxies, and quasars. Orbital elements and trend parameters are provided for some $800\,000$ astrometric, spectroscopic and eclipsing binaries. More than $150\,000$ Solar System objects, including new discoveries, with preliminary orbital solutions and individual epoch observations are part of this release. Reflectance spectra derived from the epoch BPRP spectral data are published for about 60\,000 asteroids. Finally, an additional data set is provided, namely the Gaia Andromeda Photometric Survey (abridged) △ Less

Submitted 30 July, 2022; originally announced August 2022.

Comments: 23 pages, 2 figures

Unbiasing the density of TTV-characterised sub-Neptunes: Update of the mass-radius relationship of 34 Kepler planets

Authors: A. Leleu, J. -B. Delisle, S. Udry, R. Mardling, M. Turbet, J. A. Egger, Y. Alibert, G. Chatel, P. Eggenberger, M. Stalport

Abstract: 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… ▽ More 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. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Gaia Data Release 3: Reflectance spectra of Solar System small bodies

Authors: Gaia Collaboration, L. Galluccio, M. Delbo, F. De Angeli, T. Pauwels, P. Tanga, F. Mignard, A. Cellino, A. G. A. Brown, K. Muinonen, A. Penttila, S. Jordan, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, L. Eyer, R. Guerra, A. Hutton, C. Jordi , et al. (422 additional authors not shown)

Abstract: The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was deriv… ▽ More The Gaia mission of the European Space Agency (ESA) has been routinely observing Solar System objects (SSOs) since the beginning of its operations in August 2014. The Gaia data release three (DR3) includes, for the first time, the mean reflectance spectra of a selected sample of 60 518 SSOs, primarily asteroids, observed between August 5, 2014, and May 28, 2017. Each reflectance spectrum was derived from measurements obtained by means of the Blue and Red photometers (BP/RP), which were binned in 16 discrete wavelength bands. We describe the processing of the Gaia spectral data of SSOs, explaining both the criteria used to select the subset of asteroid spectra published in Gaia DR3, and the different steps of our internal validation procedures. In order to further assess the quality of Gaia SSO reflectance spectra, we carried out external validation against SSO reflectance spectra obtained from ground-based and space-borne telescopes and available in the literature. For each selected SSO, an epoch reflectance was computed by dividing the calibrated spectrum observed by the BP/RP at each transit on the focal plane by the mean spectrum of a solar analogue. The latter was obtained by averaging the Gaia spectral measurements of a selected sample of stars known to have very similar spectra to that of the Sun. Finally, a mean of the epoch reflectance spectra was calculated in 16 spectral bands for each SSO. The agreement between Gaia mean reflectance spectra and those available in the literature is good for bright SSOs, regardless of their taxonomic spectral class. We identify an increase in the spectral slope of S-type SSOs with increasing phase angle. Moreover, we show that the spectral slope increases and the depth of the 1 um absorption band decreases for increasing ages of S-type asteroid families. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 30 pages, 26 figures

Gaia Data Release 3: Mapping the asymmetric disc of the Milky Way

Authors: Gaia Collaboration, R. Drimmel, M. Romero-Gomez, L. Chemin, P. Ramos, E. Poggio, V. Ripepi, R. Andrae, R. Blomme, T. Cantat-Gaudin, A. Castro-Ginard, G. Clementini, F. Figueras, M. Fouesneau, Y. Fremat, K. Jardine, S. Khanna, A. Lobel, D. J. Marshall, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou , et al. (431 additional authors not shown)

Abstract: With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provid… ▽ More With the most recent Gaia data release the number of sources with complete 6D phase space information (position and velocity) has increased to well over 33 million stars, while stellar astrophysical parameters are provided for more than 470 million sources, in addition to the identification of over 11 million variable stars. Using the astrophysical parameters and variability classifications provided in Gaia DR3, we select various stellar populations to explore and identify non-axisymmetric features in the disc of the Milky Way in both configuration and velocity space. Using more about 580 thousand sources identified as hot OB stars, together with 988 known open clusters younger than 100 million years, we map the spiral structure associated with star formation 4-5 kpc from the Sun. We select over 2800 Classical Cepheids younger than 200 million years, which show spiral features extending as far as 10 kpc from the Sun in the outer disc. We also identify more than 8.7 million sources on the red giant branch (RGB), of which 5.7 million have line-of-sight velocities, allowing the velocity field of the Milky Way to be mapped as far as 8 kpc from the Sun, including the inner disc. The spiral structure revealed by the young populations is consistent with recent results using Gaia EDR3 astrometry and source lists based on near infrared photometry, showing the Local (Orion) arm to be at least 8 kpc long, and an outer arm consistent with what is seen in HI surveys, which seems to be a continuation of the Perseus arm into the third quadrant. Meanwhile, the subset of RGB stars with velocities clearly reveals the large scale kinematic signature of the bar in the inner disc, as well as evidence of streaming motions in the outer disc that might be associated with spiral arms or bar resonances. (abridged) △ Less

Submitted 5 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 35 pages, 27 figures, accepted for publication in A&A special Gaia DR3 issue. V2: abstract completed. V3: complete author list and link to data: https://drive.google.com/drive/u/1/folders/1yOJPjYmM7QK5XVsqaiSOTuwDQNti2LlZ

Journal ref: A&A 674, A37 (2023)

Gaia Data Release 3: Pulsations in main sequence OBAF-type stars

Authors: Gaia Collaboration, J. De Ridder, V. Ripepi, C. Aerts, L. Palaversa, L. Eyer, B. Holl, M. Audard, L. Rimoldini, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux, M. Biermann, O. L. Creevey, C. Ducourant, D. W. Evans, R. Guerra, A. Hutton, C. Jordi, S. A. Klioner, U. L. Lammers, L. Lindegren , et al. (423 additional authors not shown)

Abstract: The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), del… ▽ More The third Gaia data release provides photometric time series covering 34 months for about 10 million stars. For many of those stars, a characterisation in Fourier space and their variability classification are also provided. This paper focuses on intermediate- to high-mass (IHM) main sequence pulsators M >= 1.3 Msun) of spectral types O, B, A, or F, known as beta Cep, slowly pulsating B (SPB), delta Sct, and gamma Dor stars. These stars are often multi-periodic and display low amplitudes, making them challenging targets to analyse with sparse time series. All datasets used in this analysis are part of the Gaia DR3 data release. The photometric time series were used to perform a Fourier analysis, while the global astrophysical parameters necessary for the empirical instability strips were taken from the Gaia DR3 gspphot tables, and the vsini data were taken from the Gaia DR3 esphs tables. We show that for nearby OBAF-type pulsators, the Gaia DR3 data are precise and accurate enough to pinpoint them in the Hertzsprung-Russell diagram. We find empirical instability strips covering broader regions than theoretically predicted. In particular, our study reveals the presence of fast rotating gravity-mode pulsators outside the strips, as well as the co-existence of rotationally modulated variables inside the strips as reported before in the literature. We derive an extensive period-luminosity relation for delta Sct stars and provide evidence that the relation features different regimes depending on the oscillation period. Finally, we demonstrate how stellar rotation attenuates the amplitude of the dominant oscillation mode of delta Sct stars. △ Less

Submitted 16 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Journal ref: A&A 674, A36 (2023)

Gaia Data Release 3: A Golden Sample of Astrophysical Parameters

Authors: Gaia Collaboration, O. L. Creevey, L. M. Sarro, A. Lobel, E. Pancino, R. Andrae, R. L. Smart, G. Clementini, U. Heiter, A. J. Korn, M. Fouesneau, Y. Frémat, F. De Angeli, A. Vallenari, D. L. Harrison, F. Thévenin, C. Reylé, R. Sordo, A. Garofalo, A. G. A. Brown, L. Eyer, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux , et al. (423 additional authors not shown)

Abstract: Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples… ▽ More Gaia Data Release 3 (DR3) provides a wealth of new data products for the astronomical community to exploit, including astrophysical parameters for a half billion stars. In this work we demonstrate the high quality of these data products and illustrate their use in different astrophysical contexts. We query the astrophysical parameter tables along with other tables in Gaia DR3 to derive the samples of the stars of interest. We validate our results by using the Gaia catalogue itself and by comparison with external data. We have produced six homogeneous samples of stars with high quality astrophysical parameters across the HR diagram for the community to exploit. We first focus on three samples that span a large parameter space: young massive disk stars (~3M), FGKM spectral type stars (~3M), and UCDs (~20K). We provide these sources along with additional information (either a flag or complementary parameters) as tables that are made available in the Gaia archive. We furthermore identify 15740 bone fide carbon stars, 5863 solar-analogues, and provide the first homogeneous set of stellar parameters of the Spectro Photometric Standard Stars. We use a subset of the OBA sample to illustrate its usefulness to analyse the Milky Way rotation curve. We then use the properties of the FGKM stars to analyse known exoplanet systems. We also analyse the ages of some unseen UCD-companions to the FGKM stars. We additionally predict the colours of the Sun in various passbands (Gaia, 2MASS, WISE) using the solar-analogue sample. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 35 pages, (incl 6 pages references, acknowledgements, affiliations), 37 figures, A&A accepted

Journal ref: A&A 674, A39 (2023)

Gaia Data Release 3: The extragalactic content

Authors: Gaia Collaboration, C. A. L. Bailer-Jones, D. Teyssier, L. Delchambre, C. Ducourant, D. Garabato, D. Hatzidimitriou, S. A. Klioner, L. Rimoldini, I. Bellas-Velidis, R. Carballo, M. I. Carnerero, C. Diener, M. Fouesneau, L. Galluccio, P. Gavras, A. Krone-Martins, C. M. Raiteri, R. Teixeira, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, F. Arenou, C. Babusiaux , et al. (422 additional authors not shown)

Abstract: The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G~21 mag. Due to the nature of the Gaia onboard selection algorithms, these are mostly point-source-like objects. Using data prov… ▽ More The Gaia Galactic survey mission is designed and optimized to obtain astrometry, photometry, and spectroscopy of nearly two billion stars in our Galaxy. Yet as an all-sky multi-epoch survey, Gaia also observes several million extragalactic objects down to a magnitude of G~21 mag. Due to the nature of the Gaia onboard selection algorithms, these are mostly point-source-like objects. Using data provided by the satellite, we have identified quasar and galaxy candidates via supervised machine learning methods, and estimate their redshifts using the low resolution BP/RP spectra. We further characterise the surface brightness profiles of host galaxies of quasars and of galaxies from pre-defined input lists. Here we give an overview of the processing of extragalactic objects, describe the data products in Gaia DR3, and analyse their properties. Two integrated tables contain the main results for a high completeness, but low purity (50-70%), set of 6.6 million candidate quasars and 4.8 million candidate galaxies. We provide queries that select purer sub-samples of these containing 1.9 million probable quasars and 2.9 million probable galaxies (both 95% purity). We also use high quality BP/RP spectra of 43 thousand high probability quasars over the redshift range 0.05-4.36 to construct a composite quasar spectrum spanning restframe wavelengths from 72-100 nm. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: Accepted to A&A

Gaia Data Release 3: Stellar multiplicity, a teaser for the hidden treasure

Authors: Gaia Collaboration, F. Arenou, C. Babusiaux, M. A. Barstow, S. Faigler, A. Jorissen, P. Kervella, T. Mazeh, N. Mowlavi, P. Panuzzo, J. Sahlmann, S. Shahaf, A. Sozzetti, N. Bauchet, Y. Damerdji, P. Gavras, P. Giacobbe, E. Gosset, J. -L. Halbwachs, B. Holl, M. G. Lattanzi, N. Leclerc, T. Morel, D. Pourbaix, P. Re Fiorentin , et al. (425 additional authors not shown)

Abstract: The Gaia DR3 Catalogue contains for the first time about eight hundred thousand solutions with either orbital elements or trend parameters for astrometric, spectroscopic and eclipsing binaries, and combinations of them. This paper aims to illustrate the huge potential of this large non-single star catalogue. Using the orbital solutions together with models of the binaries, a catalogue of tens of t… ▽ More The Gaia DR3 Catalogue contains for the first time about eight hundred thousand solutions with either orbital elements or trend parameters for astrometric, spectroscopic and eclipsing binaries, and combinations of them. This paper aims to illustrate the huge potential of this large non-single star catalogue. Using the orbital solutions together with models of the binaries, a catalogue of tens of thousands of stellar masses, or lower limits, partly together with consistent flux ratios, has been built. Properties concerning the completeness of the binary catalogues are discussed, statistical features of the orbital elements are explained and a comparison with other catalogues is performed. Illustrative applications are proposed for binaries across the H-R diagram. The binarity is studied in the RGB/AGB and a search for genuine SB1 among long-period variables is performed. The discovery of new EL CVn systems illustrates the potential of combining variability and binarity catalogues. Potential compact object companions are presented, mainly white dwarf companions or double degenerates, but one candidate neutron star is also presented. Towards the bottom of the main sequence, the orbits of previously-suspected binary ultracool dwarfs are determined and new candidate binaries are discovered. The long awaited contribution of Gaia to the analysis of the substellar regime shows the brown dwarf desert around solar-type stars using true, rather than minimum, masses, and provides new important constraints on the occurrence rates of substellar companions to M dwarfs. Several dozen new exoplanets are proposed, including two with validated orbital solutions and one super-Jupiter orbiting a white dwarf, all being candidates requiring confirmation. Beside binarity, higher order multiple systems are also found. △ Less

Submitted 11 June, 2022; originally announced June 2022.

Comments: 60 pages, 60 figures. Accepted for publication in Astronomy & Astrophysics (2022-06-09). The catalogue of binary masses is available for download from the ESA Gaia DR3 Archive and will be available from the CDS/VizieR service

Journal ref: A&A 674, A34 (2023)

Gaia Data Release 3: Chemical cartography of the Milky Way

Authors: Gaia Collaboration, A. Recio-Blanco, G. Kordopatis, P. de Laverny, P. A. Palicio, A. Spagna, L. Spina, D. Katz, P. Re Fiorentin, E. Poggio, P. J. McMillan, A. Vallenari, M. G. Lattanzi, G. M. Seabroke, L. Casamiquela, A. Bragaglia, T. Antoja, C. A. L. Bailer-Jones, R. Andrae, M. Fouesneau, M. Cropper, T. Cantat-Gaudin, U. Heiter, A. Bijaoui, A. G. A. Brown , et al. (425 additional authors not shown)

Abstract: Gaia DR3 opens a new era of all-sky spectral analysis of stellar populations thanks to the nearly 5.6 million stars observed by the RVS and parametrised by the GSP-spec module. The all-sky Gaia chemical cartography allows a powerful and precise chemo-dynamical view of the Milky Way with unprecedented spatial coverage and statistical robustness. First, it reveals the strong vertical symmetry of the… ▽ More Gaia DR3 opens a new era of all-sky spectral analysis of stellar populations thanks to the nearly 5.6 million stars observed by the RVS and parametrised by the GSP-spec module. The all-sky Gaia chemical cartography allows a powerful and precise chemo-dynamical view of the Milky Way with unprecedented spatial coverage and statistical robustness. First, it reveals the strong vertical symmetry of the Galaxy and the flared structure of the disc. Second, the observed kinematic disturbances of the disc -- seen as phase space correlations -- and kinematic or orbital substructures are associated with chemical patterns that favour stars with enhanced metallicities and lower [alpha/Fe] abundance ratios compared to the median values in the radial distributions. This is detected both for young objects that trace the spiral arms and older populations. Several alpha, iron-peak elements and at least one heavy element trace the thin and thick disc properties in the solar cylinder. Third, young disc stars show a recent chemical impoverishment in several elements. Fourth, the largest chemo-dynamical sample of open clusters analysed so far shows a steepening of the radial metallicity gradient with age, which is also observed in the young field population. Finally, the Gaia chemical data have the required coverage and precision to unveil galaxy accretion debris and heated disc stars on halo orbits through their [alpha/Fe] ratio, and to allow the study of the chemo-dynamical properties of globular clusters. Gaia DR3 chemo-dynamical diagnostics open new horizons before the era of ground-based wide-field spectroscopic surveys. They unveil a complex Milky Way that is the outcome of an eventful evolution, shaping it to the present day (abridged). △ Less

Submitted 11 June, 2022; originally announced June 2022.

Comments: Astronomy and Astrophysics (accepted, in press)

Journal ref: A&A 674, A38 (2023)

Gaia DR3 astrometric orbit determination with Markov Chain Monte Carlo and Genetic Algorithms. Systems with stellar, substellar, and planetary mass companions

Authors: B. Holl, A. Sozzetti, J. Sahlmann, P. Giacobbe, D. Ségransan, N. Unger, J. B. Delisle, D. Barbato, M. G. Lattanzi, R. Morbidelli, D. Sosnowska

Abstract: Astrometric discovery of sub-stellar mass companions orbiting stars is exceedingly hard due to the required sub-milliarcsecond precision, limiting the application of this technique to only a few instruments on a target-per-target basis as well as the global astrometry space missions Hipparcos and Gaia. The third Gaia data release includes the first Gaia astrometric orbital solutions, whose sensiti… ▽ More Astrometric discovery of sub-stellar mass companions orbiting stars is exceedingly hard due to the required sub-milliarcsecond precision, limiting the application of this technique to only a few instruments on a target-per-target basis as well as the global astrometry space missions Hipparcos and Gaia. The third Gaia data release includes the first Gaia astrometric orbital solutions, whose sensitivity in terms of estimated companion mass extends down into the planetary-mass regime. We present the contribution of the `exoplanet pipeline' to the Gaia DR3 sample of astrometric orbital solutions by describing the methods used for fitting the orbits, the identification of significant solutions, and their validation. We then present an overview of the statistical properties of the solution parameters. Using both a Markov Chain Monte Carlo and Genetic Algorithm we fit the 34 months of Gaia DR3 astrometric time series with a single Keplerian astrometric-orbit model. Verification and validation steps are taken using significance tests, internal consistency checks using the Gaia radial velocity measurements (when available), as well as literature radial velocity and astrometric data, leading to a subset of candidates that are labelled as 'validated'. We determined astrometric-orbit solutions for 1162 sources and 198 solutions have been assigned the 'validated' label. Precise companion mass estimates are presented elsewhere. From internal and external verification and validation we estimate the level of spurious/incorrect solutions in our sample to be of the order of ~5-10% in our non-'validated' candidate samples. We demonstrate that Gaia is able to confirm and sometimes refine known orbital companion orbits as well as identify new candidates, providing us with a positive outlook of the expected harvest from the full mission data in future data releases. △ Less

Submitted 30 September, 2022; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: 21 Figures, 6 Tables, submitted to A&A (v3: first revision with improvements in the text in response to the anonymous referee report, and fixed incorrect T0 quoted in text for HIP66074 noted by Marcus Marcussen)

Journal ref: A&A 674, A10 (2023)

A general stability-driven approach for the refinement of multi-planet systems

Authors: M. Stalport, J. -B. Delisle, S. Udry, E. C. Matthews, V. Bourrier, A. Leleu

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: 14 pages, 14 figures. Accepted for publication in A&A

Journal ref: A&A 664, A53 (2022)

Gaia Early Data Release 3: The celestial reference frame (Gaia-CRF3)

Authors: Gaia Collaboration, S. A. Klioner, L. Lindegren, F. Mignard, J. Hernández, M. Ramos-Lerate, U. Bastian, M. Biermann, A. Bombrun, A. de Torres, E. Gerlach, R. Geyer, T. Hilger, D. Hobbs, U. L. Lammers, P. J. McMillan, H. Steidelmüller, D. Teyssier, C. M. Raiteri, S. Bartolomé, M. Bernet, J. Castañeda, M. Clotet, M. Davidson, C. Fabricius , et al. (426 additional authors not shown)

Abstract: Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. We describe the c… ▽ More Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. We describe the construction of Gaia-CRF3, and its properties in terms of the distributions in magnitude, colour, and astrometric quality. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasars (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio-optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). The Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13 to 21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 $μ$as yr${}^{-1}$ on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but exceeds 4 mas in either coordinate for 127 sources. We outline the future of the Gaia-CRF in the next Gaia data releases. △ Less

Submitted 30 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Journal ref: A&A 667, A148 (2022)

A continuous multiple hypothesis testing framework for optimal exoplanet detection

Authors: Nathan C. Hara, Thibault de Poyferré, Jean-Baptiste Delisle, Marc Hoffmann

Abstract: When searching for exoplanets, one wants to count how many planets orbit a given star, and to determine what their orbital parameters are. If the estimated orbital elements are too far from those of a planet truly present, this should be considered as a false detection. This setting is a particular instance of a general one: aiming to retrieve which parametric components are in a dataset corrupted… ▽ More When searching for exoplanets, one wants to count how many planets orbit a given star, and to determine what their orbital parameters are. If the estimated orbital elements are too far from those of a planet truly present, this should be considered as a false detection. This setting is a particular instance of a general one: aiming to retrieve which parametric components are in a dataset corrupted by nuisance signals, with a certain accuracy on their parameters. We exhibit a detection criterion minimizing false and missed detections, either as a function of their relative cost, or when the expected number of false detections is bounded. If the components can be separated in a technical sense discussed in detail, the optimal detection criterion is a posterior probability obtained as a by-product of Bayesian evidence calculations, and we discuss how it can be calibrated. We show on two simulations emulating exoplanet searches that it can significantly outperform other criteria. Finally, we show that our framework offers solutions for the identification of components of mixture models, and Bayesian false discovery rate control when hypotheses are not discrete. △ Less

Submitted 22 August, 2023; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Annals of Applied Statistics

Journal ref: Ann. Appl. Stat. 18(1): 749-769 (March 2024)

Alleviating the Transit Timing Variations bias in transit surveys. II. RIVERS: Twin resonant Earth-sized planets around Kepler-1972 recovered from Kepler's false positive

Authors: A. Leleu, J. -B. Delisle, R. Mardling, S. Udry, G. Chatel, Y. Alibert, P. Eggenberger

Abstract: 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 Kepl… ▽ More 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. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: arXiv admin note: text overlap with arXiv:2111.06825

Efficient modeling of correlated noise. III. Scalable methods for jointly modeling several observables' time series with Gaussian processes

Authors: J. -B. Delisle, N. Unger, N. C. Hara, D. Ségransan

Abstract: The radial velocity method is a very productive technique used to detect and confirm extrasolar planets. The most recent spectrographs, such as ESPRESSO or EXPRES, have the potential to detect Earth-like planets around Sun-like stars. However, stellar activity can induce radial velocity variations that dilute or even mimic the signature of a planet. A widely recognized method for disentangling the… ▽ More The radial velocity method is a very productive technique used to detect and confirm extrasolar planets. The most recent spectrographs, such as ESPRESSO or EXPRES, have the potential to detect Earth-like planets around Sun-like stars. However, stellar activity can induce radial velocity variations that dilute or even mimic the signature of a planet. A widely recognized method for disentangling these signals is to model the radial velocity time series, jointly with stellar activity indicators, using Gaussian processes and their derivatives. However, such modeling is prohibitive in terms of computational resources for large data sets, as the cost typically scales as the total number of measurements cubed. Here, we present S+LEAF 2, a Gaussian process framework that can be used to jointly model several time series, with a computational cost that scales linearly with the data set size. This framework thus provides a state-of-the-art Gaussian process model, with tractable computations even for large data sets. We illustrate the power of this framework by reanalyzing the 246 HARPS radial velocity measurements of the nearby K2 dwarf HD 138038, together with two activity indicators. We reproduce the results of a previous analysis of these data, but with a strongly decreased computational cost (more than two order of magnitude). The gain would be even greater for larger data sets. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 659, A182 (2022)

Spin-orbit coupling for close-in planets

Authors: Alexandre C. M. Correia, Jean-Baptiste Delisle

Abstract: We study the spin evolution of close-in planets in multi-body systems and present a very general formulation of the spin-orbit problem. This includes a simple way to probe the spin dynamics from the orbital perturbations, a new method for computing forced librations and tidal deformation, and general expressions for the tidal torque and capture probabilities in resonance. We show that planet-plane… ▽ More We study the spin evolution of close-in planets in multi-body systems and present a very general formulation of the spin-orbit problem. This includes a simple way to probe the spin dynamics from the orbital perturbations, a new method for computing forced librations and tidal deformation, and general expressions for the tidal torque and capture probabilities in resonance. We show that planet-planet perturbations can drive the spin of Earth-size planets into asynchronous or chaotic states, even for nearly circular orbits. We apply our results to Mercury and to the KOI-1599 system of two super-Earths in a 3/2 mean motion resonance. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: 12 pages, 9 figures

Journal ref: Astronomy & Astrophysics, vol. 630, A102 (2019)

Alleviating the transit timing variation bias in transit surveys. I. RIVERS: Method and detection of a pair of resonant super-Earths around Kepler-1705

Authors: A. Leleu, G. Chatel, S. Udry, Y. Alibert, J. -B. Delisle, R. Mardling

Abstract: 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 detect… ▽ More 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. △ Less

Submitted 12 November, 2021; originally announced November 2021.

The Rossiter-McLaughlin effect Revolutions: An ultra-short period planet and a warm mini-Neptune on perpendicular orbits

Authors: V. Bourrier, C. Lovis, M. Cretignier, R. Allart, X. Dumusque, J. -B. Delisle, A. Deline, S. G. Sousa, V. Adibekyan, Y. Alibert, S. C. C. Barros, F. Borsa, S. Cristiani, O. Demangeon, D. Ehrenreich, P. Figueira, J. I. González Hernández, M. Lendl, J. Lillo-Box, G. Lo Curto, P. Di Marcantonio, C. J. A. P. Martins, D. Mégevand, A. Mehner, G. Micela , et al. (12 additional authors not shown)

Abstract: Comparisons of the alignment of exoplanets with a common host star can be used to distinguish among concurrent evolution scenarios. However, multi-planet systems usually host mini-Neptunes and super-Earths, whose size make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions technique, which can access spin-orbit angles of small planets by exploiti… ▽ More Comparisons of the alignment of exoplanets with a common host star can be used to distinguish among concurrent evolution scenarios. However, multi-planet systems usually host mini-Neptunes and super-Earths, whose size make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions technique, which can access spin-orbit angles of small planets by exploiting the full information contained in spectral transit time series. We validated the technique on published HARPS-N data of the mini-Neptune HD3167c, refining its high sky-projected spin-orbit angle (-108.9+5.4-5.5 deg), and we applied it to new ESPRESSO observations of the super-Earth HD3167b, revealing an aligned orbit (-6.6+6.6-7.9 deg). Surprisingly different variations in the contrast of the stellar lines occulted by the planets can be reconciled with a latitudinal dependence of the stellar line shape. In this scenario, a joint fit to both datasets constrains the inclination of the star (111.6+3.1-3.3 deg) and the 3D spin-orbit angles of HD3167b (29.5+7.2-9.4 deg) and HD3167c (107.7+5.1-4.9 deg). The projected spin-orbit angles do not depend on the model for the line contrast variations, and so, with a mutual inclination of 102.3+7.4-8.0 deg, we conclude that the two planets are on perpendicular orbits. This could be explained by HD3167b being strongly coupled to the star and retaining its primordial alignment, whereas HD3167c would have been brought to a nearly polar orbit via secular gravitational interactions with an outer companion. Follow-up observations and dynamical evolution simulations are required to search for this companion and explore this scenario. HD3167b is the smallest exoplanet with a confirmed spectroscopic Rossiter-McLaughlin signal. Our new technique opens the way to determining the orbital architectures of the super-Earth and Earth-sized planet populations. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 20 pages, 19 figures, accepted for publication in A&A (14 September 2021)

The HARPS search for southern extra-solar planets. XLVI: 12 super-Earths around the solar type stars HD39194, HD93385, HD96700, HD154088, and HD189567

Authors: N. Unger, D. Ségransan, D. Queloz, S. Udry, C. Lovis, C. Mordasini, E. Ahrer, W. Benz, F. Bouchy, J. -B. Delisle, R. F. Díaz, X. Dumusque, G. Lo Curto, M. Marmier, M. Mayor, F. Pepe, N. C. Santos, M. Stalport, R. Alonso, A. Collier Cameron, M. Deleuil, P. Figueira, M. Gillon, C. Moutou, D. Pollacco , et al. (1 additional authors not shown)

Abstract: Context. We present precise radial-velocity measurements of five solar-type stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). With a time span of more than 10 years and a fairly dense sampling, the survey is sensitive to low mass planets down to super-Earths on orbital periods up to 100 days. Aims. Our goal was to search for planetary compa… ▽ More Context. We present precise radial-velocity measurements of five solar-type stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). With a time span of more than 10 years and a fairly dense sampling, the survey is sensitive to low mass planets down to super-Earths on orbital periods up to 100 days. Aims. Our goal was to search for planetary companions around the stars HD39194, HD93385, HD96700, HD154088, and HD189567 and use Bayesian model comparison to make an informed choice on the number of planets present in the systems based on the radial velocity observations. These findings will contribute to the pool of known exoplanets and better constrain their orbital parameters. Methods. A first analysis was performed using the DACE (Data & Analysis Center for Exoplanets) online tools to assess the activity level of the star and the potential planetary content of each system. We then used Bayesian model comparison on all targets to get a robust estimate of the number of planets per star. We did this using the nested sampling algorithm PolyChord. For some targets, we also compared different noise models to disentangle planetary signatures from stellar activity. Lastly, we ran an efficient MCMC (Markov chain Monte Carlo) algorithm for each target to get reliable estimates for the planets' orbital parameters. Results. We identify 12 planets within several multiplanet systems. These planets are all in the super-Earth and sub-Neptune mass regime with minimum masses ranging between 4 and 13 M$_\oplus$ and orbital periods between 5 and 103 days. Three of these planets are new, namely HD 93385 b, HD 96700 c, and HD 189567 c. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 19 pages, 17 figures

Journal ref: A&A 654, A104 (2021)

Transit detection of the long-period volatile-rich super-Earth $ν^2$ Lupi d with $CHEOPS$

Authors: Laetitia Delrez, David Ehrenreich, Yann Alibert, Andrea Bonfanti, Luca Borsato, Luca Fossati, Matthew J. Hooton, Sergio Hoyer, Francisco J. Pozuelos, Sébastien Salmon, Sophia Sulis, Thomas G. Wilson, Vardan Adibekyan, Vincent Bourrier, Alexis Brandeker, Sébastien Charnoz, Adrien Deline, Pascal Guterman, Jonas Haldemann, Nathan Hara, Mahmoudreza Oshagh, Sergio G. Sousa, Valérie Van Grootel, Roi Alonso, Guillem Anglada Escudé , et al. (53 additional authors not shown)

Abstract: Exoplanets transiting bright nearby stars are key objects for advancing our knowledge of planetary formation and evolution. The wealth of photons from the host star gives detailed access to the atmospheric, interior, and orbital properties of the planetary companions. $ν^2$ Lupi (HD 136352) is a naked-eye ($V = 5.78$) Sun-like star that was discovered to host three low-mass planets with orbital pe… ▽ More Exoplanets transiting bright nearby stars are key objects for advancing our knowledge of planetary formation and evolution. The wealth of photons from the host star gives detailed access to the atmospheric, interior, and orbital properties of the planetary companions. $ν^2$ Lupi (HD 136352) is a naked-eye ($V = 5.78$) Sun-like star that was discovered to host three low-mass planets with orbital periods of 11.6, 27.6, and 107.6 days via radial velocity monitoring (Udry et al. 2019). The two inner planets (b and c) were recently found to transit (Kane et al. 2020), prompting a photometric follow-up by the brand-new $CHaracterising\:ExOPlanets\:Satellite\:(CHEOPS)$. Here, we report that the outer planet d is also transiting, and measure its radius and mass to be $2.56\pm0.09$ $R_{\oplus}$ and $8.82\pm0.94$ $M_{\oplus}$, respectively. With its bright Sun-like star, long period, and mild irradiation ($\sim$5.7 times the irradiation of Earth), $ν^2$ Lupi d unlocks a completely new region in the parameter space of exoplanets amenable to detailed characterization. We refine the properties of all three planets: planet b likely has a rocky mostly dry composition, while planets c and d seem to have retained small hydrogen-helium envelopes and a possibly large water fraction. This diversity of planetary compositions makes the $ν^2$ Lupi system an excellent laboratory for testing formation and evolution models of low-mass planets. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: Published in Nature Astronomy. 60 pages, 18 Figures, 6 Tables. This is the authors' version of the manuscript. The final authenticated version is available online at https://www.nature.com/articles/s41550-021-01381-5

Journal ref: Nature Astronomy, Volume 5, Pages 775-787, June 2021

Testing whether a signal is strictly periodic. Application to disentangling planets and stellar activity in radial velocities

Authors: Nathan C. Hara, Jean-Baptiste Delisle, Nicolas Unger, Xavier Dumusque

Abstract: Searches for periodicity in time series are often done with models of periodic signals, whose statistical significance is assessed via false alarm probabilities or Bayes factors. However, a statistically significant periodic model might not originate from a strictly periodic source. In astronomy in particular, one expects transient signals that show periodicity for a certain amount of time before… ▽ More Searches for periodicity in time series are often done with models of periodic signals, whose statistical significance is assessed via false alarm probabilities or Bayes factors. However, a statistically significant periodic model might not originate from a strictly periodic source. In astronomy in particular, one expects transient signals that show periodicity for a certain amount of time before vanishing. This situation is encountered for instance in the search for planets in radial velocity data. While planetary signals are expected to have a stable phase, amplitude and frequency - except when strong planet-planet interactions are present - signals induced by stellar activity will typically not exhibit the same stability. In the present article, we explore the use of periodic functions multiplied by time windows to diagnose whether an apparently periodic signal is truly so. We suggest diagnostics to check whether a signal is consistently present in the time series, and has a stable phase, amplitude and period. The tests are expressed both in a periodogram and Bayesian framework. Our methods are applied to the Solar HARPS-N data as well as HD 215152, HD 69830 and HD 13808. We find that (i) the HARPS-N Solar data exhibits signals at the Solar rotation period and its first harmonic ($\sim$ 13.4 days). The frequency and phase of the 13.4 days signal appear constant within the estimation uncertainties, but its amplitude presents significant variations which can be mapped to activity levels. (ii) as previously reported, we find four, three and two planets orbiting HD 215152, HD 69830 and HD 13808. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

MSC Class: 62P35 ACM Class: G.3

Journal ref: Astronomy & Astrophysics, Volume 658, A177, 16 pp., February 2022

Improving exoplanet detection capabilities with the false inclusion probability. Comparison with other detection criteria in the context of radial velocities

Authors: Nathan C. Hara, Nicolas Unger, Jean-Baptiste Delisle, Rodrigo Díaz, Damien Ségransan

Abstract: Context. In exoplanet searches with radial velocity data, the most common statistical significance metrics are the Bayes factor and the false alarm probability (FAP). Both have proved useful, but do not directly address whether an exoplanet detection should be claimed. Furthermore, it is unclear which detection threshold should be taken and how robust the detections are to model misspecification.… ▽ More Context. In exoplanet searches with radial velocity data, the most common statistical significance metrics are the Bayes factor and the false alarm probability (FAP). Both have proved useful, but do not directly address whether an exoplanet detection should be claimed. Furthermore, it is unclear which detection threshold should be taken and how robust the detections are to model misspecification. Aims. The present work aims at defining a detection criterion which conveys as precisely as possible the information needed to claim an exoplanet detection. We compare this new criterion to existing ones in terms of sensitivity and robustness. Methods. We define a significance metric called the false inclusion probability (FIP) based on the posterior probability of presence of a planet. Posterior distributions are computed with the nested sampling package Polychord. We show that for FIP and Bayes factor calculations, defining priors on linear parameters as Gaussian mixture models allows to significantly speed up computations. The performances of the FAP, Bayes factor and FIP are studied with simulations as well as analytical arguments. We compare the methods assuming the model is correct, then evaluate their sensitivity to the prior and likelihood choices. Results. Among other properties, the FIP offers ways to test the reliability of the significance levels, it is particularly efficient to account for aliasing and allows to exclude the presence of planets with a certain confidence. We find that, in our simulations, the FIP outperforms existing detection metrics. We show that planet detections are sensitive to priors on period and semi-amplitude and that letting free the noise parameters offers better performances than fixing a noise model based on a fit to ancillary indicators. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

MSC Class: 62P35; 62-08 ACM Class: G.3

Journal ref: A&A 663, A14 (2022)

Six transiting planets and a chain of Laplace resonances in TOI-178

Authors: A. Leleu, Y. Alibert, N. C. Hara, M. J. Hooton, T. G. Wilson, P. Robutel, J. -B. Delisle, J. Laskar, S. Hoyer, C. Lovis, E. M. Bryant, E. Ducrot, J. Cabrera, L. Delrez, J. S. Acton, V. Adibekyan, R. Allart, C. Allende Prieto, R. Alonso, D. Alves, D. R. Anderson, D. Angerhausen, G. Anglada Escudé, J. Asquier, D. Barrado , et al. (130 additional authors not shown)

Abstract: Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this cont… ▽ More Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Gaia Early Data Release 3: The Galactic anticentre

Authors: Gaia Collaboration, T. Antoja, P. McMillan, G. Kordopatis, P. Ramos, A. Helmi, E. Balbinot, T. Cantat-Gaudin, L. Chemin, F. Figueras, C. Jordi, S. Khanna, M. Romero-Gomez, G. Seabroke, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen , et al. (395 additional authors not shown)

Abstract: We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current d… ▽ More We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in-situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. We find that: i) the dynamics of the Galactic disc are very complex with vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and v) the open clusters Berkeley~29 and Saurer~1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits. We demonstrate how, once again, the Gaia are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history. △ Less

Submitted 26 April, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: Gaia EDR3 performance verification paper, version 2 closer to published version in A&A, complete list of authors

Journal ref: A&A 649, A8 (2021)

Gaia Early Data Release 3: The Gaia Catalogue of Nearby Stars

Authors: Gaia Collaboration, R. L. Smart, L. M. Sarro, J. Rybizki, C. Reylé, A. C. Robin, N. C. Hambly, U. Abbas, M. A. Barstow, J. H. J. de Bruijne, B. Bucciarelli, J. M. Carrasco, W. J. Cooper, S. T. Hodgkin, E. Masana, D. Michalik, J. Sahlmann, A. Sozzetti, A. G. A. Brown, A. Vallenari, T. Prusti, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (398 additional authors not shown)

Abstract: We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of obj… ▽ More We produce a clean and well-characterised catalogue of objects within 100\,pc of the Sun from the \G\ Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. The selection of objects within 100\,pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100\,pc is included in the catalogue. We have produced a catalogue of \NFINAL\ objects that we estimate contains at least 92\% of stars of stellar type M9 within 100\,pc of the Sun. We estimate that 9\% of the stars in this catalogue probably lie outside 100\,pc, but when the distance probability function is used, a correct treatment of this contamination is possible. We produced luminosity functions with a high signal-to-noise ratio for the main-sequence stars, giants, and white dwarfs. We examined in detail the Hyades cluster, the white dwarf population, and wide-binary systems and produced candidate lists for all three samples. We detected local manifestations of several streams, superclusters, and halo objects, in which we identified 12 members of \G\ Enceladus. We present the first direct parallaxes of five objects in multiple systems within 10\,pc of the Sun. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 45 Pages, 39 figures in main part and 18 in appendix, tables on CDS

Journal ref: A&A 649, A6 (2021)

Gaia Early Data Release 3: Acceleration of the solar system from Gaia astrometry

Authors: Gaia Collaboration, S. A. Klioner, F. Mignard, L. Lindegren, U. Bastian, P. J. McMillan, J. Hernández, D. Hobbs, M. Ramos-Lerate, M. Biermann, A. Bombrun, A. de Torres, E. Gerlach, R. Geyer, T. Hilger, U. Lammers, H. Steidelmüller, C. A. Stephenson, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, O. L. Creevey, D. W. Evans , et al. (392 additional authors not shown)

Abstract: Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the… ▽ More Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar system barycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. The effect of the acceleration is obtained as a part of the general expansion of the vector field of proper motions in Vector Spherical Harmonics (VSH). Various versions of the VSH fit and various subsets of the sources are tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution are used to get a better idea on possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is $(2.32 \pm 0.16) \times 10^{-10}$ m s${}^{-2}$ (or $7.33 \pm 0.51$ km s$^{-1}$ Myr${}^{-1}$) towards $α= 269.1^\circ \pm 5.4^\circ$, $δ= -31.6^\circ \pm 4.1^\circ$, corresponding to a proper motion amplitude of $5.05 \pm 0.35$ $μ$as yr${}^{-1}$. This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 $μ$as yr${}^{-1}$. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: A&A, accepted

Journal ref: A&A 649, A9 (2021)

Gaia Early Data Release 3: Structure and properties of the Magellanic Clouds

Authors: Gaia Collaboration, X. Luri, L. Chemin, G. Clementini, H. E. Delgado, P. J. McMillan, M. Romero-Gómez, E. Balbinot, A. Castro-Ginard, R. Mor, V. Ripepi, L. M. Sarro, M. -R. L. Cioni, C. Fabricius, A. Garofalo, A. Helmi, T. Muraveva, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans , et al. (395 additional authors not shown)

Abstract: We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasib… ▽ More We compare the Gaia DR2 and Gaia EDR3 performances in the study of the Magellanic Clouds and show the clear improvements in precision and accuracy in the new release. We also show that the systematics still present in the data make the determination of the 3D geometry of the LMC a difficult endeavour; this is at the very limit of the usefulness of the Gaia EDR3 astrometry, but it may become feasible with the use of additional external data. We derive radial and tangential velocity maps and global profiles for the LMC for the several subsamples we defined. To our knowledge, this is the first time that the two planar components of the ordered and random motions are derived for multiple stellar evolutionary phases in a galactic disc outside the Milky Way, showing the differences between younger and older phases. We also analyse the spatial structure and motions in the central region, the bar, and the disc, providing new insights into features and kinematics. Finally, we show that the Gaia EDR3 data allows clearly resolving the Magellanic Bridge, and we trace the density and velocity flow of the stars from the SMC towards the LMC not only globally, but also separately for young and evolved populations. This allows us to confirm an evolved population in the Bridge that is slightly shift from the younger population. Additionally, we were able to study the outskirts of both Magellanic Clouds, in which we detected some well-known features and indications of new ones. △ Less

Submitted 4 January, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: This paper is part of the "demonstration papers" released with Gaia EDR3: https://www.cosmos.esa.int/web/gaia/earlydr3

Journal ref: A&A 649, A7 (2021)

Gaia Early Data Release 3: Summary of the contents and survey properties

Authors: Gaia Collaboration, A. G. A Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, M. Biermann, O. L. Creevey, D. W. Evans, L. Eyer, A. Hutton, F. Jansen, C. Jordi, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, C. Soubiran, N. A. Walton, F. Arenou , et al. (401 additional authors not shown)

Abstract: We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motio… ▽ More We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the (G_BP-G_RP) colour are also available. The passbands for G, G_BP, and G_RP are provided as part of the release. For ease of use, the 7 million radial velocities from Gaia DR2 are included in this release, after the removal of a small number of spurious values. New radial velocities will appear as part of Gaia DR3. Finally, Gaia EDR3 represents an updated materialisation of the celestial reference frame (CRF) in the optical, the Gaia-CRF3, which is based solely on extragalactic sources. The creation of the source list for Gaia EDR3 includes enhancements that make it more robust with respect to high proper motion stars, and the disturbing effects of spurious and partially resolved sources. The source list is largely the same as that for Gaia DR2, but it does feature new sources and there are some notable changes. The source list will not change for Gaia DR3. Gaia EDR3 represents a significant advance over Gaia DR2, with parallax precisions increased by 30 percent, proper motion precisions increased by a factor of 2, and the systematic errors in the astrometry suppressed by 30--40 percent for the parallaxes and by a factor ~2.5 for the proper motions. The photometry also features increased precision, but above all much better homogeneity across colour, magnitude, and celestial position. A single passband for G, G_BP, and G_RP is valid over the entire magnitude and colour range, with no systematics above the 1 percent level. △ Less

Submitted 9 June, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Comments: Accepted for A&A Special Issue on Gaia EDR3, 21 pages, 2 figures. This version includes the updates in the erratum (https://doi.org/10.1051/0004-6361/202039657e)

Journal ref: A&A 650, C3 (2021)

Why do warm Neptunes present nonzero eccentricity?

Authors: A. C. M. Correia, V. Bourrier, J. -B. Delisle

Abstract: Most Neptune-mass planets in close-in orbits (orbital periods less than a few days) present nonzero eccentricity, typically around 0.15. This is somehow unexpected, as these planets undergo strong tidal dissipation that should circularize their orbits in a time-scale shorter than the age of the system. In this paper we discuss some mechanisms that can oppose to bodily tides, namely, thermal atmosp… ▽ More Most Neptune-mass planets in close-in orbits (orbital periods less than a few days) present nonzero eccentricity, typically around 0.15. This is somehow unexpected, as these planets undergo strong tidal dissipation that should circularize their orbits in a time-scale shorter than the age of the system. In this paper we discuss some mechanisms that can oppose to bodily tides, namely, thermal atmospheric tides, evaporation of the atmosphere, and excitation from a distant companion. In the first two cases, the eccentricity can increase consistently, while in the last one, the eccentricity can only be excited for a limited amount of time (that may nevertheless exceed the age of the system). We show the limitations of these different mechanisms and how some of them could, depending on specific properties of the observed planetary systems, account for their presently observed eccentricities. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 8 pages, 2 figures

Journal ref: Astronomy & Astrophysics, vol. 635, A37 (2020)

Efficient modeling of correlated noise II. A flexible noise model with fast and scalable methods

Authors: J. -B. Delisle, N. Hara, D. Ségransan

Abstract: Correlated noise affects most astronomical datasets and to neglect accounting for it can lead to spurious signal detections, especially in low signal-to-noise conditions, which is often the context in which new discoveries are pursued. For instance, in the realm of exoplanet detection with radial velocity time series, stellar variability can induce false detections. However, a white noise approxim… ▽ More Correlated noise affects most astronomical datasets and to neglect accounting for it can lead to spurious signal detections, especially in low signal-to-noise conditions, which is often the context in which new discoveries are pursued. For instance, in the realm of exoplanet detection with radial velocity time series, stellar variability can induce false detections. However, a white noise approximation is often used because accounting for correlated noise when analyzing data implies a more complex analysis. Moreover, the computational cost can be prohibitive as it typically scales as the cube of the dataset size. For some restricted classes of correlated noise models, there are specific algorithms that can be used to help bring down the computational cost. This improvement in speed is particularly useful in the context of Gaussian process regression, however, it comes at the expense of the generality of the noise model. Here, we present the S+LEAF noise model, which allows us to account for a large class of correlated noises with a linear scaling of the computational cost with respect to the size of the dataset. The S+LEAF model includes, in particular, mixtures of quasiperiodic kernels and calibration noise. This efficient modeling is made possible by a sparse representation of the covariance matrix of the noise and the use of dedicated algorithms for matrix inversion, solving, determinant computation, etc. We applied the S+LEAF model to reanalyze the HARPS radial velocity time series of HD 136352. We illustrate the flexibility of the S+LEAF model in handling various sources of noise. We demonstrate the importance of taking correlated noise into account, and especially calibration noise, to correctly assess the significance of detected signals. We provide an open-source implementation of the S+LEAF model, available at https://gitlab.unige.ch/jean-baptiste.delisle/spleaf. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: Accepted in A&A

Journal ref: A&A 638, A95 (2020)

SPHERE+: Imaging young Jupiters down to the snowline

Authors: A. Boccaletti, G. Chauvin, D. Mouillet, O. Absil, F. Allard, S. Antoniucci, J. -C. Augereau, P. Barge, A. Baruffolo, J. -L. Baudino, P. Baudoz, M. Beaulieu, M. Benisty, J. -L. Beuzit, A. Bianco, B. Biller, B. Bonavita, M. Bonnefoy, S. Bos, J. -C. Bouret, W. Brandner, N. Buchschache, B. Carry, F. Cantalloube, E. Cascone , et al. (108 additional authors not shown)

Abstract: SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with S… ▽ More SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei. △ Less

Submitted 13 March, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: White paper submitted to ESO on Feb. 20th, 2020

Spectral and atmospheric characterisation of a new benchmark brown dwarf HD13724B

Authors: E. L. Rickman, D. Ségransan, J. Hagelberg, J. -L. Beuzit, A. Cheetham, J. -B. Delisle, T. Forveille, S. Udry

Abstract: Context. HD13724 is a nearby solar-type star at 43.48 $\pm$ 0.06 pc hosting a long-period low-mass brown dwarf detected with the CORALIE echelle spectrograph as part of the historical CORALIE radial-velocity search for extra-solar planets. The companion has a minimum mass of $26.77^{+4.4}_{-2.2} M_{\mathrm{Jup}}$ and an expected semi-major axis of $\sim$ 240 mas making it a suitable target for fur… ▽ More Context. HD13724 is a nearby solar-type star at 43.48 $\pm$ 0.06 pc hosting a long-period low-mass brown dwarf detected with the CORALIE echelle spectrograph as part of the historical CORALIE radial-velocity search for extra-solar planets. The companion has a minimum mass of $26.77^{+4.4}_{-2.2} M_{\mathrm{Jup}}$ and an expected semi-major axis of $\sim$ 240 mas making it a suitable target for further characterisation with high-contrast imaging, in particular to measure its inclination, mass, and spectrum and thus establish its substellar nature. Aims. Using high-contrast imaging with the SPHERE instrument on the Very Large Telescope (VLT), we are able to directly image a brown dwarf companion to HD13724 and obtain a low-resolution spectrum. Methods. We combine the radial-velocity measurements of CORALIE and HARPS taken over two decades and high contrast imaging from SPHERE to obtain a dynamical mass estimate. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to J band, as well as photometric measurements from IRDIS in the J, H and K bands. Results. Using high-contrast imaging with the SPHERE instrument at the VLT, we report the first images of a brown dwarf companion to the host star HD13724. It has an angular separation of 175.6 $\pm$ 4.5 mas and H-band contrast of $10.61\pm0.16$ mag and, using the age estimate of the star to be $\sim$1 Gyr, gives an isochronal mass estimate of $\sim$44 $M_{\mathrm{Jup}}$. By combining radial-velocity and imaging data we also obtain a dynamical mass of $50.5^{+3.3}_{-3.5} M_{\mathrm{Jup}}$. Through fitting an atmospheric model, we estimate a surface gravity of $\log g = 5.5$ and an effective temperature of 1000K. A comparison of its spectrum with observed T dwarfs estimates a spectral type of T4 or T4.5, with a T4 object providing the best fit. △ Less

Submitted 24 February, 2020; v1 submitted 18 February, 2020; originally announced February 2020.

Comments: 11 pages, Accepted to A&A, updated with language edits

Journal ref: A&A 635, A203 (2020)

Efficient modeling of correlated noise I. Statistical significance of periodogram peaks

Authors: J. -B. Delisle, N. Hara, D. Ségransan

Abstract: Periodograms are common tools used to search for periodic signals in unevenly spaced time series. The significance of periodogram peaks is often assessed using false alarm probability (FAP), which in most studies assumes uncorrelated noise and is computed using numerical methods such as bootstrapping or Monte Carlo. These methods have a high computational cost, especially for low FAP levels, which… ▽ More Periodograms are common tools used to search for periodic signals in unevenly spaced time series. The significance of periodogram peaks is often assessed using false alarm probability (FAP), which in most studies assumes uncorrelated noise and is computed using numerical methods such as bootstrapping or Monte Carlo. These methods have a high computational cost, especially for low FAP levels, which are of most interest. We present an analytical estimate of the FAP of the periodogram in the presence of correlated noise, which is fundamental to analyze astronomical time series correctly. The analytical estimate that we derive provides a very good approximation of the FAP at a much lower cost than numerical methods. We validate our analytical approach by comparing it with Monte Carlo simulations. Finally, we discuss the sensitivity of the method to different assumptions in the modeling of the noise. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: Accepted in A&A

Journal ref: A&A 635, A83 (2020)

The SOPHIE search for northern extrasolar planets. XVI. HD 158259: A compact planetary system in a near-3:2 mean motion resonance chain

Authors: N. C. Hara, F. Bouchy, M. Stalport, I. Boisse, J. Rodrigues, J-. B. Delisle, A. Santerne, G. W. Henry, L. Arnold, N. Astudillo-Defru, S. Borgniet, X. Bonfils, V. Bourrier, B. Brugger, B. Courcol, S. Dalal, M. Deleuil, X. Delfosse, O. Demangeon, R. F. D'iaz, X. Dumusque, T. Forveille, G. Hébrard, M. Hobson, F. Kiefer , et al. (10 additional authors not shown)

Abstract: Since 2011, the SOPHIE spectrograph has been used to search for Neptunes and super-Earths in the Northern Hemisphere. As part of this observational program, 290 radial velocity measurements of the 6.4 V magnitude star HD 158259 were obtained. Additionally, TESS photometric measurements of this target are available. We present an analysis of the SOPHIE data and compare our results with the output o… ▽ More Since 2011, the SOPHIE spectrograph has been used to search for Neptunes and super-Earths in the Northern Hemisphere. As part of this observational program, 290 radial velocity measurements of the 6.4 V magnitude star HD 158259 were obtained. Additionally, TESS photometric measurements of this target are available. We present an analysis of the SOPHIE data and compare our results with the output of the TESS pipeline. The radial velocity data, ancillary spectroscopic indices, and ground-based photometric measurements were analyzed with classical and $\ell_1$ periodograms. The stellar activity was modeled as a correlated Gaussian noise and its impact on the planet detection was measured with a new technique. The SOPHIE data support the detection of five planets, each with $m \sin i \approx 6 M_\oplus$, orbiting HD 158259 in 3.4, 5.2, 7.9, 12, and 17.4 days. Though a planetary origin is strongly favored, the 17.4 d signal is classified as a planet candidate due to a slightly lower statistical significance and to its proximity to the expected stellar rotation period. The data also present low frequency variations, most likely originating from a magnetic cycle and instrument systematics. Furthermore, the TESS pipeline reports a significant signal at 2.17 days corresponding to a planet of radius $\approx 1.2 R_\oplus$. A compatible signal is seen in the radial velocities, which confirms the detection of an additional planet and yields a $\approx 2 M_\oplus$ mass estimate. We find a system of five planets and a strong candidate near a 3:2 mean motion resonance chain orbiting HD 158259. The planets are found to be outside of the two and three body resonances. △ Less

Submitted 26 March, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Bias and robustness of eccentricity estimates from radial velocity data

Authors: N. C. Hara, G. Boué, J. Laskar, J. B. Delisle, N. Unger

Abstract: Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present inter… ▽ More Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present interesting features of classical estimators. Secondly, we study whether the eccentricity estimates are to be trusted when the data contain incorrectly modelled signals, such as missed planetary companions, non Gaussian noises, correlated noises with unknown covariance, etc. Our main conclusion is that data analysis via posterior distributions, with a model including a free error term gives reliable results provided two conditions. First, convergence of the numerical methods needs to be ascertained. Secondly, the noise power spectrum should not have a particularly strong peak at the semi period of the planet of interest. As a consequence, it is difficult to determine if the signal of an apparently eccentric planet might be due to another inner companion in 2:1 mean motion resonance. We study the use of Bayes factors to disentangle these cases. Finally, we suggest methods to check if there are hints of an incorrect model in the residuals. We show on simulated data the performance of our methods and comment on the eccentricities of Proxima b and 55 Cnc f. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS