Source author record

C. Lovis

C. Lovis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

142works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

142 published item(s)

preprint2025arXiv

The Phase Induced Amplitude Apodizer and Nuller -- High transmission, high dispersion coronagraphy at 2λ/D

Context. Proxima Cen b is the prime target for the search of life around a nearby exoplanet by characterizing its atmosphere in reflected light. Due to the very high star/companion contrast (<1E-6), High Dispersion Coronagraphy is the most promising technique to perform such a characterization. Aims. With a maximum separation of 37 mas, Proxima b can be observed with a VLT in the visible. It requires a coronagraph providing high contrast (< 1E-4 ) very close from the star (< 2 λ/D ) over a broad spectral range (~30%), with a high transmission of the companion (> 50%). We look for an optimal solution that takes benefit of the properties of single-mode fibers. Methods. We introduce the Phase Induced Amplitude Apodizer and Nuller (PIAAN), a coronagraphic integral field unit, designed to feed a diffraction limited spectrograph. It uses a pupil remapping optics with moderate apodization, combined to a single mode fiber integral field unit. It exploits the properties of single mode fibers to null the star light without reducing the companion coupling. The study focuses on a proper tolerance analysis and proposes a wavefront optimization strategy. A prototype is built to demonstrate its performance. Results. We show that the PIAAN can theoretically provide contrasts of 7E-7 and a transmission of 72% at 2 λ/D over a bandwidth of 30%. A prototype is built and characterized and the proposed wavefront control strategy is also demonstrated in the lab. We reach contrast levels of 3E-5 over the full bandwidth, as expected from the tolerance analysis. Conclusions. We demonstrated the potential of the new PIAAN coronagraph, from simulations to prototype. Its performance will eventually be limited by the XAO capabilities. It is the main coronagraph candidate for the RISTRETTO instrument to observe Proxima Cen b on the VLT and its first technology milestone.

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.

preprint2023arXiv

The pristine nature of SMSS 1605$-$1443 revealed by ESPRESSO

SMSS J160540.18$-$144323.1 is the carbon-enhanced metal-poor (CEMP) star with the lowest iron abundance ever measured, [Fe/H]=-6.2, which was first reported with the SkyMapper telescope. The carbon abundance is A(C)~6.1 in the low-C band, as the majority of the stars in this metallicity range. Yet, constraining the isotopic ratio of key species, such as carbon, sheds light on the properties and origin of these elusive stars. We performed high-resolution observations of SMSS1605$-$1443 with the ESPRESSO spectrograph to look for variations in the radial velocity ($v_{rad}$) with time. These data have been combined with older MIKE and UVES archival observations to enlarge the temporal baseline. The $^{12}$C/$^{13}$C isotopic ratio is also studied to explore the possibility of mass transfer from a binary companion. A cross-correlation function against a natural template was applied to detect $v_{rad}$ variability and a spectral synthesis technique was used to derive $^{12}$C/$^{13}$C in the stellar atmosphere. We confirm previous indications of binarity in SMSS1605$-$1443 and measured a lower limit $^{12}$C/$^{13}$C$>60$ at more than a 3$σ$ confidence level, proving that this system is chemically unmixed and that no mass transfer from the unseen companion has happened so far. Thus, we confirm the CEMP-no nature of SMSS1605$-$1443 and show that the pristine chemical composition of the cloud from which it formed is currently imprinted in its stellar atmosphere free of contamination.

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'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'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's activity using a Gaussian process (GP). We use the CCF'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'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'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'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'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

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

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'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'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 "end-to-end" 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

RISTRETTO: coronagraph and AO designs enabling High Dispersion Coronagraphy at 2 lambda/D

RISTRETTO is the evolution of the original idea of coupling the VLT instruments SPHERE and ESPRESSO, aiming at High Dispersion Coronagraphy. RISTRETTO is a visitor instrument that should enable the characterization of the atmospheres of nearby exoplanets in reflected light, by using the technique of high-contrast, high-resolution spectroscopy. Its goal is to observe Prox Cen b and other planets placed at about 35mas from their star, i.e. 2lambda/D at lambda=750nm. The instrument is composed of an extreme adaptive optics, a coronagraphic Integral Field Unit, and a diffraction-limited spectrograph (R=140.000, lambda=620-840 nm). We present the status of our studies regarding the coronagraphic IFU and the XAO system. The first in particular is based on a modified version of the PIAA apodizer, allowing nulling on the first diffraction ring. Our proposed design has the potential to reach > 50% coupling and <1E-4 contrast at 2lambda/D in median seeing conditions.

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 "ultra-hot Jupiters" 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 Hot Neptune WASP-166 b with ESPRESSO II: Confirmation of atmospheric sodium

The hot Neptune desert, a distinct lack of highly irradiated planets in the size range of Neptune, remains one of the most intriguing results of exoplanet population studies. A deeper understanding of the atmosphere of exoplanets sitting at the edge or even within the Neptune desert will allow us to better understand if planetary formation or evolution processes are at the origin of the desert. A detection of sodium in WASP-166b was presented previously with tentative line broadening at the 3.4 sigma with the HARPS spectrograph. We update this result with two transits observed with the ESPRESSO spectrograph, confirming the detection in each night and the broadened character of the line. This result marks the first confirmed resolved sodium detection within the Neptune desert. In this work, we additionally highlight the importance of treating low-SNR spectral regions, particularly where absorption lines of stellar sodium and planetary sodium overlap at mid-transit - an important caveat for future observations of the system.

preprint2022arXiv

The Hot Neptune WASP-166~b with ESPRESSO I: Refining the Planetary Architecture and Stellar Variability

In this paper, we present high-resolution spectroscopic transit observations from ESPRESSO of the super-Neptune WASP-166~b. In addition to spectroscopic ESPRESSO data, we analyse photometric data from {\sl TESS} of six WASP-166~b transits along with simultaneous NGTS observations of the ESPRESSO runs. These observations were used to fit for the planetary parameters as well as assessing the level of stellar activity (e.g. spot crossings, flares) present during the ESPRESSO observations. We utilise the Reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterising the centre-to-limb convection-induced variations, and to refine the star-planet obliquity. We find WASP-166~b has a projected obliquity of $λ= -15.52^{+2.85}_{-2.76}$$^{\circ}$ and $v\sin(i) = 4.97 \pm 0.09$~kms$^{-1}$ which is consistent with the literature. We were able to characterise centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few kms$^{-1}$ for the first time. We modelled the centre-to-limb convective variations using a linear, quadratic and cubic model with the cubic being preferred. In addition, by modelling the differential rotation and centre-to-limb convective variations simultaneously we were able to retrieve a potential anti-solar differential rotational shear ($α\sim$ -0.5) and stellar inclination ($i_*$ either 42.03$^{+9.13}_{-9.60}$$^{\circ}$ or 133.64$^{+8.42}_{-7.98}$$^{\circ}$ if the star is pointing towards or away from us). Finally, we investigate how the shape of the cross-correlation functions change as a function of limb angle and compare our results to magnetohydrodynamic simulations.

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.

preprint2022arXiv

Transmission spectroscopy of MASCARA-1b with ESPRESSO: Challenges of overlapping orbital and Doppler tracks

Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$α$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b.

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

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

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

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

Detection of Na, K, and H$α$ absorption in the atmosphere of WASP-52b using ESPRESSO

WASP-52b is a low density hot Jupiter orbiting a moderately active K2V star. Previous low-resolution studies have revealed a cloudy atmosphere and found atomic Na above the cloud deck. Here we report the detection of excess absorption at the Na doublet, the H$α$ line, and the K D$_1$ line. We derived a high-resolution transmission spectrum based on three transits of WASP-52b, observed with the ultra stable, high-resolution spectrograph ESPRESSO at the VLT. We measure a line contrast of $1.09\pm 0.16$% for Na D$_1$, $1.31\pm 0.13$% for Na D$_2$, $0.86\pm 0.13$% for H$α$, and $0.46\pm 0.13$% for K D$_1$, with a line FWHM range of 11-22 km s$^{-1}$. We also find that the velocity shift of these detected lines during the transit is consistent with the planet orbital motion, thus confirming their planetary origin. We do not observe any significant net blueshift or redshift that can be attributed to planetary winds. We use activity indicator lines as control but find no excess absorption. However, we do notice signatures arising from the Center-to-Limb variation (CLV) and the Rossiter-McLaughlin (RM) effect at these control lines. This highlights the importance of the CLV+RM correction in correctly deriving the transmission spectrum, which, if not corrected, could resemble or cancel out planetary absorption in certain cases. WASP-52b is the second non-ultra-hot Jupiter to show excess H$α$ absorption, after HD 189733b. Future observations targeting non-ultra-hot Jupiters that show H$α$ could help reveal the relation between stellar activity and the heating processes in the planetary upper atmosphere.

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

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

Stellar chromospheric activity of 1,674 FGK stars from the AMBRE-HARPS sample I. A catalogue of homogeneous chromospheric activity

We present a catalogue of homogeneous determined chromospheric emission (CE), stellar atmospheric parameters and ages for 1,674 FGK main sequence (MS), subgiant, and giant stars. The analysis of CE level and variability is also performed. We measured CE in the CaII lines using more than 180,000 high-resolution spectra from the HARPS spectrograph, as compiled in the AMBRE project, obtained between 2003 and 2019. We converted the fluxes to bolometric and photospheric corrected chromospheric emission ratio, $R'_\text{HK}$. Stellar atmospheric parameters $T_\text{eff}$, $\log g$, and [Fe/H] were retrieved from the literature or determined using an homogeneous method. $M_\star$, $R_\star$, and ages were determined from isochrone fitting. We analysed the CE distribution for the different luminosity classes and spectral types and confirmed the existence of the very inactive stars (VIS) and very active stars (VAS) populations at $\log R'_\text{HK}< -5.1$ and $> -4.2$ dex, respectively. We found indications that the VIS population is composed mainly of subgiant and giant stars and that $\log R'_\text{HK}= -5.1$ dex marks a transition in stellar evolution. There appears to be at least three regimes of variability, for inactive, active and very active stars, with the inactive and active regimes separated by a diagonal Vaughan-Preston gap. We show that stars with low activity levels do not necessarily have low variability. In the case of K dwarfs which show high CE variability, inactive and active stars have similar levels of activity variability. This means that activity levels alone are not enough to infer about the activity variability of a star. We also explained the shape of the VP gap observed in the distribution of CE by using the CE variability-level diagram. In the CE variability-level diagram, the Sun is located in the high variability region of the inactive MS stars zone. (Abridged)

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Radial velocities and activity indicators from cross-correlation functions with weighted binary masks

For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology to first build weighted binary masks and, second, to compute the CCF of spectral observations with these masks from which we derive radial velocities and activity indicators. These methods are implemented in a python code that is publicly available. To build the masks, we selected a large number of sharp absorption lines based on the profile of the minima present in high signal-to-noise ratio (S/N) spectrum templates built from observations of reference stars. We computed the CCFs of observed spectra and derived RVs and the following three standard activity indicators: full-width-at-half-maximum as well as contrast and bisector inverse slope.We applied our methodology to CARMENES high-resolution spectra and obtain RV and activity indicator time series of more than 300 M dwarf stars observed for the main CARMENES survey. Compared with the standard CARMENES template matching pipeline, in general we obtain more precise RVs in the cases where the template used in the standard pipeline did not have enough S/N. We also show the behaviour of the three activity indicators for the active star YZ CMi and estimate the absolute RV of the M dwarfs analysed using the CCF RVs.

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

Measuring precise radial velocities on individual spectral lines : II. Dependance of stellar activity signal on line depth

Context. Although the new generation of radial-velocity (RV) instruments such as ESPRESSO are expected to reach the long-term precision required to find other earths, the RV measurements are contaminated by some signal from stellar activity. This makes these detections hard. Aims. On real observations, we here demonstrate for the first time the effect of stellar activity on the RV of individual spectral lines. Recent studies have shown that this is probably the key for mitigating this perturbing signal. Results. We estimate that at least 89% of the lines that appear in the spectrum of alpha Cen B for which we measure a reliable RV are correlated with the stellar activity signal. This can be interpreted as those lines being sensitive to the inhibition of the convective blueshift observed in active regions. Because the velocity of the convective blueshift increases with physical depth inside the stellar atmosphere, we find that the effect induced by stellar activity on the RV of individual spectral lines is inversely proportional to the line depth. The stellar activity signal can be mitigated down to 0.8-0.9 m/s, either by selecting lines that are less sensitive to activity or by using the difference between the RV of the spectral lines that are formed at different depths in the stellar atmosphere as an activity proxy. Conclusions. This paper shows for the first time that based on real observations of solar-type stars, it is possible to measure the RV effect of stellar activity on the RV of individual spectral lines. Our results are very promising and demonstrate that analysing the RV of individual spectral lines is probably one of the solutions to mitigate stellar activity signal in RV measurements down to a level enabling the detection of other earths.

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.

preprint2016arXiv

An extreme planetary system around HD219828. One long-period super Jupiter to a hot-neptune host star

With about 2000 extrasolar planets confirmed, the results show that planetary systems have a whole range of unexpected properties. We present a full investigation of the HD219828 system, a bright metal-rich star for which a hot neptune has previously been detected. We used a set of HARPS, SOPHIE, and ELODIE radial velocities to search for the existence of orbiting companions to HD219828. A dynamical analysis is also performed to study the stability of the system and to constrain the orbital parameters and planet masses. We announce the discovery of a long period (P=13.1years) massive (msini=15.1MJup) companion (HD219828c) in a very eccentric orbit (e=0.81). The same data confirms the existence of a hot-neptune, HD219828b, with a minimum mass of 21 MEarth and a period of 3.83days. The dynamical analysis shows that the system is stable. The HD219828 system is extreme and unique in several aspects. First, among all known exoplanet systems it presents an unusually high mass ratio. We also show that systems like HD219828, with a hot neptune and a long-period massive companion are more frequent than similar systems with a hot jupiter instead. This suggests that the formation of hot neptunes follows a different path than the formation of their hot jovian counterparts. The high mass, long period, and eccentricity of HD219828c also make it a good target for Gaia astrometry as well as a potential target for atmospheric characterisation, using direct imaging or high-resolution spectroscopy. Astrometric observations will allow us to derive its real mass and orbital configuration. If a transit of HD219828b is detected, we will be able to fully characterise the system, including the relative orbital inclinations. With a clearly known mass, HD219828c may become a benchmark object for the range in between giant planets and brown dwarfs.

preprint2016arXiv

EELT-HIRES the high-resolution spectrograph for the E-ELT

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

preprint2016arXiv

Is the activity level of HD 80606 influenced by its eccentric planet?

Aims: Several studies suggest that the activity level of a planet-host star can be influenced by the presence of a close-by orbiting planet. Moreover, the interaction mechanisms that have been proposed, magnetic interaction and tidal interaction, exhibit a very different dependence on orbital separation between the star and the planet. A detection of activity enhancement and characterization of its dependence on planetary orbital distance can, in principle, allow us to characterize the physical mechanism behind the activity enhancement. Methods: We used the HARPS-N spectrograph to measure the stellar activity level of HD 80606 during the planetary periastron passage and compared the activity measured to that close to apastron. Being characterized by an eccentricity of 0.93 and an orbital period of 111 days, the system's extreme variation in orbital separation makes it a perfect target to test our hypothesis. Results: We find no evidence for a variation in the activity level of the star as a function of planetary orbital distance, as measured by all activity indicators employed log($R'_{HK}$), H$_α$, NaI, and HeI. None of the models employed, whether magnetic interaction or tidal interaction, provides a good description of the data. Conclusions: We find no evidence for star-planet interaction in HD\,80606 at the moment of the periastron passage of its very eccentric planet. The straightforward explanation for the non-detection is the absence of interaction as a result of a low magnetic field strength on either the planet or the star and of the low level of tidal interaction between the two. However, we cannot exclude two scenarios: i) the interaction can be instantaneous and of magnetic origin, being concentrated on the substellar point and its surrounding area, and ii) the interaction can lead to a delayed activity enhancement. (abridged)

preprint2016arXiv

Rossiter--McLaughlin models and their effect on estimates of stellar rotation, illustrated using six WASP systems

We present new measurements of the projected spin--orbit angle $λ$ for six WASP hot Jupiters, four of which are new to the literature (WASP-61, -62, -76, and -78), and two of which are new analyses of previously measured systems using new data (WASP-71, and -79). We use three different models based on two different techniques: radial velocity measurements of the Rossiter--McLaughlin effect, and Doppler tomography. Our comparison of the different models reveals that they produce projected stellar rotation velocities ($v \sin I_{\rm s}$) measurements often in disagreement with each other and with estimates obtained from spectral line broadening. The Boué model for the Rossiter--McLaughlin effect consistently underestimates the value of $v\sin I_{\rm s}$ compared to the Hirano model. Although $v \sin I_s$ differed, the effect on $λ$ was small for our sample, with all three methods producing values in agreement with each other. Using Doppler tomography, we find that WASP-61\,b ($λ=4^\circ.0^{+17.1}_{-18.4}$), WASP-71\,b ($λ=-1^\circ.9^{+7.1}_{-7.5}$), and WASP-78\,b ($λ=-6^\circ.4\pm5.9$) are aligned. WASP-62\,b ($λ=19^\circ.4^{+5.1}_{-4.9}$) is found to be slightly misaligned, while WASP-79\,b ($λ=-95^\circ.2^{+0.9}_{-1.0}$) is confirmed to be strongly misaligned and has a retrograde orbit. We explore a range of possibilities for the orbit of WASP-76\,b, finding that the orbit is likely to be strongly misaligned in the positive $λ$ direction.

preprint2016arXiv

Search for giant planets in M67 III: excess of hot Jupiters in dense open clusters

Since 2008 we used high-precision radial velocity (RV) measurements obtained with different telescopes to detect signatures of massive planets around main-sequence and evolved stars of the open cluster (OC) M67. We aimed to perform a long-term study on giant planet formation in open clusters and determine how this formation depends on stellar mass and chemical composition. A new hot Jupiter (HJ) around the main-sequence star YBP401 is reported in this work. An update of the RV measurements for the two HJ host-stars YBP1194 and YBP1514 is also discussed. Our sample of 66 main-sequence and turnoff stars includes 3 HJs, which indicates a high rate of HJs in this cluster (~5.6% for single stars and ~4.5% for the full sample ). This rate is much higher than what has been discovered in the field, either with RV surveys or by transits. High metallicity is not a cause for the excess of HJs in M67, nor can the excess be attributed to high stellar masses. When combining this rate with the non-zero eccentricity of the orbits, our results are qualitatively consistent with a HJ formation scenario dominated by strong encounters with other stars or binary companions and subsequent planet-planet scattering, as predicted by N-body simulations.

preprint2016arXiv

Searching for Li-rich giants in a sample of 12 open clusters: Li enhancement in two stars with substellar companions

The aim of this work is to search for Li-rich giants in a sample of clusters where planets have been searched, thus we can study the planet engulfment scenario to explain Li replenishment using a proper comparison sample of stars without detected giant planets. We derived Li abundances for a sample of 67 red giant stars in 12 different open clusters using standard spectral synthesis techniques and high resolution spectra (from HARPS and UVES). We also determined masses, ages and radius from PARSEC stellar isochrones to constrain the evolutionary stage of these stars. We find three stars in different clusters with clearly enhanced Li abundances compared to other stars within the cluster. Interestingly, the only two stars with a detected substellar companion in our sample belong to that group. One of the planet hosts, NGC2423No3, might lie close to the luminosity bump on the HR diagram, a phase where Li production by the Cameron-Fowler process is supported by extra-mixing to bring fresh Li up to the surface. On the other hand, NGC4349No127 is a more massive and more evolved giant that does not seem to be in the evolutionary phase where other Li-rich stars are found. We discuss the possibility that the Li enhancement of this star is triggered by the engulfment of a planet, considering that close-in planets hardly survive the RGB-tip and the early-AGB phases.

preprint2016arXiv

The HARPS search for southern extra-solar planets XL. Searching for Neptunes around metal-poor stars

Stellar metallicity -- as a probe of the metallicity of proto-planetary disks -- is an important ingredient for giant planet formation, likely through its effect on the timescales in which rocky/icy planet cores can form. Giant planets have been found to be more frequent around metal-rich stars, in agreement with predictions based on the core-accretion theory. In the metal-poor regime, however, the frequency of planets, especially low-mass planets, and how it depends on metallicity are still largely unknown. As part of a planet search programme focused on metal-poor stars, we study the targets from this survey that were observed with HARPS on more than 75 nights. The main goals are to assess the presence of low-mass planets and provide a first estimate of the frequency of Neptunes and super-Earths around metal-poor stars. We perform a systematic search for planetary companions, both by analysing the periodograms of the radial-velocities and by comparing, in a statistically-meaningful way, models with an increasing number of Keplerians. A first constraint on the frequency of planets in our metal-poor sample is calculated considering the previous detection (in our sample) of a Neptune-sized planet around HD175607 and one candidate planet (with an orbital period of 68.42d and minimum mass $M_p \sin i = 11.14 \pm 2.47 M_{\oplus}$) for HD87838, announced in the present study. This frequency is determined to be close to 13% and is compared with results for solar-metallicity stars.

preprint2016arXiv

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system's dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the in- from the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (> 99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (lambda ~ -0.4 +/- 0.2 degrees) and true 3D obliquity (psi ~ 7^+12_-4 degrees). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.

preprint2016arXiv

The SOPHIE search for northern extrasolar planets. X. Detection and characterization of giant planets by the dozen

We present new radial velocity measurements of eight stars secured with the spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory allowing the detection and characterization of new giant extrasolar planets. The host stars are dwarfs of spectral types between F5 and K0 and magnitudes between 6.7 and 9.6; the planets have minimum masses M_p sin i between 0.4 to 3.8 M_Jup and orbital periods of several days to several months. The data allow only single planets to be discovered around the first six stars (HD143105, HIP109600, HD35759, HIP109384, HD220842, and HD12484), but one of them shows the signature of an additional substellar companion in the system. The seventh star, HIP65407, allows the discovery of two giant planets, just outside the 12:5 resonance in weak mutual interaction. The last star, HD141399, was already known to host a four-planetary system; our additional data and analyses allow new constraints to be put on it. We present Keplerian orbits of all systems, together with dynamical analyses of the two multi-planetary systems. HD143105 is one of the brightest stars known to host a hot Jupiter, which could allow numerous follow-up studies to be conducted despite this is not a transiting system. The giant planets HIP109600b, HIP109384b, and HD141399c are located in the habitable zone of their host star.

preprint2016arXiv

The Sun as a planet-host star: Proxies from SDO images for HARPS radial-velocity variations

The Sun is the only star whose surface can be directly resolved at high resolution, and therefore constitutes an excellent test case to explore the physical origin of stellar radial-velocity (RV) variability. We present HARPS observations of sunlight scattered off the bright asteroid 4/Vesta, from which we deduced the Sun's activity-driven RV variations. In parallel, the HMI instrument onboard the Solar Dynamics Observatory provided us with simultaneous high spatial resolution magnetograms, Dopplergrams, and continuum images of the Sun in the Fe I 6173A line. We determine the RV modulation arising from the suppression of granular blueshift in magnetised regions and the flux imbalance induced by dark spots and bright faculae. The rms velocity amplitudes of these contributions are 2.40 m/s and 0.41 m/s, respectively, which confirms that the inhibition of convection is the dominant source of activity-induced RV variations at play, in accordance with previous studies. We find the Doppler imbalances of spot and plage regions to be only weakly anticorrelated. Lightcurves can thus only give incomplete predictions of convective blueshift suppression. We must instead seek proxies that track the plage coverage on the visible stellar hemisphere directly. The chromospheric flux index R'_HK derived from the HARPS spectra performs poorly in this respect, possibly because of the differences in limb brightening/darkening in the chromosphere and photosphere. We also find that the activity-driven RV variations of the Sun are strongly correlated with its full-disc magnetic flux density, which may become a useful proxy for activity-related RV noise.

preprint2015arXiv

Chemical abundances and kinematics of 257 G-, K-type field giants. Setting a base for further analysis of giant-planet properties orbiting evolved stars

We performed a uniform and detailed abundance analysis of 12 refractory elements (Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) for a sample of 257 G- and K-type evolved stars from the CORALIE planet search program. To date, only one of these stars is known to harbor a planetary companion. We aimed to characterize this large sample of evolved stars in terms of chemical abundances and kinematics, thus setting a solid base for further analysis of planetary properties around giant stars. This sample, being homogeneously analyzed, can be used as a comparison sample for other planet-related studies, as well as for different type of studies related to stellar and Galaxy astrophysics. The abundances of the chemical elements were determined using an LTE abundance analysis relative to the Sun, with the spectral synthesis code MOOG and a grid of Kurucz ATLAS9 atmospheres. To separate the Galactic stellar populations both a purely kinematical approach and a chemical method were applied. We confirm the overabundance of Na in giant stars compared to the field FGK dwarfs. This enhancement might have a stellar evolutionary character, but departures from LTE may also produce a similar enhancement. Our chemical separation of stellar populations also suggests a "gap" in metallicity between the thick-disk and high-alpha metal-rich stars, as previously observed in dwarfs sample from HARPS. The present sample, as most of the giant star samples, also suffers from the B - V colour cut-off, which excludes low-log g stars with high metallicities, and high-logg star with low-[Fe/H]. For future studies of planet occurrence dependence on stellar metallicity around these evolved stars we suggest to use a sub-sample of stars in a "cut-rectangle" in the logg - [Fe/H] diagram to overcome the aforementioned issue.

preprint2015arXiv

Detecting ring systems around exoplanets using high resolution spectroscopy: the case of 51Pegb

In this paper we explore the possibility that the recently detected reflected light signal of 51\,Peg\,b could be caused by a ring system around the planet. We use a simple model to compare the observed signal with the expected signal from a short-period giant planet with rings. We also use simple dynamical arguments to understand the possible geometry of such a system. We provide evidence that, to a good approximation, the observations are compatible with the signal expected from a ringed planet, assuming that the rings are non-coplanar with the orbital plane. However, based on dynamical arguments, we also show that this configuration is unlikely. In the case of coplanar rings we then demonstrate that the incident flux on the ring surface is about 2\% the value received by the planet, a value that renders the ring explanation unlikely. The results suggest that the signal observed cannot in principle be explained by a planet+ring system. We discuss, however, the possibility of using reflected light spectra to detect and characterize the presence of rings around short-period planets. Finally, we show that ring systems could have already been detected by photometric transit campaigns, but their signal could have been easily misinterpreted by the expected light curve of an eclipsing binary.

preprint2015arXiv

Determination of the spectroscopic stellar parameters for 257 field giant stars

The study of stellar parameters of planet-hosting stars, such as metallicity and chemical abundances, help us to understand the theory of planet formation and stellar evolution. Here, we present a catalogue of accurate stellar atmospheric parameters and iron abundances for a sample of 257 K and G field evolved stars that are being surveyed for planets using precise radial--velocity measurements as part of the CORALIE programme to search for planets around giants. The analysis was done using a set of high--resolution and high--signal-to-noise Ultraviolet and Visible Echelle Spectrograph spectra. The stellar parameters were derived using Fe I and II ionization and excitation equilibrium methods. To take into account possible effects related to the choice of the lines on the derived parameters, we used three different iron line-list sets in our analysis, and the results differ among themselves by a small factor for most of stars. {For those stars with previous literature parameter estimates, we found very good agreement with our own values.} In the present catalogue we are providing new precise spectroscopic measurements of effective temperature, surface gravity, microturbulence, and metallicity for 190 stars for which it has not been found or published in previous articles.

preprint2015arXiv

Evidence for a spectroscopic direct detection of reflected light from 51 Peg b

The detection of reflected light from an exoplanet is a difficult technical challenge at optical wavelengths. Even though this signal is expected to replicate the stellar signal, not only is it several orders of magnitude fainter, but it is also hidden among the stellar noise. We apply a variant of the cross-correlation technique to HARPS observations of 51 Peg to detect the reflected signal from planet 51 Peg b. Our method makes use of the cross-correlation function of a binary mask with high-resolution spectra to amplify the minute planetary signal that is present in the spectra by a factor proportional to the number of spectral lines when performing the cross correlation. The resulting cross-correlation functions are then normalized by a stellar template to remove the stellar signal. Carefully selected sections of the resulting normalized CCFs are stacked to increase the planetary signal further. The recovered signal allows probing several of the planetary properties, including its real mass and albedo. We detect evidence for the reflected signal from planet 51 Peg b at a significance of 3σ_noise. The detection of the signal permits us to infer a real mass of 0.46^+0.06_-0.01 M_Jup (assuming a stellar mass of 1.04\;M_Sun) for the planet and an orbital inclination of 80^+10_-19 degrees. The analysis of the data also allows us to infer a tentative value for the (radius-dependent) geometric albedo of the planet. The results suggest that 51Peg b may be an inflated hot Jupiter with a high albedo (e.g., an albedo of 0.5 yields a radius of 1.9 \pm 0.3 R_Jup for a signal amplitude of 6.0\pm0.4 x 10^-5). We confirm that the method we perfected can be used to retrieve an exoplanet's reflected signal, even with current observing facilities. The advent of next generation of observing facilities will yield new opportunities for this type of technique to probe deeper into exoplanets.

preprint2015arXiv

Radial velocity information content of M dwarf spectra in the near-infrared

Aims: We evaluate the radial velocity (RV) information content and achievable precision on M0-M9 spectra covering the ZYJHK bands. We do so while considering both a perfect atmospheric transmission correction and discarding areas polluted by deep telluric features, as done in previous works. Methods: To simulate the M-dwarf spectra, PHOENIX-ACES model spectra were employed; they were convolved with rotational kernels and instrumental profiles to reproduce stars with a $v.sin{i}$ of 1.0, 5.0, and 10.0 km/s when observed at resolutions of 60 000, 80 000, and 100 000. We considered the RV precision as calculated on the whole spectra, after discarding strongly polluted areas, and after applying a perfect telluric correction. In our simulations we paid particular attention to the details of the convolution and sampling of the spectra, and we discuss their impact on the final spectra. Results: Our simulations show that the most important parameter ruling the difference in attainable precision between the considered bands is the spectral type. For M0-M3 stars, the bands that deliver the most precise RV measurements are the Z, Y, and H band, with relative merits depending on the parameters of the simulation. For M6-M9 stars, the bands show a difference in precision that is within a factor of $\sim$2 and does not clearly depend on the band; this difference is reduced to a factor smaller than $\sim$1.5 if we consider a non-rotating star seen at high resolution. We also show that an M6-M9 spectrum will deliver a precision about two times better as an M0-M3 spectra with the same signal-to-noise ratio. Finally, we note that the details of modelling the Earth atmosphere and interpreting the results have a significant impact on which wavelength regions are discarded when setting a limit threshold at 2-3%. (abridged)

preprint2015arXiv

Spectrally resolved detection of sodium in the atmosphere of HD189733b with the HARPS spectrograph

Atmospheric properties of exoplanets can be constrained with transit spectroscopy. The signature of atomic sodium NaI, known to be present above the clouds, is a powerful probe of the upper atmosphere, where it can be best detected and characterized at high spectral resolution. Our goal is to obtain a high-resolution transit spectrum of HD189733b in the region around the resonance doublet of NaI at 589 nm, to characterize the absorption signature previously detected from space at low resolution. We analyze archival transit data of HD189733b obtained with the HARPS spectrograph. We retrieve the transit spectrum and light curve of the planet, implementing corrections for telluric contamination and planetary orbital motion. We spectrally resolve the NaI D doublet and measure line contrasts of $0.64\pm0.07\%$ (D2) and $0.40\pm0.07\%$ (D1) and FWHMs of $0.52\pm0.08~Å$. This corresponds to a detection at the 10-$σ$ level of excess of absorption of $0.32\pm0.03\%$ in a passband of $2\times0.75\ Å$ centered on each line. We derive temperatures of $2\,600\pm600$ K and $3270\pm330$ K at altitudes of $9\,800\pm2\,800$ km and $12\,700\pm2\,600$ km in the NaI D1 and D2 line cores, respectively. We measure a temperature gradient of $\sim0.2$ K km$^{-1}$ from comparison with theoretical models. We also detect a blueshift of $0.16\pm0.04\ Å$ (4 $σ$) in the line positions. This blueshift may be due to winds blowing at $8\pm2$ km s$^{-1}$ in the upper layers of the atmosphere. We demonstrate the relevance of studying exoplanet atmospheres with high-resolution spectrographs mounted on 4-meter-class telescopes. Our results pave the way towards in-depth characterization of physical conditions in the atmospheres of many exoplanetary systems with future spectrographs such as ESPRESSO on the VLT or HiReS and METIS on the E-ELT.

preprint2015arXiv

The HARPS search for southern extra-solar planets. XXXIX. HD175607 b, the most metal-poor G dwarf with an orbiting sub-Neptune

Context. The presence of a small-mass planet (M$_p<$0.1\,M$_{Jup}$) seems, to date, not to depend on metallicity, however, theoretical simulations have shown that stars with subsolar metallicities may be favoured for harbouring smaller planets. A large, dedicated survey of metal-poor stars with the HARPS spectrograph has thus been carried out to search for Neptunes and super-Earths. Aims. In this paper, we present the analysis of \object{HD175607}, an old G6 star with metallicity [Fe/H] = -0.62. We gathered 119 radial velocity measurements in 110 nights over a time span of more than nine years. Methods. The radial velocities were analysed using Lomb-Scargle periodograms, a genetic algorithm, a Markov chain Monte Carlo analysis, and a Gaussian processes analysis. The spectra were also used to derive stellar properties. Several activity indicators were analysed to study the effect of stellar activity on the radial velocities. Results. We find evidence for the presence of a small Neptune-mass planet (M$_{p}\sin i = 8.98\pm1.10$\,M$_{\oplus}$) orbiting this star with an orbital period $P = 29.01\pm0.02$\, days in a slightly eccentric orbit ($e=0.11\pm0.08$). The period of this Neptune is close to the estimated rotational period of the star. However, from a detailed analysis of the radial velocities together with the stellar activity, we conclude that the best explanation of the signal is indeed the presence of a planetary companion rather than stellar related. An additional longer period signal ($P\sim 1400$\,d) is present in the data, for which more measurements are needed to constrain its nature and its properties. Conclusions. HD\,175607 is the most metal-poor FGK dwarf with a detected low-mass planet amongst the currently known planet hosts. This discovery may thus have important consequences for planet formation and evolution theories.

preprint2015arXiv

The HARPS search for southern extra-solar planets. XXXVIII. Bayesian re-analysis of three systems. New super-Earths, unconfirmed signals, and magnetic cycles

We present the analysis of the entire HARPS observations of three stars that host planetary systems: HD1461, HD40307, and HD204313. The data set spans eight years and contains more than 200 nightly averaged velocity measurements for each star. This means that it is sensitive to both long-period and low-mass planets and also to the effects induced by stellar activity cycles. We modelled the data using Keplerian functions that correspond to planetary candidates and included the short- and long-term effects of magnetic activity. A Bayesian approach was taken both for the data modelling, which allowed us to include information from activity proxies such as $\log{(R'_{\rm HK})}$ in the velocity modelling, and for the model selection, which permitted determining the number of significant signals in the system. The Bayesian model comparison overcomes the limitations inherent to the traditional periodogram analysis. We report an additional super-Earth planet in the HD1461 system. Four out of the six planets previously reported for HD40307 are confirmed and characterised. We discuss the remaining two proposed signals. In particular, we show that when the systematic uncertainty associated with the techniques for estimating model probabilities are taken into account, the current data are not conclusive concerning the existence of the habitable-zone candidate HD40307 g. We also fully characterise the Neptune-mass planet that orbits HD204313 in 34.9 days.

preprint2015arXiv

The HARPS-N Rocky Planet Search I. HD219134b: A transiting rocky planet in a multi-planet system at 6.5 pc from the Sun

We present here the detection of a system of four low-mass planets around the bright (V=5.5) and close-by (6.5 pc) star HD219134. This is the first result of the Rocky Planet Search program with HARPS-N on the TNG in La Palma. The inner planet orbits the star in 3.0937 +/-0.0004 days, on a quasi-circular orbit with a semi-major axis of 0.0382 +/- 0.0003 AU. Spitzer observations allowed us to detect the transit of the planet in front of the star making HD219134b the nearest known transiting planet to date. From the amplitude of the radial-velocity variation (2.33 +/- 0.24 m/s) and observed depth of the transit (359 +/- 38 ppm), the planet mass and radius are estimated to be 4.46 +/- 0.47 M_{\oplus} and 1.606 +/- 0.086 R_{\oplus} leading to a mean density of 5.89 +/- 1.17 g/cc, suggesting a rocky composition. One additional planet with minimum mass of 2.67 +/- 0.59 M_{\oplus} moves on a close-in, quasi-circular orbit with a period of 6.765 +/- 0.005 days. The third planet in the system has a period of 46.78 +/- 0.16 days and a minimum mass of 8.7 +/- 1.1 M{\oplus}, at 0.234 +/- 0.002 AU from the star. Its eccentricity is 0.32 +/- 0.14. The period of this planet is close to the rotational period of the star estimated from variations of activity indicators (42.3 +/- 0.1 days). The planetary origin of the signal is, however, the preferred solution as no indication of variation at the corresponding frequency is observed for activity-sensitive parameters. Finally, a fourth additional longer-period planet of mass of 62 +/- 6 M_{\oplus} orbits the star in 1190 days, on an eccentric orbit (e=0.27 +/- 0.11) at a distance of 2.14 +/- 0.27 AU.

preprint2014arXiv

A global analysis of Spitzer and new HARPS data confirms the loneliness and metal-richness of GJ 436 b

Context. GJ 436b is one of the few transiting warm Neptunes for which a detailed characterisation of the atmosphere is possible, whereas its non-negligible orbital eccentricity calls for further investigation. Independent analyses of several individual datasets obtained with Spitzer have led to contradicting results attributed to the different techniques used to treat the instrumental effects. Aims. We aim at investigating these previous controversial results and developing our knowledge of the system based on the full Spitzer photometry dataset combined with new Doppler measurements obtained with the HARPS spectrograph. We also want to search for additional planets. Methods. We optimise aperture photometry techniques and the photometric deconvolution algorithm DECPHOT to improve the data reduction of the Spitzer photometry spanning wavelengths from 3-24 μm. Adding the high precision HARPS radial velocity data, we undertake a Bayesian global analysis of the system considering both instrumental and stellar effects on the flux variation. Results. We present a refined radius estimate of RP=4.10 +/- 0.16 R_Earth, mass MP=25.4 +/- 2.1 M_Earth and eccentricity e= 0.162 +/- 0.004 for GJ 436b. Our measured transit depths remain constant in time and wavelength, in disagreement with the results of previous studies. In addition, we find that the post-occultation flare-like structure at 3.6 μm that led to divergent results on the occultation depth measurement is spurious. We obtain occultation depths at 3.6, 5.8, and 8.0 μm that are shallower than in previous works, in particular at 3.6 μm. However, these depths still appear consistent with a metal-rich atmosphere depleted in methane and enhanced in CO/CO2, although perhaps less than previously thought. We find no evidence for a potential planetary companion, stellar activity, nor for a stellar spin-orbit misalignment. [ABRIDGED]

preprint2014arXiv

Characterization of the Kepler-101 planetary system with HARPS-N. A hot super-Neptune with an Earth-sized low-mass companion

We report on the characterization of the Kepler-101 planetary system, thanks to a combined DE-MCMC analysis of Kepler data and forty radial velocities obtained with the HARPS-N spectrograph. This system was previously validated by Rowe et al. (2014) and is composed of a hot super-Neptune, Kepler-101b, and an Earth-sized planet, Kepler-101c. These two planets orbit the slightly evolved and metal-rich G-type star in 3.49 and 6.03 days, respectively. With mass $M_{\rm p}=51.1_{-4.7}^{+5.1}~M_{\oplus}$, radius $R_{\rm p}=5.77_{-0.79}^{+0.85}~R_{\oplus}$, and density $ρ_{\rm p}=1.45_{-0.48}^{+0.83} \rm g\;cm^{-3}$, Kepler-101b is the first fully-characterized super-Neptune, and its density suggests that heavy elements make up a significant fraction of its interior; more than $60\%$ of its total mass. Kepler-101c has a radius of $1.25_{-0.17}^{+0.19}~R_{\oplus}$, which implies the absence of any H/He envelope, but its mass could not be determined due to the relative faintness of the parent star for highly precise radial-velocity measurements ($K_{\rm p}=13.8$) and the limited number of radial velocities. The $1~σ$ upper limit, $M_{\rm p} < 3.8~M_{\oplus}$, excludes a pure iron composition with a $68.3\%$ probability. The architecture of the Kepler-101 planetary system - containing a close-in giant planet and an outer Earth-sized planet with a period ratio slightly larger than the 3:2 resonance - is certainly of interest for planet formation and evolution scenarios. This system does not follow the trend, seen by Ciardi et al. (2013), that in the majority of Kepler systems of planet pairs with at least one Neptune-size or larger planet, the larger planet has the longer period.

preprint2014arXiv

CoRoT-22 b: a validated 4.9 RE exoplanet in 10-day orbit

The CoRoT satellite has provided high-precision photometric light curves for more than 163,000 stars and found several hundreds of transiting systems compatible with a planetary scenario. If ground-based velocimetric observations are the best way to identify the actual planets among many possible configurations of eclipsing binary systems, recent transit surveys have shown that it is not always within reach of the radial-velocity detection limits. In this paper, we present a transiting exoplanet candidate discovered by CoRoT whose nature cannot be established from ground-based observations, and where extensive analyses are used to validate the planet scenario. They are based on observing constraints from radial-velocity spectroscopy, adaptive optics imaging and the CoRoT transit shape, as well as from priors on stellar populations, planet and multiple stellar systems frequency. We use the fully Bayesian approach developed in the PASTIS analysis software, and conclude that the planet scenario is at least 1400 times more probable than any other false positive scenario. The primary star is a metallic solar-like dwarf, with Ms = 1.099+-0.049 Msun and Rs = 1.136 (+0.038,-0.090) Rsun . The validated planet has a radius of Rp = 4.88 (+0.17,-0.39) RE and mass less than 49 ME. Its mean density is smaller than 2.56 g/cm^3 and orbital period is 9.7566+-0.0012 days. This object, called CoRoT-22 b, adds to a large number of validated Kepler planets. These planets do not have a proper measurement of the mass but allow statistical characterization of the exoplanet population.

preprint2014arXiv

ESPRESSO: The next European exoplanet hunter

The acronym ESPRESSO stems for Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations; this instrument will be the next VLT high resolution spectrograph. The spectrograph will be installed at the Combined-Coudé Laboratory of the VLT and linked to the four 8.2 m Unit Telescopes (UT) through four optical Coudé trains. ESPRESSO will combine efficiency and extreme spectroscopic precision. ESPRESSO is foreseen to achieve a gain of two magnitudes with respect to its predecessor HARPS, and to improve the instrumental radial-velocity precision to reach the 10 cm/s level. It can be operated either with a single UT or with up to four UTs, enabling an additional gain in the latter mode. The incoherent combination of four telescopes and the extreme precision requirements called for many innovative design solutions while ensuring the technical heritage of the successful HARPS experience. ESPRESSO will allow to explore new frontiers in most domains of astrophysics that require precision and sensitivity. The main scientific drivers are the search and characterization of rocky exoplanets in the habitable zone of quiet, nearby G to M-dwarfs and the analysis of the variability of fundamental physical constants. The project passed the final design review in May 2013 and entered the manufacturing phase. ESPRESSO will be installed at the Paranal Observatory in 2016 and its operation is planned to start by the end of the same year.

preprint2014arXiv

Exoplanet Science with the European Extremely Large Telescope. The Case for Visible and Near-IR Spectroscopy at High Resolution

Exoplanet science is booming. In 20 years our knowledge has expanded considerably, from the first discovery of a Hot Jupiter, to the detection of a large population of Neptunes and super-Earths, to the first steps toward the characterization of exoplanet atmospheres. Between today and 2025, the field will evolve at an even faster pace with the advent of several space-based transit search missions, ground-based spectrographs, high-contrast imaging facilities, and the James Webb Space Telescope. Especially the ESA M-class PLATO mission will be a game changer in the field. From 2024 onwards, PLATO will find transiting terrestrial planets orbiting within the habitable zones of nearby, bright stars. These objects will require the power of Extremely Large Telescopes (ELTs) to be characterized further. The technique of ground-based high-resolution spectroscopy is establishing itself as a crucial pathway to measure chemical composition, atmospheric structure and atmospheric circulation in transiting exoplanets. A high-resolution spectrograph covering the visible and near-IR domains, mounted on the European ELT, will be able to detect molecules such as water vapour, carbon dioxide and oxygen in the atmospheres of habitable planets under favourable circumstances. E-ELT HiRES is the perfect ground-based match to the PLATO space mission and represents a unique opportunity for Europe to lead the world into the era of exploration of exoplanets with habitable conditions. HiRES will also be extremely complementary to other E-ELT planned instruments specialising in different kinds of planets, such as METIS and EPICS.

preprint2014arXiv

Impact of micro-telluric lines on precise radial velocities and its correction

Context: In the near future, new instruments such as ESPRESSO will arrive, allowing us to reach a precision in radial-velocity measurements on the order of 10 cm/s. At this level of precision, several noise sources that until now have been outweighed by photon noise will start to contribute significantly to the error budget. The telluric lines that are not neglected by the masks for the radial velocity computation, here called micro-telluric lines, are one such noise source. Aims: In this work we investigate the impact of micro-telluric lines in the radial velocities calculations. We also investigate how to correct the effect of these atmospheric lines on radial velocities. Methods: The work presented here follows two parallel lines. First, we calculated the impact of the micro-telluric lines by multiplying a synthetic solar-like stellar spectrum by synthetic atmospheric spectra and evaluated the effect created by the presence of the telluric lines. Then, we divided HARPS spectra by synthetic atmospheric spectra to correct for its presence on real data and calculated the radial velocity on the corrected spectra. When doing so, one considers two atmospheric models for the synthetic atmospheric spectra: the LBLRTM and TAPAS. Results: We find that the micro-telluric lines can induce an impact on the radial velocities calculation that can already be close to the current precision achieved with HARPS, and so its effect should not be neglected, especially for future instruments such as ESPRESSO. Moreover, we find that the micro-telluric lines' impact depends on factors, such as the radial velocity of the star, airmass, relative humidity, and the barycentric Earth radial velocity projected along the line of sight at the time of the observation.

preprint2014arXiv

Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system

Since the discovery of the transiting super-Earth CoRoT-7b, several investigations have yielded different results for the number and masses of planets present in the system, mainly owing to the star's high level of activity. We re-observed CoRoT-7 in January 2012 with both HARPS and CoRoT, so that we now have the benefit of simultaneous radial-velocity and photometric data. This allows us to use the off-transit variations in the star's light curve to estimate the radial-velocity variations induced by the suppression of convective blueshift and the flux blocked by starspots. To account for activity-related effects in the radial-velocities which do not have a photometric signature, we also include an additional activity term in the radial-velocity model, which we treat as a Gaussian process with the same covariance properties (and hence the same frequency structure) as the light curve. Our model was incorporated into a Monte Carlo Markov Chain in order to make a precise determination of the orbits of CoRoT-7b and CoRoT-7c. We measure the masses of planets b and c to be 4.73 +/- 0.95 Mearth and 13.56 +/- 1.08 Mearth, respectively. The density of CoRoT-7b is (6.61 +/- 1.72)(Rp/1.58 Rearth)^(-3) g.cm^(-3), which is compatible with a rocky composition. We search for evidence of an additional planet d, identified by previous authors with a period close to 9 days. We are not able to confirm the existence of a planet with this orbital period, which is close to the second harmonic of the stellar rotation at around 7.9 days. Using Bayesian model selection we find that a model with two planets plus activity-induced variations is most favoured.

preprint2014arXiv

The GAPS programme with HARPS-N@TNG IV: A planetary system around XO-2S

We performed an intensive radial velocity monitoring of XO-2S, the wide companion of the transiting planet-host XO-2N, using HARPS-N at TNG in the framework of the GAPS programme. The radial velocity measurements indicate the presence of a new planetary system formed by a planet that is slightly more massive than Jupiter at 0.48 au and a Saturn-mass planet at 0.13 au. Both planetary orbits are moderately eccentric and were found to be dynamically stable. There are also indications of a long-term trend in the radial velocities. This is the first confirmed case of a wide binary whose components both host planets, one of which is transiting, which makes the XO-2 system a unique laboratory for understanding the diversity of planetary systems.

preprint2014arXiv

The HARPS search for southern extra-solar planets XXXV. Planetary systems and stellar activity of the M dwarfs GJ 3293, GJ 3341, and GJ 3543

Context. Planetary companions of a fixed mass induce larger amplitude reflex motions around lower-mass stars, which helps make M dwarfs excellent targets for extra-solar planet searches. State of the art velocimeters with $\sim$1m/s stability can detect very low-mass planets out to the habitable zone of these stars. Low-mass, small, planets are abundant around M dwarfs, and most known potentially habitable planets orbit one of these cool stars. Aims. Our M-dwarf radial velocity monitoring with HARPS on the ESO 3.6m telescope at La Silla observatory makes a major contribution to this sample. Methods. We present here dense radial velocity (RV) time series for three M dwarfs observed over $\sim5$ years: GJ 3293 (0.42M$_\odot$), GJ 3341 (0.47M$_\odot$), and GJ 3543 (0.45M$_\odot$). We extract those RVs through minimum $χ^2$ matching of each spectrum against a high S/N ratio stack of all observed spectra for the same star. We then vet potential orbital signals against several stellar activity indicators, to disentangle the Keplerian variations induced by planets from the spurious signals which result from rotational modulation of stellar surface inhomogeneities and from activity cycles. Results. Two Neptune-mass planets - $msin(i)=1.4\pm0.1$ and $1.3\pm0.1M_{nept}$ - orbit GJ 3293 with periods $P=30.60\pm0.02$ d and $P=123.98\pm0.38$ d, possibly together with a super-Earth - $msin(i)\sim7.9\pm1.4M_\oplus$ - with period $P=48.14\pm0.12\;d$. A super-Earth - $msin(i)\sim6.1M_\oplus$ - orbits GJ 3341 with $P=14.207\pm0.007\;d$. The RV variations of GJ 3543, on the other hand, reflect its stellar activity rather than planetary signals.

preprint2014arXiv

The HARPS search for southern extra-solar planets XXXV. The interesting case of HD41248: stellar activity, no planets?

The search for planets orbiting metal-poor stars is of uttermost importance for our understanding of the planet formation models. However, no dedicated searches have been conducted so far for very low mass planets orbiting such objects. Only a few cases of low mass planets orbiting metal-poor stars are thus known. Amongst these, HD41248 is a metal-poor, solar-type star on which a resonant pair of super-Earth like planets has In the present paper we present a new planet search program that is using the HARPS spectrograph to search for Neptunes and Super-Earths orbiting a sample of metal-poor FGK dwarfs. We then present a detailed analysis of an additional 162 radial velocity measurements of HD41248, obtained within this program, with the goal of confirming the existence of the proposed planetary system. We analyzed the precise radial velocities, obtained with the HARPS spectrograph, together with several stellar activity diagnostics and line profile indicators. A careful analysis shows no evidence for the planetary system previously announced. One of the signals, with a period of about 25 days, is shown to be related to the rotational period of the star, and is clearly seen in some of the activity proxies. The remaining signal (P~18 days) could not be convincingly retrieved in the new data set. We discuss possible causes for the complex (evolving) signals observed in the data of HD41248, proposing that they may be explained by the appearance and disappearance of active regions on the surface of a star with strong differential rotation, or by a combination of the sparse data sampling and active region evolution.

preprint2014arXiv

The HARPS search for southern extra-solar planets XXXVI: Five new long-period giant planets and a system update

We describe radial-velocity time series obtained by HARPS on the 3.60 m telescope in La Silla (ESO, Chile) over ten years and report the discovery of five new giant exoplanets in distant orbits; these new planets orbit the stars HD 564, HD 30669, HD 108341, and BD-114672. Their periods range from 492 to 1684 days, semi-major axes range from 1.2 to 2.69 AU, and eccentricities range from 0 to 0.85. Their minimum mass ranges from 0.33 to 3.5 Mjup. We also refine the parameters of two planets announced previously around HD 113538, based on a longer series of measurements. The planets have a period of 663+-8 and 1818+-25 days, orbital eccentricities of 0.14+-0.08 and 0.20+-0.04, and minimum masses of 0.36+-0.04 and 0.93+-0.06 Mjup. Finally, we report the discovery of a new hot-Jupiter planet around an active star, HD 103720; the planet has a period of 4.5557+-0.0001 days and a minimum mass of 0.62+-0.025 Mjup. We discuss the fundamental parameters of these systems and limitations due to stellar activity in quiet stars with typical 2m/s radial velocity precision.

preprint2014arXiv

The well-aligned orbit of WASP-84b: evidence for disc migration

We report the sky-projected orbital obliquity (spin-orbit angle) of WASP-84b, a 0.70-$M_{\rm Jup}$ planet in a 8.52-day orbit around a G9V/K0V star, to be $λ= 0.3 \pm 1.7^\circ$. We obtain a true obliquity of $ψ= 14.8 \pm 8.0^\circ$ from a measurement of the inclination of the stellar spin axis with respect to the sky plane. Due to the young age and the weak tidal forcing of the system, we suggest that the orbit of WASP-84b is unlikely to have both realigned and circularised from the misaligned and/or eccentric orbit likely to have arisen from high-eccentricity migration. Therefore we conclude that the planet probably migrated via interaction with the protoplanetary disc. This would make it the first short-orbit, giant planet to have been shown to have migrated via this pathway. Further, we argue that the distribution of obliquities for planets orbiting cool stars ($T_{\rm eff}$ < 6250 K) suggests that high-eccentricity migration is an important pathway for the formation of short-orbit, giant planets.

preprint2014arXiv

Three planetary companions around M67 stars

For the past six years we have carried out a search for massive planets around main sequence and evolved stars in the open cluster (OC) M67, using radial velocity (RV) measurements obtained with HARPS at ESO (La Silla), SOPHIE at OHP and HRS at HET. Additional RV data come from CORALIE at the Euler Swiss Telescope. We aim to perform a long-term study on giant planet formation in open clusters and determine how it depends on stellar mass and chemical composition. We report the detection of three new extrasolar planets: two in orbit around the two G dwarfs YBP1194 and YBP1514, and one around the evolved star S364. The orbital solution for YBP1194 yields a period of 6.9 days, an eccentricity of 0.24, and a minimum mass of 0.34 Mj. YBP1514 shows periodic RV variations of 5.1 days, a minimum mass of 0.40 Mj, and an eccentricity of 0.39. The best Keplerian solution for S364 yields a period of 121.7 days, an eccentricity of 0.35 and a minimum mass of 1.54 Mj. An analysis of H_alpha core flux measurements as well as of the line bisectors spans revealed no correlation with the RV periods, indicating that the RV variations are best explained by the presence of a planetary companion. Remarkably, YBP1194 is one of the best solar twins identified so far, and YBP1194b is the first planet found around a solar twin that belongs to a stellar cluster. In contrast with early reports and in agreement with recent findings, our results show that massive planets around stars of open clusters are as frequent as those around field stars.

preprint2014arXiv

Transit confirmation and improved stellar and planet parameters for the super-Earth HD 97658 b and its host star

Super-Earths transiting nearby bright stars are key objects that simultaneously allow for accurate measurements of both their mass and radius, providing essential constraints on their internal composition. We present here the confirmation, based on Spitzer transit observations, that the super-Earth HD 97658 b transits its host star. HD 97658 is a low-mass ($M_*=0.77\pm0.05\,M_{\odot}$) K1 dwarf, as determined from the Hipparcos parallax and stellar evolution modeling. To constrain the planet parameters, we carry out Bayesian global analyses of Keck-HIRES radial velocities, and MOST and Spitzer photometry. HD 97658 b is a massive ($M_P=7.55^{+0.83}_{-0.79} M_{\oplus}$) and large ($R_{P} = 2.247^{+0.098}_{-0.095} R_{\oplus}$ at 4.5 $μ$m) super-Earth. We investigate the possible internal compositions for HD 97658 b. Our results indicate a large rocky component, by at least 60% by mass, and very little H-He components, at most 2% by mass. We also discuss how future asteroseismic observations can improve the knowledge of the HD 97658 system, in particular by constraining its age. Orbiting a bright host star, HD 97658 b will be a key target for coming space missions TESS, CHEOPS, PLATO, and also JWST, to characterize thoroughly its structure and atmosphere.

preprint2014arXiv

Transiting exoplanets from the CoRoT space mission: XXVI. CoRoT-24: A transiting multi-planet system

We present the discovery of a candidate multiply-transiting system, the first one found in the CoRoT mission. Two transit-like features with periods of 5.11 and 11.76d are detected in the CoRoT light curve, around a main sequence K1V star of r=15.1. If the features are due to transiting planets around the same star, these would correspond to objects of 3.7$\pm$0.4 and 5.0$\pm$0.5 R_earth respectively. Several radial velocities serve to provide an upper limit of 5.7 M_earth for the 5.11~d signal, and to tentatively measure a mass of 28$^{+11}_{-11}$ M_earth for the object transiting with a 11.76~d period. These measurements imply low density objects, with a significant gaseous envelope. The detailed analysis of the photometric and spectroscopic data serve to estimate the probability that the observations are caused by transiting Neptune-sized planets as $>$26$\times$ higher than a blend scenario involving only one transiting planet, and $>$900$\times$ higher than a scenario involving two blends and no planets. The radial velocities show a long term modulation that might be attributed to a 1.5 M_jup planet orbiting at 1.8~A.U. from the host, but more data are required to determine the precise orbital parameters of this companion.

preprint2013arXiv

A Community Science Case for E-ELT HIRES

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

preprint2013arXiv

Line-profile variations in radial-velocity measurements: Two alternative indicators for planetary searches

Aims. We introduce two methods to identify false-positive planetary signals in the context of radial-velocity exoplanet searches. The first is the bi-Gaussian cross-correlation function fitting, and the second is the measurement of asymmetry in radial-velocity spectral line information content, Vasy. Methods. We make a systematic analysis of the most used common line profile diagnosis, Bisector Inverse Slope and Velocity Span, along with the two proposed ones. We evaluate all these diagnosis methods following a set of well-defined common criteria and using both simulated and real data. We apply them to simulated cross-correlation functions created with the program SOAP and which are affected by the presence of stellar spots, and to real cross-correlation functions, calculated from HARPS spectra, for stars with a signal originating both in activity and created by a planet. Results. We demonstrate that the bi-Gaussian method allows a more precise characterization of the deformation of line profiles than the standard bisector inverse slope. The calculation of the deformation indicator is simpler and its interpretation more straightforward. More importantly, its amplitude can be up to 30% larger than that of the bisector span, allowing the detection of smaller-amplitude correlations with radial-velocity variations. However, a particular parametrization of the bisector inverse slope is shown to be more efficient on high-signal-to-noise data than both the standard bisector and the bi-Gaussian. The results of the Vasy method show that this indicator is more effective than any of the previous ones, being correlated with the radial-velocity with more significance for signals resulting from a line deformation. Moreover, it provides a qualitative advantage over the bisector, showing significant correlations with RV for active stars for which bisector analysis is inconclusive. (abridged)

preprint2013arXiv

On the long-term correlation between the flux in the Ca II H & K and Halpha lines for FGK stars

The re-emission in the cores of the Ca II H & K and H$α$ lines, are well known proxies of stellar activity. However, these activity indices probe different activity phenomena, the first being more sensitive to plage variation, while the other one being more sensitive to filaments. In this paper we study the long-term correlation between $\log R'_{HK}$ and $\log I_{Hα}$, two indices based on the Ca II H & K and H$α$ lines respectively, for a sample of 271 FGK stars using measurements obtained over a $\sim$9 year time span. Because stellar activity is one of the main obstacles to the detection of low-mass and long-period planets, understanding further this activity index correlation can give us some hints about the optimal target to focus on, and ways to correct for these activity effects. We found a great variety of long-term correlations between $\log R'_{HK}$ and $\log I_{Hα}$. Around 20% of our sample has strong positive correlation between the indices while about 3% show strong negative correlation. These fractions are compatible with those found for the case of early-M dwarfs. Stars exhibiting a positive correlation have a tendency to be more active when compared to the median of the sample, while stars showing a negative correlation are more present among higher metallicity stars. There is also a tendency for the positively correlated stars to be more present among the coolest stars, a result which is probably due to the activity level effect on the correlation. Activity level and metallicity seem therefore to be playing a role on the correlation between $\log R'_{HK}$ and $\log I_{Hα}$. Possible explanations based on the influence of filaments for the diversity in the correlations between these indices are discussed in this paper.

preprint2013arXiv

The CORALIE survey for southern extrasolar planets XVII. New and updated long period and massive planets

Since 1998, a planet-search program around main sequence stars within 50 pc in the southern hemisphere, is carried out with the CORALIE echelle spectrograph at La Silla Observatory. With an observing time span of more than 14 years, the CORALIE survey is now able to unveil Jovian planets on Jupiter's period domain. This growing period-interval coverage is important regarding to formation and migration models since observational constraints are still weak for periods beyond the ice line. Long-term precise Doppler measurements with the CORALIE echelle spectrograph, together with a few additional observations made with the HARPS spectrograph on the ESO 3.6m telescope, reveal radial velocity signatures of massive planetary companions in long period orbits. In this paper we present seven new planets orbiting HD27631, HD98649, HD106515A, HD166724, HD196067, HD219077, and HD220689 together with the CORALIE orbital parameters for three already known planets around HD10647, HD30562, and HD86226. The period range of the new planetary companions goes from 2200 to 5500 days and covers a mass domain between 1 and 10.5 MJup. Surprisingly, five of them present quite high eccentricities above e>0.57. A pumping scenario by Kozai mechanism may be invoked for HD106515Ab and HD196067b which are both orbiting stars in multiple systems. As the presence of a third massive body can't be inferred from the data of HD98649b, HD166724b, and HD219077b, the origin of the eccentricity of these systems remains unknown. Except for HD10647b, no constraint on the upper mass of the planets is provided by Hipparcos astrometric data. Finally it is interesting to note that the hosts of these long period planets show no metallicity excess.

preprint2013arXiv

The GAPS programme with HARPS-N at TNG. I: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting system Qatar-1

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

preprint2013arXiv

The GAPS Programme with HARPS-N@TNG II: No giant planets around the metal-poor star HIP 11952

In the context of the program Global Architecture of Planetary Systems (GAPS), we have performed radial velocity monitoring of the metal-poor star HIP 11952 on 35 nights over about 150 days using the newly installed high resolution spectrograph HARPS-N at the TNG and HARPS at ESO 3.6m telescope. The radial velocities show a scatter of 7 m/s, compatible with the measurement errors for such a moderately warm metal-poor star (Teff = 6040+-120K; [Fe/H] =-1.9+-0.1). We then exclude the presence of the two giant planets with periods of 6.95+-0.01 d and 290.0+-16.2 d and radial velocity semi-amplitudes of 100.3+-19.4 m/s and 105.2+-14.7 m/s, respectively, which had recently been announced. This result is important considering that HIP 11952 was thought to be the most metal-poor star hosting a planetary system with giant planets, thus challenging some models of planet formation.

preprint2013arXiv

The HARPS search for southern extra-solar planets. XXXIV. A planetary system around the nearby M dwarf GJ163, with a super-Earth possibly in the habitable zone

The meter-per-second precision achieved by today velocimeters enables the search for 1-10 M_Earth planets in the habitable zone of cool stars. This paper reports on the detection of 3 planets orbiting GJ163 (HIP19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ163 over a period of 8 years. We searched the RV time series for coherent signals and found 5 distinct periodic variabilities. We investigated the stellar activity and casted doubts on the planetary interpretation for 2 signals. Before more data can be acquired we concluded that at least 3 planets are orbiting GJ163. They have orbital periods of P_b=8.632+-0.002, P_c=25.63+-0.03 and P_d=604+-8 days and minimum masses msini = 10.6+-0.6, 6.8+-0.9, and 29+-3 M_Earth, respectively. We hold our interpretations for the 2 additional signals with periods P_(e)=19.4 and P_(f)=108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. GJ163c, in particular, is a super-Earth with an equilibrium temperature of T_eq = (302+-10) (1-A)^(1/4) K and may lie in the so called habitable zone for albedo values (A=0.34-0.89) moderately higher than that of Earth (A_Earth=0.2-0.3).

preprint2013arXiv

The HARPS search for southern extrasolar planets: XXXIII. New multi-planet systems in the HARPS volume limited sample: a super-Earth and a Neptune in the habitable zone

The vast diversity of planetary systems detected to date is defying our capability of understanding their formation and evolution. Well-defined volume-limited surveys are the best tool at our disposal to tackle the problem, via the acquisition of robust statistics of the orbital elements. We are using the HARPS spectrograph to conduct our survey of ~850 nearby solar-type stars, and in the course of the past nine years we have monitored the radial velocity of HD103774, HD109271, and BD-061339. In this work we present the detection of five planets orbiting these stars, with m*sin(i) between 0.6 and 7 Neptune masses, four of which are in two multiple systems, comprising one super-Earth and one planet within the habitable zone of a late-type dwarf. Although for strategic reasons we chose efficiency over precision in this survey, we have the capability to detect planets down to the Neptune and super-Earth mass range, as well as multiple systems, provided that enough data points are made available.

preprint2013arXiv

WASP-80b: a gas giant transiting a cool dwarf

We report the discovery of a planet transiting the star WASP-80 (1SWASP J201240.26-020838.2; 2MASS J20124017-0208391; TYC 5165-481-1; BPM 80815; V=11.9, K=8.4). Our analysis shows this is a 0.55 +/- 0.04 Mjup, 0.95 +/- 0.03 Rjup gas giant on a circular 3.07 day orbit around a star with a spectral type between K7V and M0V. This system produces one of the largest transit depths so far reported, making it a worthwhile target for transmission spectroscopy. We find a large discrepancy between the v sin i inferred from stellar line broadening and the observed amplitude of the Rossiter-McLaughlin effect. This can be understood either by an orbital plane nearly perpendicular to the stellar spin or by an additional, unaccounted for source of broadening.

preprint2012arXiv

A hot Uranus transiting the nearby M dwarf GJ3470. Detected with HARPS velocimetry. Captured in transit with TRAPPIST photometry

We report on the discovery of GJ3470b, a transiting hot Uranus of mass m_p = 14.0+-1.8 Mearth, radius R_p = 4.2+-0.6 Rearth and period P=3.3371+-0.0002 day. Its host star is a nearby (d=25.2+-2.9pc) M1.5 dwarf of mass M_s=0.54+-0.07 Msol and radius R_s=0.50+-0.06 Rsol. The detection originates from a radial-velocity campaign with HARPS that focused on the search for short-period planets orbiting M dwarfs. Once the planet was discovered and the transit-search window narrowed to about 10% of an orbital period, a photometric search started with TRAPPIST and quickly detected the ingress of the planet. Additional observations with TRAPPIST, EulerCam and NITES definitely confirmed the transiting nature of GJ3470b and allow for the determination of its true mass and radius. The star's visible or infrared brightness (V=12.3, K=8.0 mag), together with a large eclipse depth D=0.57+-0.05%, ranks GJ3470b among the most favorable planets for follow-up characterizations.

preprint2012arXiv

Comparing HARPS and Kepler surveys: The alignment of multiple-planet systems

Aims. We study a subset of the planetary population characterized both by HARPS and Kepler surveys. We compare the statistical properties of planets in systems with m.sin i >5-10 M_Earth and R>2 R_Earth. If we assume that the underlying population has the same characteristics, the different detection sensitivity to the orbital inclination relative to the line of sight allows us to probe the planets' mutual inclination. Methods. We considered the frequency of systems with one, two and three planets as dictated by HARPS data. We used Kepler's planetary period and host mass and radii distributions (corrected from detection bias) to model planetary systems in a simple yet physically plausible way. We then varied the mutual inclination between planets in a system according to different prescriptions (completely aligned, Rayleigh distributions and isotropic) and compared the transit frequencies with one, two or three planets with those measured by Kepler. Results. The results show that the two datasets are compatible, a remarkable result especially because there are no tunable knobs other than the assumed inclination distribution. For m.sin i cutoffs of 7-10 M_Earth, which are those expected to correspond to the radius cutoff of 2 R_Earth, we conclude that the results are better described by a Rayleigh distribution with mode of 1 deg or smaller. We show that the best-fit scenario only becomes a Rayleigh distribution with mode of 5 deg if we assume a rather extreme mass-radius relationship for the planetary population. Conclusions. These results have important consequences for our understanding of the role of several proposed formation and evolution mechanisms. They confirm that planets are likely to have been formed in a disk and show that most planetary systems evolve quietly without strong angular momentum exchanges (abridged).

preprint2012arXiv

Improved precision on the radius of the nearby super-Earth 55 Cnc e

We report on new transit photometry for the super-Earth 55 Cnc e obtained with Warm Spitzer/IRAC at 4.5 microns. An individual analysis of these new data leads to a planet radius of 2.21-0.16+0.15 Rearth, in good agreement with the values previously derived from the MOST and Spitzer transit discovery data. A global analysis of both Spitzer transit time-series improves the precision on the radius of the planet at 4.5 microns to 2.20+-0.12 Rearth. We also performed an independent analysis of the MOST data, paying particular attention to the influence of the systematic effects of instrumental origin on the derived parameters and errors by including them in a global model instead of performing a preliminary detrending-filtering processing. We deduce an optical planet radius of 2.04+0.15 Rearth from this reanalysis of MOST data, which is consistent with the previous MOST result and with our Spitzer infrared radius. Assuming the achromaticity of the transit depth, we performed a global analysis combining Spitzer and MOST data that results in a planet radius of 2.17+-0.10 Rearth (13,820+-620 km). These results point to 55 Cnc e having a gaseous envelope overlying a rocky nucleus, in agreement with previous works. A plausible composition for the envelope is water which would be in super-critical form given the equilibrium temperature of the planet.

preprint2012arXiv

Long-term magnetic activity of a sample of M-dwarf stars from the HARPS program II. Activity and radial velocity

Due to their low mass and luminosity, M dwarfs are ideal targets if one hopes to find low-mass planets similar to Earth by using the radial velocity (RV) method. However, stellar magnetic cycles could add noise or even mimic the RV signal of a long-period companion. Following our previous work that studied the correlation between activity cycles and long-term RV variations for K dwarfs we now expand that research to the lower-end of the main sequence. Our objective is to detect any correlations between long-term activity variations and the observed RV of a sample of M dwarfs. We used a sample of 27 M-dwarfs with a median observational timespan of 5.9 years. The cross-correlation function (CCF) with its parameters RV, bisector inverse slope (BIS), full-width-at-half- maximum (FWHM) and contrast have been computed from the HARPS spectrum. The activity index have been derived using the Na I D doublet. These parameters were compared with the activity level of the stars to search for correlations. We detected RV variations up to ~5 m/s that we can attribute to activity cycle effects. However, only 36% of the stars with long-term activity variability appear to have their RV affected by magnetic cycles, on the typical timescale of ~6 years. Therefore, we suggest a careful analysis of activity data when searching for extrasolar planets using long-timespan RV data.

preprint2012arXiv

Metallicities for 6 nearby open clusters from high-resolution spectra of giant stars. [Fe/H] values for a planet search sample

We present a study of the stellar parameters and iron abundances of 18 giant stars in 6 open clusters. The analysis was based on high-resolution and high-S/N spectra obtained with the UVES spectrograph (VLT-UT2). The results complement our previous study where 13 clusters were already analyzed. The total sample of 18 clusters is part of a program to search for planets around giant stars. The results show that the 18 clusters cover a metallicity range between -0.23 and +0.23 dex. Together with the derivation of the stellar masses, these metallicities will allow the metallicity and mass effects to be disentangled when analyzing the frequency of planets as a function of these stellar parameters.

preprint2012arXiv

Overabundance of alpha-elements in exoplanet host stars

We present the results for a chemical abundance analysis between planet-hosting and stars without planets for 12 refractory elements for a total of 1111 nearby FGK dwarf stars observed within the context of the HARPS GTO programs. Of these stars, 109 are known to harbour high-mass planetary companions and 26 stars are hosting exclusively Neptunians and super-Earths. We found that the [X/Fe] ratios for Mg, Al, Si, Sc, and Ti both for giant and low-mass planet hosts are systematically higher than those of comparison stars at low metallicities ([Fe/H] < from -0.2 to 0.1 dex depending on the element). The most evident discrepancy between planet-hosting and stars without planets is observed for Mg. Our data suggest that the planet incidence is greater among the thick disk population than among the thin disk for mettallicities bellow -0.3 dex. After examining the [alpha/Fe] trends of the planet host and non-host samples we conclude that a certain chemical composition, and not the Galactic birth place of the stars, is the determinating factor for that. The inspection of the Galactic orbital parameters and kinematics of the planet-hosting stars shows that Neptunian hosts tend to belong to the "thicker" disk compared to their high-mass planet-hosting counterparts.We also found that Neptunian hosts follow the distribution of high-alpha stars in the UW vs V velocities space, but they are more enhanced in Mg than high-alpha stars without planetary companions. Our results indicate that some metals other than iron may also have an important contribution to planet formation if the amount of iron is low. These results may provide strong constraints for the models of planet formation, especially for planets with low mass.

preprint2012arXiv

Search for giant planets in M67 I. Overview

Precise stellar radial velocities are used to search for massive (Jupiter masses or higher) exoplanets around the stars of the open cluster M67. We aim to obtain a census of massive exoplanets in a cluster of solar metallicity and age in order to study the dependence of planet formation on stellar mass and to compare in detail the chemical composition of stars with and without planets. This first work presents the sample and the observations, discusses the cluster characteristics and the radial velocity (RV) distribution of the stars, and individuates the most likely planetary host candidates. We observed a total of 88 main-sequence stars, subgiants, and giants all highly probable members of M67, using four telescopes and instrument combinations. We investigate whether exoplanets are present by obtaining radial velocities with precisions as good as 10 m/s. To date, we have performed 680 single observations (Dec. 2011) and a preliminary analysis of data, spanning a period of up to eight years. Although the sample was pre-selected to avoid the inclusion of binaries, we identify 11 previously unknown binary candidates. Eleven stars clearly displayed larger RV variability and these are candidates to host long-term substellar companions. The average RV is also independent of the stellar magnitude and evolutionary status, confirming that the difference in gravitational redshift between giants and dwarfs is almost cancelled by the atmospheric motions. We use the subsample of solar-type stars to derive a precise true RV for this cluster. We finally create a catalog of binaries and use it to clean the color magnitude diagram (CMD). As conclusion, by pushing the search for planets to the faintest possible magnitudes, it is possible to observe solar analogues in open clusters, and we propose 11 candidates to host substellar companions.

preprint2012arXiv

The HARPS search for southern extra-solar planets XXXV. Super-Earths around the M-dwarf neighbors Gl433 and Gl667C

M dwarfs have been found to often have super-Earth planets with short orbital periods. Such stars are thus preferential targets in searches for rocky or ocean planets in the solar neighbourhood. In a recent paper (Bonfils et al. 2011), we announced the discovery of respectively 1 and 2 low mass planets around the M1.5V stars Gl433 and Gl667C. We found those planets with the HARPS spectrograph on the ESO~3.6-m telescope at La Silla Observatory, from observations obtained during the Guaranteed Time Observing program of that instrument. We have obtained additional HARPS observations of those two stars, for a total of respectively 67 and 179 Radial Velocity measurements for Gl433 and Gl667C, and present here an orbital analysis of those extended data sets and our main conclusion about both planetary systems. One of the three planets, Gl667Cc, has a mass of only M2.sin(i)~4.25 M_earth and orbits in the central habitable zone of its host star. It receives just 10% less stellar energy from Gl667C than the Earth receives from the Sun. However planet evolution in habitable zone can be very different if the host star is a M dwarf or a solar-like star, without necessarily questioning the presence of water. The two other planets, Gl433b and Gl667Cb, both have M2.sin(i) of ~5.5 M_earth and periods of ~7 days. The Radial Velocity measurements of both stars contain longer time scale signals, which we fit as longer period Keplerians. For Gl433 that signal probably originates in a Magnetic Cycle, while a longer time span will be needed to conclude for Gl667C. The metallicity of Gl433 is close to solar, while Gl667C is metal poor with [Fe/H] ~ -0.6. This reinforces the recent conclusion that the occurence of Super-Earth planets does not strongly correlate with stellar metallicity.

preprint2012arXiv

The SOPHIE search for northern extrasolar planets. V. Follow-up of ELODIE candidates: Jupiter-analogs around Sun-like stars

We present radial-velocity measurements obtained in a programs underway to search for extrasolar planets with the spectrograph SOPHIE at the 1.93-m telescope of the Haute-Provence Observatory. Targets were selected from catalogs observed with ELODIE, mounted previously at the telescope, in order to detect long-period planets with an extended database close to 15 years. Two new Jupiter-analog candidates are reported to orbit the bright stars HD150706 and HD222155 in 16.1 and 10.9 yr at 6.7 (+4.0,-1.4) and 5.1(+0.6,-0.7) AU and to have minimum masses of 2.71 (+1.44,-0.66) and 1.90 (+0.67,-0.53) M_Jup, respectively. Using the measurements from ELODIE and SOPHIE, we refine the parameters of the long-period planets HD154345b and HD89307b, and publish the first reliable orbit for HD24040b. This last companion has a minimum mass of 4.01 +/- 0.49 M_Jup orbiting its star in 10.0 yr at 4.92 +/- 0.38 AU. Moreover, the data provide evidence of a third bound object in the HD24040 system. With a surrounding dust debris disk, HD150706 is an active G0 dwarf for which we partially corrected the effect of the stellar spot on the SOPHIE radial-velocities. HD222155 is an inactive G2V star. On the basis of the previous findings of Lovis and collaborators and since no significant correlation between the radial-velocity variations and the activity index are found in the SOPHIE data, these variations are not expected to be only due to stellar magnetic cycles. Finally, we discuss the main properties of this new population of long-period Jupiter-mass planets, which for the moment, consists of fewer than 20 candidates. These stars are preferential targets either for direct-imaging or astrometry follow-up to constrain the system parameters and for higher precision radial-velocity to search for lower mass planets, aiming to find a Solar System twin.

preprint2011arXiv

A short-period super-Earth orbiting the M2.5 dwarf GJ3634. Detection with Harps velocimetry and transit search with Spitzer photometry

We report on the detection of GJ3634b, a super-Earth of mass m sin i = 7.0 +/-0.9 Mearth and period P = 2.64561 +/- 0.00066 day. Its host star is a M2.5 dwarf, has a mass of 0.45+/-0.05 Msun, a radius of 0.43+/-0.03 Rsun and lies 19.8+/-0.6 pc away from our Sun. The planet is detected after a radial-velocity campaign using the ESO/Harps spectrograph. GJ3634b had an a priori geometric probability to undergo transit of ~7% and, if telluric in composition, a non-grazing transit would produce a photometric dip of <~0.1%. We therefore followed-up upon the RV detection with photometric observations using the 4.5-micron band of the IRAC imager onboard Spitzer. Our six-hour long light curve excludes that a transit occurs for 2 sigma of the probable transit window, decreasing the probability that GJ3634b undergoes transit to ~0.5%.

preprint2011arXiv

Comparing radial velocities of atmospheric lines with radiosonde measurements

The precision of radial velocity (RV) measurements depends on the precision attained on the wavelength calibration. One of the available options is using atmospheric lines as a natural, freely available wavelength reference. Figueira et al. (2010) measured the RV of O2 lines using HARPS and showed that the scatter was only of ~10 m/s over a timescale of 6 yr. Using a simple but physically motivated empirical model, they demonstrated a precision of 2 m/s, roughly twice the average photon noise contribution. In this paper we take advantage of a unique opportunity to confirm the sensitivity of the telluric absorption lines RV to different atmospheric and observing conditions: by means of contemporaneous in-situ wind measurements by radiosondes. The RV model fitting yielded similar results to that of Figueira et al. (2010), with lower wind magnitude values and varied wind direction. The probes confirmed the average low wind magnitude and suggested that the average wind direction is a function of time as well. The two approaches deliver the same results in what concerns wind magnitude and agree on wind direction when fitting is done in segments of a couple of hours. Statistical tests show that the model provides a good description of the data on all timescales, being always preferable to not fitting any atmospheric variation. The smaller the timescale on which the fitting can be performed (down to a couple of hours), the better the description of the real physical parameters. We conclude then that the two methods deliver compatible results, down to better than 5 m/s and less than twice the estimated photon noise contribution on O2 lines RV measurement. However, we cannot rule out that parameters alpha and gamma (dependence on airmass and zero-point, respectively) have a dependence on time or exhibit some cross-talk with other parameters (abridged).

preprint2011arXiv

Detection of a transit of the super-Earth 55 Cnc e with Warm Spitzer

We report on the detection of a transit of the super-Earth 55 Cnc e with warm Spitzer in IRAC's 4.5-micron band. Our MCMC analysis includes an extensive modeling of the systematic effects affecting warm Spitzer photometry, and yields a transit depth of 410 +- 63 ppm, which translates to a planetary radius of 2.08 +- 0.16 R_Earth as measured in IRAC 4.5-micron channel. A planetary mass of 7.81 +- 0.58 M_Earth is derived from an extensive set of radial-velocity data, yielding a mean planetary density of 4.8 +- 1.3 g cm-3. Thanks to the brightness of its host star (V = 6, K = 4), 55 Cnc e is a unique target for the thorough characterization of a super-Earth orbiting around a solar-type star.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

preprint2011arXiv

Planetary detection limits taking into account stellar noise. II. Effect of stellar spot groups on radial-velocities

The detection of small mass planets with the radial-velocity technique is now confronted with the interference of stellar noise. HARPS can now reach a precision below the meter-per-second, which corresponds to the amplitudes of different stellar perturbations, such as oscillation, granulation, and activity. Solar spot groups induced by activity produce a radial-velocity noise of a few meter-per-second. The aim of this paper is to simulate this activity and calculate detection limits according to different observational strategies. Based on Sun observations, we reproduce the evolution of spot groups on the surface of a rotating star. We then calculate the radial-velocity effect induced by these spot groups as a function of time. Taking into account oscillation, granulation, activity, and a HARPS instrumental error of 80 cm/s, we simulate the effect of different observational strategies in order to efficiently reduce all sources of noise. Applying three measurements per night of 10 minutes every three days, 10 nights a month seems the best tested strategy. Depending on the level of activity considered, from log(R'_HK)= -5 to -4.75, this strategy would allow us to find planets of 2.5 to 3.5 Earth masses in the habitable zone of a K1V dwarf. Using Bern's model of planetary formation, we estimate that for the same range of activity level, 15 to 35 % of the planets between 1 and 5 Earth masses and with a period between 100 and 200 days should be found with HARPS. A comparison between the performance of HARPS and ESPRESSO is also emphasized by our simulations. Using the same optimized strategy, ESPRESSO could find 1.3 Earth mass planets in the habitable zone of early-K dwarfs. In addition, 80 % of planets with mass between 1 and 5 Earth masses and with a period between 100 and 200 days could be detected.

preprint2011arXiv

Planetary transit candidates in the CoRoT LRa01 field

Context: CoRoT is a pioneering space mission whose primary goals are stellar seismology and extrasolar planets search. Its surveys of large stellar fields generate numerous planetary candidates whose lightcurves have transit-like features. An extensive analytical and observational follow-up effort is undertaken to classify these candidates. Aims: The list of planetary transit candidates from the CoRoT LRa01 star field in the Monoceros constellation towards the Galactic anti-center is presented. The CoRoT observations of LRa01 lasted from 24 October 2007 to 3 March 2008. Methods: 7470 chromatic and 3938 monochromatic lightcurves were acquired and analysed. Instrumental noise and stellar variability were treated with several filtering tools by different teams from the CoRoT community. Different transit search algorithms were applied to the lightcurves. Results: Fifty-one stars were classified as planetary transit candidates in LRa01. Thirty-seven (i.e., 73 % of all candidates) are "good" planetary candidates based on photometric analysis only. Thirty-two (i.e., 87 % of the "good" candidates) have been followed-up. At the time of this writing twenty-two cases have been solved and five planets have been discovered: three transiting hot-Jupiters (CoRoT-5b, CoRoT-12b, and CoRoT-21b), the first terrestrial transiting planet (CoRoT-7b), and another planet in the same system (CoRoT-7c, detected by radial velocity survey only). Evidences of another non-transiting planet in the CoRoT-7 system, namely CoRoT-7d, have been recently found.

preprint2011arXiv

Report of IAU Commission 30 on Radial Velocities (2009-2012)

Brief summaries are given of the following subjects of interest to IAU Commission 30: Large-scale radial-velocity surveys; The role of radial-velocity measurements in studies of stellar angular momentum evolution and stellar age; Radial velocities in open clusters; Toward higher radial-velocity precision; High-precision radial velocities applied to studies of binary stars; Doppler boosting effect; Working groups (Stellar radial velocity bibliography; Radial velocity standards; Catalogue of orbital elements of spectroscopic binaries [SB9]).

preprint2011arXiv

Spectroscopic characterization of a sample of metal-poor solar-type stars from the HARPS planet search program. Precise spectroscopic parameters and mass estimation

Stellar metallicity strongly correlates with the presence of planets and their properties. To check for new correlations between stars and the existence of an orbiting planet, we determine precise stellar parameters for a sample of metal-poor solar-type stars. This sample was observed with the HARPS spectrograph and is part of a program to search for new extrasolar planets. The stellar parameters were determined using an LTE analysis based on equivalent widths (EW) of iron lines and by imposing excitation and ionization equilibrium. The ARES code was used to allow automatic and systematic derivation of the stellar parameters. Precise stellar parameters and metallicities were obtained for 97 low metal-content stars. We also present the derived masses, luminosities, and new parallaxes estimations based on the derived parameters, and compare our spectroscopic parameters with an infra-red flux method calibration to check the consistency of our method in metal poor stars. Both methods seems to give the same effective temperature scale. Finally we present a new calibration for the temperature as a function of \textit{B-V} and [Fe/H]. This was obtained by adding these new metal poor stars in order to increase the range in metallicity for the calibration. The standard deviation of this new calibration is $\sim$ 50 K.

preprint2011arXiv

The HARPS search for Earth-like planets in the habitable zone: I -- Very low-mass planets around HD20794, HD85512 and HD192310

In 2009 we started an intense radial-velocity monitoring of a few nearby, slowly-rotating and quiet solar-type stars within the dedicated HARPS-Upgrade GTO program. The goal of this campaign is to gather very-precise radial-velocity data with high cadence and continuity to detect tiny signatures of very-low-mass stars that are potentially present in the habitable zone of their parent stars. Ten stars were selected among the most stable stars of the original HARPS high-precision program that are uniformly spread in hour angle, such that three to four of them are observable at any time of the year. For each star we recorded 50 data points spread over the observing season. The data points consist of three nightly observations with a total integration time of 10 minutes each and are separated by two hours. This is an observational strategy adopted to minimize stellar pulsation and granulation noise. We present the first results of this ambitious program. The radial-velocity data and the orbital parameters of five new and one confirmed low-mass planets around the stars HD20794, HD85512, and HD192310 are reported and discussed, among which is a system of three super-Earths and one that harbors a 3.6 Earth-mass planet at the inner edge of the habitable zone. This result already confirms previous indications that low-mass planets seem to be very frequent around solar-type stars and that this may occur with a frequency higher than 30%

preprint2011arXiv

The HARPS search for southern extra-solar planets XXXI. The M-dwarf sample

(Abridged) Searching for planets around stars with different masses probes the outcome of planetary formation for different initial conditions. This drives observations of a sample of 102 southern nearby M dwarfs, using a fraction of our guaranteed time on the ESO/HARPS spectrograph (Feb. 11th, 2003 to Apr. 1st 2009). This paper makes available the sample's time series, presents their precision and variability. We apply systematic searches and diagnostics to discriminate whether the observed Doppler shifts are caused by stellar surface inhomogeneities or by the radial pull of orbiting planets. We recover the planetary signals corresponding to 9 planets already announced by our group (Gl176b, Gl581b, c, d & e, Gl674b, Gl433b, Gl 667Cb and c). We present radial velocities that confirm GJ 849 hosts a Jupiter-mass planet, plus a long-term radial-velocity variation. We also present RVs that precise the planetary mass and period of Gl 832b. We detect long-term RV changes for Gl 367, Gl 680 and Gl 880 betraying yet unknown long-period companions. We identify candidate signals in the radial-velocity time series and demonstrate they are most probably caused by stellar surface inhomogeneities. Finally, we derive a first estimate of the occurrence of M-dwarf planets as a function of their minimum mass and orbital period. In particular, we find that giant planets (m sin i = 100-1,000 Mearth) have a low frequency (e.g. f<1% for P=1-10 d and f=0.02^{+0.03}_{-0.01} for P=10-100 d), whereas super-Earths (m sin i = 1-10 Mearth) are likely very abundant (f=0.36^{+0.25}_{-0.10} for P=1-10 d and f=0.35^{+0.45}_{-0.11} for P=10-100 d). We also obtained eta_earth=0.41^{+0.54}_{-0.13}, the frequency of habitable planets orbiting M dwarfs (1<m sin i<10 Mearth). For the first time, eta_earth is a direct measure and not a number extrapolated from the statistic of more massive and/or shorter-period planets.

preprint2011arXiv

The HARPS search for southern extra-solar planets XXXII. Only 4 planets in the Gl~581 system

The Gl 581 planetary system has generated wide interest, because its 4 planets include both the lowest mass planet known around a main sequence star other than the Sun and the first super-Earth planet in the habitable zone of its star. A recent paper announced the possible discovery of two additional super-Earth planets in that system, one of which would be in the middle of the habitable zone of Gl 581. The statistical significance of those two discoveries has, however, been questioned. We have obtained 121 new radial velocity measurements of Gl 581 with the HARPS spectrograph on the ESO 3.6 m telescope, and analyse those together with our previous 119 measurements of that star to examine these potential additional planets. We find that neither is likely to exist with their proposed parameters. We also obtained photometric observations with the 2.5 m Isaac Newton Telescope during a potential transit of the inner planet, Gl 581e, which had a 5% geometric transit probability. Those observations exclude transits for planet densities under 4 times the Earth density within -0.2 sigma to +2.7 sigma of the predicted transit center.

preprint2011arXiv

The HARPS search for southern extra-solar planets XXXIV. Occurrence, mass distribution and orbital properties of super-Earths and Neptune-mass planets

We report on the results of an 8-year survey carried out at the La Silla Observatory with the HARPS spectrograph to detect and characterize planets in the super-Earth and Neptune mass regime. The size of our star sample and the precision achieved with HARPS have allowed the detection of a sufficiently large number of low-mass planets to study the statistical properties of their orbital elements, the correlation of the host-star metallicity with the planet masses, as well as the occurrence rate of planetary systems around solar-type stars. A robust estimate of the frequency of systems shows that more than 50% of solar-type stars harbor at least one planet of any mass and with period up to 100 days. Different properties are observed for the population of planets less massive than about 30M-Earth compared to the population of gaseous giant planets. The mass distribution of Super-Earths and Neptune-mass planets (SEN) is strongly increasing between 30 and 15M-Earth. The SEN occurence rate does not exhibit a preference for metal rich stars. Most of the SEN planets belong to multi-planetary systems. The orbital eccentricities of the SEN planets seems limited to 0.45. At the opposite, the occurence rate of gaseous giant planets is growing with the logarithm of the period, and is strongly increasing with the host-star metallicity. About 14% of solar-type stars have a planetary companion more massive than 50M-Earth? on an orbit with a period shorter than 10 years. Orbital eccentricities of giant planets are observed up to 0.9 and beyond. The precision of HARPS-type spectrographs opens the possibility to detect planets in the habitable zone of solar-type stars. Identification of a significant number of super-Earths orbiting solar-type of the Sun vicinity is achieved by Doppler spectroscopy. 37 newly discovered planets are announced in the Appendix of this paper, among which 15 Super-Earths.

preprint2011arXiv

The HARPS search for southern extra-solar planets. XXX. Planetary systems around stars with solar-like magnetic cycles and short-term activity variation

We present the discovery of four new long-period planets within the HARPS high-precision sample: \object{HD137388}b ($M\sin{i}$ = 0.22 $M_J$), \object{HD204941}b ($M\sin{i}$ = 0.27 $M_J$), \object{HD7199}b ($M\sin{i}$ = 0.29 $M_J$), \object{HD7449}b ($M\sin{i}$ = 1.04 $M_J$). A long-period companion, probably a second planet, is also found orbiting HD7449. Planets around HD137388, HD204941, and HD7199 have rather low eccentricities (less than 0.4) relative to the 0.82 eccentricity of HD7449b. {All these planets were discovered even though their hosting stars have clear signs of activity. Solar-like magnetic cycles, characterized by long-term activity variations, can be seen for HD137388, HD204941 and HD7199, whereas the measurements of HD7449 reveal a short-term activity variation, most probably induced by magnetic features on the stellar surface. We confirm that magnetic cycles induce a long-term radial velocity variation and propose a method to reduce considerably the associated noise.} The procedure consists of fitting the activity index and applying the same solution to the radial velocities because a linear correlation between the activity index and the radial velocity is found. Tested on HD137388, HD204941, and HD7199, this correction reduces considerably the stellar noise induced by magnetic cycles and allows us to derive precisely the orbital parameters of planetary companions.

preprint2011arXiv

The HARPS search for southern extra-solar planets. XXXI. Magnetic activity cycles in solar-type stars: statistics and impact on precise radial velocities

Context. Searching for extrasolar planets through radial velocity measurements relies on the stability of stellar photospheres. Aims. We aim at characterizing the statistical properties of magnetic activity cycles, and studying their impact on spectroscopic measurements such as radial velocities, line bisectors and line shapes. Methods. We use data from the HARPS high-precision planet-search sample comprising 304 FGK stars followed over about 7 years. We obtain high-precision Ca II H&K chromospheric activity measurements and convert them to R'HK indices using an updated calibration taking into account stellar metallicity. We study R'HK variability as a function of time and search for possible correlations with radial velocities and line shape parameters. Results. We classify stars according to the magnitude and timescale of the Ca II H&K variability, and identify activity cycles whenever possible. We find that 39+/-8% of old solar-type stars in the solar neighborhood do not show any activity cycles (or only very weak ones), while 61+/-8% do have one. Non-cycling stars are almost only found among G dwarfs and at mean activity levels log R'HK < -4.95. Magnetic cycle amplitude generally decreases with decreasing activity level. A significant fraction of stars exhibit small variations in radial velocities and line shape parameters that are correlated with activity cycles. The sensitivity of radial velocities to magnetic cycles increases towards hotter stars, while late K dwarfs are almost insensitive. Conclusions. Activity cycles do induce long-period, low-amplitude radial velocity variations, at levels up to ~25 m/s. Caution is therefore mandatory when searching for long-period exoplanets. However, these effects can be corrected to high precision by detrending the radial velocity data using simultaneous measurements of Ca II H&K flux and line shape parameters. (Abridged)

preprint2011arXiv

The HARPS search for southern extrasolar planets. XXV. Four new planets in orbit around the moderatly active dwarfs HD 63765, HD 104067, HD 125595, and HIP 70849

We report the detection of four new extrasolar planets in orbit around the moderately active stars HD 63765, HD 104067, HIP 70849, and HD 125595 with the HARPS Echelle spectrograph mounted on the ESO 3.6-m telescope at La Silla. The first planet, HD 63765 b, has a minimum mass of 0.64MJup, a period of 358 days, and an eccentricity of 0.24. It orbits a G9 dwarf at 0.94 AU. The second planet, HD 104067 b, is a 3.6 Neptune-mass-planet with a 55.8-day-period. It orbits its parent K2 dwarf, in a circular orbit with a semi-major axis of a=0.26 AU. Radial velocity measurements present a \approx 500-day-oscillation that reveals significant magnetic cycles. The third planet is a 0.77 Neptune-mass-planet in circular orbit around the K4 dwarf, HD 12595, with a 9.67-day-period. Finally, HIP 7849 b is a long-period (5< P <75 years) and massive planet of m. sin i \approx 3.5-15 MJup that orbits a late K7 dwarf.

preprint2011arXiv

The retrograde orbit of the HAT-P-6b exoplanet

We observed with the SOPHIE spectrograph (OHP, France) the transit of the HAT-P-6b exoplanet across its host star. The resulting stellar radial velocities display the Rossiter-McLaughlin anomaly and reveal a retrograde orbit: the planetary orbital spin and the stellar rotational spin point towards approximately opposite directions. A fit to the anomaly measures a sky-projected angle lambda = 166 +/- 10 degrees between these two spin axes. All seven known retrograde planets are hot jupiters with masses M_p < 3 M_Jup. About two thirds of the planets in this mass range however are prograde and aligned (lambda ~ 0). By contrast, most of the more massive planets (M_p > 4 M_Jup) are prograde but misaligned. Different mechanisms may therefore be responsible for planetary obliquities above and below ~3.5 M_Jup.

preprint2011arXiv

The SOPHIE search for northern extrasolar planets. IV. Massive companions in the planet-brown dwarf boundary

The mass domain where massive extrasolar planets and brown dwarfs lay is still poorly understood. Indeed, not even a clear dividing line between massive planets and brown dwarfs has been established yet. This is partly due to the paucity of this kind of objects orbiting close to solar-type stars, the so-called brown dwarf desert, that hinders setting up a strong observational base to compare to models and theories of formation and evolution. We search to increase the current sample of massive sub-stellar objects with precise orbital parameters, and to constrain the true mass of detected sub-stellar candidates. The initial identification of sub-stellar candidates is done using precise radial velocity measurements obtained with the SOPHIE spectrograph at the 1.93-m telescope of the Haute-Provence Observatory. Subsequent characterisation of these candidates, with the principal aim of identifying stellar companions in low-inclination orbits, is done by means of different spectroscopic diagnostics, as the measurement of the bisector velocity span and the study of the correlation mask effect. With this objective, we also employed astrometric data from the Hipparcos mission and a novel method of simulating stellar cross-correlation functions. Seven new objects with minimum masses between ~ 10 Mjup and ~90 Mjup are detected. Out of these, two are identified as low-mass stars in low-inclination orbits, and two others have masses below the theoretical deuterium-burning limit, and are therefore planetary candidates. The remaining three are brown dwarf candidates; the current upper limits for their the masses do not allow us to conclude on their nature. Additionally, we have improved on the parameters of an already-known brown dwarf (HD137510b), confirmed by astrometry.

preprint2011arXiv

Transiting exoplanets from the CoRoT space mission - XIX. CoRoT-23b: a dense hot Jupiter on an eccentric orbit

We report the detection of CoRoT-23b, a hot Jupiter transiting in front of its host star with a period of 3.6314 \pm 0.0001 days. This planet was discovered thanks to photometric data secured with the CoRoT satellite, combined with spectroscopic radial velocity (RV) measurements. A photometric search for possible background eclipsing binaries conducted at CFHT and OGS concluded with a very low risk of false positives. The usual techniques of combining RV and transit data simultaneously were used to derive stellar and planetary parameters. The planet has a mass of Mp = 2.8 \pm 0.3 MJup, a radius of Rpl = 1.05 \pm 0.13 RJup, a density of \approx 3 g cm-3. RV data also clearly reveal a non zero eccentricity of e = 0.16 \pm 0.02. The planet orbits a mature G0 main sequence star of V =15.5 mag, with a mass M\star = 1.14 \pm 0.08 M\odot, a radius R\star = 1. 61 \pm 0.18 R\odot and quasi-solar abundances. The age of the system is evaluated to be 7 Gyr, not far from the transition to subgiant, in agreement with the rather large stellar radius. The two features of a significant eccentricity of the orbit and of a fairly high density are fairly uncommon for a hot Jupiter. The high density is, however, consistent with a model of contraction of a planet at this mass, given the age of the system. On the other hand, at such an age, circularization is expected to be completed. In fact, we show that for this planetary mass and orbital distance, any initial eccentricity should not totally vanish after 7 Gyr, as long as the tidal quality factor Qp is more than a few 105, a value that is the lower bound of the usually expected range. Even if Corot-23b features a density and an eccentricity that are atypical of a hot Jupiter, it is thus not an enigmatic object.

preprint2011arXiv

Transiting exoplanets from the CoRoT space mission XVII. The hot Jupiter CoRoT-17b: a very old planet

We report on the discovery of a hot Jupiter-type exoplanet, CoRoT-17b, detected by the CoRoT satellite. It has a mass of $2.43\pm0.30$\Mjup and a radius of $1.02\pm0.07$\Rjup, while its mean density is $2.82\pm0.38$ g/cm$^3$. CoRoT-17b is in a circular orbit with a period of $3.7681\pm0.0003$ days. The host star is an old ($10.7\pm1.0$ Gyr) main-sequence star, which makes it an intriguing object for planetary evolution studies. The planet's internal composition is not well constrained and can range from pure H/He to one that can contain $\sim$380 earth masses of heavier elements.

preprint2011arXiv

Transiting exoplanets from the CoRoT space mission: XIII. CoRoT-14b: an unusually dense very hot Jupiter

In this paper, the CoRoT Exoplanet Science Team announces its 14th discovery. Herein, we discuss the observations and analyses that allowed us to derive the parameters of this system: a hot Jupiter with a mass of $7.6 \pm 0.6$ Jupiter masses orbiting a solar-type star (F9V) with a period of only 1.5 d, less than 5 stellar radii from its parent star. It is unusual for such a massive planet to have such a small orbit: only one other known exoplanet with a higher mass orbits with a shorter period.

preprint2010arXiv

Astrometric search for a planet around VB 10

We observed VB 10 in August and September 2009 using the FORS2 camera of the VLT with the aim of measuring its astrometric motion and of probing the presence of the announced planet VB 10b. We used the published STEPS astrometric positions of VB 10 over a time-span of 9 years, which allowed us to compare the expected motion of VB 10 due to parallax and proper motion with the observed motion and to compute precise deviations. The achieved single-epoch precisions of our observations are about 0.1 mas and the data showed no significant residual trend, while the presence of the planet should have induced an apparent proper motion larger than 10 mas/yr. Subtraction of the predicted orbital motion from the observed data produces a large trend in position residuals of VB 10. We estimated the probability that this trend is caused by random noise. Taking all the uncertainties into account and using Monte-Carlo resampling of the data, we are able to reject the existence of VB 10b with the announced mass of 6.4 M_J with the false alarm probability of only 0.0005. A 3.2 M_J planet is also rejected with a false alarm probability of 0.023.

preprint2010arXiv

Consequences of spectrograph illumination for the accuracy of radial-velocimetry

For fiber-fed spectrographs with a stable external wavelength source, scrambling properties of optical fibers and, homogeneity and stability of the instrument illumination are important for the accuracy of radial-velocimetry. Optical cylindric fibers are known to have good azimuthal scrambling. In contrast, the radial one is not perfect. In order to improve the scrambling ability of the fiber and to stabilize the illumination, optical double scrambler are usually coupled to the fibers. Despite that, our experience on SOPHIE and HARPS has lead to identified remaining radial-velocity limitations due to the non-uniform illumination of the spectrograph. We conducted tests on SOPHIE with telescope vignetting, seeing variation and centering errors on the fiber entrance. We simulated the light path through the instrument in order to explain the radial velocity variation obtained with our tests. We then identified the illumination stability and uniformity has a critical point for the extremely high-precision radial velocity instruments (ESPRESSO@VLT, CODEX@E-ELT). Tests on square and octagonal section fibers are now under development and SOPHIE will be used as a bench test to validate these new feed optics.

preprint2010arXiv

Do stellar magnetic cycles influence the measurement of precise radial velocities?

The ever increasing level of precision achieved by present and future radial-velocity instruments is opening the way to discovering very low-mass, long-period planets (e.g. solar-system analogs). These systems will be detectable as low-amplitude signals in radial-velocity (RV). However, an important obstacle to their detection may be the existence of stellar magnetic cycles on similar timescales. Here we present the results of a long-term program to simultaneously measure radial-velocities and stellar-activity indicators (CaII, H_alpha, HeI) for a sample of stars with known activity cycles. Our results suggest that all these stellar activity indexes can be used to trace the stellar magnetic cycle in solar-type stars. Likewise, we find clear indications that different parameters of the HARPS cross-correlation function (BIS, FWHM, and contrast) are also sensitive to activity level variations. Finally, we show that, although in a few cases slight correlations or anti-correlations between radial-velocity and the activity level of the star exist, their origin is still not clear. We can, however, conclude that for our targets (early-K dwarfs) we do not find evidence of any radial-velocity variations induced by variations of the stellar magnetic cycle with amplitudes significantly above ~1 m/s.

preprint2010arXiv

Evaluating the stability of atmospheric lines with HARPS

Context: In the search for extrasolar systems by radial velocity technique, a precise wavelength calibration is necessary for high-precision measurements. The choice of the calibrator is a particularly important question in the infra-red domain, where the precision and exploits still fall behind the achievements of the optical. Aims: We investigate the long-term stability of atmospheric lines as a precise wavelength reference and analyze their sensitivity to different atmospheric and observing conditions. Methods: We use HARPS archive data on three bright stars, Tau Ceti, Mu Arae and Epsilon Eri, spanning 6 years and containing high-cadence measurements over several nights. We cross-correlate this data with an O2 mask and evaluate both radial velocity and bisector variations down to a photon noise of 1 m/s. Results: We find that the telluric lines in the three data-sets are stable down to 10 m/s (r.m.s.) over the 6 years. We also show that the radial velocity variations can be accounted for by simple atmospheric models, yielding a final precision of 1-2 m/s. Conclusions: The long-term stability of atmospheric lines was measured as being of 10 m/s over six years, in spite of atmospheric phenomena. Atmospheric lines can be used as a wavelength reference for short-time-scales programs, yielding a precision of 5 m/s "out-of-the box". A higher precision, down to 2 m/s can be reached if the atmospheric phenomena are corrected for by the simple atmospheric model described, making it a very competitive method even on long time-scales.

preprint2010arXiv

Evidence against the young hot-Jupiter around BD +20 1790

Context. The young active star BD +20 1790 is believed to host a substellar companion, revealed by radial-velocity measurements that detected the reflex motion induced on the parent star. Aims. A complete characterisation of the radial-velocity signal is necessary in order to assess its nature. Methods. We used CORALIE spectrograph to obtain precise (~10 m/s) velocity measurements on this active star, while characterizing the bisector span variations. Particular attention was given to correctly sample both the proposed planetary orbital period, of 7.8 days, and the stellar rotation period, of 2.4 days. Results. A smaller radial-velocity signal (with peak-to-peak variations <500 m/s) than had been reported previously was detected, with different amplitude on two different campaigns. A periodicity similar to the rotational period is found on the data, as well as a clear correlation between radial-velocities and bisector span. This evidence points towards a stellar origin of the radial-velocity variations of the star instead of a barycentric movement of the star, and repudiates the reported detection of a hot-Jupiter.

preprint2010arXiv

Observation of the full 12-hour-long transit of the exoplanet HD80606b. Warm-Spitzer photometry and SOPHIE spectroscopy

We present new observations of a transit of the 111-day-period exoplanet HD80606b. Using the Spitzer Space Telescope and its IRAC camera on the post-cryogenic mission, we performed a 19-hour-long photometric observation of HD80606 that covers the full transit of 13-14 January 2010. We complement this photometric data by new spectroscopic observations that we simultaneously performed with SOPHIE at Haute-Provence Observatory. This provides radial velocity measurements of the first half of the transit that was previously uncovered with spectroscopy. This new data set allows the parameters of this singular planetary system to be significantly refined. We obtained a planet-to-star radius ratio R_p/R_* = 0.1001 +/- 0.0006 that is slightly lower than the one measured from previous ground observations. We detected a feature in the Spitzer light curve that could be due to a stellar spot. We also found a transit timing about 20 minutes earlier than the ephemeris prediction; this could be caused by actual TTVs due to an additional body in the system or by underestimated systematic uncertainties. The sky-projected angle between the spin-axis of HD80606 and the normal to the planetary orbital plane is found to be lambda = 42 +/- 8 degrees thanks to the fit of the Rossiter-McLaughlin anomaly. This allows scenarios with aligned spin-orbit to be definitively rejected. Over the twenty planetary systems with measured spin-orbit angles, a few of them are misaligned; this is probably the signature of two different evolution scenarios for misaligned and aligned systems, depending if they experienced or not gravitational interaction with a third body. As in the case of HD80606b, most of the planetary systems including a massive planet are tilted; this could be the signature of a separate evolution scenario for massive planets in comparison with Jupiter-mass planets.

preprint2010arXiv

Stellar noise and planet detection. I. Oscillations, granulation and sun-like spots

Spectrographs like HARPS can now reach a sub-m/s precision in radial-velocity (RV) (Pepe & Lovis 2008). At this level of accuracy, we start to be confronted with stellar noise produced by 3 different physical phenomena: oscillations, granulation phenomena (granulation, meso- and super-granulation) and activity. On solar type stars, these 3 types of perturbation can induce m/s RV variation, but on different time scales: 3 to 15 minutes for oscillations, 15 minutes to 1.5 days for granulation phenomena and 10 to 50 days for activity. The high precision observational strategy used on HARPS, 1 measure per night of 15 minutes, on 10 consecutive days each month, is optimized, due to a long exposure time, to average out the noise coming from oscillations (Dumusque et al. 2010) but not to reduce the noise coming from granulation and activity. The smallest planets found with this strategy (Mayor et al. 2009) seems to be at the limit of the actual observational strategy and not at the limit of the instrumental precision. To be able to find Earth mass planets in the habitable zone of solar-type stars (200 days for a K0 dwarf), new observational strategies, averaging out simultaneously all type of stellar noise, are required.

preprint2010arXiv

Stellar noise and planet detection. II. Radial-velocity noise induced by magnetic cycles

For the 451 stars of the HARPS high precision program, we study correlations between the radial-velocity (RV) variation and other parameters of the Cross Correlated Function (CCF). After a careful target selection, we found a very good correlation between the slope of the RV-activity index (log(R'HK) correlation and the Teff for dwarf stars. This correlation allow us to correct RV from magnetic cycles given the activity index and the Teff.

preprint2010arXiv

The HARPS search for southern extra-solar planets XIX. Characterization and dynamics of the GJ876 planetary system

Precise radial-velocity measurements for data acquired with the HARPS spectrograph infer that three planets orbit the M4 dwarf star GJ876. In particular, we confirm the existence of planet "d", which orbits every 1.93785 days. We find that its orbit may have significant eccentricity (e=0.14), and deduce a more accurate estimate of its minimum mass of 6.3 Earth masses. Dynamical modeling of the HARPS measurements combined with literature velocities from the Keck Observatory strongly constrain the orbital inclinations of the "b" and "c" planets. We find that i_b = 48.9 degrees and i_c = 48.1 degrees, which infers the true planet masses of M_b = 2.64 Jupiter masses and M_c = 0.83 Jupiter masses, respectively. Radial velocities alone, in this favorable case, can therefore fully determine the orbital architecture of a multi-planet system, without the input from astrometry or transits. The orbits of the two giant planets are nearly coplanar, and their 2:1 mean motion resonance ensures stability over at least 5 Gyr. The libration amplitude is smaller than 2 degrees, suggesting that it was damped by some dissipative process during planet formation. The system has space for a stable fourth planet in a 4:1 mean motion resonance with planet "b", with a period around 15 days. The radial velocity measurements constrain the mass of this possible additional planet to be at most that of the Earth.

preprint2010arXiv

The HARPS search for southern extra-solar planets: XXVI: Seven new planetary systems

We are conducting a planet search survey with HARPS since seven years. The volume-limited stellar sample includes all F2 to M0 main-sequence stars within 57.5 pc, where extrasolar planetary signatures are systematically searched for with the radial-velocity technics. In this paper, we report the discovery of new substellar companions of seven main-sequence stars and one giant star, detected through multiple Doppler measurements with the instrument HARPS installed on the ESO 3.6m telescope, La Silla, Chile. These extrasolar planets orbit the stars HD 1690, HD 25171, HD 33473A, HD 89839, HD 113538, HD 167677 and HD 217686. The already-published giant planet around HD 72659 is also analysed here, and its elements are better determined by the addition of HARPS and Keck data. The other discoveries are giant planets in distant orbits, ranging from 0.3 to 29 MJup, in mass and between 0.7 and 10 years in orbital period. The low metallicity of most of these new planet-hosting stars reinforces the current trend for long-distance planets around metal-poor stars. Long-term radial-velocity surveys allow probing the outskirts of extrasolar planetary systems, although confidence in the solution may be low until more than one orbital period is fully covered by the observations. For many systems discussed in this paper, longer baselines are necessary to refine the radial-velocity fit and derive planetary parameters. The radial-velocity time series of stars BD -114672 and HIP 21934 are also analysed and their behaviour interpreted in terms of the activity cycle of the star, rather than long-period planetary companions.

preprint2010arXiv

The HARPS search for southern extra-solar planets. XXIV. Companions to HD 85390, HD 90156 and HD 103197: A Neptune analogue and two intermediate mass planets

We report the detection of three new extrasolar planets orbiting the solar type stars HD 85390, HD 90156 and HD 103197 with the HARPS spectrograph mounted on the ESO 3.6-m telescope at La Silla observatory. HD 85390 has a planetary companion with a projected intermediate mass (42.0 Earth masses) on a 788-day orbit (a=1.52 AU) with an eccentricity of 0.41, for which there is no analogue in the solar system. A drift in the data indicates the presence of another companion on a long period orbit, which is however not covered by our measurements. HD 90156 is orbited by a warm Neptune analogue with a minimum mass of 17.98 Earth masses (1.05 Neptune masses), a period of 49.8 days (a=0.25 AU) and an eccentricity of 0.31. HD 103197 has an intermediate mass planet on a circular orbit (P=47.8 d, Msini=31.2 Earth masses). We discuss the formation of planets of intermediate mass (about 30-100 Earth masses) which should be rare inside a few AU according to core accretion formation models.

preprint2010arXiv

The HARPS search for southern extra-solar planets. XXVII. Up to seven planets orbiting HD 10180: probing the architecture of low-mass planetary systems

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, opening a new window on planetary systems. Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. Methods. Progress has been possible in this field thanks in particular to the sub-m/s radial velocity precision achieved by HARPS. We present here new high-quality measurements from this instrument. Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging from 12 to 25 M_Earth, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is present at a longer period, probably due to a 65-M_Earth object. Moreover, another body with a minimum mass as low as 1.4 M_Earth may be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems exhibit a rather "packed" orbital architecture with little or no space left for additional planets. (Abridged)

preprint2010arXiv

The HARPS search for southern extrasolar planets XXI. Three new giant planets orbiting the metal-poor stars HD5388, HD181720, and HD190984

We present the discovery of three new giant planets around three metal-deficient stars: HD5388b (1.96M_Jup), HD181720b (0.37M_Jup), and HD190984b (3.1M_Jup). All the planets have moderately eccentric orbits (ranging from 0.26 to 0.57) and long orbital periods (from 777 to 4885 days). Two of the stars (HD181720 and HD190984) were part of a program searching for giant planets around a sample of ~100 moderately metal-poor stars, while HD5388 was part of the volume-limited sample of the HARPS GTO program. Our discoveries suggest that giant planets in long period orbits are not uncommon around moderately metal-poor stars.

preprint2010arXiv

The HARPS search for southern extrasolar planets XXV. Results from the metal-poor sample

Searching for extrasolar planets around stars of different metallicity may provide strong constraints to the models of planet formation and evolution. In this paper we present the overall results of a HARPS (a high-precision spectrograph mostly dedicated to deriving precise radial velocities) program to search for planets orbiting a sample of 104 metal-poor stars (selected [Fe/H] below -0.5). Radial velocity time series of each star are presented and searched for signals using several statistical diagnostics. Stars with detected signals are presented, including 3 attributed to the presence of previously announced giant planets orbiting the stars HD171028, HD181720, and HD190984. Several binary stars and at least one case of a coherent signal caused by activity-related phenomena are presented. One very promising new, possible giant planet orbiting the star HD107094 is discussed, and the results are analyzed in light of the metallicity-giant planet correlation. We conclude that the frequency of giant planets orbiting metal-poor stars may be higher than previously thought, probably reflecting the higher precision of the HARPS survey. In the metallicity domain of our sample, we also find evidence that the frequency of planets is a steeply rising function of the stellar metal content, as found for higher metallicity stars.

preprint2010arXiv

The SOPHIE search for northern extrasolar planets III. A Jupiter-mass companion around HD 109246

We report the detection of a Jupiter-mass planet discovered with the SOPHIE spectrograph mounted on the 1.93-m telescope at the Haute-Provence Observatory. The new planet orbits HD109246, a G0V star slightly more metallic than the Sun. HD109246b has a minimum mass of 0.77 MJup, an orbital period of 68 days, and an eccentricity of 0.12. It is placed in a sparsely populated region of the period distribution of extrasolar planets. We also present a correction method for the so-called seeing effect that affects the SOPHIE radial velocities. We complement this discovery announcement with a description of some calibrations that are implemented in the SOPHIE automatic reduction pipeline. These calibrations allow the derivation of the photon-noise radial velocity uncertainty and some useful stellar properties (vsini, [Fe/H], logR'HK) directly from the SOPHIE data.

preprint2010arXiv

The SOPHIE search for northern extrasolar planets. II. A multi-planet system around HD9446

We report the discovery of a planetary system around HD9446, performed from radial velocity measurements secured with the spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory during more than two years. At least two planets orbit this G5V, active star: HD9446b has a minimum mass of 0.7 M_Jup and a slightly eccentric orbit with a period of 30 days, whereas HD9446c has a minimum mass of 1.8 M_Jup and a circular orbit with a period of 193 days. As for most of the known multi-planet systems, the HD9446-system presents a hierarchical disposition, with a massive outer planet and a lighter inner planet.

preprint2010arXiv

The Spitzer search for the transits of HARPS low-mass planets - I. No transit for the super-Earth HD 40307b

We have used Spitzer and its IRAC camera to search for the transit of the super-Earth HD 40307b. The transiting nature of the planet could not be firmly discarded from our first photometric monitoring of a transit window because of the uncertainty coming from the modeling of the photometric baseline. To obtain a firm result, two more transit windows were observed and a global Bayesian analysis of the three IRAC time series and the HARPS radial velocities was performed. Unfortunately, any transit of the planet during the observed phase window is firmly discarded, while the probability that the planet transits but that the eclipse was missed by our observations is nearly negligible (0.26%).

preprint2010arXiv

Transiting exoplanets from the CoRoT space mission. XI. CoRoT-8b: a hot and dense sub-Saturn around a K1 dwarf

We report the discovery of CoRoT-8b, a dense small Saturn-class exoplanet that orbits a K1 dwarf in 6.2 days, and we derive its orbital parameters, mass, and radius. We analyzed two complementary data sets: the photometric transit curve of CoRoT-8b as measured by CoRoT and the radial velocity curve of CoRoT-8 as measured by the HARPS spectrometer. We find that CoRoT-8b is on a circular orbit with a semi-major axis of 0.063 +/- 0.001 AU. It has a radius of 0.57 +/- 0.02 RJ, a mass of 0.22 +/- 0.03 MJ, and therefore a mean density 1.6 +/- 0.1 g/cm^3. With 67 % of the size of Saturn and 72 % of its mass, CoRoT-8b has a density comparable to that of Neptune (1.76 g/cm^3). We estimate its content in heavy elements to be 47-63 Earth masses, and the mass of its hydrogen-helium envelope to be 7-23 Earth masses. At 0.063 AU, the thermal loss of hydrogen of CoRoT-8b should be no more than about 0.1 % over an assumed integrated lifetime of 3~Ga.

preprint2010arXiv

Transiting exoplanets from the CoRoT space mission. XIII. CoRoT-13b: a dense hot Jupiter in transit around a star with solar metallicity and super-solar lithium content

We announce the discovery of the transiting planet CoRoT-13b. Ground based follow-up in CFHT and IAC80 confirmed CoRoT's observations. The mass of the planet was measured with the HARPS spectrograph and the properties of the host star were obtained analyzing HIRES spectra from the Keck telescope. It is a hot Jupiter-like planet with an orbital period of 4.04 days, 1.3 Jupiter masses, 0.9 Jupiter radii, and a density of 2.34 g cm-3. It orbits a G0V star with Teff=5945K, M*=1.09 Msun, R*=1.01 Rsun, solar metallicity, a lithium content of +1.45 dex, and an estimated age between 0.12 and 3.15 Gyr. The lithium abundance of the star is consistent with its effective temperature, activity level, and age range derived from the stellar analysis. The density of the planet is extreme for its mass. It implies the existence of an amount of heavy elements with a mass between about 140 and 300 Mearth.

preprint2009arXiv

Extrasolar planets and brown dwarfs around A-F type stars - VII. Theta Cygni radial velocity variations: planets or stellar phenomenon?

(abridged) In the frame of the search for extrasolar planets and brown dwarfs around early-type main-sequence stars, we present the results obtained on the early F-type star Theta Cygni. Elodie and Sophie at OHP were used to obtain the spectra. Our dedicated radial-velocity measurement method was used to monitor the star's radial velocities over five years. We also use complementary, high angular resolution and high-contrast images taken with PUEO at CFHT. We show that Theta Cygni radial velocities are quasi-periodically variable, with a ~150-day period. These variations are not due to the ~0.35-Msun stellar companion that we detected in imaging at more than 46 AU from the star. The absence of correlation between the bisector velocity span variations and the radial velocity variations for this 7 km/s vsini star, as well as other criteria indicate that the observed radial velocity variations are not due to stellar spots. The observed amplitude of the bisector velocity span variations also seems to rule out stellar pulsations. However, we observe a peak in the bisector velocity span periodogram at the same period as the one found in the radial velocity periodogram, which indicates a probable link between these radial velocity variations and the low amplitude lineshape variations which are of stellar origin. Long-period variations are not expected from this type of star to our knowledge. If a stellar origin (hence of new type) was to be confirmed for these long-period radial velocity variations, this would have several consequences on the search for planets around main-sequence stars, both in terms of observational strategy and data analysis. An alternative explanation for these variable radial velocities is the presence of at least one planet of a few Jupiter masses orbiting at less than 1 AU. (abridged)

preprint2009arXiv

Radial Velocities with CRIRES: Pushing precision down to 5-10 m/s

With the advent of high-resolution infrared spectrographs, Radial Velocity (RV) searches enter into a new domain. As of today, the most important technical question to address is which wavelength reference is the most suitable for high-precision RV measurements. In this work we explore the usage of atmospheric absorption features. We make use of CRIRES data on two programs and three different targets. We re-analyze the data of the TW Hya campaign, reaching a dispersion of about 6 m/s on the RV standard in a time scale of roughly 1 week. We confirm the presence of a low-amplitude RV signal on TW Hya itself, roughly 3 times smaller than the one reported at visible wavelengths. We present RV measurements of Gl 86 as well, showing that our approach is capable of detecting the signal induced by a planet and correctly quantifying it. Our data show that CRIRES is capable of reaching a RV precision of less than 10 m/s in a time-scale of one week. The limitations of this particular approach are discussed, and the limiting factors on RV precision in the IR in a general way. The implications of this work on the design of future dedicated IR spectrographs are addressed as well.

preprint2009arXiv

The HARPS search for southern extra-solar planets XVIII. An Earth-mass planet in the GJ 581 planetary system

The GJ 581 planetary system was already known to harbour three planets, including two super-Earths planets which straddle its habitable zone. We report here the detection of an additional planet -- GJ 581e -- with a minimum mass of 1.9 M_earth. With a period of 3.15 days, it is the innermost planet of the system and has a ~5% transit probability. We also correct our previous confusion of the orbital period of GJ 581d (the outermost planet) with a one-year alias, thanks to an extended time span and many more measurements. The revised period is 66.8 days, and locates the semi-major axis inside the habitable zone of the low mass star. The dynamical stability of the 4-planet system imposes an upper bound on the orbital plane inclination. The planets cannot be more massive than approximately 1.6 times their minimum mass.

preprint2009arXiv

The HARPS search for southern extra-solar planets. XX. Planets around the active star BD-08:2823

We report the detection of a planetary system around BD-08:2823, that includes at least one Uranus-mass planet and one Saturn-mass planet. This discovery serendipitously originates from a search for planetary transits in the Hipparcos photometry database. This program preferentially selected active stars and did not allow the detection of new transiting planets. It allowed however the identification of the K3V star BD-08:2823 as a target harboring a multiplanet system, that we secured and characterized thanks to an intensive monitoring with the HARPS spectrograph at the 3.6-m ESO telescope in La Silla. The stellar activity level of BD-08:2823 complicates the analysis but does not prohibit the detection of two planets around this star. BD-08:2823b has a minimum mass of 14.4+/-2.1 M_Earth and an orbital period of 5.60 days, whereas BD-08:2823c has a minimum mass of 0.33+/-0.03 M_Jup and an orbital period of 237.6 days. This new system strengthens the fact that low-mass planets are preferentially found in multiplanetary systems, but not around high-metallicity stars as this is the case for massive planets. It also supports the belief that active stars should not be neglected in exoplanet searches, even when searching for low-mass planets.

preprint2007arXiv

An extrasolar planetary system with three Neptune-mass planets

Over the past two years, the search for low-mass extrasolar planets has led to the detection of seven so-called 'hot Neptunes' or 'super-Earths' around Sun-like stars. These planets have masses 5-20 times larger than the Earth and are mainly found on close-in orbits with periods of 2-15 days. Here we report a system of three Neptune-mass planets with periods of 8.67, 31.6 and 197 days, orbiting the nearby star HD 69830. This star was already known to show an infrared excess possibly caused by an asteroid belt within 1 AU (the Sun-Earth distance). Simulations show that the system is in a dynamically stable configuration. Theoretical calculations favour a mainly rocky composition for both inner planets, while the outer planet probably has a significant gaseous envelope surrounding its rocky/icy core; the outer planet orbits within the habitable zone of this star.

preprint2007arXiv

The CODEX-ESPRESSO experiment: cosmic dynamics, fundamental physics, planets and much more..

CODEX, a high resolution, super-stable spectrograph to be fed by the E-ELT, the most powerful telescope ever conceived, will for the first time provide the possibility of directly measuring the change of the expansion rate of the Universe with time and much more, from the variability of fundamental constants to the search for other earths. A study for the implementation at the VLT of a precursor of CODEX, dubbed ESPRESSO, is presently carried out by a collaboration including ESO, IAC, INAF, IoA Cambridge and Observatoire de Geneve. The present talk is focused on the cosmological aspects of the experiment.