Researcher profile

S. Udry

S. Udry contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
58works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects' orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects' bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects' atmospheres. It also constrains these objects' bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2025arXiv

Improved characterization of the TOI-2141 system: a dense sub-Neptune with non-transiting inner and outer companions

We aim to refine the fundamental parameters of the TOI-2141 planetary system, which includes a transiting sub-Neptune orbiting a Sun-like star in a relatively long orbit of 18.26 days, by combining new photometric and spectroscopic observations. We analyze new space-based photometry from TESS and CHEOPS as well as 61 radial velocity measurements from HARPS-N. We perform individual and joint photometric and RV analyses using several modeling tools within a Bayesian model comparison framework. We refine the radius and mass of the transiting planet TOI-2141 b to 3.15 $\pm$ 0.04 $R_\oplus$ and 20.1 $\pm$ 1.6 $M_\oplus$, respectively, five and two times more precise than the previously reported values. Our radial velocity analysis reveals two additional non-transiting companions with orbital periods of 5.46 and 60.45 days. Despite the innermost planet's high geometric transit probability, we find no evidence for transits in the photometric data. The bulk properties of TOI-2141 b suggest a significant volatile envelope atop an Earth-like core, with modeling indicating a hydrogen-rich atmosphere that may have experienced mild photoevaporation over the system's history. Planets b and c must exhibit a modest mutual inclination of at least 2.4 degrees.

preprint2023arXiv

DREAM I. Orbital architecture orrery

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

preprint2022arXiv

55 Cancri e's occultation captured with CHEOPS

Past occultation and phase-curve observations of the ultra-short period super-Earth 55 Cnc e obtained at visible and infrared wavelengths have been challenging to reconcile with a planetary reflection and emission model. In this study, we analyse a set of 41 occultations obtained over a two-year timespan with the CHEOPS satellite. We report the detection of 55 Cnc e&#39;s occultation with an average depth of $12\pm3$ ppm. We derive a corresponding 2-$σ$ upper limit on the geometric albedo of $A_g < 0.55$ once decontaminated from the thermal emission measured by Spitzer at 4.5$μ$m. CHEOPS&#39;s photometric performance enables, for the first time, the detection of individual occultations of this super-Earth in the visible and identifies short-timescale photometric corrugations likely induced by stellar granulation. We also find a clear 47.3-day sinusoidal pattern in the time-dependent occultation depths that we are unable to relate to stellar noise, nor instrumental systematics, but whose planetary origin could be tested with upcoming JWST occultation observations of this iconic super-Earth.

preprint2022arXiv

A candidate short-period sub-Earth orbiting Proxima Centauri

Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough evaluation of the presence of a third low-mass planetary companion, which started emerging during a previous campaign. Radial velocities (RVs) were calculated using both a cross-correlation function (CCF) and a template matching approach. The RV analysis includes a component to model Proxima&#39;s activity using a Gaussian process (GP). We use the CCF&#39;s full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit.

preprint2022arXiv

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

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.

preprint2022arXiv

A novel framework for semi-Bayesian radial velocities through template matching

The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alternative for M-dwarf stars. In this paper, we describe a new implementation of template matching within a semi-Bayesian framework, providing a more statistically principled characterization of the RV measurements. In this context, a common RV shift is used to describe the difference between each spectral order of a given stellar spectrum and a template built from the available observations. Posterior probability distributions are obtained for the relative RV associated with each spectrum, after marginalizing with respect to the continuum. This methodology was named S-BART: Semi-Bayesian Approach for RVs with Template-matching, and it can be applied to HARPS and ESPRESSO. The application of our method to HARPS archival observations of Barnard&#39;s star allowed us to validate our implementation against HARPS-TERRA and SERVAL. Then, we applied it to 33 ESPRESSO targets, evaluating its performance and comparing it with the CCF method. We found a decrease in the median RV scatter of \sim 10\% and \sim 4\% for M- and K-type stars, respectively. S-BART yields more precise RV estimates than the CCF method, particularly in the case of M-type stars where a median uncertainty of \sim 15 cm/s is achieved over 309 observations. Further, we estimated the nightly zero point (NZP) of ESPRESSO, finding a weighted NZP scatter below \sim 0.7 m/s. As this includes stellar variability, photon noise, and potential planetary signals, it should be taken as an upper limit of the RV precision attainable with ESPRESSO data.

preprint2022arXiv

A Second Planet Transiting LTT 1445A and a Determination of the Masses of Both Worlds

LTT 1445 is a hierarchical triple M-dwarf star system located at a distance of 6.86 parsecs. The primary star LTT 1445A (0.257 M_Sun) is known to host the transiting planet LTT 1445Ab with an orbital period of 5.4 days, making it the second closest known transiting exoplanet system, and the closest one for which the host is an M dwarf. Using TESS data, we present the discovery of a second planet in the LTT 1445 system, with an orbital period of 3.1 days. We combine radial velocity measurements obtained from the five spectrographs ESPRESSO, HARPS, HIRES, MAROON-X, and PFS to establish that the new world also orbits LTT 1445A. We determine the mass and radius of LTT 1445Ab to be 2.87+/-0.25 M_Earth and 1.304^{+0.067}_{-0.060} R_Earth, consistent with an Earth-like composition. For the newly discovered LTT 1445Ac, we measure a mass of 1.54^{+0.20}_{-0.19} M_Earth and a minimum radius of 1.15 R_Earth, but we cannot determine the radius directly as the signal-to-noise of our light curve permits both grazing and non-grazing configurations. Using MEarth photometry and ground-based spectroscopy, we establish that star C (0.161 M_Sun) is likely the source of the 1.4-day rotation period, and star B (0.215 M_Sun) has a likely rotation period of 6.7 days. We estimate a probable rotation period of 85 days for LTT 1445A. Thus, this triple M-dwarf system appears to be in a special evolutionary stage where the most massive M dwarf has spun down, the intermediate mass M dwarf is in the process of spinning down, while the least massive stellar component has not yet begun to spin down.

preprint2022arXiv

A stellar occultation by the transneptunian object (50000) Quaoar observed by CHEOPS

Stellar occultation is a powerful technique that allows the determination of some physical parameters of the occulting object. The result depends on the photometric accuracy, the temporal resolution, and the number of chords obtained. Space telescopes can achieve high photometric accuracy as they are not affected by atmospheric scintillation. Using ESA&#39;s CHEOPS space telescope, we observed a stellar occultation by the Transneptunian object (50000) Quaoar. We compare the obtained chord with previous occultations by this object and determine its astrometry with sub-milliarcsecond precision. Also, we determine upper limits to the presence of a global methane atmosphere on the occulting body. We predicted and observed a stellar occultation by Quaoar using the CHEOPS space telescope. We measured the occultation light curve from this data-set and determined the dis- and re-appearance of the star behind the occulting body. Furthermore, a ground-based telescope in Australia was used to constrain Quaoar&#39;s limb. Combined with results from previous works, these measurements allow us to obtain a precise position of Quaoar at the occultation time. We present results obtained from the first stellar occultation by a Transneptunian object (TNO) using space telescope orbiting Earth. It was the occultation by Quaoar observed on 2020 June 11. We used the CHEOPS light curve to obtain a surface pressure upper limit of 85 nbar for the detection of a global methane atmosphere. Also, combining this observation with a ground-based observation we fit Quaoar&#39;s limb to determine its astrometric position with an uncertainty below 1.0 mas. This observation is a first of its kind, and it shall be considered as a proof of concept of stellar occultation observations of Transneptunian objects with space telescopes orbiting Earth. Moreover, it shows significant prospects for the James Webb Space Telescope.

preprint2022arXiv

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

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 &#34;not transit-like&#34; 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.

preprint2022arXiv

Analysis of Early Science observations with the CHaracterising ExOPlanets Satellite (CHEOPS) using pycheops

CHEOPS(CHaracterising ExOPlanet Satellite) is an ESA S-class mission that observes bright stars at high cadence from low-Earth orbit. The main aim of the mission is to characterize exoplanets that transit nearby stars using ultrahigh precision photometry. Here we report the analysis of transits observed by CHEOPS during its Early Science observing programme for four well-known exoplanets: GJ436b, HD106315b, HD97658b and GJ1132b. The analysis is done using pycheops, an open-source software package we have developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques that are fully described herein. We show that the precision of the transit parameters measured using CHEOPS is comparable to that from larger space telescopes such as Spitzer Space Telescope and Kepler. We use the updated planet parameters from our analysis to derive new constraints on the internal structure of these four exoplanets.

preprint2022arXiv

CaRM: Exploring the chromatic Rossiter-McLaughlin effect. The cases of HD 189733b and WASP-127b

In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO.

preprint2022arXiv

CHEOPS geometric albedo of the hot Jupiter HD 209458b

We report the detection of the secondary eclipse of the hot Jupiter HD 209458b in optical/visible light using the CHEOPS space telescope. Our measurement of 20.4 +/- 3.3 ppm translates into a geometric albedo of A_g = 0.096 +/- 0.016. The previously estimated dayside temperature of about 1500 K implies that our geometric albedo measurement consists predominantly of reflected starlight and is largely uncontaminated by thermal emission. This makes the present result one of the most robust measurements of A_g for any exoplanet. Our calculations of the bandpass-integrated geometric albedo demonstrate that the measured value of A_g is consistent with a cloud-free atmosphere, where starlight is reflected via Rayleigh scattering by hydrogen molecules, and the water and sodium abundances are consistent with stellar metallicity. We predict that the bandpass-integrated TESS geometric albedo is too faint to detect and that a phase curve of HD 209458b observed by CHEOPS would have a distinct shape associated with Rayleigh scattering if the atmosphere is indeed cloud free.

preprint2022arXiv

Coralie radial velocity search for companions around evolved stars (CASCADES) II. Seismic masses for three red giants orbited by long-period massive planets

The advent of asteroseismology as the golden path to precisely characterize single stars naturally led to synergies with the field of exoplanetology. Today, the precise determination of stellar masses, radii and ages for exoplanet-host stars is a driving force in the development of dedicated software and techniques to achieve this goal. However, as various approaches exist, it is clear that they all have advantages and inconveniences and that there is a trade-off between accuracy, efficiency, and robustness of the techniques. We aim to compare and discuss various modelling techniques for exoplanet-host red giant stars for which TESS data are available. The results of the seismic modelling are then used to study the dynamical evolution and atmospheric evaporation of the planetary systems. We study, in detail, the robustness, accuracy and precision of various seismic modelling techniques when applied to four exoplanet-host red giants observed by TESS. We discuss the use of global seismic indexes, the use of individual radial frequencies and that of non-radial oscillations. In each case, we discuss the advantages and inconveniences of the modelling technique. We determine precise and accurate masses of exoplanet-host red giant stars orbited by long-period Jupiter-like planets using various modelling techniques. For each target, we also provide a model-independent estimate of the mass from a mean density inversion combined with radii values from Gaia and spectroscopic data. We show that no engulfment or migration is observed for these targets, even if their evolution is extended beyond their estimated seismic ages up the red giant branch.

preprint2022arXiv

Coralie radial-velocity search for companions around evolved stars (CASCADES) III. A new Jupiter host-star: in-depth analysis of HD 29399 using TESS data

Context. Increasing the number of detected exoplanets is far from anecdotal, especially for long-period planets that require a long duration of observation. More detections imply a better understanding of the statistical properties of exoplanet populations, and detailed modelling of their host stars also enables thorough discussions of star-planet interactions and orbital evolution of planetary systems. Aims. In the context of the discovery of a new planetary system, we aim to perform a complete study of HD 29399 and its companion by means of radial-velocity measurements, seismic characterisation of the host-star, and modelling of the orbital evolution of the system. Methods. High-resolution spectra of HD 29399 were acquired with the CORALIE spectrograph mounted on the 1.2- m Swiss telescope located at La Silla Observatory (Chile) as part of the CASCADES survey. We used the moments of the cross-correlation function profile as well as the photometric variability of the star as diagnostics to distinguish between stellar and planetary-induced signals. To model the host star we combined forward modelling with global and local minimisation approaches and inversion techniques. We also studied the orbital history of the system under the effects of both dynamical and equilibrium tides. Results. We present the detection of a long-period giant planet. Combining these measurements with photometric observations by TESS, we are able to thoroughly model the host star and study the orbital evolution of the system. We derive stellar and planetary masses of $1.17 \pm 0.10~ M_{\odot}$ and $1.59 \pm 0.08 ~M_{Jup}$, respectively, and an age for the system of 6.2 Gyr. We show that neither dynamical nor equilibrium tides have been able to affect the orbital evolution of the planet. Moreover, no engulfment is predicted for the future evolution of the system.

preprint2022arXiv

CORALIE radial-velocity search for companions around evolved stars (CASCADES). I. Sample definition and first results: Three new planets orbiting giant stars

Context. Following the first discovery of a planet orbiting a giant star in 2002, we started the CORALIE radial-velocity search for companions around evolved stars (CASCADES). We present the observations of three stars conducted at the 1.2 m Leonard Euler Swiss telescope at La Silla Observatory, Chile, using the CORALIE spectrograph. Aims. We aim to detect planetary companions to intermediate-mass G- and K- type evolved stars and perform a statistical analysis of this population. We searched for new planetary systems around the stars HD22532 (TIC 200851704), HD64121 (TIC 264770836), and HD69123 (TIC 146264536). Methods. We have followed a volume-limited sample of 641 red giants since 2006 to obtain high-precision radial-velocity measurements. We used the Data & Analysis Center for Exoplanets (DACE) platform to perform a radial-velocity analysis to search for periodic signals in the line profile and activity indices, to distinguish between planetary-induced radial-velocity variations and stellar photospheric jitter, and to search for significant signals in the radial-velocity time series to fit a corresponding Keplerian model. Results. In this paper, we present the survey in detail, and we report on the discovery of the first three planets of the sample around the giant stars HD22532, HD64121, and HD69123.

preprint2022arXiv

Detailed stellar activity analysis and modelling of GJ 832: Reassessment of the putative habitable zone planet GJ 832c

Context. Gliese 832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 $\pm$ 9.3 days) is close to the orbital period of putative planet c (35.68 $\pm$ 0.03 days). Aims. We aim to confirm or dismiss the planetary nature of the RV signature attributed to GJ 832c, by adding 119 new RV data points, new photometric data, and an analysis of the spectroscopic stellar activity indicators. Additionally, we update the orbital parameters of the planetary system and search for additional signals. Methods. We performed a frequency content analysis of the RVs to search for periodic and stable signals. Radial velocity time series were modelled with Keplerians and Gaussian process (GP) regressions alongside activity indicators to subsequently compare them within a Bayesian framework. Results. We updated the stellar rotational period of GJ 832 from activity indicators, obtaining $37.5^{+1.4}_{-1.5}$ days, improving the precision by a factor of 6. The new photometric data are in agreement with this value. We detected an RV signal near 18 days (FAP < 4.6%), which is half of the stellar rotation period. Two Keplerians alone fail at modelling GJ 832b and a second planet with a 35-day orbital period. Moreover, the Bayesian evidence from the GP analysis of the RV data with simultaneous activity indices prefers a model without a second Keplerian, therefore negating the existence of planet c.

preprint2022arXiv

Fundamental effective temperature measurements for eclipsing binary stars -- III. SPIRou near-infrared spectroscopy and CHEOPS photometry of the benchmark G0V star EBLM J0113+31

EBLM J0113+31 is moderately bright (V=10.1), metal-poor ([Fe/H]$\approx-0.3$) G0V star with a much fainter M dwarf companion on a wide, eccentric orbit (=14.3 d). We have used near-infrared spectroscopy obtained with the SPIRou spectrograph to measure the semi-amplitude of the M dwarf&#39;s spectroscopic orbit, and high-precision photometry of the eclipse and transit from the CHEOPS and TESS space missions to measure the geometry of this binary system. From the combined analysis of these data together with previously published observations we obtain the following model-independent masses and radii: $M_1 = 1.029 \pm 0.025 M_{\odot}$, $M_2 = 0.197 \pm 0.003 M_{\odot}$, $R_1 = 1.417 \pm 0.014 R_{\odot}$, $R_2 = 0.215 \pm 0.002 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3 we find that this star&#39;s angular diameter is $θ= 0.0745 \pm 0.0007$ mas. The apparent bolometric flux of the G0V star corrected for both extinction and the contribution from the M dwarf ($<0.2$ per cent) is ${\mathcal F}_{\oplus,0} = (2.62\pm 0.05)\times10^{-9}$ erg.cm$^{-2}$.s$^{-1}$. Hence, this G0V star has an effective temperature $T_{\rm eff,1} = 6124{\rm\,K} \pm 40{\rm \,K\,(rnd.)} \pm 10 {\rm \,K\,(sys.)}$. EBLM J0113+31 is an ideal benchmark star that can be used for &#34;end-to-end&#34; tests of the stellar parameters measured by large-scale spectroscopic surveys, or stellar parameters derived from asteroseismology with PLATO. The techniques developed here can be applied to many other eclipsing binaries in order to create a network of such benchmark stars.

preprint2022arXiv

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N and TESS

We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($ρ_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies.

preprint2022arXiv

New binaries from the SHINE survey

We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose.

preprint2022arXiv

Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection

The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $Ψ$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $μ$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $Ψ=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations, we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we could separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $Ψ$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved JWST phase curves.

preprint2022arXiv

The atmosphere and architecture of WASP-189 b probed by its CHEOPS phase curve

Gas giants orbiting close to hot and massive early-type stars can reach dayside temperatures that are comparable to those of the coldest stars. These &#34;ultra-hot Jupiters&#34; have atmospheres made of ions and atomic species from molecular dissociation and feature strong day-to-night temperature gradients. Photometric observations at different orbital phases provide insights on the planet atmospheric properties. We analyse the photometric observations of WASP-189 acquired with the instrument CHEOPS to derive constraints on the system architecture and the planetary atmosphere. We implement a light curve model suited for asymmetric transit shape caused by the gravity-darkened photosphere of the fast-rotating host star. We also model the reflective and thermal components of the planetary flux, the effect of stellar oblateness and light-travel time on transit-eclipse timings, the stellar activity and CHEOPS systematics. From the asymmetric transit, we measure the size of the ultra-hot Jupiter WASP-189 b, $R_p=1.600^{+0.017}_{-0.016}\,R_J$, with a precision of 1%, and the true orbital obliquity of the planetary system $Ψ_p=89.6\pm1.2°$ (polar orbit). We detect no significant hotspot offset from the phase curve and obtain an eclipse depth $δ_\text{ecl}=96.5^{+4.5}_{-5.0}\,\text{ppm}$, from which we derive an upper limit on the geometric albedo: $A_g<0.48$. We also find that the eclipse depth can only be explained by thermal emission alone in the case of extremely inefficient energy redistribution. Finally, we attribute the photometric variability to the stellar rotation, either through superficial inhomogeneities or resonance couplings between the convective core and the radiative envelope.

preprint2022arXiv

The young HD 73583 (TOI-560) planetary system: Two 10-M$_\oplus$ mini-Neptunes transiting a 500-Myr-old, bright, and active K dwarf

We present the discovery and characterisation of two transiting planets observed by \textit{TESS} in the light curves of the young and bright (V=9.67) star HD73583 (TOI-560). We perform an intensive spectroscopic and photometric space- and ground-based follow-up in order to confirm and characterise the system. We found that HD73583 is a young ($\sim 500$~Myr) active star with a rotational period of $12.08 \pm 0.11 $\,d, and a mass and radius of $ 0.73 \pm 0.02 M_\odot$ and $0.65 \pm 0.02 R_\odot$, respectively. HD73583 b ($P_b=6.3980420 _{ - 0.0000062 }^{+0.0000067}$ d) has a mass and radius of $10.2 _{-3.1}^{+3.4} M_\oplus$ and$2.79 \pm 0.10 R_\oplus$, respectively, that gives a density of $2.58 _{-0.81}^{ 0.95} {\rm g\,cm^{-3}}$. HD73583 c ($P_c= 18.87974 _{-0.00074 }^{+0.00086}$) has a mass and radius of $9.7_{-1.7} ^ {+1.8} M_\oplus$ and $2.39_{-0.09}^{+0.10} R_\oplus$, respectively, this translates to a density of $3.88 _{-0.80}^{+0.91} {\rm g\,cm^{-3}}$. Both planets are consistent with worlds made of a solid core surrounded by a volatile envelope. Because of their youth and host star brightness, they both are excellent candidates to perform transmission spectroscopy studies. We expect ongoing atmospheric mass-loss for both planets caused by stellar irradiation. We estimate that the detection of evaporating signatures on H and He would be challenging, but doable with present and future instruments.

preprint2022arXiv

Transit timing variations of AU Microscopii b and c

Here we report large-amplitude transit timing variations (TTVs) for AU\,Microcopii b and c as detected in combined TESS (2018, 2020) and CHEOPS (2020, 2021) transit observations. AU Mic is a young planetary system with a debris disk and two transiting warm Neptunes. A TTV on the order of several minutes was previously reported for AU Mic b, which was suggested to be an outcome of mutual perturbations between the planets in the system. In 2021, we observed AU Mic b (five transits) and c (three transits) with the CHEOPS space telescope to follow-up the TTV of AU Mic b and possibly detect a TTV for AU Mic c. When analyzing TESS and CHEOPS 2020--2021 measurements together, we find that a prominent TTV emerges with a full span of ~23 minutes between the two TTV extrema. Assuming that the period change results from a periodic process - such as mutual perturbations - we demonstrate that the times of transits in the summer of 2022 are expected to be 30--85 minutes later than predicted by the available linear ephemeris.

preprint2021arXiv

A hot mini-Neptune in the radius valley orbiting solar analogue HD 110113

We report the discovery of HD 110113 b (TOI-755.01), a transiting mini-Neptune exoplanet on a 2.5-day orbit around the solar-analogue HD 110113 (Teff = 5730K). Using TESS photometry and HARPS radial velocities gathered by the NCORES program, we find HD 110113 b has a radius of $2.05\pm0.12$ $R_\oplus$ and a mass of $4.55\pm0.62$ $M_\oplus$. The resulting density of $2.90^{+0.75}_{-0.59}$ g cm^{-3} is significantly lower than would be expected from a pure-rock world; therefore, HD 110113 b must be a mini-Neptune with a significant volatile atmosphere. The high incident flux places it within the so-called radius valley; however, HD 110113 b was able to hold onto a substantial (0.1-1\%) H-He atmosphere over its $\sim4$ Gyr lifetime. Through a novel simultaneous gaussian process fit to multiple activity indicators, we were also able to fit for the strong stellar rotation signal with period $20.8\pm1.2$ d from the RVs and confirm an additional non-transiting planet with a mass of $10.5\pm1.2$ $M_\oplus$ and a period of $6.744^{+0.008}_{-0.009}$ d.

preprint2021arXiv

CHEOPS observations of the HD 108236 planetary system: A fifth planet, improved ephemerides, and planetary radii

The detection of a super-Earth and three mini-Neptunes transiting the bright ($V$ = 9.2 mag) star HD 108236 (also known as TOI-1233) was recently reported on the basis of TESS and ground-based light curves. We perform a first characterisation of the HD 108236 planetary system through high-precision CHEOPS photometry and improve the transit ephemerides and system parameters. We characterise the host star through spectroscopic analysis and derive the radius with the infrared flux method. We constrain the stellar mass and age by combining the results obtained from two sets of stellar evolutionary tracks. We analyse the available TESS light curves and one CHEOPS transit light curve for each known planet in the system. We find that HD 108236 is a Sun-like star with $R_{\star}=0.877\pm0.008 R_{\odot}$, $M_{\star}=0.869^{+0.050}_{-0.048} M_{\odot}$, and an age of $6.7_{-5.1}^{+4.0}$ Gyr. We report the serendipitous detection of an additional planet, HD 108236 f, in one of the CHEOPS light curves. For this planet, the combined analysis of the TESS and CHEOPS light curves leads to a tentative orbital period of about 29.5 days. From the light curve analysis, we obtain radii of $1.615\pm0.051$, $2.071\pm0.052$, $2.539_{-0.065}^{+0.062}$, $3.083\pm0.052$, and $2.017_{-0.057}^{+0.052}$ $R_{\oplus}$ for planets HD 108236 b to HD 108236 f, respectively. These values are in agreement with previous TESS-based estimates, but with an improved precision of about a factor of two. We perform a stability analysis of the system, concluding that the planetary orbits most likely have eccentricities smaller than 0.1. We also employ a planetary atmospheric evolution framework to constrain the masses of the five planets, concluding that HD 108236 b and HD 108236 c should have an Earth-like density, while the outer planets should host a low mean molecular weight envelope.

preprint2021arXiv

First light of the VLT planet finder SPHERE. III. New spectrophotometry and astrometry of the HR8799 exoplanetary system

The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $μ$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $μ$m, 1.667 $μ$m), $K1K2$ (2.110 $μ$m, 2.251 $μ$m), and broadband $J$ (1.245 $μ$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $μ$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Loève image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

preprint2021arXiv

HD207897 b: A dense sub-Neptune transiting a nearby and bright K-type star

We present the discovery and characterization of a transiting sub-Neptune orbiting with a 16.20 day period around a nearby (28 pc) and bright(V=8.37) K0V star HD207897 (TOI-1611). This discovery is based on photometric measurements from the Transiting Exoplanet Survey Satellite(TESS) mission and radial velocity (RV) observations from the SOPHIE, Automated Planet Finder (APF) and HIRES high precision spectrographs. We used EXOFASTv2 for simultaneously modeling the parameters of the planet and its host star, combining photometric and RV data to determine the planetary system parameters. We show that the planet has a radius of 2.50+/-0.08 RE and a mass of either 14.4+/-1.6 ME or 15.9+/-1.6 ME with nearly equal probability; the two solutions correspond to two possibilities for the stellar activity period. Hence, the density is either 5.1+/-0.7 g cm^-3 or 5.5^{+0.8}_{-0.7} g cm^-3, making it one of the relatively rare dense sub-Neptunes. The existence of such a dense planet at only 0.12 AU from its host star is unusual in the currently observed sub-Neptune (2<RE<4) population. The most likely scenario is that this planet has migrated to its current position.

preprint2021arXiv

NGTS-14Ab: a Neptune-sized transiting planet in the desert

Context: The sub-Jovian or Neptunian desert is a previously-identified region of parameter space where there is a relative dearth of intermediate-mass planets at short orbital periods. Aims: We present the discovery of a new transiting planetary system within the Neptunian desert, NGTS-14. Methods: Transits of NGTS-14Ab were discovered in photometry from the Next Generation Transit Survey (NGTS). Follow-up transit photometry was conducted from several ground-based facilities, as well as extracted from TESS full-frame images. We combine radial velocities from the HARPS spectrograph with the photometry in a global analysis to determine the system parameters. Results: NGTS-14Ab has a radius about 30 per cent larger than that of Neptune ($0.444\pm0.030~\mathrm{R_{Jup}}$), and is around 70 per cent more massive than Neptune ($0.092 \pm 0.012~\mathrm{M_{Jup}}$). It transits the main-sequence K1 star, NGTS-14A, with a period of 3.54 days, just far enough to have maintained at least some of its primordial atmosphere. We have also identified a possible long-period stellar mass companion to the system, NGTS-14B, and we investigate the binarity of exoplanet host stars inside and outside the Neptunian desert using Gaia.

preprint2021arXiv

Retrieving the transmission spectrum of HD 209458b using CHOCOLATE: A new chromatic Doppler tomography technique

Multiband photometric transit observations or low-resolution spectroscopy (spectro-photometry) are normally used to retrieve the broadband transmission spectra of transiting exoplanets in order to assess the chemical composition of their atmospheres. In this paper, we present an alternative approach for recovering the broadband transmission spectra using chromatic Doppler tomography based on physical modeling through the SOAP tool: CHOCOLATE (CHrOmatiC line prOfiLe tomogrAphy TEchnique). To validate the method and examine its performance, we use observational data recently obtained with the ESPRESSO instrument to retrieve the transmission spectra of the archetypal hot Jupiter HD 209458b. Our findings indicate that the recovered transmission spectrum is in good agreement with the results presented in previous studies, which used different methodologies to extract the spectrum, achieving similar precision. We explored several atmospheric models and inferred from spectral retrieval that a model containing H2O and NH3 is the preferred scenario. The CHOCOLATE methodology is particularly interesting for future studies of exoplanets around young and active stars or moderate to fast rotating stars, considering SOAP&#39;s ability to model stellar active regions and the fact that the rotational broadening of spectral lines favors its application. Furthermore, CHOCOLATE will allow the broad transmission spectrum of a planet to be retrieved using high signal-to-noise ratio, high-resolution spectroscopy with the next generation of Extremely Large Telescopes (ELTs), where low-resolution spectroscopy will not always be accessible.

preprint2021arXiv

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

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&#39;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.

preprint2021arXiv

The atmosphere of HD 209458b seen with ESPRESSO. No detectable planetary absorptions at high resolution

We observed two transits of the iconic gas giant HD 209458b between 380 and 780 nm, using the high-resolution ESPRESSO spectrograph. The derived planetary transmission spectrum exhibits features at all wavelengths where the parent star shows strong absorption lines, for example, NaI, MgI, FeI, FeII, CaI, VI, H$α$, and KI. We interpreted these features as the signature of the deformation of the stellar line profiles due to the Rossiter-McLaughlin effect, combined with the centre-to-limb effects on the stellar surface, which is in agreement with similar reports recently presented in the literature. We also searched for species that might be present in the planetary atmosphere but not in the stellar spectra, such as TiO and VO, and obtained a negative result. Thus, we find no evidence of any planetary absorption, including previously reported NaI, in the atmosphere of HD 209458b. The high signal-to-noise ratio in the transmission spectrum allows us to compare the modelled deformation of the stellar lines in assuming different one-dimensional stellar atmospheric models. We conclude that the differences among various models and observations remain within the precision limits of the data. However, the transmission light curves are better explained when the centre-to-limb variation is not included in the computation and only the Rossiter-McLaughlin deformation is considered. This demonstrates that ESPRESSO is currently the best facility for spatially resolving the stellar surface spectrum in the optical range using transit observations and carrying out empirical validations of stellar models.

preprint2021arXiv

The HARPS search for southern extra-solar planets XLV. Two Neptune mass planets orbiting HD 13808: a study of stellar activity modelling&#39;s impact on planet detection

We present a comprehensive analysis of 10 years of HARPS radial velocities of the K2V dwarf star HD 13808, which has previously been reported to host two unconfirmed planet candidates. We use the state-of-the-art nested sampling algorithm PolyChord to compare a wide variety of stellar activity models, including simple models exploiting linear correlations between RVs and stellar activity indicators, harmonic models for the activity signals, and a more sophisticated Gaussian process regression model. We show that the use of overly-simplistic stellar activity models that are not well-motivated physically can lead to spurious `detections&#39; of planetary signals that are almost certainly not real. We also reveal some difficulties inherent in parameter and model inference in cases where multiple planetary signals may be present. Our study thus underlines the importance both of exploring a variety of competing models and of understanding the limitations and precision settings of one&#39;s sampling algorithm. We also show that at least in the case of HD 13808, we always arrive at consistent conclusions about two particular signals present in the RV, regardless of the stellar activity model we adopt; these two signals correspond to the previously-reported though unconfirmed planet candidate signals. Given the robustness and precision with which we can characterize these two signals, we deem them secure planet detections. In particular, we find two planets orbiting HD 13808 at distances of 0.11, 0.26 AU with periods of 14.2, 53.8 d, and minimum masses of 11, 10 Earth masses.

preprint2021arXiv

The HD 137496 system: A dense, hot super-Mercury and a cold Jupiter

Most of the currently known planets are small worlds with radii between that of the Earth and that of Neptune. The characterization of planets in this regime shows a large diversity in compositions and system architectures, with distributions hinting at a multitude of formation and evolution scenarios. Using photometry from the K2 satellite and radial velocities measured with the HARPS and CORALIE spectrographs, we searched for planets around the bright and slightly evolved Sun-like star HD 137496. We precisely estimated the stellar parameters, $M_*$ = 1.035 +/- 0.022 $M_\odot$, $R_*$ = 1.587 +/- 0.028 $R_\odot$, $T_\text{eff}$ = 5799 +/- 61 K, together with the chemical composition of the slightly evolved star. We detect two planets orbiting HD 137496. The inner planet, HD 137496 b, is a super-Mercury (an Earth-sized planet with the density of Mercury) with a mass of $M_b$ = 4.04 +/- 0.55 $M_\oplus$, a radius of $R_b = 1.31_{-0.05}^{+0.06} R_\oplus,$ and a density of $ρ_b = 10.49_{-1.82}^{+2.08}$ $\mathrm{g cm^{-3}}$. With an interior modeling analysis, we find that the planet is composed mainly of iron, with the core representing over 70% of the planet&#39;s mass ($M_{core}/M_{total} = 0.73^{+0.11}_{-0.12}$). The outer planet, HD 137496 c, is an eccentric ($e$ = 0.477 +/- 0.004), long period ($P$ = $479.9_{-1.1}^{+1.0}$ days) giant planet ($M_c\sin i_c$ = 7.66 +/- 0.11 $M_{Jup}$) for which we do not detect a transit. HD 137496 b is one of the few super-Mercuries detected to date. The accurate characterization reported here enhances its role as a key target to better understand the formation and evolution of planetary systems. The detection of an eccentric long period giant companion also reinforces the link between the presence of small transiting inner planets and long period gas giants.

preprint2021arXiv

The PLATO field selection process I. Identification and content of the long-pointing fields

PLATO (PLAnetary Transits and Oscillations of stars) is an ESA M-class satellite planned for launch by end 2026 and dedicated to the wide-field search of transiting planets around bright and nearby stars, with a strong focus on discovering habitable rocky planets hosted by solar-like stars. The choice of the fields to be pointed at is a crucial task since it has a direct impact on the scientific return of the mission. In this paper we describe and discuss the formal requirements and the key scientific prioritization criteria that have to be taken into account in the Long-duration Observation Phase (LOP) field selection, and apply a quantitative metric to guide us in this complex optimization process. We identify two provisional LOP fields, one for each hemisphere (LOPS1, LOPN1), and discuss their properties and stellar content. While additional fine-tuning shall be applied to LOP selection before the definitive choice (to be made two years before launch), we expect their position will not move by more than a few degrees with respect to what is proposed in this paper.

preprint2020arXiv

A hot terrestrial planet orbiting the bright M dwarf L 168-9 unveiled by TESS

We report the detection of a transiting super-Earth-sized planet (R=1.39+-0.09 Rearth) in a 1.4-day orbit around L 168-9 (TOI-134),a bright M1V dwarf (V=11, K=7.1) located at 25.15+-0.02 pc. The host star was observed in the first sector of the Transiting Exoplanet Survey Satellite (TESS) mission and, for confirmation and planet mass measurement, was followed up with ground-based photometry, seeing-limited and high-resolution imaging, and precise radial velocity (PRV) observations using the HARPS and PFS spectrographs. Combining the TESS data and PRV observations, we find the mass of L168-9 b to be 4.60+-0.56 Mearth, and thus the bulk density to be 1.74+0.44-0.33 times larger than that of the Earth. The orbital eccentricity is smaller than 0.21 (95% confidence). This planet is a Level One Candidate for the TESS Mission&#39;s scientific objective - to measure the masses of 50 small planets - and is one of the most observationally accessible terrestrial planets for future atmospheric characterization.

preprint2020arXiv

An ultra-short period rocky super-Earth orbiting the G2-star HD 80653

Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise masses and radii, which also includes estimates of the level of irradiation and information about possible companions. Here we report a detailed characterization of a USP planet around the solar-type star HD 80653 $\equiv$ EP 251279430 using the K2 light curve and 108 precise radial velocities obtained with the HARPS-N spectrograph, installed on the Telescopio Nazionale Galileo. From the K2 C16 data, we found one super-Earth planet ($R_{b}=1.613\pm0.071 R_{\oplus}$) transiting the star on a short-period orbit ($P_{\rm b}=0.719573\pm0.000021$ d). From our radial velocity measurements, we constrained the mass of HD 80653 b to $M_{b}=5.60\pm0.43 M_{\oplus}$. We also detected a clear long-term trend in the radial velocity data. We derived the fundamental stellar parameters and determined a radius of $R_{\star}=1.22\pm0.01 R_{\odot}$ and mass of $M_{\star}=1.18\pm0.04 M_{\odot}$, suggesting that HD 80653, has an age of $2.7\pm1.2$ Gyr. The bulk density ($ρ_{b} = 7.4 \pm 1.1$ g cm$^{-3}$) of the planet is consistent with an Earth-like composition of rock and iron with no thick atmosphere. Our analysis of the K2 photometry also suggests hints of a shallow secondary eclipse with a depth of 8.1$\pm$3.7 ppm. Flux variations along the orbital phase are consistent with zero. The most important contribution might come from the day-side thermal emission from the surface of the planet at $T\sim3480$ K.

preprint2020arXiv

Atmospheric Rossiter-McLaughlin effect and transmission spectroscopy of WASP-121b with ESPRESSO

WASP-121b is one of the most studied Ultra-hot Jupiters: many recent analyses of its atmosphere report interesting features at different wavelength ranges. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in 1-telescope mode, and one partial transit taken during the commissioning of the instrument in 4-telescope mode. We investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. By analysing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter-McLaughlin effect. We measured the height of the planetary atmospheric layer that correlates with the stellar mask (mainly Fe) to be 1.052$\pm$0.015 Rp and we also confirmed the blueshift of the planetary atmosphere. By examining the planetary absorption signal on the stellar cross-correlation functions we confirmed the presence of a temporal variation of its blueshift during transit, which could be investigated spectrum-by-spectrum. We detected significant absorption in the transmission spectrum for Na, H, K, Li, CaII, Mg, and we certified their planetary nature by using the 2D tomographic technique. Particularly remarkable is the detection of Li, with a line contrast of $\sim$0.2% detected at the 6$σ$ level. With the cross-correlation technique we confirmed the presence of FeI, FeII, CrI and VI. H$α$ and CaII are present up to very high altitudes in the atmosphere ($\sim$1.44 Rp and $\sim$2 Rp, respectively), and also extend beyond the transit-equivalent Roche lobe radius of the planet. These layers of the atmosphere have a large line broadening that is not compatible with being caused by the tidally-locked rotation of the planet alone, and could arise from vertical winds or high-altitude jets in the evaporating atmosphere.

preprint2020arXiv

Broadband transmission spectroscopy of HD209458b with ESPRESSO: Evidence for Na, TiO, or both

The detection and characterization of exoplanet atmospheres is currently one of the main drivers pushing the development of new observing facilities. In this context, high-resolution spectrographs are {proving} their potential and showing that high-resolution spectroscopy will be paramount in this field. We aim to make use of ESPRESSO high-resolution spectra, which cover two transits of HD209458b, to probe the broadband transmission optical spectrum of the planet. We applied the chromatic Rossiter-McLaughin method to derive the transmission spectrum of HD209458b. We compared the results with previous HST observations and with synthetic spectra. We recover a transmission spectrum of HD209458b similar to the one obtained with HST data. The models suggest that the observed signal can be explained by only Na, only TiO, or both Na and TiO, even though none is fully capable of explaining our observed transmission spectrum. Extra absorbers may be needed to explain the full dataset, though modeling approximations and observational errors can also be responsible for the observed mismatch. Using the chromatic Rossiter-McLaughlin technique, ESPRESSO is able to provide broadband transmission spectra of exoplanets from the ground, in conjunction with space-based facilities, opening good perspectives for similar studies of other planets.

preprint2020arXiv

Discovery and characterization of the exoplanets WASP-148b and c. A transiting system with two interacting giant planets

We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric survey, then characterized it with the SOPHIE spectrograph. Our radial velocity measurements subsequently revealed a second planet in the system, WASP-148c, with an orbital period of 34.5 days and a minimum mass of 0.40 M_Jup. No transits of this outer planet were detected. The orbits of both planets are eccentric and fall near the 4:1 mean-motion resonances. This configuration is stable on long timescales, but induces dynamical interactions so that the orbits differ slightly from purely Keplerian orbits. In particular, WASP-148b shows transit-timing variations of typically 15 minutes, making it the first interacting system with transit-timing variations that is detected on ground-based light curves. We establish that the mutual inclination of the orbital plane of the two planets cannot be higher than 35 degrees, and the true mass of WASP-148c is below 0.60 M_Jup. We present photometric and spectroscopic observations of this system that cover a time span of ten years. We also provide their Keplerian and Newtonian analyses; these analyses should be significantly improved through future TESS~observations.

preprint2020arXiv

ESPRESSO high resolution transmission spectroscopy of WASP-76b

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b.

preprint2020arXiv

ESPRESSO@VLT -- On-sky performance and first results

ESPRESSO is the new high-resolution spectrograph of ESO&#39;s Very-Large Telescope (VLT). It was designed for ultra-high radial-velocity precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UT) of the VLT at a spectral resolving power of 140,000 or 190,000 over the 378.2 to 788.7 nm wavelength range, or with all UTs together, turning the VLT into a 16-m diameter equivalent telescope in terms of collecting area, while still providing a resolving power of 70,000. We provide a general description of the ESPRESSO instrument, report on the actual on-sky performance, and present our Guaranteed-Time Observation (GTO) program with its first results. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1st, 2018, but improvements to the instrument and re-commissioning runs were conducted until July 2019. The measured overall optical throughput of ESPRESSO at 550 nm and a seeing of 0.65 arcsec exceeds the 10% mark under nominal astro-climatic conditions. We demonstrate a radial-velocity precision of better than 25 cm/s during one night and 50 cm/s over several months. These values being limited by photon noise and stellar jitter show that the performanceis compatible with an instrumental precision of 10 cm/s. No difference has been measured across the UTs neither in throughput nor RV precision. The combination of the large collecting telescope area with the efficiency and the exquisite spectral fidelity of ESPRESSO opens a new parameter space in RV measurements, the study of planetary atmospheres, fundamental constants, stellar characterisation and many other fields.

preprint2020arXiv

HATS-47b, HATS-48Ab, HATS-49b and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of $0.369_{-0.021}^{+0.031}$ $M_{J}$, $0.243_{-0.030}^{+0.022}$ $M_{J}$, $0.353_{-0.027}^{+0.038}$ $M_{J}$ and $0.1254\pm0.0039$ $M_{J}$, respectively, and radii of $1.117\pm0.014$ $R_{J}$, $0.800\pm0.015$ $R_{J}$, $0.765\pm0.013$ $R_{J}$, and $0.7224\pm0.0032$ $R_{J}$, respectively. The planets orbit close to their host stars with orbital periods of $3.9228$ d, $3.1317$ d, $4.1480$ d and $7.3279$ d, respectively. The hosts are main sequence K dwarfs with masses of $0.674_{-0.012}^{+0.016}$ $M_{\odot}$, $0.7279\pm0.0066$ $M_{\odot}$, $0.7133\pm0.0075$ $M_{\odot}$, and $0.7311\pm0.0028$ $M_{\odot}$ and with $V$-band magnitudes of $V = 14.829\pm0.010$, $14.35\pm0.11$, $14.998\pm0.040$ and $12.469\pm0.010$. The Super-Neptune HATS-72b (a.k.a.\ WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS and Magellan/PFS. We also incorporate TESS observations of the warm Saturn-hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved $5.\!\!^{\prime\prime}4$ neighbor in Gaia~DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of $0.22$ $M_{\odot}$ and a current projected physical separation of $\sim$1,400 au.

preprint2020arXiv

Hot Exoplanet Atmospheres Resolved with Transit Spectroscopy (HEARTS) III. Atmospheric structure of the misaligned ultra-hot Jupiter WASP-121b

Ultra-hot Jupiters offer interesting prospects for expanding our theories on dynamical evolution and the properties of extremely irradiated atmospheres. In this context, we present the analysis of new optical spectroscopy for the transiting ultra-hot Jupiter WASP-121b. We first refine the orbital properties of WASP-121b, which is on a nearly polar (obliquity $ψ^{\rm North}$=88.1$\pm$0.25$^{\circ}$ or $ψ^{\rm South}$=91.11$\pm$0.20$^{\circ}$) orbit, and exclude a high differential rotation for its fast-rotating (P$<$1.13 days), highly inclined ($i_\mathrm{\star}^{\rm North}$=8.1$\stackrel{+3.0}{_{-2.6}}^{\circ}$ or $i_\mathrm{\star}^{\rm South}$=171.9$\stackrel{+2.5}{_{-3.4}}^{\circ}$) star. We then present a new method that exploits the reloaded Rossiter-McLaughlin technique to separate the contribution of the planetary atmosphere and of the spectrum of the stellar surface along the transit chord. Its application to HARPS transit spectroscopy of WASP-121b reveals the absorption signature from metals, likely atomic iron, in the planet atmospheric limb. The width of the signal (14.3$\pm$1.2 km/s) can be explained by the rotation of the tidally locked planet. Its blueshift (-5.2$\pm$0.5 km/s) could trace strong winds from the dayside to the nightside, or the anisotropic expansion of the planetary thermosphere.

preprint2020arXiv

Mass determinations of the three mini-Neptunes transiting TOI-125

The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission&#39;s primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS&#39;s observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2σ$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system.

preprint2020arXiv

Masses for the seven planets in K2-32 and K2-233. Four diverse planets in resonant chain and the first young rocky worlds

High-precision planetary densities are key to derive robust atmospheric properties for extrasolar planets. Measuring precise masses is the most challenging part, especially in multi-planetary systems. We measure the masses and densities of a four-planet near resonant chain system (K2-32), and a young ($\sim400$ Myr old) planetary system consisting of three close-in small planets (K2-233). We obtained 199 new HARPS observations for K2-32 and 124 for K2-233 covering a more than three year baseline. We find that K2-32 is a compact scaled-down version of the Solar System&#39;s architecture, with a small rocky inner planet (M$_e=2.1^{+1.3}_{-1.1}$~M$_{\oplus}$, P$_e\sim4.35$~days) followed by an inflated Neptune-mass planet (M$_b=15.0^{+1.8}_{-1.7}$~M$_{\oplus}$, P$_b\sim8.99$~days) and two external sub-Neptunes (M$_c=8.1\pm2.4$~M$_{\oplus}$, P$_c\sim20.66$~days; M$_d=6.7\pm2.5$~M$_{\oplus}$, P$_d\sim31.72$~days). K2-32 becomes one of the few multi-planetary systems with four or more planets known with measured masses and radii. Additionally, we constrain the masses of the three planets in K2-233. For the two inner Earth-size planets we constrain their masses to be smaller than M$_b<11.3$ M$_{\oplus}$ (P$_b\sim2.47$~days), M$_c<12.8$ M$_{\oplus}$ (P$_c\sim7.06$~days). The outer planet is a sub-Neptune size planet with an inferred mass of M$_d=8.3^{+5.2}_{-4.7}$ M$_{\oplus}$ (M$_d<21.1$ M$_{\oplus}$, P$_d\sim24.36$~days). Our observations of these two planetary systems confirm for the first time the rocky nature of two planets orbiting a young star, with relatively short orbital periods ($<7$ days). They provide key information for planet formation and evolution models of telluric planets. Additionally, the Neptune-like derived masses of the three planets K2-32 b, c, d puts them in a relatively unexplored regime of incident flux and planet mass, key for transmission spectroscopy studies.

preprint2020arXiv

Orbital and spectral characterization of the benchmark T-type brown dwarf HD 19467B

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

preprint2020arXiv

Revisiting Proxima with ESPRESSO

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

preprint2020arXiv

Spectral and atmospheric characterisation of a new benchmark brown dwarf HD13724B

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.

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

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&#39;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.

preprint2020arXiv

SPIRou: nIR velocimetry & spectropolarimetry at the CFHT

This paper presents an overview of SPIRou, the new-generation near-infrared spectropolarimeter / precision velocimeter recently installed on the 3.6-m Canada-France-Hawaii Telescope (CFHT). Starting from the two main science goals, namely the quest for planetary systems around nearby M dwarfs and the study of magnetized star / planet formation, we outline the instrument concept that was designed to efficiently address these forefront topics, and detail the in-lab and on-sky instrument performances measured throughout the intensive testing phase that SPIRou was submitted to before passing the final acceptance review in early 2019 and initiating science observations. With a central position among the newly started programmes, the SPIRou Legacy Survey (SLS) Large Programme was allocated 300 CFHT nights until at least mid 2022. We also briefly describe a few of the first results obtained in the various science topics that SPIRou started investigating, focusing in particular on planetary systems of nearby M dwarfs, transiting exoplanets and their atmospheres, magnetic fields of young stars, but also on alternate science goals like the atmospheres of M dwarfs and the Earth&#39;s atmosphere. We finally conclude on the essential role that SPIRou and the CFHT can play in coordination with forthcoming major facilities like the JWST, the ELTs, PLATO and ARIEL over the decade.

preprint2020arXiv

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

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.

preprint2020arXiv

Three Short Period Jupiters from TESS

We report the confirmation and mass determination of three hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HIP 65Ab (TOI-129, TIC-201248411) is an ultra-short-period Jupiter orbiting a bright (V=11.1 mag) K4-dwarf every 0.98 days. It is a massive 3.213 +/- 0.078 Mjup planet in a grazing transit configuration with an impact parameter of b = 1.17 +0.10/-0.08. As a result the radius is poorly constrained, 2.03 +0.61/-0.49 Rjup. The planet&#39;s distance to its host star is less than twice the separation at which it would be destroyed by Roche lobe overflow. It is expected to spiral into HIP 65A on a timescale ranging from 80 Myr to a few gigayears, assuming a reduced tidal dissipation quality factor of Qs&#39; = 10^7 - 10^9. We performed a full phase-curve analysis of the TESS data and detected both illumination- and ellipsoidal variations as well as Doppler boosting. HIP 65A is part of a binary stellar system, with HIP 65B separated by 269 AU (3.95 arcsec on sky). TOI-157b (TIC 140691463) is a typical hot Jupiter with a mass of 1.18 +/- 0.13 Mjup and a radius of 1.29 +/- 0.02 Rjup. It has a period of 2.08 days, which corresponds to a separation of just 0.03 AU. This makes TOI-157 an interesting system, as the host star is an evolved G9 sub-giant star (V=12.7). TOI-169b (TIC 183120439) is a bloated Jupiter orbiting a V=12.4 G-type star. It has a mass of 0.79 +/- 0.06 Mjup and a radius of 1.09 +0.08/-0.05 Rjup. Despite having the longest orbital period (P = 2.26 days) of the three planets, TOI-169b receives the most irradiation and is situated on the edge of the Neptune desert. All three host stars are metal rich with [Fe/H] ranging from 0.18 - 0.24.

preprint2020arXiv

WASP-127b: A misaligned planet with a partly cloudy atmosphere and tenuous sodium signature seen by ESPRESSO

The study of exoplanet atmospheres is essential to understand the formation, evolution and composition of exoplanets. The transmission spectroscopy technique is playing a significant role in this domain. In particular, the combination of state-of-the-art spectrographs at low- and high-spectral resolution is key to our understanding of atmospheric structure and composition. Two transits of the close-in sub Saturn-mass planet,WASP-127b, have been observed with ESPRESSO in the frame of the Guaranteed Time Observations Consortium. Transit observations allow us to study simultaneously the system architecture and the exoplanet atmosphere. We found that this planet is orbiting its slowly rotating host star (veq sin(i)=0.53+/-0.07 km/s) on a retrograde misaligned orbit (lambda=-128.41+/-5.60 deg). We detected the sodium line core at the 9-sigma confidence level with an excess absorption of 0.3+/-0.04%, a blueshift of 2.7+/-0.79 km/s and a FWHM of 15.18+/-1.75 km/s. However, we did not detect the presence of other atomic species but set upper-limits of only few scale heights. Finally, we put a 3-sigma upper limit, to the average depth of the 1600 strongest water lines at equilibrium temperature in the visible band, of 38 ppm. This constrains the cloud-deck pressure between 0.3 and 0.5 mbar by combining our data with low-resolution data in the near-infrared and models computed for this planet. To conclude, WASP-127b, with an age of about 10 Gyr, is an unexpected exoplanet by its orbital architecture but also by the small extension of its sodium atmosphere (~7 scale heights). ESPRESSO allows us to take a step forward in the detection of weak signals, thus bringing strong constraints on the presence of clouds in exoplanet atmospheres. The framework proposed in this work can be applied to search for molecular species and study cloud-decks in other exoplanets.

preprint2019arXiv

Optical phase curve of the ultra-hot Jupiter WASP-121b

We present the analysis of TESS optical photometry of WASP-121b, which reveal the phase curve of this transiting ultra-hot Jupiter. Its hotspot is located at the substellar point, showing inefficient heat transport from the dayside (2870 K) to the nightside ($<$ 2200 K) at the altitudes probed by TESS. The TESS eclipse depth, measured at the shortest wavelength to date for WASP-121b, confirms the strong deviation from blackbody planetary emission. Our atmospheric retrieval on the complete emission spectrum supports the presence of a temperature inversion, which can be explained by the presence of VO and possibly TiO and FeH. The strong planetary emission at short wavelengths could arise from an H$^{-}$ continuum.

preprint2019arXiv

RefPlanets: Search for reflected light from extra-solar planets with SPHERE/ZIMPOL

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6&#34;. We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13&#34;, 16.3-18.3 mag at 0.5&#34; and 18.8-20.4 mag at 1.5&#34;. The contrast limits closer in (<0.6&#34;) depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6&#34;), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS I: Description, data reduction and observing strategy

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS II: Characterization and correction of instrumental polarization effects

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE&#39;s internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE&#39;s first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.