Source author record

P. Figueira

P. Figueira 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

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

91 published item(s)

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

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

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

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

preprint2022arXiv

Stellar signal components seen in HARPS and HARPS-N solar radial velocities

Context. Radial velocity (RV) measurements induced by the presence of planets around late-type stars are contaminated by stellar signals that are of the order of a few meters per second in amplitude, even for the quietest stars. Those signals are induced by acoustic oscillations, convective granulation patterns, active regions co-rotating with the stellar surface, and magnetic activity cycles. Aims. This study investigates the properties of all coherent stellar signals seen on the Sun on timescales up to its sidereal rotational period. By combining HARPS and HARPS-N solar data spanning several years, we are able to clearly resolve signals on timescales from minutes to several months. Methods. We use a Markov Chain Monte Carlo (MCMC) mixture model to determine the quality of the solar data based on the expected airmass-magnitude extinction law. We then fit the velocity power spectrum of the cleaned and heliocentric RVs with all known variability sources, to recreate the RV contribution of each component. Results. After rejecting variations caused by poor weather conditions, we are able to improve the average intra-day root mean square (RMS) value by a factor of ~1.8. On sub-rotational timescales, we are able to fully recreate the observed RMS of the RV variations. In order to also include rotational components and their strong alias peaks introduced by nightly sampling gaps, the alias powers are accounted for by being redistributed to the central frequencies of the rotational harmonics. Conclusions. In order to enable a better understanding and mitigation of stellar activity sources, their respective impact on the total RV must be well-measured and characterized. We are able to recreate RV components up to rotational timescales, which can be further used to analyse the impact of each individual source of stellar signals on the detectability of exoplanets.

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

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

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

preprint2022arXiv

TOI-1759 b: a transiting sub-Neptune around a low mass star characterized with SPIRou and TESS

We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial velocities were extracted from the SPIRou spectra using both the CCF and the LBL methods, optimizing the velocity measurements in the near infrared domain. We analyzed the broadband SED of the star and the high-resolution SPIRou spectra to constrain the stellar parameters and thus improve the accuracy of the derived planet parameters. A LSD analysis of the SPIRou Stokes $V$ polarized spectra detects Zeeman signatures in TOI-1759. We model the rotational modulation of the magnetic stellar activity using a GP regression with a quasi-periodic covariance function, and find a rotation period of $35.65^{+0.17}_{-0.15}$ d. We reconstruct the large-scale surface magnetic field of the star using ZDI, which gives a predominantly poloidal field with a mean strength of $18\pm4$ G. Finally, we perform a joint Bayesian MCMC analysis of the TESS photometry and SPIRou RVs to optimally constrain the system parameters. At $0.1176\pm0.0013$ au from the star, the planet receives $6.4$ times the bolometric flux incident on Earth, and its equilibrium temperature is estimated at $433\pm14$ K. TOI-1759 b is a likely gas-dominated sub-Neptune with an expected high rate of photoevaporation. Therefore, it is an interesting target to search for neutral hydrogen escape, which may provide important constraints on the planetary formation mechanisms responsible for the observed sub-Neptune radius desert.

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

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth'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 HD 137496 system: A dense, hot super-Mercury and a cold Jupiter

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

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

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

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

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

preprint2020arXiv

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

SPIRou: nIR velocimetry & spectropolarimetry at the CFHT

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

preprint2020arXiv

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 correlation between photometric variability and radial velocity jitter, based on TESS and HARPS observations

The current and upcoming high precision photometric surveys such as TESS, CHEOPS, and PLATO will provide the community with thousands of new exoplanet candidates. As a consequence, the presence of such a correlation is crucial in selecting the targets with the lowest RV jitter for efficient RV follow-up of exoplanetary candidates. Studies of this type are also crucial to design optimized observational strategies to mitigate RV jitter when searching for Earth-mass exoplanets. Our goal is to assess the correlation between high-precision photometric variability measurements and high-precision RV jitter over different time scales. We analyze 171 G, K, and M stars with available TESS high precision photometric time-series and HARPS precise RVs. We derived the stellar parameters for the stars in our sample and measured the RV jitter and photometric variability. We also estimated chromospheric Ca II H $\&$ K activity indicator $log(R' _{HK})$, $\textit{v sin i}$, and the stellar rotational period. Finally, we evaluate how different stellar parameters and a RV sampling subset can have an impact on the potential correlations. We find a varying correlation between the photometric variability and RV jitter as function of time intervals between the TESS photometric observation and HARPS RV. As the time intervals of the observations considered for the analysis increases, the correlation value and significance becomes smaller and weaker, to the point that it becomes negligible. We also find that for stars with a photometric variability above 6.5 ppt the correlation is significantly stronger. We show that such a result can be due to the transition between the spot-dominated and the faculae-dominated regime. We quantified the correlations and updated the relationship between chromospheric Ca II H $\&$ K activity indicator $log(R' _{HK})$ and RV jitter.

preprint2020arXiv

Three planets transiting the evolved star EPIC 249893012: a hot 8.8-M$_\oplus$ super-Earth and two warm 14.7 and 10.2-M$_\oplus$ sub-Neptunes

We report the discovery of a new planetary system with three transiting planets, one super-Earth and two sub-Neptunes, that orbit EPIC\,249893012, a G8\,IV-V evolved star ($M_\star$\,=\,1.05\,$\pm$\,0.05\,$M_\odot$, $R_\star$\,=\,1.71\,$\pm$\,0.04\,$R_\odot$, $T_\mathrm{eff}$\,=5430\,$\pm$\,85\,K). The star is just leaving the main sequence. We combined \ktwo \ photometry with IRCS adaptive-optics imaging and HARPS, HARPS-N, and CARMENES high-precision radial velocity measurements to confirm the planetary system, determine the stellar parameters, and measure radii, masses, and densities of the three planets. With an orbital period of $3.5949^{+0.0007}_{-0.0007}$ days, a mass of $8.75^{+1.09}_{-1.08}\ M_{\oplus}$ , and a radius of $1.95^{+0.09}_{-0.08}\ R_{\oplus}$, the inner planet b is compatible with nickel-iron core and a silicate mantle ($ρ_b= 6.39^{+1.19}_{-1.04}$ g cm$^{-3}$). Planets c and d with orbital periods of $15.624^{+0.001}_{-0.001}$ and $35.747^{+0.005}_{-0.005}$ days, respectively, have masses and radii of $14.67^{+1,84}_{-1.89}\ M_{\oplus}$ and $3.67^{+0.17}_{-0.14}\ R_{\oplus}$ and $10.18^{+2.46}_{-2.42}\ M_{\oplus}$ and $3.94^{+0.13}_{-0.12}\ R_{\oplus}$, respectively, yielding a mean density of $1.62^{+0.30}_{-0.29}$ and $0.91^{+0.25}_{-0.23}$ g cm$^{-3}$, respectively. The radius of planet b lies in the transition region between rocky and gaseous planets, but its density is consistent with a rocky composition. Its semimajor axis and the corresponding photoevaporation levels to which the planet has been exposed might explain its measured density today. In contrast, the densities and semimajor axes of planets c and d suggest a very thick atmosphere. The singularity of this system, which orbits a slightly evolved star that is just leaving the main sequence, makes it a good candidate for a deeper study from a dynamical point of view.

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

Decoding the radial velocity variations of HD41248 with ESPRESSO

Twenty-four years after the discoveries of the first exoplanets, the radial-velocity (RV) method is still one of the most productive techniques to detect and confirm exoplanets. But stellar magnetic activity can induce RV variations large enough to make it difficult to disentangle planet signals from the stellar noise. In this context, HD41248 is an interesting planet-host candidate, with RV observations plagued by activity-induced signals. We report on ESPRESSO observations of HD41248 and analyse them together with previous observations from HARPS with the goal of evaluating the presence of orbiting planets. Using different noise models within a general Bayesian framework designed for planet detection in RV data, we test the significance of the various signals present in the HD41248 dataset. We use Gaussian processes as well as a first-order moving average component to try to correct for activity-induced signals. At the same time, we analyse photometry from the TESS mission, searching for transits and rotational modulation in the light curve. The number of significantly detected Keplerian signals depends on the noise model employed, which can range from 0 with the Gaussian process model to 3 with a white noise model. We find that the Gaussian process alone can explain the RV data while allowing for the stellar rotation period and active region evolution timescale to be constrained. The rotation period estimated from the RVs agrees with the value determined from the TESS light curve. Based on the data that is currently available, we conclude that the RV variations of HD41248 can be explained by stellar activity (using the Gaussian process model) in line with the evidence from activity indicators and the TESS photometry.

preprint2016arXiv

$ζ^2$ Ret, its debris disk, and its lonely stellar companion $ζ^1$ Ret. Different $T_{\mathrm{c}}$ trends for different spectra

Several studies have reported a correlation between the chemical abundances of stars and condensation temperature (known as Tc trend). Very recently, a strong Tc trend was reported for the $ζ$ Reticuli binary system, which consists of two solar analogs. The observed trend in $ζ^2$ Ret relative to its companion was explained by the presence of a debris disk around $ζ^2$ Ret. Our goal is to re-evaluate the presence and variability of the Tc trend in the $ζ$ Reticuli system and to understand the impact of the presence of the debris disk on a star. We used very high-quality spectra of the two stars retrieved from the HARPS archive to derive very precise stellar parameters and chemical abundances. We derived the stellar parameters with the classical (nondifferential) method, while we applied a differential line-by-line analysis to achieve the highest possible precision in abundances, which are fundamental to explore for very tiny differences in the abundances between the stars. We confirm that the abundance difference between $ζ^2$ Ret and $ζ^1$ Ret shows a significant ($\sim$ 2 $σ$) correlation with Tc. However, we also find that the Tc trends depend on the individual spectrum used (even if always of very high quality). In particular, we find significant but varying differences in the abundances of the same star from different individual high-quality spectra. Our results for the $ζ$ Reticuli system show, for example, that nonphysical factors, such as the quality of spectra employed and errors that are not accounted for, can be at the root of the Tc trends for the case of individual spectra.

preprint2016arXiv

A pragmatic Bayesian perspective on correlation analysis: The exoplanetary gravity - stellar activity case

We apply the Bayesian framework to assess the presence of a correlation between two quantities. To do so, we estimate the probability distribution of the parameter of interest, $ρ$, characterizing the strength of the correlation. We provide an implementation of these ideas and concepts using python programming language and the pyMC module in a very short ($\sim$130 lines of code, heavily commented) and user-friendly program. We used this tool to assess the presence and properties of the correlation between planetary surface gravity and stellar activity level as measured by the log($R'_{\mathrm{HK}}$) indicator. The results of the Bayesian analysis are qualitatively similar to those obtained via p-value analysis, and support the presence of a correlation in the data. The results are more robust in their derivation and more informative, revealing interesting features such as asymmetric posterior distributions or markedly different credible intervals, and allowing for a deeper exploration. We encourage the reader interested in this kind of problem to apply our code to his/her own scientific problems. The full understanding of what the Bayesian framework is can only be gained through the insight that comes by handling priors, assessing the convergence of Monte Carlo runs, and a multitude of other practical problems. We hope to contribute so that Bayesian analysis becomes a tool in the toolkit of researchers, and they understand by experience its advantages and limitations.

preprint2016arXiv

Abundance trend with condensation temperature for stars with different Galactic birth places

During the past decade, several studies reported a correlation between chemical abundances of stars and condensation temperature (also known as Tc trend). However, the real astrophysical nature of this correlation is still debated. The main goal of this work is to explore the possible dependence of the Tc trend on stellar Galactocentric distances, Rmean. We used high-quality spectra of about 40 stars observed with the HARPS and UVES spectrographs to derive precise stellar parameters, chemical abundances, and stellar ages. A differential line-by-line analysis was applied to achieve the highest possible precision in the chemical abundances. We confirm previous results that [X/Fe] abundance ratios depend on stellar age and that for a given age, some elements also show a dependence on Rmean. When using the whole sample of stars, we observe a weak hint that the Tc trend depends on Rmean. The observed dependence is very complex and disappears when only stars with similar ages are considered. To conclude on the possible dependence of the Tc trend on the formation place of stars, a larger sample of stars with very similar atmospheric parameters and stellar ages observed at different Galactocentric distances is needed

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

Can stellar activity make a planet seem misaligned?

Several studies have shown that the occultation of stellar active regions by the transiting planet can generate anomalies in the high-precision transit light curves, and these anomalies may lead to an inaccurate estimate of the planetary parameters (e.g., the planet radius). Since the physics and geometry behind the transit light curve and the Rossiter- McLaughlin effect (spectroscopic transit) are the same, the Rossiter-McLaughlin observations are expected to be affected by the occultation of stellar active regions in a similar way. In this paper we perform a fundamental test on the spin-orbit angles as derived by Rossiter-McLaughlin measurements, and we examine the impact of the occultation of stellar active regions by the transiting planet on the spin-orbit angle estimations. Our results show that the inaccurate estimation on the spin-orbit angle due to stellar activity can be quite significant (up to 30 degrees), particularly for the edge-on, aligned, and small transiting planets. Therefore, our results suggest that the aligned transiting planets are the ones that can be easily misinterpreted as misaligned owing to the stellar activity. In other words, the biases introduced by ignoring stellar activity are unlikely to be the culprit for the highly misaligned systems.

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

Modelling the Rossiter-McLaughlin Effect: Impact of the Convective Centre-to-Limb Variations in the Stellar Photosphere

Observations of the Rossiter-McLaughlin (RM) effect provide information on star-planet alignments, which can inform planetary migration and evolution theories. Here, we go beyond the classical RM modelling and explore the impact of a convective blueshift that varies across the stellar disc and non-Gaussian stellar photospheric profiles. We simulated an aligned hot Jupiter with a 4 d orbit about a Sun-like star and injected centre-to-limb velocity (and profile shape) variations based on radiative 3D magnetohydrodynamic simulations of solar surface convection. The residuals between our modelling and classical RM modelling were dependent on the intrinsic profile width and v sin i; the amplitude of the residuals increased with increasing v sin i, and with decreasing intrinsic profile width. For slowly rotating stars the centre-to-limb convective variation dominated the residuals (with amplitudes of 10s of cm/s to ~1 m/s); however, for faster rotating stars the dominant residual signature was due a non-Gaussian intrinsic profile (with amplitudes from 0.5-9 m/s). When the impact factor was 0, neglecting to account for the convective centre-to-limb variation led to an uncertainty in the obliquity of ~10-20 degrees, even though the true v sin i was known. Additionally, neglecting to properly model an asymmetric intrinsic profile had a greater impact for more rapidly rotating stars (e.g. v sin i = 6 km/s), and caused systematic errors on the order of ~20 degrees in the measured obliquities. Hence, neglecting the impact of stellar surface convection may bias star-planet alignment measurements and consequently also theories on planetary migration and evolution.

preprint2016arXiv

Spectroscopic characterisation of microlensing events Towards a new interpretation of OGLE-2011-BLG-0417

The microlensing event OGLE-2011-BLG-0417 is an exceptionally bright lens binary that was predicted to present radial velocity variation at the level of several km/s. Pioneer radial velocity follow-up observations with the UVES spectrograph at the ESO - VLT of this system clearly ruled out the large radial velocity variation, leaving a discrepancy between the observation and the prediction. In this paper, we further characterise the microlensing system by analysing its spectral energy distribution (SED) derived using the UVES spectrum and new observations with the ARCoIRIS (CTIO) near-infrared spectrograph and the Keck adaptive optics instrument NIRC2 in the J, H, and Ks bands. We determine the mass and distance of the stars independently from the microlensing modelling. We find that the SED is compatible with a giant star in the Galactic bulge and a foreground star with a mass of 0.94+/-0.09Msun at a distance of 1.07+/-0.24kpc. We find that this foreground star is likely the lens. Its parameters are not compatible with the ones previously reported in the literature (0.52+/-0.04Msun at 0.95+/-0.06kpc), based on the microlensing light curve. A thoughtful re-analysis of the microlensing event is mandatory to fully understand the reason of this new discrepancy. More importantly, this paper demonstrates that spectroscopic follow-up observations of microlensing events are possible and provide independent constraints on the parameters of the lens and source stars, hence breaking some degeneracies in the analysis. UV-to-NIR low-resolution spectrographs like X-SHOOTER (ESO - VLT) could substantially contribute to this follow-up efforts, with magnitude limits above all microlensing events detected so far.

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

Uncovering the planets and stellar activity of CoRoT-7 using only radial velocities

Stellar activity can induce signals in the radial velocities of stars, complicating the detection of orbiting low-mass planets. We present a method to determine the number of planetary signals present in radial-velocity datasets of active stars, using only radial-velocity observations. Instead of considering separate fits with different number of planets, we use a birth-death Markov chain Monte Carlo algorithm to infer the posterior distribution for the number of planets in a single run. In a natural way, the marginal distributions for the orbital parameters of all planets are also inferred. This method is applied to HARPS data of CoRoT-7. We confidently recover both CoRoT-7b and CoRoT-7c although the data show evidence for additional signals.

preprint2015arXiv

Constraining planet structure from stellar chemistry: the cases of CoRoT-7, Kepler-10, and Kepler-93

We explore the possibility that the stellar relative abundances of different species can be used to constrain the bulk abundances of known transiting rocky planets. We use high resolution spectra to derive stellar parameters and chemical abundances for Fe, Si, Mg, O, and C in three stars hosting low mass, rocky planets: CoRoT-7, Kepler-10, and Kepler-93. These planets follow the same line along the mass-radius diagram, pointing toward a similar composition. The derived abundance ratios are compared with the solar values. With a simple stoichiometric model, we estimate the iron mass fraction in each planet, assuming stellar composition. We show that in all cases, the iron mass fraction inferred from the mass-radius relationship seems to be in good agreement with the iron abundance derived from the host star's photospheric composition. The results suggest that stellar abundances can be used to add constraints on the composition of orbiting rocky planets.

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

Eclipsing binaries and fast rotators in the Kepler sample. Characterization via radial velocity analysis from Calar Alto

The Kepler mission has provided high-accurate photometric data in a long time span for more than two hundred thousands stars, looking for planetary transits. Among the detected candidates, the planetary nature of around 15% has been established or validated by different techniques. But additional data is needed to characterize the rest of the candidates and reject other possible configurations. We started a follow-up program to validate, confirm, and characterize some of the planet candidates. In this paper we present the radial velocity analysis (RV) of those presenting large variations, compatible with being eclipsing binaries. We also study those showing large rotational velocities, which prevents us from obtaining the necessary precision to detect planetary-like objects. We present new RV results for 13 Kepler objects of interest (KOIs) obtained with the CAFE spectrograph at the Calar Alto Observatory, and analyze their high-spatial resolution images and the Kepler light curves of some interesting cases. We have found five spectroscopic and eclipsing binaries. Among them, the case of KOI-3853 is of particular interest. This system is a new example of the so-called heartbeat stars, showing dynamic tidal distortions in the Kepler light curve. We have also detected duration and depth variations of the eclipse. We suggest possible scenarios to explain such effect, including the presence of a third substellar body possibly detected in our RV analysis. We also provide upper mass limits to the transiting companions of other six KOIs with large rotational velocities. This property prevents the RV method to obtain the necessary precision to detect planetary-like masses. Finally, we analyze the large RV variations of other two KOIs, incompatible with the presence of planetary-mass objects. These objects are likely to be stellar binaries but a longer timespan is still needed.

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

From stellar to planetary composition: Galactic chemical evolution of Mg/Si mineralogical ratio

The main goal of this work is to study element ratios that are important for the formation of planets of different masses. We study potential correlations between the existence of planetary companions and the relative elemental abundances of their host stars. We use a large sample of FGK-type dwarf stars for which precise Mg, Si, and Fe abundances have been derived using HARPS high-resolution and high-quality data. A first analysis of the data suggests that low-mass planet host stars show higher [Mg/Si] ratios, while giant planet hosts present [Mg/Si] that is lower than field stars. However, we found that the [Mg/Si] ratio significantly depends on metallicity through Galactic chemical evolution. After removing the Galactic evolution trend only the difference in the [Mg/Si] elemental ratio between low-mass planet hosts and non-hosts was present in a significant way. These results suggests that low-mass planets are more prevalent around stars with high [Mg/Si]. Our results demonstrate the importance of Galactic chemical evolution and indicate that it may play an important role in the planetary internal structure and composition.

preprint2015arXiv

Further constraints on the optical transmission spectrum of HAT-P-1b

We report on novel observations of HAT-P-1 aimed at constraining the optical transmission spectrum of the atmosphere of its transiting Hot-Jupiter exoplanet. Ground-based differential spectrophotometry was performed over two transit windows using the DOLORES spectrograph at the Telescopio Nazionale Galileo (TNG). Our measurements imply an average planet to star radius ratio equal to $\rm R_p/R_{\star}$=(0.1159$\pm$0.0005). This result is consistent with the value obtained from recent near infrared measurements of this object but differs from previously reported optical measurements being lower by around 4.4 exoplanet scale heights. Analyzing the data over 5 different spectral bins 600Å$\,$ wide we observed a single peaked spectrum (3.7 $\rmσ$ level) with a blue cut-off corresponding to the blue edge of the broad absorption wing of sodium and an increased absorption in the region in between 6180-7400Å. We also infer that the width of the broad absorption wings due to alkali metals is likely narrower than the one implied by solar abundance clear atmospheric models. We interpret the result as evidence that HAT-P-1b has a partially clear atmosphere at optical wavelengths with a more modest contribution from an optical absorber than previously reported.

preprint2015arXiv

Identifying the best iron-peak and $α$-capture elements for chemical tagging: The impact of the number of lines on measured scatter

The main goal of this work is to explore which elements carry the most information about the birth origin of stars and as such that are best suited for chemical tagging. We explored different techniques to minimize the effect of outlier value lines in the abundances by using Ni abundances derived for 1111 FGK type stars.We evaluated how the limited number of spectral lines can affect the final chemical abundance. Then we were able to make an efficient even footing comparison of the [X/Fe] scatter between the elements that have different number of observable spectral lines in the studied spectra. We found that the most efficient way of calculating the average abundance of elements when several spectral lines are available is to use a weighted mean (WM) where as a weight we considered the distance from the median abundance. This method can be effectively used without removing suspected outlier lines.We showed that when the same number of lines is used to determine chemical abundances, the [X/Fe] star-to-star scatter for iron group and $α$-capture elements is almost the same. On top of this, but at a lower level the largest scatter was observed for Al and the smallest for Cr and Ni. We recommend caution when comparing [X/Fe] scatters among elements that have a different number of spectral lines available. A meaningful comparison is necessary to identify elements that show the largest intrinsic scatter and can be thus used for chemical tagging.

preprint2015arXiv

Kepler-447b: a hot-Jupiter with an extremely grazing transit

We present the radial velocity confirmation of the extrasolar planet Kepler-447b, initially detected as a candidate by the Kepler mission. In this work, we analyze its transit signal and the radial velocity data obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). By simultaneously modeling both datasets, we obtain the orbital and physical properties of the system. According to our results, Kepler-447b is a Jupiter-mass planet ($M_p=1.37^{+0.48}_{-0.46}~M_{\rm Jup}$), with an estimated radius of $R_p=1.65^{+0.59}_{-0.56}~R_{\rm Jup}$ (uncertainties provided in this work are $3σ$ unless specified). This translates into a sub-Jupiter density. The planet revolves every $\sim7.8$ days in a slightly eccentric orbit ($e=0.123^{+0.037}_{-0.036}$) around a G8V star with detected activity in the Kepler light curve. Kepler-447b transits its host with a large impact parameter ($b=1.076^{+0.112}_{-0.086}$), being one of the few planetary grazing transits confirmed so far and the first in the Kepler large crop of exoplanets. We estimate that only around 20% of the projected planet disk occults the stellar disk. The relatively large uncertainties in the planet radius are due to the large impact parameter and short duration of the transit. Planets with such an extremely large impact parameter can be used to detect and analyze interesting configurations such as additional perturbing bodies, stellar pulsations, rotation of a non-spherical planet, or polar spot-crossing events. All these scenarios would periodically modify the transit properties (depth, duration, and time of mid-transit), what could be detectable with sufficient accurate photometry. Short-cadence photometric data (at the 1 minute level) would help in the search for these exotic configurations in grazing planetary transits like that of Kepler-447b.

preprint2015arXiv

PASTIS: Bayesian extrasolar planet validation II. Constraining exoplanet blend scenarios using spectroscopic diagnoses

The statistical validation of transiting exoplanets proved to be an efficient technique to secure the nature of small exoplanet signals which cannot be established by purely spectroscopic means. However, the spectroscopic diagnoses are providing us with useful constraints on the presence of blended stellar contaminants. In this paper, we present how a contaminating star affects the measurements of the various spectroscopic diagnoses as function of the parameters of the target and contaminating stars using the model implemented into the PASTIS planet-validation software. We find particular cases for which a blend might produce a large radial velocity signal but no bisector variation. It might also produce a bisector variation anti-correlated with the radial velocity one, as in the case of stellar spots. In those cases, the full width half maximum variation provides complementary constraints. These results can be used to constrain blend scenarios for transiting planet candidates or radial velocity planets. We review all the spectroscopic diagnoses reported in the literature so far, especially the ones to monitor the line asymmetry. We estimate their uncertainty and compare their sensitivity to blends. Based on that, we recommend the use of BiGauss which is the most sensitive diagnosis to monitor line-profile asymmetry. In this paper, we also investigate the sensitivity of the radial velocities to constrain blend scenarios and develop a formalism to estimate the level of dilution of a blended signal. Finally, we apply our blend model to re-analyse the spectroscopic diagnoses of HD16702, an unresolved face-on binary which exhibits bisector variations.

preprint2015arXiv

Polar stellar-spots and grazing planetary transits: possible explanation for the low number of discovered grazing planets

We assess a physically feasible explanation for the low number of discovered (near-)grazing planetary transits through all ground and space based transit surveys. We performed simulations to generate the synthetic distribution of detectable planets based on their impact parameter, and found that a larger number of (near-)grazing planets should have been detected than have been detected. Our explanation for the insufficient number of (near-)grazing planets is based on a simple assumption that a large number of (near-)grazing planets transit host stars which harbor dark giant polar spot, and thus the transit light-curve vanishes due to the occultation of grazing planet and the polar spot. We conclude by evaluating the properties required of polar spots in order to make disappear the grazing transit light-curve, and we conclude that their properties are compatible with the expected properties from observations.

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

Solar analogs with and without planets: T$_c$ trends and galactic evolution

We explore a sample of 148 solar-like stars to search for a possible correlation between the slopes of the abundance trends versus condensation temperature (known as the Tc slope) both with stellar parameters and Galactic orbital parameters in order to understand the nature of the peculiar chemical signatures of these stars and the possible connection with planet formation. We find that the Tc slope correlates at a significant level with the stellar age and the stellar surface gravity. We also find tentative evidence that the Tc slope correlates with the mean galactocentric distance of the stars (Rmean), suggesting that stars that originated in the inner Galaxy have fewer refractory elements relative to the volatile ones. We found that the chemical peculiarities (small refractory-to-volatile ratio) of planet-hosting stars is probably a reflection of their older age and their inner Galaxy origin. We conclude that the stellar age and probably Galactic birth place are key to establish the abundances of some specific elements.

preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

preprint2015arXiv

The first radial velocity measurements of a microlensing event: no evidence for the predicted binary

The gravitational microlensing technique allows the discovery of exoplanets around stars distributed in the disk of the galaxy towards the bulge. However, the alignment of two stars that led to the discovery is unique over the timescale of a human life and cannot be re-observed. Moreover, the target host is often very faint and located in a crowded region. These difficulties hamper and often make impossible the follow-up of the target and study of its possible companions. Gould et al. (2013) predicted the radial-velocity curve of a binary system, OGLE-2011-BLG-0417, discovered and characterised from a microlensing event by Shin et al. (2012). We used the UVES spectrograph mounted at the VLT, ESO to derive precise radial-velocity measurements of OGLE-2011-BLG-0417. To gather high-precision on faint targets of microlensing events, we proposed to use the source star as a reference to measure the lens radial velocities. We obtained ten radial velocities on the putative V=18 lens with a dispersion of ~100 m/s, spread over one year. Our measurements do not confirm the microlensing prediction for this binary system. The most likely scenario is that the assumed V=18 mag lens is actually a blend and not the primary lens that is 2 magnitude fainter. Further observations and analyses are needed to understand the microlensing observation and infer on the nature and characteristics of the lens itself.

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

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

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 hosts reveal lithium depletion: Results from a homogeneous statistical analysis

Aims. We study the impact of the presence of planets on the lithium abundance of host stars and evaluate the previous claim that planet hosts exhibit lithium depletion when compared to their non-host counterparts. Methods. Using previously published lithium abundances, we remove the confounding effect of the different fundamental stellar parameters by applying a multivariable regression on our dataset. In doing so, we explicitly make an assumption made implicitly by different authors: that lithium abundance depends linearly on fundamental stellar parameters. Using a moderator variable to distinguish stars with planets from those without, we evaluate the existence of an offset in lithium abundances between the two groups. We perform this analysis first for stars that present a clear lithium detection exclusively and include in a second analysis upper lithium measurements. Results. Our analysis shows that under the above-mentioned assumption of linearity, a statistically significant negative offset in lithium abundance between planet hosts and non-hosts is recovered. We concluded that an inflation on the lithium uncertainty estimations by a factor of larger than 5 is required to render the measured offset compatible with zero at less than 3-4 σ and make it non-significant. We demonstrated that the offset as delivered by our method depends on the different nature of the stars in the two samples. We did so by showing that the offset is reduced down to zero if the planet-host stars are replaced by comparison stars. Moreover, the measured depletion is still significant when one imposes different constraints on the dataset, such as a limit in planetary mass or constrain the host temperature to around solar value. We conclude then that planet-host stars exhibit enhanced lithium depletion when compared with non-host stars. (abridged)

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

Impact of occultations of stellar active regions on transmission spectra: Can occultation of a plage mimic the signature of a blue sky?

Transmission spectroscopy during planetary transits, which is based on the measurements of the variations of planet-to-star radius ratio as a function of wavelength, is a powerful technique to study the atmospheric properties of transiting planets. One of the main limitation of this technique is the effects of stellar activity, which up until now, have been taken into account only by assessing the effect of non-occulted stellar spots on the estimates of planet-to-star radius ratio. In this paper, we study, for the first time, the impact of the occultation of a stellar spot and plage on the transmission spectra of transiting exoplanets. We simulated this effect by generating a large number of transit light curves for different transiting planets, stellar spectral types, and for different wavelengths. Results of our simulations indicate that the anomalies inside the transit light curve can lead to a significant underestimation or overestimation of the planet-to-star radius ratio as a function of wavelength. At short wavelengths, the effect can reach to a difference of up to 10% in the planet-to-star radius ratio, mimicking the signature of light scattering in the planetary atmosphere. Atmospheric scattering has been proposed to interpret the increasing slopes of transmission spectra toward blue for exoplanets HD 189733b and GJ 3470b. Here we show that these signatures can be alternatively interpreted by the occultation of stellar plages. Results also suggest that the best strategy to identify and quantify the effects of stellar activities on the transmission spectrum of a planet is to perform several observations during the transit epoch at the same wavelength. This will allow for identifying the possible variations in transit depth as a function of time due to stellar activity variability.

preprint2014arXiv

Li abundances in F stars: planets, rotation and galactic evolution

We find that hot jupiter host stars within the T$_\mathrm{eff}$ range 5900-6300K show lower Li abundances, by 0.14 dex, than stars without detected planets. This offset has a significance at the level 7$σ$, pointing to a stronger effect of planet formation on Li abundances when the planets are more massive and migrate close to the star. However, we also find that the average v \textit{sin}i of (a fraction of) stars with hot jupiters is higher on average than for single stars in the same T$_\mathrm{eff}$ region, suggesting that rotationally-induced mixing (and not the presence of planets) might be the cause for a greater depletion of Li. We confirm that the mass-metallicity dependence of the Li dip is extended towards [Fe/H] $\sim$ 0.4 dex (beginning at [Fe/H] $\sim$ -0.4 dex for our stars) and that probably reflects the mass-metallicity correlation of stars of the same T$_\mathrm{eff}$ on the Main Sequence. We find that for the youngest stars ($<$ 1.5 Gyr) around the Li dip, the depletion of Li increases with v \textit{sin}i values, as proposed by rotationally-induced depletion models. This suggests that the Li dip consists of fast rotators at young ages whereas the most Li-depleted old stars show lower rotation rates (probably caused by the spin-down during their long lifes). We have also explored the Li evolution with [Fe/H] taking advantage of the metal-rich stars included in our sample. We find that Li abundance reaches its maximum around solar metallicity but decreases in the most metal-rich stars, as predicted by some models of Li Galactic production.

preprint2014arXiv

On the origin of stars with and without planets. Tc trends and clues to Galactic evolution

We explore a sample of 148 solar-like stars to search for a possible correlation between the slopes of the abundance trends versus condensation temperature (known as the Tc slope) with stellar parameters and Galactic orbital parameters in order to understand the nature of the peculiar chemical signatures of these stars and the possible connection with planet formation. We find that the Tc slope significantly correlates (at more than 4sigma) with the stellar age and the stellar surface gravity. We also find tentative evidence that the Tc slope correlates with the mean galactocentric distance of the stars (Rmean), suggesting that those stars that originated in the inner Galaxy have fewer refractory elements relative to the volatiles. While the average Tc slope for planet-hosting solar analogs is steeper than that of their counterparts without planets, this difference probably reflects the difference in their age and Rmean. We conclude that the age and probably the Galactic birth place are determinant to establish the star's chemical properties. Old stars (and stars with inner disk origin) have a lower refractory-to-volatile ratio.

preprint2014arXiv

Radial velocity confirmation of Kepler-91 b. Additional evidence of its planetary nature using the Calar Alto/CAFE instrument

The object transiting the star Kepler-91 was recently assessed as being of planetary nature. The confirmation was achieved by analysing the light-curve modulations observed in the Kepler data. However, quasi-simultaneous studies claimed a self-luminous nature for this object, thus rejecting it as a planet. In this work, we apply an {independent} approach to confirm the planetary mass of Kepler-91b by using multi-epoch high-resolution spectroscopy obtained with the Calar Alto Fiber-fed Echelle spectrograph (CAFE). We obtain the physical and orbital parameters with the radial velocity technique. In particular, we derive a value of $1.09 \pm 0.20\,M_{\mathrm{Jup}}$ for the mass of Kepler-91b, in excellent agreement with our previous estimate that was based on the orbital brightness modulation.

preprint2014arXiv

Revisiting the correlation between stellar activity and planetary surface gravity

Aims: We re-evaluate the correlation between planetary surface gravity and stellar host activity as measured by the index log($R'_{HK}$). This correlation, previously identified by Hartman (2010), is now analyzed in light of an extended measurements dataset, roughly 3 times larger than the original one. Methods: We calculated the Spearman's rank correlation coefficient between the two quantities and its associated p-value. The correlation coefficient was calculated for both the full dataset and the star-planet pairs that follow the conditions proposed by Hartman (2010). In order to do so, we considered effective temperatures both as collected from the literature and from the SWEET-Cat catalog, which provides a more homogeneous and accurate effective temperature determination. Results: The analysis delivers significant correlation coefficients, but with a lower value than those obtained by Hartman (2010). Yet, the two datasets are compatible, and we show that a correlation coefficient as large as previously published can arise naturally from a small-number statistics analysis of the current dataset. The correlation is recovered for star-planet pairs selected using the different conditions proposed by Hartman (2010). Remarkably, the usage of SWEET-Cat temperatures leads to larger correlation coefficient values. We highlight and discuss the role of the correlation betwen different parameters such as effective temperature and activity index. Several additional effects on top of those discussed previously were considered, but none fully explains the detected correlation. In light of the complex issue discussed here, we encourage the different follow-up teams to publish their activity index values in the form of log($R'_{HK}$) index so that a comparison across stars and instruments can be pursued.

preprint2014arXiv

Stars with and without planets: Where do they come from?

A long and thorough investigation of chemical abundances of planet-hosting stars that lasted for more than a decade has finally beared fruit. We explore a sample of 148 solar-like stars to search for a possible correlation between the slopes of the abundance trends versus condensation temperature (known as the Tc slope) both with stellar parameters and Galactic orbital parameters in order to understand the nature of the peculiar chemical signatures of these stars and the possible connection with planet formation. We find that the Tc slope correlates at a significant level (at more than 4sigma) with the stellar age and the stellar surface gravity. We also find tentative evidence that the Tc slope correlates with the mean galactocentric distance of the stars (Rmean), suggesting that stars that originated in the inner Galaxy have fewer refractory elements relative to the volatile ones. We found that the chemical peculiarities (small refractory-to-volatile ratio) of planet-hosting stars is merely a reflection of their older age and their inner Galaxy origin. We conclude that the stellar age and probably Galactic birth place are key to establish the abundances of some specific elements.

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 PLATO 2.0 Mission

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

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

Is there plenty of metal-poor stars with planets in the Galactic thick disk?

We performed an uniform spectroscopic analysis of 1111 FGK dwarfs observed as part of the HARPS GTO planet search program. We applied a purely chemical approach, based on [alpha/Fe] ratio, to distinguish the various stellar components in the Galaxy. Apart from the well known Galactic thick and thin disks, we separated an alpha-enhanced stellar family at super-solar metallicities. The metal-rich high-alpha stars have orbits similar to the thin disk stars, but they are similar to thick disk stars in terms of age. Our data indicate that the incidence of stars with planets are greater among the chemically separated thick disk stars with [Fe/H] < -0.3 dex than they are among thin disk stars in the same [Fe/H] interval. Our results allow us to suppose that a certain chemical composition, and not the Galactic birth place of the stars, is the causative factor for that.

preprint2013arXiv

Kinematics and chemical properties of the Galactic stellar populations. The HARPS FGK dwarfs sample

(Abridged) We analyze chemical and kinematical properties of about 850 FGK solar neighborhood long-lived dwarfs observed with the HARPS high-resolution spectrograph. The stars in the sample have logg > 4 dex, 5000 < Teff < 6500 K, and -1.39 < [Fe/H] < 0.55 dex. We apply a purely chemical analysis approach based on the [alpha/Fe] vs. [Fe/H] plot to separate Galactic stellar populations into the thin disk, thick disk and high-alpha metal-rich (hamr). Our analysis shows a negative gradient of the rotational velocity of the thin disk stars with [Fe/H] (-17 km s^-1 dex^-1), and a steep positive gradient for both the thick disk and hamr stars with the same magnitude of about +42 km s^-1 dex^-1. For the thin disk stars we observed no correlation between orbital eccentricities and metallicity, but observed a steep negative gradient for the thick disk and hamr stars with practically the same magnitude (about -0.18 dex^-1). Our results suggest that radial migration played an important role in the formation and evolution of the thin disk. For the thick disk stars it is not possible to reach a firm conclusion about their origin. Based on the eccentricity distribution of the thick disk stars only their accretion origin can be ruled out, and the heating and migration scenario could explain the positive steep gradient of V_phi with [Fe/H]. Analyzing the hamr stellar population we found that they share properties of both the thin and thick disk population. A comparison of the properties of the hamr stars with that of the subsample of stars from the N-body/SPH simulation using radial migration suggest that they may have originated from the inner Galaxy. Further detailed investigations would help to clarify their exact nature and origin.

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

Orbital and physical properties of planets and their hosts: new insights on planet formation and evolution

We explore the relations between physical and orbital properties of planets and properties of their host stars to identify the main observable signatures of the formation and evolution processes of planetary systems. We use a large sample of FGK dwarf planet hosts with stellar parameters derived in a homogeneous way from the SWEET-Cat database to study the relation between stellar metallicity and position of planets in the period-mass diagram. In the second part we use all the RV-detected planets orbiting FGK stars to explore the role of planet-disk and planet-planet interaction on the evolution of orbital properties of planets with masses above 1MJup. We show that planets orbiting metal-poor stars have longer periods than those in metal-rich systems. This trend is valid for masses at least from 10MEarth to 4MJup. Earth-like planets orbiting metal-rich stars always show shorter periods (fewer than 20 days) than those orbiting metal-poor stars. We also found statistically significant evidence that very high mass giants have on average more eccentric orbits than giant planets with lower mass.Finally, we show that the eccentricity of planets with masses higher than 4MJup tends to be lower for planets with shorter periods. Our results suggest that the planets in the P-MP diagram are evolving differently because of a mechanism that operates over a wide range of planetary masses. This mechanism is stronger or weaker depending on the metallicity of the respective system. One possibility is that planets in metal-poor disks form farther out from their central star and/or they form later and do not have time to migrate as far as the planets in metal-rich systems. The trends and dependencies obtained for very high mass planetary systems suggest that planet-disk interaction is a very important and orbit-shaping mechanism for planets in the high-mass domain. Shortened.

preprint2013arXiv

Probing the effect of gravitational microlensing on the measurements of the Rossiter-McLaughlin effect

In general, in the studies of transit light-curves and the Rossiter-McLaughlin (RM), the contribution of the planet's gravitational microlensing is neglected. Theoretical studies, have, however shown that the planet's microlensing can affect the transit light-curve and in some extreme cases cause the transit depth to vanish. In this letter, we present the results of our quantitative analysis of microlening on the RM effect. Results indicate that for massive planets in on long period orbits, the planet's microlensing will have considerable contribution to the star's RV measurements. We present the details of our study, and discuss our analysis and results.

preprint2013arXiv

The contribution of secondary eclipses as astrophysical false positives to exoplanet transit surveys

We investigate in this paper the astrophysical false-positive configuration in exoplanet-transit surveys that involves eclipsing binaries and giant planets which present only a secondary eclipse, as seen from the Earth. To test how an eclipsing binary configuration can mimic a planetary transit, we generate synthetic light curve of three examples of secondary-only eclipsing binary systems that we fit with a circular planetary model. Then, to evaluate its occurrence we model a population of binaries in double and triple system based on binary statistics and occurrence. We find that 0.061% +/- 0.017% of main-sequence binary stars are secondary-only eclipsing binaries mimicking a planetary transit candidate down to the size of the Earth. We then evaluate the occurrence that an occulting-only giant planet can mimic an Earth-like planet or even smaller planet. We find that 0.009% +/- 0.002% of stars harbor a giant planet that present only the secondary transit. Occulting-only giant planets mimic planets smaller than the Earth that are in the scope of space missions like Kepler and PLATO. We estimate that up to 43.1 +/- 5.6 Kepler Objects of Interest can be mimicked by this new configuration of false positives, re-evaluating the global false-positive rate of the Kepler mission from 9.4% +/- 0.9% to 11.3% +/- 1.1%. We note however that this new false-positive scenario occurs at relatively long orbital period compared with the median period of Kepler candidates.

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

World-leading science with SPIRou - the nIR spectropolarimeter / high-precision velocimeter for CFHT

SPIRou is a near-infrared (nIR) spectropolarimeter / velocimeter proposed as a new-generation instrument for CFHT. SPIRou aims in particular at becoming world-leader on two forefront science topics, (i) the quest for habitable Earth-like planets around very- low-mass stars, and (ii) the study of low-mass star and planet formation in the presence of magnetic fields. In addition to these two main goals, SPIRou will be able to tackle many key programs, from weather patterns on brown dwarf to solar-system planet atmospheres, to dynamo processes in fully-convective bodies and planet habitability. The science programs that SPIRou proposes to tackle are forefront (identified as first priorities by most research agencies worldwide), ambitious (competitive and complementary with science programs carried out on much larger facilities, such as ALMA and JWST) and timely (ideally phased with complementary space missions like TESS and CHEOPS). SPIRou is designed to carry out its science mission with maximum efficiency and optimum precision. More specifically, SPIRou will be able to cover a very wide single-shot nIR spectral domain (0.98-2.35 μm) at a resolving power of 73.5K, providing unpolarized and polarized spectra of low-mass stars with a ~15% average throughput and a radial velocity (RV) precision of 1 m/s.

preprint2012arXiv

A new analysis of the WASP-3 system: no evidence for an additional companion

In this work we investigate the problem concerning the presence of additional bodies gravitationally bounded with the WASP-3 system. We present eight new transits of this planet and analyse all the photometric and radial velocity data published so far. We did not observe significant periodicities in the Fourier spectrum of the observed minus calculated (O-C) transit timing and radial velocity diagrams (the highest peak having false-alarm probabilities of 56 per cent and 31 per cent, respectively) or long-term trends. Combining all the available information, we conclude that the radial velocity and transit timing techniques exclude, at 99 per cent confidence limit, any perturber more massive than M \gtrsim 100 M_Earth with periods up to 10 times the period of the inner planet. We also investigate the possible presence of an exomoon on this system and determined that considering the scatter of the O-C transit timing residuals a coplanar exomoon would likely produce detectable transits. This hypothesis is however apparently ruled out by observations conducted by other researchers. In case the orbit of the moon is not coplanar the accuracy of our transit timing and transit duration measurements prevents any significant statement. Interestingly, on the basis of our reanalysis of SOPHIE data we noted that WASP-3 passed from a less active (log R'_hk=-4.95) to a more active (log R'_hk=-4.8) state during the 3 yr monitoring period spanned by the observations. Despite no clear spot crossing has been reported for this system, this analysis claims for a more intensive monitoring of the activity level of this star in order to understand its impact on photometric and radial velocity measurements.

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

Transit-timing measurements with the model-independent barycenter method: Application to the LHS 6343 system

We present a model-independent technique for calculating the time of mid-transits. This technique, named "barycenter method", uses the light-curve's symmetry to determine the transit timing by calculating the transit light-curve barycenter. Unlike the other methods of calculating mid-transit timing, this technique does not depend on the parameters of the system and central star. We demonstrate the capabilities of the barycenter method by applying this technique to some known transiting systems including several \emph{Kepler} confirmed planets. Results indicate that for complete and symmetric transit lightcurves, the barycenter method achieves the same precision as other techniques, but with fewer assumptions and much faster. Among the transiting systems studied with the barycenter method, we focus in particular on LHS 6343C, a brown dwarf that transits a member of an M+M binary system, LHS 6343AB. We present the results of our analysis, which can be used to set an upper limit on the period and mass of a possible second small perturber.

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

Exoplanets transmission spectroscopy: accounting for eccentricity and longitude of periastron. Superwinds in the upper atmosphere of HD209458b?

Context: Several studies have so far placed useful constraints on planetary atmospheric properties using transmission spectrsocopy, and in the case of HD209458b even the radial velocity of the planet during the transit event has been reconstructed opening a new range of possibilities. AIMS. In this contribution we highlight the importance to account for the orbital eccentricity and longitude of periastron of the planetary orbit to accurately interpret the measured planetary radial velocity during the transit. Methods: We calculate the radial velocity of a transiting planet in an eccentric orbit. Given the larger orbital speed of planets with respect to their stellar companions even small eccentricities can result in detectable blue or redshift radial velocity offsets during the transit with respect to the systemic velocity, the exact value depending also on the longitude of the periastron of the planetary orbit. For an hot-jupiter planet, an eccentricity of only e=0.01 can produce a radial velocity offset of the order of the km/s. Conclusions: We propose an alternative interpretation of the recently claimed radial velocity blueshift (~2 km/s) of the planetary spectral lines of HD209458b which implies that the orbit of this system is not exactly circular. In this case, the longitude of the periastron of the stellar orbit is most likely confined in the first quadrant (and that one of the planet in the third quadrant). We highlight that transmission spectroscopy allows not only to study the compositional properties of planetary atmospheres, but also to refine their orbital parameters and that any conclusion regarding the presence of windflows on planetary surfaces coming from transmission spectroscopy measurements requires precise known orbital parameters from RV.

preprint2011arXiv

Finding proto-spectroscopic binaries: Precise multi-epoch radial velocities of 7 protostars in rho-Ophiuchus

The formation of spectroscopic binaries (SB) may be a natural byproduct of star formation. The early dynamical evolution of multiple stellar systems after the initial fragmentation of molecular clouds leaves characteristic imprints on the properties of young, multiple stars. The discovery and the characterization of the youngest SB will allow us to infer the mechanisms and timescales involved in their formation. Our work aims to find spectroscopic companions around young stellar objects (YSO). We present a near-IR high-resolution (R ~ 60000) multi-epoch radial velocity survey of 7 YSO in the star forming region (SFR) rho-Ophiuchus. The radial velocities of each source were derived using a two-dimensional cross-correlation function, using the zero-point established by the Earth's atmosphere as reference. More than 14 spectral lines in the CO (0-2) bandhead window were used in the cross-correlation against LTE atmospheric models to compute the final results. We found that the spectra of the protostars in our sample agree well with the predicted stellar photospheric profiles, indicating that the radial velocities derived are indeed of stellar nature. Three of the targets analyzed exhibit large radial velocity variations during the three observation epochs. These objects - pending further confirmation and orbital characteristics - may become the first evidence for proto-spectroscopic binaries, and will provide important constraints on their formation. Our preliminary binary fraction (BF) of ~71% (when merging our results with those of previous studies) is in line with the notion that multiplicity is very high at young ages and therefore a byproduct of star formation

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

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

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.

preprint2009arXiv

Bulk composition of the transiting hot Neptune around GJ 436

The hot Neptune orbiting around GJ 436 is a unique example of an intermediate mass planet. Its close-in orbit suggests that the planet has undergone migration and its study is fundamental to understanding planet formation and evolution. As it transits its parent star, it is the only Neptune-mass extrasolar planet of known mass and radius, being slightly larger and more massive than Neptune (M=22.6 M_Earth, R=4.19R_Earth). In this regime, several bulk compositions are possible: from an Earth-like core with a thick hydrogen envelope to a water-rich planet with a thin hydrogen envelope comprising a Neptune-like structure. We combine planet-structure modeling with an advanced planet-formation model to assess the likelihood of the different possible bulk compositions of GJ 436 b. We find that both an envelope-free water planet ("Ocean planet") as well as a diminute version of a gaseous giant planet are excluded. Consisting of a rocky core with a thick hydrogen/helium envelope, a "dry" composition produces not only too small a radius but is also a very unlikely outcome of planet formation around such a low-mass star. We conclude that GJ 436 b is probably of much higher rock content than Neptune (more than 45% in mass), with a small H-He envelope (10 - 20% in mass). This is the expected outcome of the gathering of materials during the migration process in the inner disk, creating a population of which the hot Neptune is representative.

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 CORALIE survey for southern extrasolar planets. XVI. Discovery of a planetary system around HD 147018 and of two long period and massive planets orbiting HD 171238 and HD 204313

We report the detection of a double planetary system around HD 140718 as well as the discovery of two long period and massive planets orbiting HD 171238 and HD 204313. Those discoveries were made with the CORALIE Echelle spectrograph mounted on the 1.2-m Euler Swiss telescope located at La Silla Observatory, Chile. The planetary system orbiting the nearby G9 dwarf HD 147018 is composed of an eccentric inner planet (e=0.47) with twice the mass of Jupiter (2.1 MJup ) and with an orbital period of 44.24 days. The outer planet is even more massive (6.6 MJup) with a slightly eccentric orbit (e=0.13) and a period of 1008 days. The planet orbiting HD 171238 has a minimum mass of 2.6 MJup, a period of 1523 days and an eccentricity of 0.40. It orbits a G8 dwarfs at 2.5 AU. The last planet, HD 204313 b, is a 4.0 MJup -planet with a period of 5.3 years and has a low eccentricity (e = 0.13). It orbits a G5 dwarfs at 3.1 AU. The three parent stars are metal rich, which further strengthened the case that massive planets tend to form around metal rich stars.