sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Faria, J. P.

Normalized to: Faria, J.

45 article(s) in total. 286 co-authors, from 1 to 40 common article(s). Median position in authors list is 6,0.

[1]  oai:arXiv.org:2007.06410  [pdf] - 2132073
A precise architecture characterization of the $\pi$ Men planetary system
Comments: Accepted for publication on A&A
Submitted: 2020-07-13
The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.
[2]  oai:arXiv.org:2005.14008  [pdf] - 2103526
Scheduling ESPRESSO follow-up of TESS Targets. I. Myopic versus non-myopic sampling
Comments: 15 pages, 10 figures, submitted to MNRAS
Submitted: 2020-05-28
Radial-velocity follow-up of stars harbouring transiting planets detected by TESS is expected to require very large amounts of expensive telescope time in the next few years. Therefore, scheduling strategies should be implemented to maximize the amount of information gathered about the target planetary systems. We present one random scheduler and two types of uniform-in-phase schedulers: one myopic, which selects targets one-at-a-time, and one non-myopic that efficiently explores all the possible combinations between stars to be observed and available time slots. We compare these strategies with respect to the bias, accuracy and precision achieved in recovering the mass and orbital parameters of transiting and non-transiting planets from the mock radial-velocity follow-up of a sample of 50 TESS target stars, with simulated planetary systems containing at least one transiting planet with a radius below $4R_{\oplus}$. For each system and strategy, 10 radial-velocity datasets were generated and analysed using a fully Bayesian methodology. We find the myopic strategies lead to a biased estimation of the order of 5\% of the mass of the transiting exoplanets while the non-myopic scheduler is able to provide an unbiased (<1\%) measurement of the masses while keeping the relative accuracy and precision around 16\% and 23\% respectively. The number of non-transiting planets detected is similar for all the strategies considered, although the random scheduler leads to less biased and more accurate estimates for their mass and orbital parameters, possibly due to a higher mean number of scheduled radial-velocities for the datasets associated with non-transiting planets detections.
[3]  oai:arXiv.org:2005.12114  [pdf] - 2101585
Revisiting Proxima with ESPRESSO
Comments: 25 pages, 26 figures
Submitted: 2020-05-25, last modified: 2020-05-26
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.
[4]  oai:arXiv.org:2005.10105  [pdf] - 2129363
The correlation between photometric variability and radial velocity jitter, based on TESS and HARPS observations
Comments: Accepted for publication in section 10. Planets and planetary systems of A&A
Submitted: 2020-05-20, last modified: 2020-05-21
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.
[5]  oai:arXiv.org:2002.11422  [pdf] - 2061775
Can planetary rings explain the extremely low density of HIP 41378f?
Comments: 9 figures, 8 pages, accepted by A&A
Submitted: 2020-02-26
The presence of rings around a transiting planet can cause its radius to be overestimated and lead to an underestimation of its density if the mass is known. We employ a Bayesian framework to show that the anomalously low density ($\sim$0.09 g cm${^{-3}}$) of the transiting long-period planet HIP$\,$41378$\,f$ might be due to the presence of opaque circum-planetary rings. Given our adopted model priors and data from the K2 mission, we find the statistical evidence for the ringed planet scenario to be comparable to that of the planet-only scenario. The ringed planet solution suggests a larger planetary density of $\sim$1.23$\,$g$\,$cm$^{-3}$ similar to Uranus. The associated ring extends from 1.05 to 2.59 times the planetary radius and is inclined away from the sky-plane by $\sim$25$^\mathrm{o}$. Future high-precision transit observations of HIP$\,$41378$\,f$ would be necessary to confirm/dismiss the presence of planetary rings.
[6]  oai:arXiv.org:2002.08227  [pdf] - 2065432
Can we detect the stellar differential rotation of WASP-7 through theRossiter-McLaughlin observations?
Comments: 12 pages plus appendix, 9 figures
Submitted: 2020-02-19
The Rossiter-McLaughlin (RM) effect is the radial velocity signal generated when an object transits a rotating star. Stars rotate differentially and this affects the shape and amplitude of this signal, on a level that can no longer be ignored with precise spectrographs. Highly misaligned planets provide a unique opportunity to probe stellar differential rotation via the RM effect, as they cross several stellar latitudes. In this sense, WASP-7, and its hot Jupiter with a projected misalignment of \sim 90{\deg}, is one of the most promising targets. The aim of this work is to understand if the stellar differential rotation is measurable through the RM signal for systems with a geometry similar to WASP-7. In this sense, we use a modified version of SOAP3.0 to explore the main hurdles that prevented the precise determination of the differential rotation of WASP-7. We also investigate whether the adoption of the next generation spectrographs, like ESPRESSO, would solve these issues. Additionally, we assess how instrumental and stellar noise influence this effect and the derived geometry of the system. We found that, for WASP-7, the white noise represents an important hurdle in the detection of the stellar differential rotation, and that a precision of at least 2m/s or better is essential.
[7]  oai:arXiv.org:2001.07975  [pdf] - 2046489
Improving transit characterisation with Gaussian process modelling of stellar variability
Comments: 13 pages, 5 figures, 11 tables accepted for publication in A&A
Submitted: 2020-01-22
New photometric space missions to detect and characterise transiting exoplanets are focusing on bright stars to obtain high cadence, high signal-to-noise light curves. Since these missions will be sensitive to stellar oscillations and granulation even for dwarf stars, they will be limited by stellar variability. We tested the performance of Gaussian process (GP) regression on the characterisation of transiting planets, and in particular to determine how many components of variability are necessary to describe high cadence, high signal-to-noise light curves expected from CHEOPS and PLATO. We found that the best GP stellar variability model contains four to five variability components: one stellar oscillation component, two to four granulation components, and/or one rotational modulation component. This high number of components is in contrast with the one-component GP model (1GP) commonly used in the literature for transit characterisation. Therefore, we compared the performance of the best multi-component GP model with the 1GP model in the derivation of transit parameters of simulated transits. We found that for Jupiter- and Neptune-size planets the best multi-component GP model is slightly better than the 1GP model, and much better than the non-GP model that gives biased results. For Earth-size planets, the 1GP model fails to retrieve the transit because it is a poor description of stellar activity. The non-GP model gives some biased results and the best multi-component GP is capable of retrieving the correct transit model parameters. We conclude that when characterising transiting exoplanets with high signal-to-noise ratios and high cadence light curves, we need models that couple the description of stellar variability with the transits analysis, like GPs. Moreover, for Earth-like exoplanets a better description of stellar variability improves the planetary characterisation.
[8]  oai:arXiv.org:1911.11714  [pdf] - 2057802
Decoding the radial velocity variations of HD41248 with ESPRESSO
Comments: Accepted for publication in A&A
Submitted: 2019-11-26
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.
[9]  oai:arXiv.org:1806.04683  [pdf] - 2057659
Quantifying the Bayesian Evidence for a Planet in Radial Velocity Data
Comments: Accepted in AJ. 33 pages, 8 figures, 4 tables. Evidence Challenge repository can be found at https://github.com/EPRV3EvidenceChallenge/Inputs
Submitted: 2018-06-12, last modified: 2019-10-27
We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.
[10]  oai:arXiv.org:1909.13527  [pdf] - 1987468
Exoplanet characterisation in the longest known resonant chain: the K2-138 system seen by HARPS
Comments:
Submitted: 2019-09-30
The detection of low-mass transiting exoplanets in multiple systems brings new constraints to planetary formation and evolution processes and challenges the current planet formation theories. Nevertheless, only a mere fraction of the small planets detected by Kepler and K2 have precise mass measurements, which are mandatory to constrain their composition. We aim to characterise the planets that orbit the relatively bright star K2-138. This system is dynamically particular as it presents the longest chain known to date of planets close to the 3:2 resonance. We obtained 215 HARPS spectra from which we derived the radial-velocity variations of K2-138. Via a joint Bayesian analysis of both the K2 photometry and HARPS radial-velocities (RVs), we constrained the parameters of the six planets in orbit. The masses of the four inner planets, from b to e, are 3.1, 6.3, 7.9, and 13.0 $\mathrm{M}_{\oplus}$ with a precision of 34%, 20%, 18%, and 15%, respectively. The bulk densities are 4.9, 2.8, 3.2, and 1.8 g cm$^{-3}$, ranging from Earth to Neptune-like values. For planets f and g, we report upper limits. Finally, we predict transit timing variations of the order two to six minutes from the masses derived. Given its peculiar dynamics, K2-138 is an ideal target for transit timing variation (TTV) measurements from space with the upcoming CHaracterizing ExOPlanet Satellite (CHEOPS) to study this highly-packed system and compare TTV and RV masses.
[11]  oai:arXiv.org:1909.07392  [pdf] - 1975457
Beyond the exoplanet mass-radius relation
Comments: Comments: 9 pages, 7 figures, accepted for publication in A&A, the code BEM used to predict exoplanet radii is available online, see https://github.com/soleneulmer/bem
Submitted: 2019-09-16
Mass and radius are two fundamental properties for characterising exoplanets, but only for a relatively small fraction of exoplanets are they both available. Mass is often derived from radial velocity measurements, while the radius is almost always measured using the transit method. For a large number of exoplanets, either the radius or the mass is unknown, while the host star has been characterised. Several mass-radius relations that are dependent on the planet's type have been published that often allow us to predict the radius. Our goal is to derive the radius of exoplanets using only observables extracted from spectra used primarily to determine radial velocities and spectral parameters. Our objective is to obtain a mass-radius relation independent of the planet's type. We worked with a database of confirmed exoplanets with known radii and masses, as well as the planets from our Solar System. Using random forests, a machine learning algorithm, we computed the radius of exoplanets and compared the results to the published radii. Our code, Bem, is available online. The estimated radii reproduces the spread in radius found for high mass planets better than previous mass-radius relations. The average radius error is 1.8 $R_{\oplus}$ across the whole range of radii from 1-22 $R_{\oplus}$. We find that a random forest algorithm is able to derive reliable radii, especially for planets between 4 $R_{\oplus}$ and 20 $R_{\oplus}$ for which the error is under 25$\%$. The algorithm has a low bias yet a high variance, which could be reduced by limiting the growth of the forest, or adding more data. The random forest algorithm is a promising method for deriving exoplanet properties. We show that the exoplanet's mass and equilibrium temperature are the relevant properties that constrain the radius, and do so with higher accuracy than the previous methods.
[12]  oai:arXiv.org:1908.10662  [pdf] - 1961679
Gaussian Process modelling of granulation and oscillations in red-giant stars
Comments: Accepted for publication in MNRAS; 12 pages, 10 figures, 2 tables
Submitted: 2019-08-28
The analysis of photometric time series in the context of transiting planet surveys suffers from the presence of stellar signals, often dubbed "stellar noise". These signals, caused by stellar oscillations and granulation, can usually be disregarded for main-sequence stars, as the stellar contributions average out when phase-folding the light curve. For evolved stars, however, the amplitudes of such signals are larger and the timescales similar to the transit duration of short-period planets, requiring that they be modeled alongside the transit. With the promise of TESS delivering on the order of $\sim\!10^5$ light curves for stars along the red-giant branch, there is a need for a method capable of describing the "stellar noise" while simultaneously modelling an exoplanet's transit. In this work, a Gaussian Process regression framework is used to model stellar light curves and the method validated by applying it to TESS-like artificial data. Furthermore, the method is used to characterize the stellar oscillations and granulation of a sample of well-studied \textit{Kepler} low-luminosity red-giant branch stars. The parameters determined are compared to equivalent ones obtained by modelling the power spectrum of the light curve. Results show that the method presented is capable of describing the stellar signals in the time domain and can also return an accurate and precise measurement of $\nu_\text{max}$, i.e., the frequency of maximum oscillation amplitude. Preliminary results show that using the method in transit modelling improves the precision and accuracy of the ratio between the planetary and stellar radius, $R_p/R_\star$. The method's implementation is publicly available.
[13]  oai:arXiv.org:1903.04937  [pdf] - 1847099
The Metallicity-Period-Mass Diagram of low-mass exoplanets
Comments: 12 pages, 9 Figures, Accepted for MNRAS
Submitted: 2019-03-12
The number of exoplanet detections continues to grow following the development of better instruments and missions. Key steps for the understanding of these worlds comes from their characterization and its statistical studies. We explore the metallicity-period-mass diagram for known exoplanets by using an updated version of The Stellar parameters for stars With ExoplanETs CATalog (SWEET-Cat), a unique compilation of precise stellar parameters for planet-host stars provided for the exoplanet community. Here we focus on the planets with minimum mass below 30 M$_{\oplus}$ which seems to present a possible correlation in the metallicity-period-mass diagram where the mass of the planet increases with both metallicity and period. Our analysis suggests that the general observed correlation may be not fully explained by observational biases. Additional precise data will be fundamental to confirm or deny this possible correlation.
[14]  oai:arXiv.org:1811.11172  [pdf] - 1791658
ACTIN: A tool to calculate stellar activity indices
Comments: 2 pages; published in the Journal of Open Source Software; all code available at https://github.com/gomesdasilva/ACTIN
Submitted: 2018-11-22
Magnetic activity in the atmospheres of stars produces a number of spectroscopic signatures that are visible in the shape and strength of spectral lines. These signatures can be used to access, among other things, the variability of the magnetic activity, or its infuence on other parameters such as the measured radial velocity (RV). This latter is of utmost importance for the detection and characterization of planets orbiting other stars. ACTIN is a Python program to calculate stellar activity indices. The program reads input data either from .fits files returned by the pipelines of spectrographs, or from .rdb tables. It extracts automatically the spectral data required to calculate spectral activity indices. The output is an .rdb table, with the calculated stellar activity indices for each date (Julian Date), as well as the RV and Cross-Correlation Function profile parameters, if available. It also outputs timeseries plots of the activity indices and plots the spectral lines used to compute the indices.
[15]  oai:arXiv.org:1811.07776  [pdf] - 1815070
The GAPS Programme with HARPS-N at TNG XVIII. Two new giant planets around the metal-poor stars HD 220197 and HD 233832
Comments: 15 pages, 8 figures, 5 tables, accepted for publication in A&A
Submitted: 2018-11-19
Statistical studies of exoplanets have shown that giant planets are more commonly hosted by metal-rich dwarf stars than low-metallicity ones, while such a correlation is not evident for lower-mass planets. The search for giant planets around metal-poor stars and the estimate of their occurrence $f_p$ is an important element in providing support to models of planet formation. We present results from the HARPS-N search for giant planets orbiting metal-poor ($-1.0\leq[Fe/H]\leq-0.5$ dex) stars in the northern hemisphere complementing a previous HARPS survey on southern stars in order to update the estimate of $f_p$. High-precision HARPS-N observations of 42 metal-poor stars are used to search for planetary signals to be fitted using differential evolution MCMC single-Keplerian models. We then join our detections to the results of the previous HARPS survey on 88 metal-poor stars to provide a preliminar estimate of the two-hemisphere $f_p$. We report the detection of two new giant planets around HD 220197 and HD 233832. The first companion has M$\sin{i}=0.20_{-0.04}^{+0.07}$ M$_{\rm Jup}$ and orbital period of $1728_{-80}^{+162}$ days, and for the second companion we find two solutions of equal statistical weight having periods $2058_{-40}^{+47}$ and $4047_{-117}^{+91}$ days and minimum masses of $1.78_{-0.06}^{+0.08}$ and $2.72_{-0.23}^{+0.23}$ M$_{\rm Jup}$, respectively. Joining our two detections with the three from the southern survey we obtain a preliminary and conservative estimate of global frequency of $f_p=3.84_{-1.06}^{+2.45}\%$ for giant planets around metal-poor stars. The two new giant planets orbit dwarf stars at the metal-rich end of the HARPS-N metal-poor sample, corroborating previous results suggesting that giant planet frequency still is a rising function of host star [Fe/H]. We also note that all detections in the overall sample are giant long-period planets.
[16]  oai:arXiv.org:1810.08108  [pdf] - 1791208
SWEET-Cat updated. New homogenous spectroscopic parameters
Comments: 13 pages, 4 figures, 3 tables. accepted for A&A; Updated sweet-cat: https://www.astro.up.pt/resources/sweet-cat/
Submitted: 2018-10-18
Context: Exoplanets have now been proven to be very common. The number of its detections continues to grow following the development of better instruments and missions. One key step for the understanding of these worlds is their characterization, which mostly depend on their host stars. Aims:We perform a significant update of the Stars With ExoplanETs CATalog (SWEET-Cat), a unique compilation of precise stellar parameters for planet-host stars provided for the exoplanet community. Methods: We made use of high-resolution spectra for planet-host stars, either observed by our team or found in several public archives. The new spectroscopic parameters were derived for the spectra following the same homogeneous process (ARES+MOOG). The host star parameters were then merged together with the planet properties listed in exoplanet.eu to perform simple data analysis. Results: We present new spectroscopic homogeneous parameters for 106 planet-host stars. Sixty-three planet hosts are also reviewed with new parameters. We also show that there is a good agreement between stellar parameters derived for the same star but using spectra obtained from different spectrographs. The planet-metallicity correlation is reviewed showing that the metallicity distribution of stars hosting low-mass planets (below 30 M$_{\oplus}$) is indistinguishable from that from the solar neighborhood sample in terms of metallicity distribution.
[17]  oai:arXiv.org:1809.01027  [pdf] - 1789494
Activity induced variation in spin-orbit angles as derived from Rossiter-McLaughlin measurements
Comments: 18 pages, 7 figures, 6 tables, accepted for publication in A&A
Submitted: 2018-09-04
One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling.
[18]  oai:arXiv.org:1805.08021  [pdf] - 1727202
Recovering the colour-dependent albedo of exoplanets with high-resolution spectroscopy: from ESPRESSO to the ELT
Comments: Published in MNRAS; 27 pages; 16 figures; 11 tables
Submitted: 2018-05-21, last modified: 2018-08-03
The characterization of planetary atmospheres is a daunting task, pushing current observing facilities to their limits. The next generation of high-resolution spectrographs mounted on large telescopes -- such as ESPRESSO@VLT and HIRES@ELT -- will allow us to probe and characterize exoplanetary atmospheres in greater detail than possible to this point. We present a method that permits the recovery of the colour-dependent reflectivity of exoplanets from high-resolution spectroscopic observations. Determining the wavelength-dependent albedo will provide insight into the chemical properties and weather of the exoplanet atmospheres. For this work, we simulated ESPRESSO@VLT and HIRES@ELT high-resolution observations of known planetary systems with several albedo configurations. We demonstrate how the cross correlation technique applied to theses simulated observations can be used to successfully recover the geometric albedo of exoplanets over a range of wavelengths. In all cases, we were able to recover the wavelength dependent albedo of the simulated exoplanets and distinguish between several atmospheric models representing different atmospheric configurations. In brief, we demonstrate that the cross correlation technique allows for the recovery of exoplanetary albedo functions from optical observations with the next generation of high-resolution spectrographs that will be mounted on large telescopes with reasonable exposure times. Its recovery will permit the characterization of exoplanetary atmospheres in terms of composition and dynamics and consolidates the cross correlation technique as a powerful tool for exoplanet characterization.
[19]  oai:arXiv.org:1807.09608  [pdf] - 1775607
Planets around evolved intermediate-mass stars in open clusters II. Are there really planets around IC4651No9122, NGC2423No3 and NGC4349No127?
Comments: 15 pages, 18 figures, accepted by A&A
Submitted: 2018-07-25
(shorter version)The aim of this work is to search for planets around intermediate-mass stars in open clusters by using RV data obtained with HARPS from an extensive survey with more than 15 years of observations for a sample of 142 giant stars in 17 open clusters. We present the discovery of a periodic RV signal compatible with the presence of a planet candidate in the 1.15 Gyr open cluster IC4651 orbiting the 2.06 M$_\odot$ star No. 9122. If confirmed, the planet candidate would have a minimum mass of 7.2 M$_{J}$ and a period of 747 days. However, we also find that the FWHM of the CCF varies with a period close to the RV, casting doubts on the planetary nature of the signal. We also provide refined parameters for the previously discovered planet around NGC2423 No. 3 but show evidence that the BIS of the CCF is correlated with the RV during some of the observing periods. This fact advises us that this might not be a real planet and that the RV variations could be caused by stellar activity and/or pulsations. Finally, we show that the previously reported signal by a brown dwarf around NGC4349 No. 127 is presumably produced by stellar activity modulation. The long-term monitoring of several red giants in open clusters has allowed us to find periodic RV variations in several stars. However, we also show that the follow-up of this kind of stars should last more than one orbital period to detect long-term signals of stellar origin. This work warns that although it is possible to detect planets around red giants, large-amplitude, long-period RV modulations do exist in such stars that can mimic the presence of an orbiting planetary body. Therefore, we need to better understand how such RV modulations behave as stars evolve along the RGB and perform a detailed study of all the possible stellar-induced signals (e.g. spots, pulsations, granulation) to comprehend the origin of RV variations.
[20]  oai:arXiv.org:1807.00773  [pdf] - 1767496
The TROY project: II. Multi-technique constraints on exotrojans in nine planetary systems
Comments: Accepted for publication in A&A, 18 pages, 13 figures, 12 tables. TROY project website: http://www.sc.eso.org/~jlillobo/troy/
Submitted: 2018-07-02
Co-orbital bodies are the byproduct of planet formation and evolution, as we know from the Solar System. Although planet-size co-orbitals do not exists in our planetary system, dynamical studies show that they can remain stable for long periods of time in the gravitational well of massive planets. Should they exist, their detection is feasible with the current instrumentation. In this paper, we present new ground-based observations searching for these bodies co-orbiting with nine close-in (P<5 days) planets, using different observing techniques. The combination of all of them allows us to restrict the parameter space of any possible trojan in the system. We use multi-technique observations (radial velocity, precision photometry and transit timing variations), both newly acquired in the context of the TROY project and publicly available, to constrain the presence of planet-size trojans in the Lagrangian points of nine known exoplanets. We find no clear evidence of trojans in these nine systems through any of the techniques used down to the precision of the observations. However, this allows us to constrain the presence of any potential trojan in the system, specially in the trojan mass/radius versus libration amplitude plane. In particular, we can set upper mass limits in the super-Earth mass regime for six of the studied systems.
[21]  oai:arXiv.org:1806.08305  [pdf] - 1703383
kima: Exoplanet detection in radial velocities
Comments: 3 pages; published in the Journal of Open Source Software; all code available at https://github.com/j-faria/kima
Submitted: 2018-06-21
The radial-velocity (RV) method is one of the most successful in the detection of exoplanets, but is hindered by the intrinsic RV variations of the star, which can easily mimic or hide true planetary signals. kima is a package for the detection and characterization of exoplanets using RV data. It fits a sum of Keplerian curves to a timeseries of RV measurements and calculates the evidence for models with a fixed number Np of Keplerian signals, or after marginalising over Np. Moreover, kima can use a GP with a quasi-periodic kernel as a noise model, to deal with activity-induced signals. The hyperparameters of the GP are inferred together with the orbital parameters. The code is written in C++, but includes a helper Python package, pykima, which facilitates the analysis of the results.
[22]  oai:arXiv.org:1806.06483  [pdf] - 1729054
Chemical abundances of neutron capture elements in exoplanet-hosting stars
Comments: Accepted for publication in PASP
Submitted: 2018-06-17
To understand the formation and composition of planetary systems it is important to study their host stars composition since both are formed in the same stellar nebula. In this work we analyze the behaviour of chemical abundances of Cu, Zn, Sr, Y, Zr, Ba, Ce, Nd and Eu in the large and homogeneous HARPS-GTO planet search sample ($R \sim$ 115000). This sample is composed of 120 stars hosting high-mass planets, 29 stars hosting exclusively Neptunians and Super-Earths and 910 stars without detected giant planets. We compare the [X/Fe] ratios of such elements in different metallicity bins and we find that planet hosts present higher abundances of Zn for [Fe/H]$<$--0.1 dex. On the other hand, Ba, Sr, Ce and Zr abundances are underabundant in stars with planets, with a bigger difference for stars only hosting low-mass planets. However, most of the offsets found can be explained by differences in stellar parameters and by the fact that planet hosts at low metallicity mostly belong to the Galactic thick disk. Only in the case of Ba we find a statistically significant (3$\sigma$) underabundance of 0.03 dex for low-mass planet hosts. The origin of these elements is quite complex due to their evolution during the history of the Galaxy. Therefore, it is necessary to understand and characterize the stellar populations to which planet hosts belong in order to do a fair comparison with stars without detected planets. This work demonstrates that the effects of Galactic chemical evolution and not the presence of planets mostly account for the differences we find.
[23]  oai:arXiv.org:1805.08405  [pdf] - 1686983
An Earth-sized exoplanet with a Mercury-like composition
Comments: Accepted preprint in Nature Astronomy. Publisher-edited version available at http://rdcu.be/JRE7 Supplement materials available at https://www.nature.com/articles/s41550-018-0420-5
Submitted: 2018-05-22
The Earth, Venus, Mars, and some extrasolar terrestrial planets have a mass and radius that is consistent with a mass fraction of about 30% metallic core and 70% silicate mantle. At the inner frontier of the solar system, Mercury has a completely different composition, with a mass fraction of about 70% metallic core and 30% silicate mantle. Several formation or evolution scenarios are proposed to explain this metal-rich composition, such as a giant impact, mantle evaporation, or the depletion of silicate at the inner-edge of the proto-planetary disk. These scenarios are still strongly debated. Here we report the discovery of a multiple transiting planetary system (K2-229), in which the inner planet has a radius of 1.165+/-0.066 Rearth and a mass of 2.59+/-0.43 Mearth. This Earth-sized planet thus has a core-mass fraction that is compatible with that of Mercury, while it was expected to be similar to that of the Earth based on host-star chemistry. This larger Mercury analogue either formed with a very peculiar composition or it has evolved since, e.g. by losing part of its mantle. Further characterisation of Mercury-like exoplanets like K2-229 b will help putting the detailed in-situ observations of Mercury (with Messenger and BepiColombo) into the global context of the formation and evolution of solar and extrasolar terrestrial planets.
[24]  oai:arXiv.org:1712.01805  [pdf] - 1656212
Distinguishing the albedo of exoplanets from stellar activity
Comments: Accepted for publication in A&A, 14 pages, 12 figures
Submitted: 2017-12-05
Light curves show the flux variation from the target star and its orbiting planets as a function of time. In addition to the transit features created by the planets, the flux also includes the reflected light component of each planet, which depends on the planetary albedo. This signal is typically referred to as phase curve and could be easily identified if there were no additional noise. As well as instrumental noise, stellar activity, such as spots, can create a modulation in the data, which may be very difficult to distinguish from the planetary signal. We analyze the limitations imposed by the stellar activity on the detection of the planetary albedo, considering the limitations imposed by the predicted level of instrumental noise and the short duration of the observations planned in the context of the CHEOPS mission. As initial condition, we have assumed that each star is characterized by just one orbiting planet. We built mock light curves that included a realistic stellar activity pattern, the reflected light component of the planet and an instrumental noise level, which we have chosen to be at the same level as predicted for CHEOPS. We then fit these light curves to try to recover the reflected light component, assuming the activity patterns can be modeled with a Gaussian process.We estimate that at least one full stellar rotation is necessary to obtain a reliable detection of the planetary albedo. This result is independent of the level of noise, but it depends on the limitation of the Gaussian process to describe the stellar activity when the light curve time-span is shorter than the stellar rotation. Finally, in presence of typical CHEOPS gaps in the simulations, we confirm that it is still possible to obtain a reliable albedo.
[25]  oai:arXiv.org:1711.00777  [pdf] - 1602759
Constraining planet structure and composition from stellar chemistry: trends in different stellar populations
Comments: Accepted for publication in Astronomy & Astrophysics
Submitted: 2017-11-02
The chemical composition of stars that have orbiting planets provides important clues about the frequency, architecture, and composition of exoplanet systems. We explore the possibility that stars from different galactic populations that have different intrinsic abundance ratios may produce planets with a different overall composition. We compiled abundances for Fe, O, C, Mg, and Si in a large sample of solar neighbourhood stars that belong to different galactic populations. We then used a simple stoichiometric model to predict the expected iron-to-silicate mass fraction and water mass fraction of the planet building blocks, as well as the summed mass percentage of all heavy elements in the disc. Assuming that overall the chemical composition of the planet building blocks will be reflected in the composition of the formed planets, we show that according to our model, discs around stars from different galactic populations, as well as around stars from different regions in the Galaxy, are expected to form rocky planets with significantly different iron-to-silicate mass fractions. The available water mass fraction also changes significantly from one galactic population to another. The results may be used to set constraints for models of planet formation and chemical composition. Furthermore, the results may have impact on our understanding of the frequency of planets in the Galaxy, as well as on the existence of conditions for habitability.
[26]  oai:arXiv.org:1710.01138  [pdf] - 1622283
The TROY project: Searching for co-orbital bodies to known planets. I. Project goals and first results from archival radial velocity
Comments: Accepted for publication in A&A, 18 pages, 7 figures, 6 tables. Project website: http://www.sc.eso.org/~jlillobo/troy/
Submitted: 2017-10-03
The detection of Earth-like planets, exocomets or Kuiper belts show that the different components found in the solar system should also be present in other planetary systems. Trojans are one of these components and can be considered fossils of the first stages in the life of planetary systems. Their detection in extrasolar systems would open a new scientific window to investigate formation and migration processes. In this context, the main goal of the TROY project is to detect exotrojans for the first time and to measure their occurrence rate (eta-Trojan). In this first paper, we describe the goals and methodology of the project. Additionally, we used archival radial velocity data of 46 planetary systems to place upper limits on the mass of possible trojans and investigate the presence of co-orbital planets down to several tens of Earth masses. We used archival radial velocity data of 46 close-in (P<5 days) transiting planets (without detected companions) with information from high-precision radial velocity instruments. We took advantage of the time of mid-transit and secondary eclipses (when available) to constrain the possible presence of additional objects co-orbiting the star along with the planet. This, together with a good phase coverage, breaks the degeneracy between a trojan planet signature and signals coming from additional planets or underestimated eccentricity. We identify nine systems for which the archival data provide 1-sigma evidence for a mass imbalance between L4 and L5. Two of these systems provide 2-sigma detection, but no significant detection is found among our sample. We also report upper limits to the masses at L4/L5 in all studied systems and discuss the results in the context of previous findings.
[27]  oai:arXiv.org:1709.00865  [pdf] - 1604890
Precise masses for the transiting planetary system HD 106315 with HARPS
Comments: 16 pages, 9 figures, A\&A accepted
Submitted: 2017-09-04
The multi-planetary system HD 106315 was recently found in K2 data . The planets have periods of $P_b \sim9.55$ and $P_c \sim 21.06\,$days, and radii of $ r_b = 2.44 \pm 0.17\, $ and $r_c = 4.35 \pm 0.23\, $ $R_{\oplus}$. The brightness of the host star (V=9.0 mag) makes it an excellent target for transmission spectroscopy. However, to interpret transmission spectra it is crucial to measure the planetary masses. We obtained high precision radial velocities for HD~106315 to determine the mass of the two transiting planets discovered with Kepler K2. Our successful observation strategy was carefully tailored to mitigate the effect of stellar variability. We modelled the new radial velocity data together with the K2 transit photometry and a new ground-based partial transit of HD 106315c to derive system parameters. We estimate the mass of HD 106315b to be 12.6 $\pm$ 3.2 $M_{\oplus}$ and the density to be $4.7 \pm 1.7\, g\,cm^{-3}$, while for HD 106315c we estimate a mass of 15.2 $\pm$ 3.7 $M_{\oplus}$ and a density of $1.01 \pm 0.29\, $g\,cm$^{-3}$. Hence, despite planet c having a radius almost twice as large as planet b, their masses are consistent with one another. We conclude that HD 106315c has a thick hydrogen-helium gaseous envelope. A detailed investigation of HD 106315b using a planetary interior model constrains the core mass fraction to be 5-29\%, and the water mass fraction to be 10-50\%. An alternative, not considered by our model, is that HD 106315b is composed of a large rocky core with a thick H-He envelope. Transmission spectroscopy of these planets will give insight into their atmospheric compositions and also help constrain their core compositions.
[28]  oai:arXiv.org:1707.01827  [pdf] - 1585592
Understanding stellar activity-induced radial velocity jitter using simultaneous K2 photometry and HARPS RV measurements
Comments: 13 pages, 7 figures, 2 tables, accepted for publication in A&A
Submitted: 2017-07-06
One of the best ways to improve our understanding of the stellar activity-induced signal in radial velocity (RV) measurements is through simultaneous high-precision photometric and RV observations. This is of prime importance to mitigate the RV signal induced by stellar activity and therefore unveil the presence of low-mass exoplanets. The K2 Campaign 7 and 8 field-of-views were located in the southern hemisphere, and provided a unique opportunity to gather unprecedented simultaneous high precision photometric observation with K2 and high-precision RV measurements with the HARPS spectrograph to study the relationship between photometric variability and RV jitter. We observed nine stars with different levels of activity; from quiet to very active. We probe the presence of any meaningful relation between measured RV jitter and the simultaneous photometric variation, and also other activity indicators (e.g. BIS, FWHM, $logR'_{HK}$, and F8), by evaluating the strength and significance of the correlation between RVs and each indicator. We found that for the case of very active stars, strong and significant correlations exist between almost all the observables and measured RVs; however, for lower activity levels the correlations become random. Except for the F8 which its strong correlation with RV jitter persists over a wide range of stellar activity level, and thus our result suggests that F8 might be a powerful proxy for activity induced RV jitter. Moreover, we examine the capability of two state-of-the-art modeling techniques, namely the FF' method and SOAP2.0, in accurately predicting the RV jitter amplitude using the simultaneous photometric observation. We found that for the very active stars both techniques can reasonably well predict the amplitude of the RV jitter, however, at lower activity levels the FF' method underpredicts the RV jitter amplitude.
[29]  oai:arXiv.org:1705.06090  [pdf] - 1583525
Observational evidence for two distinct giant planet populations
Comments: A&A, accepted
Submitted: 2017-05-17
Analysis of the statistical properties of exoplanets, together with those of their host stars, are providing a unique view into the process of planet formation and evolution. In this paper we explore the properties of the mass distribution of giant planet companions to solar-type stars, in a quest for clues about their formation process. With this goal in mind we studied, with the help of standard statistical tests, the mass distribution of giant planets using data from the exoplanet.eu catalog and the SWEET-Cat database of stellar parameters for stars with planets. We show that the mass distribution of giant planet companions is likely to present more than one population with a change in regime around 4\,M$_{Jup}$. Above this value host stars tend to be more metal poor and more massive and have [Fe/H] distributions that are statistically similar to those observed in field stars of similar mass. On the other hand, stars that host planets below this limit show the well-known metallicity-giant planet frequency correlation. We discuss these results in light of various planet formation models and explore the implications they may have on our understanding of the formation of giant planets. In particular, we discuss the possibility that the existence of two separate populations of giant planets indicates that two different processes of formation are at play.
[30]  oai:arXiv.org:1703.04828  [pdf] - 1581927
SIGS - Seismic Inferences for Glitches in Stars
Comments: 4 pages; 5 figures; Submitted to the Proceedings of the joint TASC2/KASC9 workshop - SPACEINN & HELAS8 conference 2016
Submitted: 2017-03-14
The increased amount of high precision seismic data for solar-like stars calls for the existence of tools that can extract information from such data. In the case of the study of acoustic glitches there are no publicly available tools and most existing ones require a deep knowledge of their implementation. In this work a tool is presented that aims to both simplify the interaction with the user and also be capable of working automatically to determine properties of acoustic glitches from seismic data of solar-like stars. This tool is shown to work with both the Sun and other solar analogs but also shows that are still severe limitations to the methods used, when considering smaller datasets.
[31]  oai:arXiv.org:1606.05549  [pdf] - 1426454
Is the activity level of HD 80606 influenced by its eccentric planet?
Comments: 9 pages, 4 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2016-06-17, last modified: 2016-06-21
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$_\alpha$, 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)
[32]  oai:arXiv.org:1606.04714  [pdf] - 1457408
Abundance trend with condensation temperature for stars with different Galactic birth places
Comments: Accepted by A&A
Submitted: 2016-06-15
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
[33]  oai:arXiv.org:1605.06942  [pdf] - 1433072
An extreme planetary system around HD219828. One long-period super Jupiter to a hot-neptune host star
Comments: Astronomy & Astrophysics, in press
Submitted: 2016-05-23
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.
[34]  oai:arXiv.org:1605.01918  [pdf] - 1418896
$\zeta^2$ Ret, its debris disk, and its lonely stellar companion $\zeta^1$ Ret. Different $T_{\mathrm{c}}$ trends for different spectra
Comments: Accepted by A&A
Submitted: 2016-05-06
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 $\zeta$ Reticuli binary system, which consists of two solar analogs. The observed trend in $\zeta^2$ Ret relative to its companion was explained by the presence of a debris disk around $\zeta^2$ Ret. Our goal is to re-evaluate the presence and variability of the Tc trend in the $\zeta$ 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 $\zeta^2$ Ret and $\zeta^1$ Ret shows a significant ($\sim$ 2 $\sigma$) 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 $\zeta$ 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.
[35]  oai:arXiv.org:1601.07495  [pdf] - 1374310
Uncovering the planets and stellar activity of CoRoT-7 using only radial velocities
Comments: 7 pages, 5 figures, 2 tables, accepted for publication in A&A, language editing revisions
Submitted: 2016-01-27, last modified: 2016-02-03
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.
[36]  oai:arXiv.org:1601.07535  [pdf] - 1388907
The HARPS search for southern extra-solar planets XL. Searching for Neptunes around metal-poor stars
Comments: 18 pages, 29 figures, 4 tables, accepted for publication in A&A
Submitted: 2016-01-27
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.
[37]  oai:arXiv.org:1601.05107  [pdf] - 1415058
A pragmatic Bayesian perspective on correlation analysis: The exoplanetary gravity - stellar activity case
Comments: 10 pages, 1 figure, accepted for publication in the peer-reviewed proceedings of the conference "Habitability in the Universe: From the Early Earth to Exoplanets"
Submitted: 2016-01-19
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, $\rho$, 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.
[38]  oai:arXiv.org:1511.03941  [pdf] - 1339020
The HARPS search for southern extra-solar planets. XXXIX. HD175607 b, the most metal-poor G dwarf with an orbiting sub-Neptune
Comments: 12 pages, 9 figures, accepted for publication in Astronomy and Astrophysics
Submitted: 2015-11-12
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.
[39]  oai:arXiv.org:1509.02419  [pdf] - 1304098
Identifying the best iron-peak and $\alpha$-capture elements for chemical tagging: The impact of the number of lines on measured scatter
Comments: Accepted by A&A
Submitted: 2015-09-08
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 $\alpha$-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.
[40]  oai:arXiv.org:1504.05962  [pdf] - 984011
Evidence for a spectroscopic direct detection of reflected light from 51 Peg b
Comments: 9 pages, 6 figures
Submitted: 2015-04-22
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\sigma_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.
[41]  oai:arXiv.org:1412.0467  [pdf] - 920448
BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram
Comments: 6 pages, 6 figures, accepted for publication in A&A, code online at https://www.astro.up.pt/exoearths/tools.html
Submitted: 2014-12-01
Context. Frequency analyses are very important in astronomy today, not least in the ever-growing field of exoplanets, where short-period signals in stellar radial velocity data are investigated. Periodograms are the main (and powerful) tools for this purpose. However, recovering the correct frequencies and assessing the probability of each frequency is not straightforward. Aims. We provide a formalism that is easy to implement in a code, to describe a Bayesian periodogram that includes weights and a constant offset in the data. The relative probability between peaks can be easily calculated with this formalism. We discuss the differences and agreements between the various periodogram formalisms with simulated examples. Methods. We used the Bayesian probability theory to describe the probability that a full sine function (including weights derived from the errors on the data values and a constant offset) with a specific frequency is present in the data. Results. From the expression for our Baysian generalised Lomb-Scargle periodogram (BGLS), we can easily recover the expression for the non-Bayesian version. In the simulated examples we show that this new formalism recovers the underlying periods better than previous versions. A Python-based code is available for the community.
[42]  oai:arXiv.org:1409.0890  [pdf] - 881958
Exoplanet hosts reveal lithium depletion: Results from a homogeneous statistical analysis
Comments: 6 pages, 2 figures, accepted for publication in A&A
Submitted: 2014-09-02
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 {\sigma} 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)
[43]  oai:arXiv.org:1405.7512  [pdf] - 862854
Asteroseismic estimate of helium abundance of a solar analog binary system
Comments: Accepted for publication in ApJ
Submitted: 2014-05-29
16 Cyg A and B are among the brightest stars observed by Kepler. What makes these stars more interesting is that they are solar analogs. 16 Cyg A and B exhibit solar-like oscillations. In this work we use oscillation frequencies obtained using 2.5 years of Kepler data to determine the current helium abundance of these stars. For this we use the fact that the helium ionization zone leaves a signature on the oscillation frequencies and that this signature can be calibrated to determine the helium abundance of that layer. By calibrating the signature of the helium ionization zone against models of known helium abundance, the helium abundance in the envelope of 16 Cyg A is found to lie in the range 0.231 to 0.251 and that of 16 Cyg B lies in the range 0.218 to 0.266.
[44]  oai:arXiv.org:1404.6135  [pdf] - 844860
The HARPS search for southern extra-solar planets XXXV. The interesting case of HD41248: stellar activity, no planets?
Comments: Accepted for publication in A&A
Submitted: 2014-04-24
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.
[45]  oai:arXiv.org:1211.6000  [pdf] - 1158039
On the possibility of using seismic probes to study the core composition in pulsating white dwarfs
Comments: 4 pages, 4 figures, The Modern Era of Helio- and Asteroseismology, to be published in Astron. Nachr. 333, No. 10, 954 - 957 (2012)
Submitted: 2012-11-26
White dwarfs correspond to the final stages of stellar evolution of solar-type stars. In these objects, production of energy by nuclear burning has ended which means that a white dwarf simply cools down over the course of the next billion years. It is now known that white dwarfs spend some of their cooling history in an instability strip. The pulsating white dwarfs with an hydrogen atmosphere (called DAV or ZZ Ceti stars) show non-radial oscillation modes with periods in the range 100 - 1200s. In this work we try to illustrate how the oscillation p-mode frequencies of idealized white dwarf models change as the result of a different chemical composition in the core, with the ultimate goal of determining the chemical stratification from seismic observations. The presence of acoustic glitches in the internal structure results in a periodic signal in the frequencies. We find that this signal depends on the chemical stratification/composition of the core in a form that can be analytically modelled.