High Accuracy Radial Velocity Planet Searcher Research Articles

Precise characterisation of HD 15337 with CHEOPS: A laboratory for planet formation and evolution

The HD 15337 (TIC 120896927, TOI-402) system was observed by the Transiting Exoplanet Survey Satellite (TESS), revealing the presence of two short-period planets situated on opposite sides of the radius gap. This offers an excellent opportunity to study theories of formation and evolution, as well as to investigate internal composition and atmospheric evaporation. We aim to constrain the internal structure and composition of two short-period planets situated on opposite sides of the radius valley:\ HD 15337 b and c. We use new transit photometry and radial velocity data. We acquired 6 new transit visits with the CHaracterising ExOPlanet Satellite (CHEOPS) and 32 new radial velocity measurements from the High Accuracy Radial Velocity Planet Searcher (HARPS) to improve the accuracy of the mass and radius estimates for both planets. We re-analysed the light curves from TESS sectors 3 and 4 and analysed new data from sector 30, correcting for long-term stellar activity. Subsequently, we performed a joint fit of the TESS and CHEOPS light curves, along with all available RV data from HARPS and the Planet Finder Spectrograph (PFS). Our model fit the planetary signals, stellar activity signal, and instrumental decorrelation model for the CHEOPS data simultaneously. The stellar activity was modelled using a Gaussian-process regression on both the RV and activity indicators. Finally, we employed a Bayesian retrieval code to determine the internal composition and structure of the planets. We derived updated and highly precise parameters for the HD 15337 system. Our improved precision on the planetary parameters makes HD 15337 b one of the most precisely characterised rocky exoplanets, with radius and mass measurements achieving a precision better than 2<!PCT!> and 7<!PCT!>, respectively. We were able to improve the precision of the radius measurement of HD 15337 c to 3<!PCT!>. Our results imply that the composition of HD 15337 b is predominantly rocky, while HD 15337 c exhibits a gas envelope with a mass of at least $0.01\ M_ oplus$. Our results lay the groundwork for future studies, which can further unravel the atmospheric evolution of these exoplanets and offer new insights into their composition and formation history as well as the causes behind the radius gap.

A.C.I.D – an improved LSD technique for accurate line profile retrieval

ABSTRACT Stellar activity and planetary effects induce radial velocity (RV) offsets and cause temporal distortions in the shape of the stellar line profile. Hence, accurately probing the stellar line profile offers a wealth of information on both the star itself and any orbiting planets. Typically, cross-correlation functions (CCFs) are used as a proxy for the stellar line profile. The shape of CCFs, however, can be distorted by line blending and aliasing limiting the stellar and planetary physics that can be probed from them. Least-squares deconvolution (LSD) offers an alternative that directly fits the mean line profile of the spectrum to produce a high-precision profile. In this paper, we introduce our novel method ACID (Accurate Continuum fItting and Deconvolution) that builds on LSD techniques by simultaneously fitting the spectral continuum and line profile as well as performing LSD in effective optical depth. Tests on model data revealed ACID can accurately identify and correct the spectral continuum to retrieve an injected line profile. ACID was also applied to archival High Accuracy Radial-velocity Planet Searcher (HARPS) data obtained during the transit of HD189733b. The application of the Reloaded Rossiter–McLaughlin technique to both ACID profiles and HARPS CCFs shows ACID residual profiles improved the out-of-line root mean square (RMS) by over 5 per cent compared to CCFs. Furthermore, ACID profiles are shown to exhibit a Voigt profile shape that better describes the expected profile shape of the stellar line profile. This improved representation shows that ACID better preserves the stellar and planetary physics encoded in the stellar line profile shape for slow rotating stars.

TOI-2498 b: a hot bloated super-Neptune within the Neptune desert

ABSTRACT We present the discovery and confirmation of a transiting hot bloated super-Neptune using photometry from the Transiting Exoplanet Survey Satellite (TESS) and the Las Cumbres Observatory Global Telescope (LCOGT) and radial velocity measurements from the High Accuracy Radial velocity Planet Searcher (HARPS). The host star TOI-2498 is a V = 11.2, G-type (Teff = 5905 ± 12 K) solar-like star with a mass of 1.12 ± 0.02 M⊙ and a radius of 1.26 ± 0.04 R⊙. The planet, TOI-2498 b, orbits the star with a period of 3.7 d, has a radius of 6.1 ± 0.3 R⊕, and a mass of 35 ± 4 M⊕. This results in a density of 0.86 ± 0.25 g cm−3. TOI-2498 b resides on the edge of the Neptune desert; a region of mass–period parameter space in which there appears to be a dearth of planets. Therefore TOI-2498 b is an interesting case to study to further understand the origins and boundaries of the Neptune desert. Through modelling the evaporation history, we determine that over its ∼3.6 Gyr lifespan, TOI-2498 b has likely reduced from a Saturn-sized planet to its current radius through photoevaporation. Moreover, TOI-2498 b is a potential candidate for future atmospheric studies searching for species like water or sodium in the optical using high resolution spectroscopy, and for carbon-based molecules in the infrared using JWST.

Photometric variations from stellar activity as an age indicator for solar-twins

ABSTRACT Stellar ages are elusive to measure, albeit being very important for understanding stellar evolution. We investigate the impact of photospheric activity on 2-min cadence light curves from the TESS/NASA mission of a selected sample of 30 solar-twins with well-determined ages. The photometric variability, $\mathcal {A}_{ {{TESS}}}$, of the light curves due to rotational modulations by the presence of active regions was estimated and correlated with chromospheric activity (Ca ii H&K lines from an extensive High Accuracy Radial velocity Planet Searcher (HARPS) at the European Southern Observatory (ESO) HARPS/ESO activity time series) and ages. Moreover, these results were compared with the total solar irradiance amplitude behaviour during the solar magnetic cycles 23 and 24, validating our findings for solar-twins. Our results show the photometric amplitude to be strongly correlated to the average level of chromospheric activity for the star sample. Also, we found a good correlation of $\mathcal {A}_{ {{TESS}}}$ with stellar age (in Gyr) described by $\log t = +12.239-0.894 \log \mathcal {A}_{ {{TESS}}}$. In conclusion, stellar photometric variability $\mathcal {A}_{ {{TESS}}}$ may be used as a simple age diagnostic for solar-twins.

High-mass eclipsing binaries: A testbed for models of interior structure and evolution

Aims. The surface chemical compositions of stars are affected by physical processes that bring the products of thermonuclear burning to the surface. Despite their potential in helping us understand the structure and evolution of stars, elemental abundances are available for only a few high-mass binary stars. We aim to enlarge this sample by determining the physical properties and photospheric abundances for four eclipsing binary systems that contain high-mass stars: V1034 Sco, GL Car, V573 Car, and V346 Cen. The components have masses of 8–17 M⊙, have effective temperatures from 22 500 to 32 200 K, and are all on the main sequence. Methods. We present new high-resolution and high signal-to-noise spectroscopy from the High Accuracy Radial velocity Planet Searcher (HARPS), which we analysed using spectral disentangling and non-local thermodynamic equilibrium spectral synthesis. We modelled existing light curves and new photometry from the Transiting Exoplanet Survey Satellite (TESS). Results. We measure the stellar masses to a 0.6–2.0% precision, radii to a 0.8–1.7% precision, effective temperatures to a 1.1–1.6% precision, and abundances of C, N, O, Mg, and Si. The abundances are similar to those found in our previous studies of high-mass eclipsing binaries; our sample now comprises 25 high-mass stars in 13 binary systems. We also find tidally excited pulsations in V346 Cen. Conclusions. These results reinforce our previous conclusions: interior chemical element transport is not as efficient in binary star components as in their single-star counterparts in the same mass regime and evolutionary stage, possibly due to the effects of tidal forces. Our ultimate goal is to provide a larger sample of OB-type stars in binaries to enable a thorough comparison to stellar evolutionary models, as well as to single high-mass stars.

VaTEST I: validation of sub-Saturn exoplanet TOI-181b in narrow orbit from its host star

ABSTRACT We present here a validation of sub-Saturn exoplanet TOI-181b orbiting a K spectral type star TOI-181 (mass: 0.822 ± 0.04 M⊙, radius: 0.745 ± 0.02 R⊙, temperature: 4994 ± 50 K) as a part of Validation of Transiting Exoplanets using Statistical Tools (VaTEST) project. TOI-181b is a planet with radius 6.95 ± 0.08 R⊕, mass 46.16 ± 2.71 M⊕, orbiting in a slightly eccentric orbit with eccentricity 0.15 ± 0.06 and semimajor axis of 0.054 ± 0.004 au, with an orbital period of 4.5320 ± 0.000002 d. The transit photometry data were collected using Transiting Exoplanet Survey Satellite (TESS) and spectroscopic data for radial velocity analysis were collected using The European Southern Observatory’s (ESO) High Accuracy Radial Velocity Planet Searcher (HARPS) telescope. Based on the radial velocity best-fit model we measured RV semi-amplitude to be 20.56 ± 2.37 m s−1. Additionally, we used VESPA and TRICERATOPS to compute the False Positive Probability (FPP), and the findings were FPP values of 1.68 × 10−14 and 3.81 × 10−04, respectively, which are significantly lower than the 1 per cent threshold. The finding of TOI-181b is significant in the perspective of future work on the formation and migration history of analogous planetary systems since warm sub-Saturns are uncommon in the known sample of exoplanets.

HARPS radial velocity search for planets in the Scorpius-Centaurus association

Context. The Scorpius-Centaurus (Sco-Cen) young and nearby massive star-forming region is particularly well suited for extrasolar planet searches with both direct imaging and radial velocity (RV) techniques. The RV search, however, is challenging, as the stars are faster rotators on average than their older stellar counterparts of similar spectral types. Moreover, the RV time series show strong signatures of stellar variability (spots and faculae) and/or stellar pulsations. Aims. Our aim is to search for giant planets (GPs) and brown dwarfs at short orbital distances around star members of the Sco-Cen association. We also aim at using these data together with others available on young stars to estimate the GP occurrence rate for young stars for periods of up to 1000 days. Methods. We used the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph on the 3.6 m telescope at the La Silla Observatory to monitor 88 A – F Sco-Cen stars. To improve our statistics and analysis, we combined this survey with two previous surveys that focused on young nearby stars (YNS) to compute companion occurrence rates from a sample of 176 young A – M stars. Results. We report the discovery of a massive hot-Jupiter candidate around HD 145467, together with the discovery of one probable short-period (P < 10 days) brown dwarf around HD 149790. In addition, we confirm the binary nature of eight single-line binaries: HD 108857, HD 108904, HD 111102, HD 114319, HD 121176, HD 126488, HD 126838, and HD 133574. From our sample, we obtain a GP (mc ∈ [1; 13] MJup) occurrence rate of 0.7−0.2+1.6% for periods between 1 and 1000 days and a brown dwarf (mc ∈ [13; 80] MJup) occurrence rate of 0.6−0.2+1.4%, in the same period range. In addition, we report a possible lack of close (P ∈ [1; 1000] days) GPs around young F-K stars compared to their older counterparts, with a confidence level of 95%.

Evolved eclipsing binary systems in the Galactic bulge: Precise physical and orbital parameters of OGLE-BLG-ECL-305487 and OGLE-BLG-ECL-116218

Aims. Our goal is to determine, with high accuracy, the physical and orbital parameters of two double-lined eclipsing binary systems, where the components are two giant stars. We also aim to study the evolutionary status of the binaries, to derive the distances towards them by using a surface brightness–colour relation, and to compare these measurements with the measurements presented by the Gaia mission. Methods. In order to measure the physical and orbital parameters of the systems, we analysed the light curves and radial-velocity curves with the Wilson–Devinney code. We used V-band and I-band photometry from the Optical Gravitational Lensing Experiment (OGLE) catalogue and near-infrared photometry obtained with the New Technology Telescope (NTT) equipped with the SOFI instrument. The spectroscopic data were collected with the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph mounted at the ESO 3.6-m telescope and the Magellan Inamori Kyocera Echelle (MIKE) spectrograph mounted at the 6.5-m Clay telescope. Results. We present the first analysis of this kind for two evolved eclipsing binary systems from the OGLE catalogue: OGLE-BLG-ECL-305487 and OGLE-BLG-ECL-116218. The masses of the components of OGLE-BLG-ECL-305487 are M1 = 1.059 ± 0.019 and M2 = 0.991 ± 0.018 M⊙, and the radii are R1 = 19.27 ± 0.28 and R2 = 29.99 ± 0.24 R⊙. For OGLE-BLG-ECL-116218, the masses are M1 = 0.969 ± 0.012 and M2 = 0.983 ± 0.012 M⊙, while the radii are R1 = 16.73 ± 0.28 and R2 = 22.06 ± 0.26 R⊙. The evolutionary status of the systems is discussed based on the PARSEC and MIST isochrones. The ages of the systems were established to be between 7.3–10.9 Gyr for OGLE-BLG-ECL-305487 and around 10 Gyr for OGLE-BLG-ECL-116218. We also measured the distances to the binaries. For OGLE-BLG-ECL-305487, d = 7.80 ± 0.18 (stat.) ± 0.19 (syst.) kpc and for OGLE-BLG-ECL-116218, d = 7.57 ± 0.28 (stat.) ± 0.19 (syst.) kpc.

Line-by-line Velocity Measurements: an Outlier-resistant Method for Precision Velocimetry

We present a new algorithm for precision radial velocity (pRV) measurements, a line-by-line (LBL) approach designed to handle outlying spectral information in a simple but efficient manner. The effectiveness of the LBL method is demonstrated on two data sets, one obtained with SPIRou on Barnard’s star, and the other with the High Accuracy Radial velocity Planet Searcher (HARPS) on Proxima Centauri. In the near-infrared, the LBL provides a framework for meters-per-second-level accuracy in pRV measurements despite the challenges associated with telluric absorption and sky emission lines. We confirm with SPIRou measurements spanning 2.7 yr that the candidate super-Earth on a 233 day orbit around Barnard’s star is an artifact due to a combination of time sampling and activity. The LBL analysis of the Proxima Centauri HARPS post-upgrade data alone easily recovers the Proxima b signal and also provides a 2σ detection of the recently confirmed 5 day Proxima d planet, but argues against the presence of the candidate Proxima c with a period of 1900 days. We provide evidence that the Proxima c signal is associated with small, unaccounted systematic effects affecting the HARPS-TERRA template-matching radial velocity extraction method for long-period signals. Finally, the LBL framework provides a very effective activity indicator, akin to the FWHM derived from the cross-correlation function, from which we infer a rotation period of days for Proxima.

Variable and Supersonic Winds in the Atmosphere of an Ultrahot Giant Planet

Abstract Hot Jupiters (HJs) receive intense irradiation from their stellar hosts. The resulting extreme environments in their atmospheres allow us to study the conditions that drive planetary atmospheric dynamics, e.g., global-scale winds. General circulation models predict day-to-nightside winds and equatorial jets with speeds of the order of a few km s−1. To test these models, we apply high-resolution transmission spectroscopy using the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) spectrograph on the Large Binocular Telescope to study the atmosphere of KELT-9 b, an ultrahot Jupiter and currently the hottest known planet. We measure ∼10 km s−1 day-to-nightside winds traced by Fe ii features in the planet’s atmosphere. This is at odds with previous literature (including data taken with PEPSI), which report no significant day-to-nightside winds on KELT-9 b. We identify the cause of this discrepancy as due to an inaccurate ephemeris for KELT-9 b in previous literature. We update the ephemeris, which shifts the midtransit time by up to 10 minutes for previous data sets, resulting in consistent detections of blueshifts in all the data sets analyzed here. Furthermore, a comparison with archival data sets from the High-accuracy Radial velocity Planet Searcher for the Northern hemisphere suggests a temporal wind variability of ∼5–8 km s−1 over timescales between weeks to years. Temporal variability of atmospheric dynamics on HJs is a phenomenon anticipated by certain general circulation models that has not been observed over these timescales until now. However, such large variability as we measure on KELT-9 b challenges general circulation models, which predict much lower amplitudes of wind variability over timescales between days to weeks.

A pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 characterized with CHEOPS

ABSTRACT We report the discovery and characterization of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in the Transiting Exoplanet Survey Satellite (TESS) photometry. To characterize the system, we performed and retrieved the CHaracterising ExOPlanets Satellite (CHEOPS), TESS, and ground-based photometry, the High Accuracy Radial velocity Planet Searcher (HARPS) high-resolution spectroscopy, and Gemini speckle imaging. We characterize the host star and determine $T_{\rm eff, \star }=4734\pm 67\,\mathrm{ K}$, $R_{\star }=0.726\pm 0.007\, \mathrm{ R}_{\odot }$, and $M_{\star }=0.748\pm 0.032\, \mathrm{ M}_{\odot }$. We present a novel detrending method based on point spread function shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of Pb = 6.44387 ± 0.00003 d, a radius of Rb = 2.59 ± 0.04 R⊕, and a mass of $M_{\rm b} = 13.5_{-1.8}^{+1.7}$ M⊕, whilst TOI-1064 c has an orbital period of $P_{\rm c} = 12.22657^{+0.00005}_{-0.00004}$ d, a radius of Rc = 2.65 ± 0.04 R⊕, and a 3σ upper mass limit of 8.5 M⊕. From the high-precision photometry we obtain radius uncertainties of ∼1.6 per cent, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterized sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further radial velocities are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass–radius space, and it allow us to identify a trend in bulk density–stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets.

Masses and compositions of three small planets orbiting the nearby M dwarf L231-32 (TOI-270) and the M dwarf radius valley

Abstract We report on precise Doppler measurements of L231-32 (TOI-270), a nearby M dwarf (d = 22 pc, M⋆ = 0.39 M⊙, R⋆ = 0.38 R⊙), which hosts three transiting planets that were recently discovered using data from the Transiting Exoplanet Survey Satellite (TESS). The three planets are 1.2, 2.4, and 2.1 times the size of Earth and have orbital periods of 3.4, 5.7, and 11.4 days. We obtained 29 high-resolution optical spectra with the newly commissioned Echelle Spectrograph for Rocky Exoplanet and Stable Spectroscopic Observations (ESPRESSO) and 58 spectra using the High Accuracy Radial velocity Planet Searcher (HARPS). From these observations, we find the masses of the planets to be 1.58 ± 0.26, 6.15 ± 0.37, and 4.78 ± 0.43 M⊕, respectively. The combination of radius and mass measurements suggests that the innermost planet has a rocky composition similar to that of Earth, while the outer two planets have lower densities. Thus, the inner planet and the outer planets are on opposite sides of the ‘radius valley’ — a region in the radius-period diagram with relatively few members, which has been interpreted as a consequence of atmospheric photo-evaporation. We place these findings into the context of other small close-in planets orbiting M dwarf stars, and use support vector machines to determine the location and slope of the M dwarf (Teff < 4000 K) radius valley as a function of orbital period. We compare the location of the M dwarf radius valley to the radius valley observed for FGK stars, and find that its location is a good match to photo-evaporation and core-powered mass loss models. Finally, we show that planets below the M dwarf radius valley have compositions consistent with stripped rocky cores, whereas most planets above have a lower density consistent with the presence of a H-He atmosphere.

An inventory of atomic species in the atmosphere of WASP-121b using UVES high-resolution spectroscopy

ABSTRACT Ultrahot Jupiters (UHJs) present excellent targets for atmospheric characterization. Their hot dayside temperatures (T ≳ 2200 K) strongly suppress the formation of condensates, leading to clear and highly inflated atmospheres extremely conducive to transmission spectroscopy. Recent studies using optical high-resolution spectra have discovered a plethora of neutral and ionized atomic species in UHJs, placing constraints on their atmospheric structure and composition. Our recent work has presented a search for molecular features and detection of Fe i in the UHJ WASP-121b using Very Large Telescope (VLT)/UV–Visual Echelle Spectrograph (UVES) transmission spectroscopy. Here, we present a systematic search for atomic species in its atmosphere using cross-correlation methods. In a single transit, we uncover potential signals of 17 atomic species that we investigate further, categorizing five as strong detections, three as tentative detections, and nine as weak signals worthy of further exploration. We confirm previous detections of Cr i, V i, Ca i, K i, and exospheric H i and Ca ii made with the High Accuracy Radial velocity Planet Searcher (HARPS) and the Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), and independently re-recover our previous detection of Fe i at 8.8σ using both the blue and red arms of the UVES data. We also add a novel detection of Sc ii at 4.2σ. Our results further demonstrate the richness of UHJs for optical high-resolution spectroscopy.

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

ABSTRACT We report the discovery of HD 110113 b (TESS object of interest-755.01), a transiting mini-Neptune exoplanet on a 2.5-d orbit around the solar-analogue HD 110113 (Teff = 5730 K). Using TESS photometry and High Accuracy Radial velocity Planet Searcher (HARPS) radial velocities gathered by the NCORES program, we find that HD 110113 b has a radius of 2.05 ± 0.12 R⊕ and a mass of 4.55 ± 0.62 M⊕. The resulting density of $2.90^{+0.75}_{-0.59}$ g cm−3 is significantly lower than would be expected from a pure-rock world; therefore HD 110113 b must be a mini-Neptune with a significant volatile atmosphere. The high incident flux places it within the so-called radius valley; however, HD 110113 b was able to hold on to a substantial (0.1–1 per cent) H–He atmosphere over its ∼4 Gyr lifetime. Through a novel simultaneous Gaussian process fit to multiple activity indicators, we were also able to fit for the strong stellar rotation signal with period 20.8 ± 1.2 d from the RVs and confirm an additional non-transiting planet, HD 110113 c, which has a mass of 10.5 ± 1.2 M⊕ and a period of $6.744^{+0.008}_{-0.009}$ d.

Measuring and characterizing the line profile of HARPS with a laser frequency comb

Aims. We study the 2D spectral line profile of the High Accuracy Radial Velocity Planet Searcher (HARPS), measuring its variation with position across the detector and with changing line intensity. The characterization of the line profile and its variations are important for achieving the precision of the wavelength scales of 10−10 or 3.0 cm s−1 necessary to detect Earth-twins in the habitable zone around solar-like stars. Methods. We used a laser frequency comb (LFC) with unresolved and unblended lines to probe the instrument line profile. We injected the LFC light – attenuated by various neutral density filters – into both the object and the reference fibres of HARPS, and we studied the variations of the line profiles with the line intensities. We applied moment analysis to measure the line positions, widths, and skewness as well as to characterize the line profile distortions induced by the spectrograph and detectors. Based on this, we established a model to correct for point spread function distortions by tracking the beam profiles in both fibres. Results. We demonstrate that the line profile varies with the position on the detector and as a function of line intensities. This is consistent with a charge transfer inefficiency effect on the HARPS detector. The estimate of the line position depends critically on the line profile, and therefore a change in the line amplitude effectively changes the measured position of the lines, affecting the stability of the wavelength scale of the instrument. We deduce and apply the correcting functions to re-calibrate and mitigate this effect, reducing it to a level consistent with photon noise.

Occurrence rates of small planets from HARPS

Context. The stars in the Milky Way thin and thick disks can be distinguished by several properties such as metallicity and kinematics. It is not clear whether the two populations also differ in the properties of planets orbiting the stars. In order to study this, a careful analysis of both the chemical composition and mass detection limits is required for a sufficiently large sample. Currently, this information is still limited only to large radial-velocity (RV) programs. Based on the recently published archival database of the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph, we present a first analysis of low-mass (small) planet occurrence rates in a sample of thin- and thick-disk stars. Aims. We aim to assess the effects of stellar properties on planet occurrence rates and to obtain first estimates of planet occurrence rates in the thin and thick disks of the Galaxy. As a baseline for comparison, we also aim to provide an updated value for the small close-in planet occurrence rate and compare it with the results of previous RV and transit (Kepler) works. Methods. We used archival HARPS RV datasets to calculate detection limits of a sample of stars that were previously analysed for their elemental abundances. For stars with known planets we first subtracted the Keplerian orbit. We then used this information to calculate planet occurrence rates according to a simplified Bayesian model in different regimes of stellar and planet properties. Results. Our results suggest that metal-poor stars and more massive stars host fewer low-mass close-in planets. We find the occurrence rates of these planets in the thin and thick disks to be comparable. In the iron-poor regimes, we find these occurrence rates to be significantly larger at the high-α region (thick-disk stars) as compared with the low-α region (thin-disk stars). In general, we find the average number of close-in small planets (2–100 days, 1–20M⊕) per star (FGK-dwarfs) to be: n¯p = 0.36 ± 0.05, while the fraction of stars with planets is Fh = 0.23−0.03+0.04. Qualitatively, our results agree well with previous estimates based on RV and Kepler surveys. Conclusions. This work provides a first estimate of the close-in small planet occurrence rates in the solar neighbourhood of the thin and thick disks of the Galaxy. It is unclear whether there are other stellar properties related to the Galactic context that affect small-planet occurrence rates, or if it is only the combined effects of stellar metal content and mass. A future larger sample of stars and planets is needed to address those questions.

Hot Exoplanet Atmospheres Resolved with Transit Spectroscopy (HEARTS)

WASP-127b is one of the puffiest exoplanets found to date, with a mass of only 3.4 Neptune masses, but a radius larger than that of Jupiter. It is located at the border of the Neptune desert, which describes the lack of highly irradiated Neptune-sized planets, and which remains poorly understood. Its large scale height and bright host star make the transiting WASP-127b a valuable target to characterise in transmission spectroscopy. We used combined EulerCam and TESS light curves to recalculate the system parameters. Additionally, we present an in-depth search for sodium in four transit observations of WASP-127b, obtained as part of the Hot Exoplanet Atmosphere Resolved with Transit Spectroscopy (HEARTS) survey with the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph. Two nights from this dataset were analysed independently by another team. The team claimed a detection of sodium that is incompatible with previous studies of data from both ground and space. We show that this strong sodium detection is due to contamination from telluric sodium emissions and the low signal-to-noise ratio in the core of the deep stellar sodium lines. When these effects are properly accounted for, the previous sodium signal is reduced to an absorption of 0.46 ± 0.20% (2.3σ), which is compatible with analyses of WASP-127b transits carried out with other instruments. We can fit a Gaussian to the D2 line, but the D1 line was not detected. This indicates an unusual line ratio if sodium exists in the atmosphere. Follow-up of WASP-127 at high resolution and with high sensitivity is required to firmly establish the presence of sodium and analyse its line shape.

Hot Exoplanet Atmospheres Resolved with Transit Spectroscopy (HEARTS)

Planet formation processes or evolution mechanisms are surmised to be at the origin of the hot Neptune desert. Studying exoplanets currently living within or at the edge of this desert could allow disentangling the respective roles of formation and evolution. We present the High Accuracy Radial velocity Planet Searcher (HARPS) transmission spectrum of the bloated super-Neptune WASP-166b, located at the outer rim of the Neptune desert. Neutral sodium is detected at the 3.4σ level (0.455 ± 0.135%), with a tentative indication of line broadening, which could be caused by winds blowing sodium farther into space, a possible manifestation of the bloated character of these highly irradiated worlds. We put this detection into context with previous work claiming a non-detection of sodium in the same observations and show that the high noise in the trace of the discarded stellar sodium lines was responsible for the non-detection. We highlight the impact of this low signal-to-noise ratio remnant on detections for exoplanets similar to WASP-166b.

Transits of Known Planets Orbiting a Naked-eye Star

Abstract Some of the most scientifically valuable transiting planets are those that were already known from radial velocity (RV) surveys. This is primarily because their orbits are well characterized and they preferentially orbit bright stars that are the targets of RV surveys. The Transiting Exoplanet Survey Satellite (TESS) provides an opportunity to survey most of the known exoplanet systems in a systematic fashion to detect possible transits of their planets. HD 136352 (Nu2 Lupi) is a naked-eye (V = 5.78) G-type main-sequence star that was discovered to host three planets with orbital periods of 11.6, 27.6, and 108.1 days via RV monitoring with the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph. We present the detection and characterization of transits for the two inner planets of the HD 136352 system, revealing radii of R ⊕ and R ⊕ for planets b and c, respectively. We combine new HARPS observations with RV data from the Keck/High Resolution Echelle Spectrometer and the Anglo-Australian Telescope, along with TESS photometry from Sector 12, to perform a complete analysis of the system parameters. The combined data analysis results in extracted bulk density values of g cm−3 and g cm−3 for planets b and c, respectively, thus placing them on either side of the radius valley. The combination of the multitransiting planet system, the bright host star, and the diversity of planetary interiors and atmospheres means this will likely become a cornerstone system for atmospheric and orbital characterization of small worlds.

Automatic Échelle Spectrograph Wavelength Calibration

Abstract Time-domain astronomy and the increasing number of exoplanet candidates call for reliable, robust, and automatic wavelength calibration. We present an algorithm for wavelength-calibrating échelle spectrographs that uses order-by-order extracted spectra and a list of laboratory wavelengths. Our approach is fully automatic and does not need the pixel locations of certain spectral features with which to anchor the wavelength solution, nor the true order number of each diffraction order. We use spectral features that are duplicated in adjacent orders to establish the scale-invariant component of the wavelength solution. We then match the central wavelengths of spectral features to laboratory wavelengths to establish the scale and higher-order components of the wavelength solution. We demonstrate our method on the four spectrographs of Las Cumbres Observatory’s Network of Robotic Échelle Spectrographs (NRES), on the High Accuracy Radial Velocity Planet Searcher (HARPS) spectrograph, and on synthetic data. We obtain a velocity-equivalent precision of ∼10 m s−1 on NRES. We achieve ∼1 m s−1 on HARPS, which agrees with the precision reported by the HARPS team. On synthetic data, we achieve the velocity precision set by Gaussian centroiding errors. Our algorithm likely holds for a wide range of spectrographs beyond the five presented here. We provide an open-source Python package, xwavecal, which outputs wavelength-calibrated spectra as well as the wavelengths of spectral features.