Radial Velocity Signal Research Articles

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.

CARMENES input catalog of M dwarfs

Context. Radial-velocity (RV) jitter caused by stellar magnetic activity is an important factor in state-of-the-art exoplanet discovery surveys such as CARMENES. Stellar rotation, along with heterogeneities in the photosphere and chromosphere caused by activity, can result in false-positive planet detections. Hence, it is necessary to determine the stellar rotation period and compare it to any putative planetary RV signature. Long-term measurements of activity indicators such as the chromospheric emission in the Ca II H&K lines (RHK′) enable the identification of magnetic activity cycles. Aims. In order to determine stellar rotation periods and study the long-term behavior of magnetic activity of the CARMENES guaranteed time observations (GTO) sample, it is advantageous to extract RHK′ time series from archival data, since the CARMENES spectrograph does not cover the blue range of the stellar spectrum containing the Ca II H&K lines. Methods. We have assembled a catalog of 11 634 archival spectra of 186 M dwarfs acquired by seven different instruments covering the Ca II H&K regime: ESPaDOnS, FEROS, HARPS, HIRES, NARVAL, TIGRE, and UVES. The relative chromospheric flux in these lines, RHK′, was directly extracted from the spectra by rectification with PHOENIX synthetic spectra via narrow passbands around the Ca II H&K line cores. Results. The combination of archival spectra from various instruments results in time series for 186 stars from the CARMENES GTO sample. As an example of the use of the catalog, we report the tentative discovery of three previously unknown activity cycles of M dwarfs. Conclusions. We conclude that the method of extracting $ R_{\mathrm{HK}}^\prime $ with the use of model spectra yields consistent results for different instruments and that the compilation of this catalog will enable the analysis of long-term activity time series for a large number of M dwarfs.

A HARPS-N mass for the elusive Kepler-37d: a case study in disentangling stellar activity and planetary signals

ABSTRACTTo date, only 18 exoplanets with radial velocity (RV) semi-amplitude <2 m s−1 have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 – known to host three transiting planets, one of which, Kepler-37d, should be on the cusp of RV detectability with modern spectrographs – as a case study in disentangling planetary and stellar activity signals. We show how two different statistical techniques – one seeking to identify activity signals in stellar spectra, and another to model activity signals in extracted RVs and activity indicators – can each enable a detection of the hitherto elusive Kepler-37d. Moreover, we show that these two approaches can be complementary, and in combination, facilitate a definitive detection and precise characterization of Kepler-37d. Its RV semi-amplitude of 1.22 ± 0.31 m s−1 (mass 5.4 ± 1.4 M⊕) is formally consistent with TOI-178b’s $1.05^{+0.25}_{-0.30}$ m s−1, the latter being the smallest detected RV signal of any transiting planet to date, though dynamical simulations suggest Kepler-37d’s mass may be on the lower end of our 1σ credible interval. Its consequent density is consistent with either a water-world or that of a gaseous envelope ($\sim 0.4{{\ \rm per\ cent}}$ by mass) surrounding a rocky core. Based on RV modelling and a re-analysis of Kepler-37 TTVs, we also suggest that the putative (non-transiting) planet Kepler-37e should be stripped of its ‘confirmed’ status.

Constraints on the mass and on the atmospheric composition and evolution of the low-density young planet DS Tucanae A b

Context. Observations of young close-in exoplanets are providing initial indications for the characteristics of the population and clues to the early stages of their evolution. Transiting planets at young ages are also key benchmarks for our understanding of planetary evolution via the verification of atmospheric escape models. Aims. We performed radial velocity (RV) monitoring of the 40 Myr old star DS Tuc A with HARPS at the ESO-3.6 m to determine the planetary mass of its 8.14-day planet, which was first revealed by the NASA TESS satellite. We also observed two planetary transits with HARPS and ESPRESSO at ESO-VLT to measure the Rossiter-McLaughlin (RM) effect and characterise the planetary atmosphere. We measured the high-energy emission of the host with XMM-Newton observations to investigate models for atmospheric evaporation. Methods. We employed a Gaussian Processes (GP) regression to model the high level of the stellar activity, which is more than 40 times larger than the expected RV planetary signal. GPs were also used to correct the stellar contribution to the RV signal of the RM effect. We extracted the transmission spectrum of DS Tuc A b from the ESPRESSO data and searched for atmospheric elements and molecules either by single-line retrieval and by performing cross-correlation with a set of theoretical templates. Through a set of simulations, we evaluated different scenarios for the atmospheric photo-evaporation of the planet induced by the strong XUV stellar irradiation. Results. While the stellar activity prevented us from obtaining a clear detection of the planetary signal from the RVs, we set a robust mass upper limit of 14.4 M⊕ for DS Tuc A b. We also confirm that the planetary system is almost (but not perfectly) aligned. The strong level of stellar activity hampers the detection of any atmospheric compounds, which is in line with other studies presented in the literature. The expected evolution of DS Tuc A b from our grid of models indicates that the planetary radius after the photo-evaporation phase will be 1.8–2.0 R⊕, falling within the Fulton gap. Conclusions. The comparison of the available parameters of known young transiting planets with the distribution of their mature counterpart confirms that the former are characterised by a low density, with DS Tuc A b being one of the less dense. A clear determination of their distribution is still affected by the lack of a robust mass measurement, particularly for planets younger than ~100 Myr.

IGRINS RV: A Precision Radial Velocity Pipeline for IGRINS Using Modified Forward Modeling in the Near-infrared*

Abstract Application of the radial velocity (RV) technique in the near-infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open-source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, an R ≡ λ/Δλ ∼ 45,000 spectrograph with simultaneous coverage of the H band (1.49–1.80 μm) and K band (1.96–2.46 μm). Using a modified forward-modeling technique, we construct high-resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 m s−1 in the K band and 31.1 m s−1 in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars, as well as the successful retrieval of planet-induced RV signals for both HD 189733 and τ Boo A; furthermore, our results affirm the presence of the Rossiter–McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use.

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.

HADES RV Programme with HARPS-N at TNG

Context. The high number of super-Earth and Earth-like planets in the habitable zone detected around M-dwarf stars in recent years has revealed these stellar objects to be the key to planetary radial velocity (RV) searches. Aims. Using the HARPS-N spectrograph within The HArps-n red Dwarf Exoplanet Survey (HADES) we have reached the precision needed to detect small planets with a few Earth masses using the spectroscopic radial velocity technique. HADES is mainly focused on the M-dwarf population of the northern hemisphere. Methods. We obtained 138 HARPS-N RV measurements between 2013 May and 2020 September of GJ 720 A, classified as an M0.5 V star located at a distance of 15.56 pc. To characterize the stellar variability and to distinguish the periodic variation due to the Keplerian signals from those related to stellar activity, the HARPS-N spectroscopic activity indicators and the simultaneous photometric observations with the APACHE and EXORAP transit surveys were analyzed. We also took advantage of TESS, MEarth, and SuperWASP photometric surveys. The combined analysis of HARPS-N RVs and activity indicators let us address the nature of the periodic signals. The final model and the orbital planetary parameters were obtained by simultaneously fitting the stellar variability and the Keplerian signal using a Gaussian process regression and following a Bayesian criterion. Results. The HARPS-N RV periodic signals around 40 days and 100 days have counterparts at the same frequencies in HARPS-N activity indicators and photometric light curves. We thus attribute these periodicities to stellar activity; the first period is likely associated with the stellar rotation. GJ 720 A shows the most significant signal at 19.466 ± 0.005 days with no counterparts in any stellar activity indices. We hence ascribe this RV signal, having a semi-amplitude of 4.72 ± 0.27 m s−1, to the presence of a sub-Neptune mass planet. The planet GJ 720 Ab has a minimum mass of 13.64 ± 0.79 M⊕, it is in circular orbit at 0.119 ± 0.002 AU from its parent star, and lies inside the inner boundary of the habitable zone around its parent star.

A super-Earth on a close-in orbit around the M1V star GJ 740

Context. M-dwarfs have proven to be ideal targets for planetary radial velocity (RV) searches due to their higher planet-star mass contrast, which favors the detection of low-mass planets. The abundance of super-Earth and Earth-like planets detected around this type of star motivates further such research on hosts without reported planetary companions. Aims. The HADES and CARMENES programs are aimed at carrying out extensive searches of exoplanetary systems around M-type stars in the northern hemisphere, allowing us to address, in a statistical sense, the properties of the planets orbiting these objects. In this work, we perform a spectroscopic and photometric study of one of the program stars (GJ 740), which exhibits a short-period RV signal that is compatible with a planetary companion. Methods. We carried out a spectroscopic analysis based on 129 HARPS-N spectra taken over a time span of 6 yr combined with 57 HARPS spectra taken over 4 yr, as well as 32 CARMENES spectra taken during more than 1 yr, resulting in a dataset with a time coverage of 10 yr. We also relied on 459 measurements from the public ASAS survey with a time-coverage of 8 yr, along with 5 yr of photometric magnitudes from the EXORAP project taken in the V, B, R, and I filters to carry out a photometric study. Both analyses were made using Markov chain Monte Carlo simulations and Gaussian process regression to model the activity of the star. Results. We present the discovery of a short-period super-Earth with an orbital period of 2.37756−0.00011+0.00013 d and a minimum mass of 2.96−0.48+0.50 M⊕. We offer an update to the previously reported characterization of the magnetic cycle and rotation period of the star, obtaining values of Prot = 35.563 ± 0.071 d and Pcycle = 2800 ± 150 d. Furthermore, the RV time series exhibits a possibly periodic long-term signal, which might be related to a Saturn-mass planet of ~100 M⊕.

Three years of HARPS-N high-resolution spectroscopy and precise radial velocity data for the Sun

Context. The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such a high-cadence, long-baseline data set is crucial for understanding spurious radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub- m s−1 systematics that needed to be corrected for. Aims. The goals of this manuscript are to (i) present a new data reduction software for HARPS-N, (ii) demonstrate the improvement brought by this new software during the first three years of the HARPS-N solar data set, and (iii) release all the obtained solar products, from extracted spectra to precise radial velocities. Methods. To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimised recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS. Results. The most significant improvement brought by the new DRS is a strong decrease in the day-to-day radial-velocity scatter, from 1.27 to 1.07 m s−1; this is thanks to a more robust method to derive wavelength solutions, but also to the use of calibrations closer in time. The newly derived solar radial-velocities are also better correlated with the chromospheric activity level of the Sun in the long term, with a Pearson correlation coefficient of 0.93 compared to 0.77 before, which is expected from our understanding of stellar signals. Finally, we also discuss how HARPS-N spectral ghosts contaminate the measurement of the calcium activity index, and we present an efficient technique to derive an index free of instrumental systematics. Conclusions. This paper presents a new data reduction software for HARPS-N and demonstrates its improvements, mainly in terms of radial-velocity precision, when applied to the first three years of the HARPS-N solar data set. Those newly reduced solar data, representing an unprecedented time series of 34 550 high-resolution spectra and precise radial velocities, are released alongside this paper. Those data are crucial to understand stellar activity signals in solar-type stars further and develop the mitigating techniques that will allow us to detect other Earths.

Tidal Effects on the Radial Velocities of V723 Mon: Additional Evidence for a Dark 3 M ⊙ Companion

Abstract Jayasinghe et al. identified a dark ≈3 M ⊙ companion on a nearly edge-on ≈60 day orbit around the red giant star V723 Monoceros as a black hole candidate in the mass gap. This scenario was shown to explain most of the data presented by Jayasinghe et al., except for periodic radial velocity (RV) residuals from the circular Keplerian model. Here we show that the RV residuals are explained by orbital phase-dependent distortion of the absorption line profile associated with changing visible fractions of the approaching and receding sides of the red giant star, whose surface is tidally deformed by and rotating synchronously with the dark companion. Our RV model constrains the companion mass M • = 2.95 ± 0.17 M ⊙ and orbital inclination (medians and 68.3% highest density intervals of the marginal posteriors) adopting the radius of the red giant 24.0 ± 0.9 R ⊙ as constrained from its SED and distance. The analysis provides independent support for the companion mass from ellipsoidal variations and the limits on the companion’s luminosity from the absence of eclipses, both derived by Jayasinghe et al. We also show that a common scheme to evaluate the tidal RV signal as the flux-weighted mean of the surface velocity field can significantly underestimate its amplitude for RVs measured with a cross-correlation technique, and present a modified prescription that directly models the distorted line profile and its effects on the measured RVs. The formulation will be useful for estimating the component masses and inclinations in other similar binaries.

Revisiting theKeplerfield withTESS: Improved ephemerides usingTESS2 min data

ABSTRACTUp to date planet ephemerides are becoming increasingly important as exoplanet science moves from detecting exoplanets to characterizing their architectures and atmospheres in depth. In this work, ephemerides are updated for 22 Kepler planets and 4 Kepler planet candidates, constituting all Kepler planets and candidates with sufficient signal to noise in the TESS 2 min data set. A purely photometric method is utilized here to allow ephemeris updates for planets even when they do not posses significant radial velocity data. The obtained ephemerides are of very high precision and at least seven years ‘fresher’ than archival ephemerides. In particular, significantly reduced period uncertainties for Kepler-411d, Kepler-538b, and the candidates K00075.01/K00076.01 are reported. O–C diagrams were generated for all objects, with the most interesting ones discussed here. Updated TTV fits of five known multiplanet systems with significant TTVs were also attempted (Kepler-18, Kepler-25, Kepler-51, Kepler-89, and Kepler-396), however these suffered from the comparative scarcity and dimness of these systems in TESS. Despite these difficulties, TESS has once again shown itself to be an incredibly powerful follow-up instrument as well as a planet-finder in its own right. Extension of the methods used in this paper to the 30 min-cadence TESS data and TESS extended mission has the potential to yield updated ephemerides of hundreds more systems in the future.

A Decade of Radial-velocity Monitoring of Vega and New Limits on the Presence of Planets

Abstract We present an analysis of 1524 spectra of Vega spanning 10 yr, in which we search for periodic radial-velocity variations. A signal with a periodicity of 0.676 day and a semi-amplitude of ∼10 m s−1 is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The activity signal appears to evolve on long timescales, which may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega’s photometric rotational modulation for the first time, with a total amplitude of only 10 ppm. A comparison of the spectroscopic and photometric amplitudes suggests that the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial-velocity injection recovery tests exclude the presence of transiting planets larger than 2 R ⊕ and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial-velocity signal with a period of 2.43 days and a semi-amplitude of 6 m s−1. If caused by an orbiting companion, its minimum mass would be ∼20 M ⊕; because of Vega’s pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet.

Validation of HD 183579b Using Archival Radial Velocities: A Warm Neptune Orbiting a Bright Solar Analog

Abstract As exoplanetary science matures into its third decade, we are increasingly offered the possibility of pre-existing, archival observations for newly detected candidates. This is particularly poignant for the TESS mission, whose survey spans bright, nearby dwarf stars in both hemispheres—precisely the types of sources targeted by previous radial velocity (RV) surveys. On this basis, we investigated whether any of the TESS Objects of Interest (TOIs) coincided with such observations, from which we find 18 single-planet candidate systems. Of these, one exhibits an RV signature that has the correct period and phase matching the transiting planetary candidates with a false-alarm probability of less than 1%. After further checks, we exploit this fact to validate HD 183579b (TOI-1055b). This planet is <4 R ⊕ and has better than 33% planetary mass measurements, thus advancing TESS’ primary objective of finding 50 such worlds. We find that this planet is among the most accessible small transiting planets for atmospheric characterization. Our work highlights that the efforts to confirm and even precisely measure the masses of new transiting planet candidates need not always depend on acquiring new observations—in some instances, these tasks can be completed with existing data.

IGRINS RV: A Precision RV Pipeline for IGRINS Using Modified Forward-Modeling in the Near-Infrared

Application of the radial velocity (RV) technique in the near infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, a R ~ 45,000 spectrograph with simultaneous coverage of the H band (1.49--1.80 $\mu$m) and K band (1.96--2.46 $\mu$m). Using a modified forward modeling technique, we construct high resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 $\rm m\,s^{-1}$ in the K band and 31.1 $\rm m\,s^{-1}$ in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars as well as the successful retrieval of planet-induced RV signals for both HD 189733 and $\tau$ Boo A; furthermore, our results affirm the presence of the Rossiter-McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use.

Radial Velocity Monitoring of the Young Star Hubble 4: Disentangling Star-spot Lifetimes from Orbital Motion*

Abstract We studied the weak-lined T Tauri star Hubble 4, a known long-period binary, and its star-spot phenomena. We used optical radial velocity (RV) data taken over a span of 14 yr (2004–2010, 2017–2019) at the McDonald Observatory 2.7 m Harlan J. Smith Telescope and single epoch imaging from the Hubble Space Telescope (HST)/Wide Field Camera 3 instrument. The observed and apparent RV variations show contributions, respectively, from the binary motion as well as from a large spot group on one of the stars, presumed to be the primary. Fitting and removing the orbital signal from the RVs, we found the lower bound on the lifetime of a previously identified large spot group on the surface of the star to be at least 5.1 yr. An ∼5 yr lower limit is a long, but not unprecedented, duration for a single spot group. The later epoch data indicate significant spot evolution has occurred, placing an upper bound on the spot group lifetime at 12 yr. We find that pre-main-sequence evolutionary models for the age of Taurus (∼2 Myr), combined with component mass estimates from the literature, permit us to reproduce the HST relative photometry and the binary-induced contribution to the apparent RV variations. The long-lived star spot we find on Hubble 4 has significant implications for dynamo models in young stars, as it adds evidence for long lifetimes of magnetic field topologies. There are also significant implications for young star exoplanet searches, as long-lived coherent RV signals may be spot induced and not the result of planetary motion.

Speckle Imaging Characterization of Radial Velocity Exoplanet Systems

Abstract We conducted speckle imaging observations of 53 stellar systems that were members of long-term radial velocity (RV) monitoring campaigns and exhibited substantial accelerations indicative of planetary or stellar companions in wide orbits. Our observations were made with blue and red filters using the Differential Speckle Survey Instrument at Gemini-South and the NN-Explore Exoplanet Stellar Speckle Imager at the WIYN telescope. The speckle imaging identifies eight luminous companions within 2″ of the primary stars. In three of these systems—HD 1388, HD 87359, and HD 104304—the properties of the imaged companion are consistent with the RV measurements, suggesting that these companions may be associated with the primary and the cause of the RV variation. For all 53 stellar systems, we derive differential magnitude limits (i.e., contrast curves) from the imaging. We extend this analysis to include upper limits on companion mass in systems without imaging detections. In 25 systems, we rule out companions with masses greater than 0.2 M ⊙, suggesting that the observed RV signals are caused by late-M dwarfs or substellar (potentially planetary) objects. On the other hand, the joint RV and imaging analysis almost entirely rules out planetary explanations of the RV signal for HD 19522 and suggests that the companion must have an angular separation below a few tenths of an arcsecond. This work highlights the importance of combined RV and imaging observations for characterizing the outer regions of nearby planetary systems.

Physical Parameters of the Multiplanet Systems HD 106315 and GJ 9827* †

Abstract HD 106315 and GJ 9827 are two bright, nearby stars that host multiple super-Earths and sub-Neptunes discovered by K2 that are well suited for atmospheric characterization. We refined the planets’ ephemerides through Spitzer transits, enabling accurate transit prediction required for future atmospheric characterization through transmission spectroscopy. Through a multiyear high-cadence observing campaign with Keck/High Resolution Echelle Spectrometer and Magellan/Planet Finder Spectrograph, we improved the planets’ mass measurements in anticipation of Hubble Space Telescope transmission spectroscopy. For GJ 9827, we modeled activity-induced radial velocity signals with a Gaussian process informed by the Calcium II H&K lines in order to more accurately model the effect of stellar noise on our data. We measured planet masses of M b = 4.87 ± 0.37 M ⊕, M c = 1.92 ± 0.49 M ⊕, and M d = 3.42 ± 0.62 M ⊕. For HD 106315, we found that such activity radial velocity decorrelation was not effective due to the reduced presence of spots and speculate that this may extend to other hot stars as well (T eff > 6200 K). We measured planet masses of M b = 10.5 ± 3.1 M ⊕ and M c = 12.0 ± 3.8 M ⊕. We investigated all of the planets’ compositions through comparison of their masses and radii to a range of interior models. GJ 9827 b and GJ 9827 c are both consistent with a 50/50 rock-iron composition, GJ 9827 d and HD 106315 b both require additional volatiles and are consistent with moderate amounts of water or hydrogen/helium, and HD 106315 c is consistent with a ∼10% hydrogen/helium envelope surrounding an Earth-like rock and iron core.

Stability constrained characterization of multiplanet systems

ABSTRACT Many discovered multiplanet systems are tightly packed. This implies that wide parameter ranges in masses and orbital elements can be dynamically unstable and ruled out. We present a case study of Kepler-23, a compact three-planet system where constraints from stability, transit timing variations (TTVs), and transit durations can be directly compared. We find that in this tightly packed system, stability can place upper limits on the masses and orbital eccentricities of the bodies that are comparable to or tighter than current state of the art methods. Specifically, stability places 68 per cent upper limits on the orbital eccentricities of 0.09, 0.04, and 0.05 for planets b, c, and d, respectively. These constraints correspond to radial velocity signals ≲ 20 cm s−1, are significantly tighter to those from transit durations, and comparable to those from TTVs. Stability also yields 68 per cent upper limits on the masses of planets b, c, and d of 2.2, 16.1, and 5.8 M⊕, respectively, which were competitive with TTV constraints for the inner and outer planets. Performing this stability constrained characterization is computationally expensive with N-body integrations. We show that SPOCK, the Stability of Planetary Orbital Configurations Klassifier, is able to faithfully approximate the N-body results over 4000 times faster. We argue that such stability constrained characterization of compact systems is a challenging ‘needle-in-a-haystack’ problem (requiring removal of 2500 unstable configurations for every stable one for our adopted priors) and we offer several practical recommendations for such stability analyses.

Investigating the young AU Mic system with SPIRou: large-scale stellar magnetic field and close-in planet mass

ABSTRACT We present a velocimetric and spectropolarimetric analysis of 27 observations of the 22-Myr M1 star AU Microscopii (AU Mic) collected with the high-resolution YJHK (0.98–2.35 μm) spectropolarimeter SPIRou from 2019 September 18 to November 14. Our radial velocity (RV) time-series exhibits activity-induced fluctuations of 45 m s−1 rms, ∼3 times smaller than those measured in the optical domain, that we filter using Gaussian Process Regression. We report a 3.9σ detection of the recently discovered 8.46 -d transiting planet AU Mic b, with an estimated mass of 17.1$^{+4.7}_{-4.5}$ M⊕ and a bulk density of 1.3 ± 0.4 g cm−3, inducing an RV signature of semi-amplitude K = 8.5$^{+2.3}_{-2.2}$ m s−1 in the spectrum of its host star. A consistent detection is independently obtained when we simultaneously image stellar surface inhomogeneities and estimate the planet parameters with Zeeman–Doppler imaging (ZDI). Using ZDI, we invert the time-series of unpolarized and circularly polarized spectra into surface brightness and large-scale magnetic maps. We find a mainly poloidal and axisymmetric field of 475 G, featuring, in particular, a dipole of 450 G tilted at 19° to the rotation axis. Moreover, we detect a strong differential rotation of dΩ = 0.167 ± 0.009 rad d−1 shearing the large-scale field, about twice stronger than that shearing the brightness distribution, suggesting that both observables probe different layers of the convective zone. Even though we caution that more RV measurements are needed to accurately pin down the planet mass, AU Mic b already appears as a prime target for constraining planet formation models, studying the interactions with the surrounding debris disc, and characterizing its atmosphere with upcoming space- and ground-based missions.

Architectures of Exoplanetary Systems. III. Eccentricity and Mutual Inclination Distributions of AMD-stable Planetary Systems

AbstractThe angular momentum deficit (AMD) of a planetary system is a measure of its orbital excitation and a predictor of long-term stability. We adopt the AMD-stability criterion to constrain the orbital architectures for exoplanetary systems. Previously, He et al. showed through forward modeling (SysSim) that the observed multiplicity distribution can be well reproduced by two populations consisting of a low and a high mutual inclination component. Here, we show that a broad distribution of mutual inclinations arising from systems at the AMD-stability limit can also match the observed Kepler population. We show that distributing a planetary system’s maximum AMD among its planets results in a multiplicity-dependent distribution of eccentricities and mutual inclinations. Systems with intrinsically more planets have lower median eccentricities and mutual inclinations, and this trend is well described by power-law functions of the intrinsic planet multiplicity (n):and, whereandare the medians of the eccentricity and inclination distributions. We also find that intrinsic single planets have higher eccentricities (σe,1 ∼ 0.25) than multiplanet systems and that the trends with multiplicity appear in the observed distributions of period-normalized transit duration ratios. We show that the observed preferences for planet-size orderings and uniform spacings are more extreme than what can be produced by the detection biases of the Kepler mission alone. Finally, we find that for systems with detected transiting planets between 5 and 10 days, there is another planet with a greater radial velocity signal ≃53% of the time.