Source author record

J. P. Faria

J. P. Faria 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

28works
4topics
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

28 published item(s)

preprint2026arXiv

ASTRA: A Python Package for Cross-Instrument Stellar and Telluric Template Construction

ASTRA is a Python package that provides a modular, instrument-independent interface for working with high-resolution stellar spectra. Designed to support data from multiple spectrographs, including ESPRESSO (Pepe et al., 2021), HARPS (Mayor et al., 2003; Pepe et al., 2002), MAROON-X (Seifahrt et al., 2022), and CARMENES (Quirrenbach et al., 2014). ASTRA offers a unified abstraction over their data formats, enabling consistent access to fluxes, wavelengths, uncertainties, and metadata across instruments. Furthermore, it applies the necessary wavelength and flux calibrations that are needed, as described by the official pipelines of each instrument.

preprint2025arXiv

Inferring planet occurrence rates from radial velocities

We introduce a new method to infer the posterior distribution for planet occurrence rates from radial-velocity (RV) observations. The approach combines posterior samples from the analysis of individual RV datasets of several stars, using importance sampling to reweight them appropriately. This eliminates the need for injection-recovery tests to compute detection limits and avoids the explicit definition of a detection threshold. We validate the method on simulated RV datasets and show that it yields unbiased estimates of the occurrence rate in different regions, with increasing precision as more stars are included in the analysis.

preprint2022arXiv

$\tt{KOBEsim}$: a Bayesian observing strategy algorithm for planet detection in radial velocity blind-search surveys

Ground-based observing time is precious in the era of exoplanet follow-up and characterization, especially in high-precision radial velocity instruments. Blind-search radial velocity surveys thus require a dedicated observational strategy in order to optimize the observing time, which is particularly crucial for the detection of small rocky worlds at large orbital periods. We develop an algorithm with the purpose of improving the efficiency of radial velocity observations in the context of exoplanet searches, and we apply it to the K-dwarfs Orbited By habitable Exoplanets (KOBE) experiment. We aim at accelerating exoplanet confirmations or, alternatively, rejecting false signals as early as possible in order to save telescope time and increase the efficiency of both blind-search surveys and follow-up of transiting candidates. Once a minimum initial number of radial velocity datapoints is reached in such a way that a periodicity starts to emerge according to generalized Lomb-Scargle (GLS) periodograms, that period is targeted with the proposed algorithm, named $\texttt{KOBEsim}$. The algorithm selects the next observing date that maximizes the Bayesian evidence for such periodicity in comparison with a model with no Keplerian orbits. By means of simulated data, we prove that the algorithm accelerates the exoplanet detection, needing $29 - 33\,\%$ less observations and $41 - 47\,\%$ less timespan of the full dataset for low-mass planets ($m_{\rm p}\,<\,10\,M_{\oplus}$) in comparison with a conventional monotonic cadence strategy. The enhancement in the number of datapoints for $20\,M_{\oplus}$ planets is also appreciable, $16\,\%$. We also test $\texttt{KOBEsim}$ with real data for a particular KOBE target, and for the confirmed planet $HD~102365\,b$. Both of them demonstrate that the strategy is capable of speeding up the detection up to a factor of $2$.

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

Modelling stellar activity with Gaussian process regression networks

Stellar photospheric activity is known to limit the detection and characterisation of extra-solar planets. In particular, the study of Earth-like planets around Sun-like stars requires data analysis methods that can accurately model the stellar activity phenomena affecting radial velocity (RV) measurements. Gaussian Process Regression Networks (GPRNs) offer a principled approach to the analysis of simultaneous time-series, combining the structural properties of Bayesian neural networks with the non-parametric flexibility of Gaussian Processes. Using HARPS-N solar spectroscopic observations encompassing three years, we demonstrate that this framework is capable of jointly modelling RV data and traditional stellar activity indicators. Although we consider only the simplest GPRN configuration, we are able to describe the behaviour of solar RV data at least as accurately as previously published methods. We confirm the correlation between the RV and stellar activity time series reaches a maximum at separations of a few days, and find evidence of non-stationary behaviour in the time series, associated with an approaching solar activity minimum.

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

Can planetary rings explain the extremely low density of HIP 41378f?

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.

preprint2020arXiv

Can we detect the stellar differential rotation of WASP-7 through theRossiter-McLaughlin observations?

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°, 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.

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

Improving transit characterisation with Gaussian process modelling of stellar variability

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.

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.

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

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

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

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

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

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.

preprint2014arXiv

BGLS: A Bayesian formalism for the generalised Lomb-Scargle periodogram

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.

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

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.

preprint2012arXiv

On the possibility of using seismic probes to study the core composition in pulsating white dwarfs

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.