Planet Search Research Articles

Searching For Transiting Planets Around Halo Stars. I. Sample Selection and Validation

Abstract By measuring the elemental abundances of a star, we can gain insight into the composition of its initial gas cloud—the formation site of the star and its planets. Planet formation requires metals, the availability of which is determined by the elemental abundance. In the case where metals are extremely deficient, planet formation can be stifled. To investigate such a scenario requires a large sample of metal-poor stars and a search for planets therein. This paper focuses on the selection and validation of a halo star sample. We select ∼17,000 metal-poor halo stars based on their Galactic kinematics, and confirm their low metallicities ([Fe/H] < −0.5), using spectroscopy from the literature. Furthermore, we perform high-resolution spectroscopic observations using LBT/PEPSI and conduct detailed metallicity ([Fe/H]) analyses on a sample of 13 previously-known halo stars that also have hot kinematics. We can use the halo star sample presented here to measure the frequency of planets and to test planet formation in extremely metal-poor environments. The result of the planet search and its implications will be presented and discussed in a companion paper by Boley et al.

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.

TOI-431/HIP 26013: a super-Earth and a sub-Neptune transiting a bright, early K dwarf, with a third RV planet

ABSTRACTWe present the bright (Vmag = 9.12), multiplanet system TOI-431, characterized with photometry and radial velocities (RVs). We estimate the stellar rotation period to be 30.5 ± 0.7 d using archival photometry and RVs. Transiting Exoplanet Survey Satellite (TESS) objects of Interest (TOI)-431 b is a super-Earth with a period of 0.49 d, a radius of 1.28 ± 0.04 R⊕, a mass of 3.07 ± 0.35 M⊕, and a density of 8.0 ± 1.0 g cm−3; TOI-431 d is a sub-Neptune with a period of 12.46 d, a radius of 3.29 ± 0.09 R⊕, a mass of $9.90^{+1.53}_{-1.49}$ M⊕, and a density of 1.36 ± 0.25 g cm−3. We find a third planet, TOI-431 c, in the High Accuracy Radial velocity Planet Searcher RV data, but it is not seen to transit in the TESS light curves. It has an Msin i of $2.83^{+0.41}_{-0.34}$ M⊕, and a period of 4.85 d. TOI-431 d likely has an extended atmosphere and is one of the most well-suited TESS discoveries for atmospheric characterization, while the super-Earth TOI-431 b may be a stripped core. These planets straddle the radius gap, presenting an interesting case-study for atmospheric evolution, and TOI-431 b is a prime TESS discovery for the study of rocky planet phase curves.

Searching For Transiting Planets Around Halo Stars. ii. Constraining the Occurrence Rate of Hot Jupiters

Abstract Jovian planet formation has been shown to be strongly correlated with host-star metallicity, which is thought to be a proxy for disk solids. Observationally, previous works have indicated that Jovian planets preferentially form around stars with solar and supersolar metallicities. Given these findings, it is challenging to form planets within metal-poor environments, particularly for hot Jupiters that are thought to form via metallicity-dependent core accretion. Although previous studies have conducted planet searches for hot Jupiters around metal-poor stars, they have been limited due to small sample sizes, which are a result of a lack of high-quality data making hot-Jupiter occurrence within the metal-poor regime difficult to constrain until now. We use a large sample of halo stars observed by TESS to constrain the upper limit of hot-Jupiter occurrence within the metal-poor regime (−2.0 ≤ [Fe/H] ≤ −0.6). Placing the most stringent upper limit on hot-Jupiter occurrence, we find the mean 1σ upper limit to be 0.18% for radii 0.8–2 R Jupiter and periods 0.5–10 days. This result is consistent with previous predictions indicating that there exists a certain metallicity below which no planets can form.

Searching for Small Circumbinary Planets. I. The STANLEY Automated Algorithm and No New Planets in Existing Systems

Abstract No circumbinary planets have been discovered smaller than 3 R ⊕, yet planets of this small size comprise over 75% of the discoveries around single stars. The observations do not prove the nonexistence of small circumbinary planets; rather, they are much harder to find than around single stars because their transit timing variations are much larger than the transit durations. We present Stanley, an automated algorithm to find small circumbinary planets. It employs custom methods to detrend eclipsing binary light curves and stack shallow transits of variable duration and interval using N-body integrations. Applied to the Kepler circumbinaries, we recover all known planets, including the three planets of Kepler-47, and constrain the absence of additional planets of similar or smaller size. We also show that we could have detected <3 R ⊕ planets in half of the known systems. Our work will ultimately be applied to a broad sample of eclipsing binaries to (hopefully) produce new discoveries and derive a circumbinary size distribution that can be compared to that for single stars.

The CARMENES search for exoplanets around M dwarfs

Context. Stellar activity poses one of the main obstacles for the detection and characterisation of small exoplanets around cool stars, as it can induce radial velocity (RV) signals that can hide or mimic the presence of planetary companions. Several indicators of stellar activity are routinely used to identify activity-related signals in RVs, but not all indicators trace exactly the same activity effects, nor are any of them always effective in all stars. Aims. We evaluate the performance of a set of spectroscopic activity indicators for M dwarf stars with different masses and activity levels with the aim of finding a relation between the indicators and stellar properties. Methods. In a sample of 98 M dwarfs observed with CARMENES, we analyse the temporal behaviour of RVs and nine spectroscopic activity indicators: cross-correlation function (CCF) full-width-at-half-maximum (FWHM), CCF contrast, CCF bisector inverse slope (BIS), RV chromatic index (CRX), differential line width (dLW), and indices of the chromospheric lines Hα and calcium infrared triplet. Results. A total of 56 stars of the initial sample show periodic signals related to activity in at least one of these ten parameters. RV is the parameter for which most of the targets show an activity-related signal. CRX and BIS are effective activity tracers for the most active stars in the sample, especially stars with a relatively high mass, while for less active stars, chromospheric lines perform best. FWHM and dLW show a similar behaviour in all mass and activity regimes, with the highest number of activity detections in the low-mass, high-activity regime. Most of the targets for which we cannot identify any activity-related signals are stars at the low-mass end of the sample (i.e. with the latest spectral types). These low-mass stars also show the lowest RV scatter, which indicates that ultracool M dwarfs could be better candidates for planet searches than earlier types, which show larger RV jitter. Conclusions. Our results show that the spectroscopic activity indicators analysed behave differently, depending on the mass and activity level of the target star. This underlines the importance of considering different indicators of stellar activity when studying the variability of RV measurements. Therefore, when assessing the origin of an RV signal, it is critical to take into account a large set of indicators, or at least the most effective ones considering the characteristics of the star, as failing to do so may lead to false planet claims.

Perturbers: SPHERE detection limits to planetary-mass companions in protoplanetary disks

The detection of a wide range of substructures such as rings, cavities, and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in the disk (perturbers), but so far the only bona fide detection has been the two giant planets carving the disk around PDS 70. Here, we present a sample of 15 protoplanetary disks showing substructures in SPHERE scattered-light images and a homogeneous derivation of planet detection limits in these systems. To obtain mass limits we rely on different post-formation luminosity models based on distinct formation conditions, which are critical in the first million years of evolution. We also estimate the mass of these perturbers through a Hill radius prescription and a comparison to ALMA data. Assuming that one single planet carves each substructure in scattered light, we find that more massive perturbers are needed to create gaps within cavities than rings, and that we might be close to a detection in the cavities of RX J1604.3-2130A, RX J1615.3-3255, Sz Cha, HD 135344B, and HD 34282. We reach typical mass limits in these cavities of 3–10 MJup. For planets in the gaps between rings, we find that the detection limits of SPHERE high-contrast imaging are about an order of magnitude away in mass, and that the gaps of PDS 66 and HD 97048 seem to be the most promising structures for planet searches. The proposed presence of massive planets causing spiral features in HD 135344B and HD 36112 are also within SPHERE’s reach assuming hot-start models. These results suggest that the current detection limits are able to detect hot-start planets in cavities, under the assumption that they are formed by a single perturber located at the centre of the cavity. More realistic planet mass constraints would help to clarify whether this is actually the case, which might indicate that perturbers are not the only way of creating substructures.

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.

Spectroscopic frequency and mode identification of γ Doradus stars HD 109799 and HD 103257

ABSTRACT Following frequency and mode identification of two candidate γ Doradus stars, HD 103257 and HD 109799, both can be categorized as bona-fide members of the class. Over 250 high-resolution spectra of the two stars were collected at University of Canterbury Mount John Observatory (UCMJO) using the High Efficiency and Resolution Canterbury University Large Echelle Spectrograph (HERCULES) spectrograph. The spectra were cross-correlated with a synthetic δ-function template to produce line-profiles and further augmented with high-quality photometry from the Transiting Exoplanet Survey Satellite (TESS) mission, the Wide Angle Search for Planets (WASP) archive, and HIPPARCOS catalogue for frequency and mode analysis. Frequency analysis was carried out using FAMIAS and SigSpec. Two pulsation frequencies were identified in the spectra for HD 103257: 1.22496 ± 0.00001 and 1.14569 ± 0.00002 d−1, explaining 58.9 ${{\ \rm per\ cent}}$ of the variation across the line profiles. Both frequencies were characterized with best-fitting modes of (ℓ, m) = (1, 1). Two pulsation frequencies were identified in the spectra for HD 109799: 1.48679 ± 0.00002 and 1.25213 ± 0.00002 d−1, explaining 32.6 ${{\ \rm per\ cent}}$ of the variation across the line profiles. Both frequencies also yielded individual mode fits of modes (ℓ, m) = (1, 1). The excellent quality photometry from TESS observations has proven complementary to the work herein, and will provide a meaningful opportunity for deeper analysis of additional stars in the HERCULES catalogue. This will be a key component in the continued development of models relating to the processes within γ Doradus stars.

A SOPHIE RV search for giant planets around young nearby stars (YNS)

Context.The search of close (a≲ 5 au) giant planet (GP) companions with radial velocity (RV) around young stars and the estimate of their occurrence rates is important to constrain the migration timescales. Furthermore, this search will allow the giant planet occurrence rates to be computed at all separations via the combination with direct imaging techniques. The RV search around young stars is a challenge as they are generally faster rotators than older stars of similar spectral types and they exhibit signatures of magnetic activity (spots) or pulsation in their RV time series. Specific analyses are necessary to characterize, and possibly correct for, this activity.Aims.Our aim is to search for planets around young nearby stars and to estimate the GP occurrence rates for periods up to 1000 days.Methods.We used the SOPHIEspectrograph on the 1.93 m telescope at the Haute-Provence Observatory to observe 63 A−Myoung (<400 Myr) stars. We used our Spectroscopic data via Analysis of the Fourier Interspectrum Radial velocities software to compute the RVs and other spectroscopic observables. We then combined this survey with the HARPSYNS survey to compute the companion occurrence rates on a total of 120 youngA−Mstars.Results.We report one new trend compatible with a planetary companion on HD 109647. We also report HD 105693 and HD 112097 as binaries, and we confirm the binarity of HD 2454, HD 13531, HD 17250 A, HD 28945, HD 39587, HD 131156, HD 142229, HD 186704 A, and HD 195943. We constrained for the first time the orbital parameters of HD 195943 B. We refute the HD 13507 single brown dwarf (BD) companion solution and propose a double BD companion solution. Two GPs were previously reported from this survey in the HD 113337 system. Based on our sample of 120 young stars, we obtain a GP occurrence rate of 1−0.3+2.2% for periods lower than 1000 days, and we obtain an upper limit on BD occurrence rate of 0.9−0.9+2% in the same period range. We report a possible lack of close (P∈ [1;1000] days) GPs around young FK stars compared to their older counterparts, with a confidence level of 90%.

Detection Limits of Low-mass, Long-period Exoplanets Using Gaussian Processes Applied to HARPS-N Solar Radial Velocities

Abstract Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s−1. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s−1) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s−1) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.

The SOPHIE search for northern extrasolar planets. XVII. A wealth of new objects: Six cool Jupiters, three brown dwarfs, and 16 low-mass binary stars

Distinguishing classes within substellar objects and understanding their formation and evolution need larger samples of substellar companions such as exoplanets, brown dwarfs, and low-mass stars. In this paper, we look for substellar companions using radial velocity surveys of FGK stars with the SOPHIE spectrograph at the Observatoire de Haute-Provence. We assign here the radial velocity variations of 27 stars to their orbital motion induced by low-mass companions. We also constrained their plane-of-the-sky motion using HIPPARCOS and Gaia Data Release 1 measurements, which constrain the true masses of some of these companions. We report the detection and characterization of six cool Jupiters, three brown dwarf candidates, and 16 low-mass stellar companions. We additionally update the orbital parameters of the low-mass star HD 8291 B, and we conclude that the radial velocity variations of HD 204277 are likely due to stellar activity despite resembling the signal of a giant planet. One of the new giant planets, BD+631405 b, adds to the population of highly eccentric cool Jupiters, and it is presently the most massive member. Two of the cool Jupiter systems also exhibit signatures of an additional outer companion. The orbital periods of the new companions span 30 days to 11.5 yr, their masses 0.72 MJ–0.61 M⊙, and their eccentricities 0.04–0.88. These discoveries probe the diversity of substellar objects and low-mass stars, which will help constrain the models of their formation and evolution.

A Stellar Activity F-statistic for Exoplanet Surveys (SAFE)

Abstract In the search for planets orbiting distant stars, the presence of stellar activity in the atmospheres of observed stars can obscure the radial velocity signal used to detect such planets. Furthermore, this stellar activity contamination is set by the star itself and cannot simply be avoided with better instrumentation. Various stellar activity indicators have been developed that may correlate with this contamination. We introduce a new stellar activity indicator called the Stellar Activity F-statistic for Exoplanet surveys (SAFE) that has higher statistical power (i.e., probability of detecting a true stellar activity signal) than many traditional stellar activity indicators in a simulation study of an active region on a Sun-like star with a moderate-to-high signal-to-noise ratio. Also through simulation, the SAFE is demonstrated to be associated with the projected area on the visible side of the star covered by active regions. We also demonstrate that the SAFE detects statistically significant stellar activity in most of the spectra for HD 22049, a star known to have high stellar variability. Additionally, the SAFE is calculated for recent observations of the three low-variability stars HD 34411, HD 10700, and HD 3651, the latter of which is known to have a planetary companion. As expected, the SAFE for these three only occasionally detects activity. Furthermore, initial exploration appears to indicate that the SAFE may be useful for disentangling stellar activity signals from planet-induced Doppler shifts.

Dry or water world? How the water contents of inner sub-Neptunes constrain giant planet formation and the location of the water ice line

In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by forming a pressure bump outside their orbit in the protoplanetary disk. Depending on the relative position of a growing planetary core relative to the water ice line, water-rich pebbles might be blocked outside or inside the water ice line. Pebbles blocked outside the water ice line do not evaporate and thus do not release their water vapor into the gas phase, resulting in a dry inner disk, while pebbles blocked inside the water ice line release their water vapor into the gas phase, resulting in water vapor diffusing into the inner disk. As a consequence, close-in sub-Neptunes that accrete some gas from the disk should be dry or wet, respectively, if outer gas giants are outside or inside the water ice line, assuming that giant planets form fast, as has been suggested for Jupiter in our Solar System. Alternatively, a sub-Neptune could form outside the water ice line, accreting a large amount of icy pebbles and then migrating inward as a very wet sub-Neptune. We suggest that the water content of inner sub-Neptunes in systems with giant planets that can efficiently block the inward drifting pebbles could constrain the formation conditions of these systems, thus making these sub-Neptunes exciting targets for detailed characterization (e.g., with JWST, ELT, or ARIEL). In addition, the search for giant planets in systems with already characterized sub-Neptunes can be used to constrain the formation conditions of giant planets as well.

Matched filtering with non-Gaussian noise for planet transit detections

ABSTRACT We develop a method for planet detection in transit data, which is based on the matched filter technique, combined with the Gaussianization of the noise outliers. The method is based on Fourier transforms and is as fast as the existing methods for planet searches. The Gaussianized matched filter (GMF) method significantly outperforms the standard baseline methods in terms of the false positive rate, enabling planet detections at up to 30 per cent lower transit amplitudes. Moreover, the method extracts all the main planet transit parameters, amplitude, period, phase, and duration. By comparison to the state-of-the-art Gaussian process methods on both simulations and real data, we show that all the transit parameters are determined with an optimal accuracy (no bias and minimum variance), meaning that the GMF method can be used for both the initial planet detection and the follow-up planet parameter analysis.

A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn’s Star

Abstract Radial–velocity (RV) planet searches are often polluted by signals caused by gas motion at the star’s surface. Stellar activity can mimic or mask changes in the RVs caused by orbiting planets, resulting in false positives or missed detections. Here we use Gaussian process regression to disentangle the contradictory reports of planets versus rotation artifacts from Kapteyn’s star. To model rotation, we use joint quasiperiodic kernels for the RV and Hα signals, requiring that their periods and correlation timescales be the same. We find that the rotation period of Kapteyn’s star is 125 days, while the characteristic active-region lifetime is 694 days. Adding a planet to the RV model produces a best-fit orbital period of 100 yr, or 10 times the observing time baseline, indicating that the observed RVs are best explained by star rotation only. We also find no significant periodic signals in residual RV data sets constructed by subtracting off realizations of the best-fit rotation model and conclude that both previously reported “planets” are artifacts of the star’s rotation and activity. Our results highlight the pitfalls of using sinusoids to model quasiperiodic rotation signals.

A search for a fifth planet around HR 8799 using the star-hopping RDI technique at VLT/SPHERE

Context.Direct imaging of extrasolar giant planets demands the highest possible contrasts (ΔH ≳ 10 mag) at the smallest angular separations (∼0.1″) from the star. We present an adaptive optics observing method, called star-hopping, recently offered as standard queue observing (service mode) for the SPHERE instrument at the VLT. The method uses reference difference imaging (RDI) but, unlike earlier RDI applications, images of a reference star for PSF subtraction are obtained within minutes of observing the target star.Aims.We aim to significantly gain in contrast beyond the conventional angular differencing imaging (ADI) method to search for a fifth planet at separations less than 10 au, interior to the four giant planets of the HR 8799 system. The most likely semimajor axes allowed for this hypothetical planet, which were estimated via dynamical simulations in earlier works, were 7.5 au and 9.7 au within a mass range of 1–8MJup.Methods.We obtained 4.5 h of simultaneous low-resolution integral field spectroscopy (R ∼ 30,Y − Hband with IFS) and dual-band imaging (K1 andK2 bands with IRDIS) of the HR 8799 system, interspersed with observations of a reference star. The reference star was observed for about one-third of the total time and generally needs to be of similar brightness (ΔR ≲ 1 mag) and separated on sky by ≲1–2°. The hops between stars were made every 6–10 min, with only 1 min gaps in on-sky integration per hop.Results.We did not detect the hypothetical fifth planet at the most plausible separations, 7.5 and 9.7 au, down to mass limits of 3.6MJupand 2.8MJup, respectively, but attained an unprecedented contrast limit of 11.2 magnitudes at 0.1″. We detected all four planets with high signal-to-noise ratios. TheYJHspectra for planetsc,dwere detected with redderH-band spectral slopes than found in earlier studies. As noted in previous works, the planet spectra are matched very closely by some red field dwarfs. Finally, comparing the current locations of the planets to orbital solutions, we found that planetseandcare most consistent with coplanar and resonant orbits. We also demonstrated that with star-hopping RDI, the contrast improvement at 0.1″ separation can be up to 2 mag.Conclusions.Since ADI, meridian transit and the concomitant sky rotation are not needed, the time of observation can be chosen from within a window that is two to three times larger. In general, star-hopping can be used for stars fainter thanR = 4 magnitudes, since for these a reference star of suitable brightness and separation is usually available.

Tracking Advanced Planetary Systems (TAPAS) with HARPS-N

Context. We present the current status of and new results from our search for exoplanets in a sample of solar-mass evolved stars observed with the HARPS-N and the 3.6 m Telescopio Nazionale Galileo (TNG), and the High-Resolution Spectrograph (HRS) and the 9.2 m Hobby-Eberly Telescope (HET). Aims. The aim of this project is to detect and characterize planetary-mass companions to solar-mass stars in a sample of 122 targets at various stages of evolution from the main sequence to the red giant branch, mostly subgiants and giants, selected from the Pennsylvania-Toruń Planet Search sample, and to use this sample to study relations between stellar properties, such as metallicity, luminosity, and the planet occurrence rate. Methods. This work is based on precise radial velocity (RV) measurements. We have observed the program stars for up to 14 yr with the HET/HRS and the TNG/HARPS-N. Results. We present the analysis of RV measurements with the HET/HRS and the TNG/HARPS-N of four solar-mass stars, HD 4760, HD 96992, BD+02 3313, and TYC 0434-04538-1. We found that HD 4760 hosts a companion with a minimum mass of 13.9 MJ (a = 1.14 au, e = 0.23); HD 96992 is a host to a m sin i = 1.14 MJ companion on an a = 1.24 au and e = 0.41 orbit, and TYC 0434-04538-1 hosts an m sin i = 6.1 MJ companion on an a = 0.66 au and e = 0.08 orbit. In the case of BD+02 3313 we found a correlation between the measured RVs and one of the stellar activity indicators, suggesting that the observed RV variations may either originate in stellar activity or be caused by the presence of an unresolved companion. We also discuss the current status of the project and a statistical analysis of the RV variations in our sample of target stars. Conclusions. In our sample of 122 solar-mass stars, 49 ± 5% of them appear to be single and 16 ± 3% spectroscopic binaries. The three giants hosting low-mass companions presented in this paper join the six previously identified giants in the sample.

The Epoch of Giant Planet Migration Planet Search Program. I. Near-infrared Radial Velocity Jitter of Young Sun-like Stars

Abstract We present early results from the Epoch of Giant Planet Migration program, a precise radial velocity (RV) survey of more than 100 intermediate-age (∼20–200 Myr) G and K dwarfs with the Habitable Zone Planet Finder spectrograph (HPF) at McDonald Observatory’s Hobby–Eberly Telescope. The goals of this program are to determine the timescale and dominant physical mechanism of giant planet migration interior to the water ice line of Sun-like stars. Here, we summarize results from the first 14 months of this program, with a focus on our custom RV pipeline for HPF, a measurement of the intrinsic near-infrared RV activity of young Solar analogs, and modeling the underlying population-level distribution of stellar jitter. We demonstrate on-sky stability at the sub-2 m s−1 level for the K2 standard HD 3765 using a least-squares matching method to extract precise RVs. Based on a subsample of 29 stars with at least three RV measurements from our program, we find a median rms level of 34 m s−1. This is nearly a factor of 2 lower than the median rms level in the optical of 60 m s−1 for a comparison sample with similar ages and spectral types as our targets. The observed near-infrared jitter measurements for this subsample are well reproduced with a log-normal parent distribution with μ = 4.15 and σ = 1.02. Finally, by compiling rms values from previous planet search programs, we show that near-infrared jitter for G and K dwarfs generally decays with age in a similar fashion to optical wavelengths, albeit with a shallower slope and lower overall values for ages ≲1 Gyr.

SPECIES

Context. As part of the search for planets around evolved stars, we can understand planet populations around significantly higher-mass stars than the Sun on the main sequence. This population is difficult to study any other way, such as using radial-velocities to measure planet masses and orbital mechanics, since the stars are too hot and rotate too fast to present the quantity of narrow stellar spectral lines that is necessary for measuring velocities at the level of a few m s−1. Aims. Our goal is to estimate stellar parameters for all of the giant stars from the EXPRESS project, which aims to detect planets orbiting evolved stars, and study their occurrence rate as a function of stellar mass. Methods. We analysed the high-resolution echelle spectra of these stars and computed their atmospheric parameters by measuring the equivalent widths for a set of iron lines, using an updated method implemented during this work. Physical parameters, such as mass and radius, were computed by interpolating through a grid of stellar evolutionary models, following a procedure that carefully takes into account the post-main sequence evolutionary phases. The atmospheric parameters, as well as the photometric and parallax data, are used as constraints during the interpolation process. The probabilities of the star being in the red giant branch (RGB) or the horizontal branch (HB) are estimated from the derived distributions. Results. We obtained atmospheric and physical stellar parameters for the whole EXPRESS sample, which comprises a total of 166 evolved stars. We find that 101 of them are most likely first ascending the RGB phase, while 65 of them have already reached the HB phase. The mean derived mass is 1.41 ± 0.46 M⊙ and 1.87 ± 0.53 M⊙ for RGB and HB stars, respectively. To validate our method, we compared our derived physical parameters with data from interferometry and asteroseismology studies. In particular, when comparing to stellar radii derived from interferometric angular diameters, we find: ΔRinter = −0.11 R⊙, which corresponds to a 1.7% difference. Similarly, when comparing with asteroseismology, we obtain the following results: Δ log g = 0.07 cgs (2.4%), ΔR = −0.12 R⊙ (1.5%), ΔM = 0.08 M⊙ (6.2%), and Δage = −0.55 Gyr (11.9%). Additionally, we compared our derived atmospheric parameters with previous spectroscopic studies. We find the following results: ΔTeff = 22 K (0.5%), Δ log g = −0.03 (1.0%) and Δ[Fe/H] = −0.04 dex (2%). We also find a mean systematic difference in the mass with respect to those presented in the EXPRESS original catalogue of ΔM = −0.28 ± 0.27 M⊙, corresponding to a systematic mean difference of 16%. For the rest of the atmospheric and physical parameters we find a good agreement between the original catalogue and the results presented here. Finally, we find excellent agreement between the spectroscopic and trigonometric log g values, showing the internal consistency and robustness of our method. Conclusions. We show that our method, which includes a re-selection of iron lines and changes in the interpolation of evolutionary models, as well as Gaia parallaxes and newer extinction maps, can greatly improve the estimates of stellar parameters for giant stars compared to those presented in our previous work. This method also results in smaller mass estimates, an issue that has been described in results for giant stars from spectroscopy studies in the literature. The results provided here will improve the physical parameter estimates of planetary companions found orbiting these stars and give us insights into their formation and the effect of stellar evolution on their survival.