Researcher profile

J. P. Faria

J. P. Faria contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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&#39;s activity using a Gaussian process (GP). We use the CCF&#39;s full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit.

preprint2022arXiv

A novel framework for semi-Bayesian radial velocities through template matching

The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alternative for M-dwarf stars. In this paper, we describe a new implementation of template matching within a semi-Bayesian framework, providing a more statistically principled characterization of the RV measurements. In this context, a common RV shift is used to describe the difference between each spectral order of a given stellar spectrum and a template built from the available observations. Posterior probability distributions are obtained for the relative RV associated with each spectrum, after marginalizing with respect to the continuum. This methodology was named S-BART: Semi-Bayesian Approach for RVs with Template-matching, and it can be applied to HARPS and ESPRESSO. The application of our method to HARPS archival observations of Barnard&#39;s star allowed us to validate our implementation against HARPS-TERRA and SERVAL. Then, we applied it to 33 ESPRESSO targets, evaluating its performance and comparing it with the CCF method. We found a decrease in the median RV scatter of \sim 10\% and \sim 4\% for M- and K-type stars, respectively. S-BART yields more precise RV estimates than the CCF method, particularly in the case of M-type stars where a median uncertainty of \sim 15 cm/s is achieved over 309 observations. Further, we estimated the nightly zero point (NZP) of ESPRESSO, finding a weighted NZP scatter below \sim 0.7 m/s. As this includes stellar variability, photon noise, and potential planetary signals, it should be taken as an upper limit of the RV precision attainable with ESPRESSO data.

preprint2022arXiv

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&#39;s Very-Large Telescope (VLT). It was designed for ultra-high radial-velocity precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UT) of the VLT at a spectral resolving power of 140,000 or 190,000 over the 378.2 to 788.7 nm wavelength range, or with all UTs together, turning the VLT into a 16-m diameter equivalent telescope in terms of collecting area, while still providing a resolving power of 70,000. We provide a general description of the ESPRESSO instrument, report on the actual on-sky performance, and present our Guaranteed-Time Observation (GTO) program with its first results. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1st, 2018, but improvements to the instrument and re-commissioning runs were conducted until July 2019. The measured overall optical throughput of ESPRESSO at 550 nm and a seeing of 0.65 arcsec exceeds the 10% mark under nominal astro-climatic conditions. We demonstrate a radial-velocity precision of better than 25 cm/s during one night and 50 cm/s over several months. These values being limited by photon noise and stellar jitter show that the performanceis compatible with an instrumental precision of 10 cm/s. No difference has been measured across the UTs neither in throughput nor RV precision. The combination of the large collecting telescope area with the efficiency and the exquisite spectral fidelity of ESPRESSO opens a new parameter space in RV measurements, the study of planetary atmospheres, fundamental constants, stellar characterisation and many other fields.

preprint2020arXiv

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.