Researcher profile

D. Ehrenreich

D. Ehrenreich contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
4topics
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

31 published item(s)

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 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 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

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

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

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

Large Interferometer For Exoplanets (LIFE): I. Improved exoplanet detection yield estimates for a large mid-infrared space-interferometer mission

One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measures the thermal emission of exoplanets. For this, we have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo simulations of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Within a 2.5-year initial search phase, an interferometer consisting of four 2 m apertures with a total instrument throughput of 5% covering a wavelength range between 4 and 18.5 $μ$m could detect up to ~550 exoplanets with radii between 0.5 and 6 R$_\oplus$ with an integrated SNR$\ge$7. At least ~160 of the detected exoplanets have radii $\le$1.5 R$_\oplus$. Depending on the observing scenario, ~25-45 rocky exoplanets (objects with radii between 0.5 and 1.5 $_{\oplus}$) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With an aperture size of 3.5 m, the total number of detections can increase to up to ~770, including ~60-80 rocky, eHZ planets. With 1 m aperture size, the maximum detection yield is ~315 exoplanets, including $\le$20 rocky, eHZ planets. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. (abridged)

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 polar orbit of the warm Neptune GJ436b seen with VLT/ESPRESSO

GJ436b might be the prototype of warm Neptunes that have undergone late migration induced by an outer companion. Precise determination of the orbital architecture of such systems is critical to constraining their dynamical history and evaluating the role of delayed migration in the exoplanet population. To this purpose we analyzed the Rossiter-McLaughlin (RM) signal of GJ436 b in two transits - recently observed with ESPRESSO - using three different techniques. The high level of precision achieved in radial velocity (RV) measurements allows us to detect the deviation from the Keplerian orbit, despite the slow rotation of the M dwarf host (vsini = 272.0+40.0-34.0 m/s), and to measure the sky-projected obliquity ($λ$ = 102.5+17.2-18.5$^{\circ}$). The Reloaded RM technique, which allows the stellar RV field along the transit chord to be analyzed, yields $λ$ = 107.5+23.6-19.3$^{\circ}$ and vsini = 292.9+41.9-49.9 m/s. The RM Revolutions technique, which allows us to fit the spectral profiles from all planet-occulted regions together, yields $λ$ = 114.1+22.8-17.8$^{\circ}$ and vsini = 300.5+45.9-57.0 m/s. The consistent results between these three techniques, and with published results from HARPS/HARPS-N data, confirm the polar orbit of GJ436b and support the hypothesis that its origin lies in Kozai migration. Results from a joint RM Revolutions analysis of the ESPRESSO, HARPS, and HARPS-N datasets ($λ$ = 113.5+23.3-17.3$^{\circ}$; vsini = 293.5+43.7-52.2 m/s) combined with a revised stellar inclination ($i_\star$ = 35.7+5.9-7.6$^{\circ}$ or 144.2+7.6-5.9$^{\circ}$) lead us to constrain the 3D obliquity $Ψ$ to 103.2+12.8-11.5$^{\circ}$.

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

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

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

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.

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

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

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 loss rate and local thermodynamic state of KELT-9 b thermosphere from the hydrogen Balmer series

KELT-9 b, the hottest known exoplanet with $T\sim4400$ K, is the archetype of a new planet class known as ultra-hot Jupiters. These exoplanets are presumed to have an atmosphere dominated by neutral and ionized atomic species. In particular, H$α$ and H$β$ Balmer lines have been detected in the KELT-9 b upper atmosphere, suggesting that hydrogen is filling the planetary Roche lobe and escaping from the planet. In this work, we detected $δ$ Scuti-type stellar pulsation (with a period $P=7.54\pm0.12$ h) and studied the Rossiter-McLaughlin effect (finding a spin-orbit angle $λ=-85.01^°\pm0.23^°$) prior to focussing on the Balmer lines (H$α$ to H$ζ$) in the optical transmission spectrum of KELT-9 b. Our HARPS-N data show significant absorption for H$α$ to H$δ$. The precise line shapes of the H$α$, H$β$, and H$γ$ absorptions allow us to put constraints on the thermospheric temperature. Moreover, the mass loss rate, and the excited hydrogen population of KELT-9 b are also constrained, thanks to a retrieval analysis performed with a new atmospheric model. We retrieved a thermospheric temperature of $T=13200^{+800}_{-720}$ K and a mass loss rate of $\dot{M}=10^{12.8\pm0.3}$ g s$^{-1}$ when the atmosphere was assumed to be in hydrodynamical expansion and in local thermodynamic equilibrium (LTE). Since the thermospheres of hot Jupiters are not expected to be in LTE, we explored atmospheric structures with non-Boltzmann equilibrium for the population of the excited hydrogen. We do not find strong statistical evidence in favor of a departure from LTE. However, our non-LTE scenario suggests that a departure from the Boltzmann equilibrium may not be sufficient to explain the retrieved low number densities of the excited hydrogen. In non-LTE, Saha equilibrium departure via photo-ionization, is also likely to be necessary to explain the data.

preprint2020arXiv

MOVES III. Simultaneous X-ray and ultraviolet observations unveiling the variable environment of the hot Jupiter HD 189733b

In this third paper of the MOVES (Multiwavelength Observations of an eVaporating Exoplanet and its Star) programme, we combine Hubble Space Telescope far-ultraviolet observations with XMM-Newton/Swift X-ray observations to measure the emission of HD 189733 in various FUV lines, and its soft X-ray spectrum. Based on these measurements we characterise the interstellar medium toward HD 189733 and derive semi-synthetic XUV spectra of the star, which are used to study the evolution of its high-energy emission at five different epochs. Two flares from HD 189733 are observed, but we propose that the long-term variations in its spectral energy distribution have the most important consequences for the environment of HD 189733b. Reduced coronal and wind activity could favour the formation of a dense population of Si$^{2+}$ atoms in a bow-shock ahead of the planet, responsible for pre- and in-transit absorption measured in the first two epochs. In-transit absorption signatures are detected in the Lyman-$α$ line in the second, third and fifth epochs, which could arise from the extended planetary thermosphere and a tail of stellar wind protons neutralised via charge-exchange with the planetary exosphere. We propose that increases in the X-ray irradiation of the planet, and decreases in its EUV irradiation causing lower photoionisation rates of neutral hydrogen, favour the detection of these signatures by sustaining larger densities of H$^{0}$ atoms in the upper atmosphere and boosting charge-exchanges with the stellar wind. Deeper and broader absorption signatures in the last epoch suggest that the planet entered a different evaporation regime, providing clues as to the link between stellar activity and the structure of the planetary environment.

preprint2020arXiv

Neutral Iron Emission Lines From The Day-side Of KELT-9b -- The GAPS Programme With HARPS-N At TNG XX

We present the first detection of atomic emission lines from the atmosphere of an exoplanet. We detect neutral iron lines from the day-side of KELT-9b (Teq $\sim$ 4, 000 K). We combined thousands of spectrally resolved lines observed during one night with the HARPS-N spectrograph (R $\sim$ 115, 000), mounted at the Telescopio Nazionale Galileo. We introduce a novel statistical approach to extract the planetary parameters from the binary mask cross-correlation analysis. We also adapt the concept of contribution function to the context of high spectral resolution observations, to identify the location in the planetary atmosphere where the detected emission originates. The average planetary line profile intersected by a stellar G2 binary mask was found in emission with a contrast of 84 $\pm$ 14 ppm relative to the planetary plus stellar continuum (40 $\pm$ 5$\%$ relative to the planetary continuum only). This result unambiguously indicates the presence of an atmospheric thermal inversion. Finally, assuming a modelled temperature profile previously published (Lothringer et al. 2018), we show that an iron abundance consistent with a few times the stellar value explains the data well. In this scenario, the iron emission originates at the $10^{-3}$-$10^{-5}$ bar level.

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

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

Wind of Change: retrieving exoplanet atmospheric winds from high-resolution spectroscopy

Context. The atmosphere of exoplanets has been studied extensively in recent years, using numerical models to retrieve chemical composition, dynamical circulation or temperature from data. One of the best observational probes in transmission is the sodium doublet, due to its large cross section. However, modelling the shape of the planetary sodium lines has proven to be challenging. Models with different assumptions regarding the atmosphere have been employed to fit the lines in the literature, yet statistically sound direct comparisons of different models are needed to paint a clear picture. Aims. We will compare different wind and temperature patterns and provide a tool to distinguish them driven by their best fit for the sodium transmission spectrum of the hot Jupiter HD 189733b. We parametrise different possible wind patterns already tested in literature and introduce the new option of an upwards driven vertical wind. Methods. We construct a forward model where the wind speed, wind geometry and temperature are injected into the calculation of the transmission spectrum. We embed this forward model in a nested sampling retrieval code to rank the models via their Bayesian evidence. Results. We retrieve a best-fit to the HD 189733b data for vertical upward winds $|\vec{v}_{\mathrm{ver}}(\mathrm{mean})|=40\pm4$ km/s at altitudes above $10^{-6}$ bar. With the current data from HARPS, we cannot distinguish wind patterns for higher pressure atmospheric layers. Conclusions. We show that vertical upwards winds in the upper atmosphere are a possible explanation for the broad sodium signature in hot Jupiters. We highlight other influences on the width of the doublet and explore strong magnetic fields acting on the lower atmosphere as one possible origin of the retrieved wind speed.