skip to main content


Title: Modelling stellar variability in archival HARPS data: I - Rotation and activity properties with multidimensional Gaussian processes
ABSTRACT

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-metre per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m s–1 in Sun-like stars. To overcome this, a multidimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S + LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to (a) measure rotation periods for 49 slow rotators in our sample, (b) explore the impact of these results on the spin-down of middle-aged late F, G, and K stars, and (c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations programme, and help fine-tune current stellar structure and evolution models.

 
more » « less
NSF-PAR ID:
10491227
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
528
Issue:
4
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 5511-5527
Size(s):
["p. 5511-5527"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Aims: We present a detailed long-term study of the single M6 III giant RZ Ari to obtain direct and simultaneous measurements of the magnetic field, activity indicators, and radial velocity in order to infer the origin of its activity. We study its magnetic activity in the context of stellar evolution, and for this purpose, we also refined its evolutionary status and Li abundance. In general, for the M giants, little is known about the properties of the magnetic activity and its causes. RZ Ari possess the strongest surface magnetic field of the known Zeeman-detected M giants and is bright enough to allow a deep study of its surface magnetic structure. The results are expected to shed light on the activity mechanism in these stars.

    Methods: We used the spectropolarimeter Narval at the Télescope Bernard Lyot (Observatoire du Pic du Midi, France) to obtain a series of Stokes I and V profiles for RZ Ari. Using the least-squares deconvolution technique, we were able to detect the Zeeman signature of the magnetic field. We measured its longitudinal component by means of the averaged Stokes I and V profiles. In addition, we also applied Zeeman-Doppler imaging (ZDI) to search for the rotation period of the star, and we constructed a tentative magnetic map. It is the first magnetic map for a star that evolved at the tip of red giant branch (RGB) or even on the asymptotic giant branch (AGB). The spectra also allowed us to monitor chromospheric emission lines, which are well-known indicators of stellar magnetic activity. From the observations obtained between September 2010 and August 2019, we studied the variability of the magnetic field of RZ Ari. We also redetermined the initial mass and evolutionary status of this star based on current stellar evolutionary tracks and on the angular diameter measured from CHARA interferometry. Results: Our results point to an initial mass of 1.5Mso that this giant is more likely an early-AGB star, but a lotaction at the tip of the RGB is not completely excluded. With a v sin i of 6.0 ±0.5 km s−1, the upper limit for the rotation period is found to be 909 days. On the basis of our dataset and AAVSO photometric data, we determined periods longer than 1100 days for the magnetic field and photometric variability, and 704 days for the spectral line activity indicators. The rotation period determined on the basis of the Stokes V profiles variability is 530 days. A similar period of 544 days is also found for the photometric data. When we take this rotation period and the convective turnover time into account, an effective action of an α-ω type dynamo seems to be unlikely, but other types of dynamo could be operating there. The star appears to lie outside the two magnetic strips on the giant branches, where the α-ω-type dynamo is expected to operate effectively, and it also has a much higher lithium content than the evolutionary model predicts. These facts suggest that a planet engulfment could speed up its rotation and trigger dynamo-driven magnetic activity. On the other hand, the period of more than 1100 days cannot be explained by rotational modulation and could be explained by the lifetime of large convective structures. The absence of linear polarization at the time the magnetic field was detected, however, suggests that a local dynamo probably does not contribute significantly to the magnetic field, at least for that time interval. 
    more » « less
  2. ABSTRACT

    We report precise radial velocity (RV) observations of HD 212657 (= K2-167), a star shown by K2 to host a transiting sub-Neptune-sized planet in a 10 d orbit. Using Transiting Exoplanet Survey Satellite (TESS) photometry, we refined the planet parameters, especially the orbital period. We collected 74 precise RVs with the HARPS-N spectrograph between August 2015 and October 2016. Although this planet was first found to transit in 2015 and validated in 2018, excess RV scatter originally limited mass measurements. Here, we measure a mass by taking advantage of reductions in scatter from updates to the HARPS-N Data Reduction System (2.3.5) and our new activity mitigation method called CCF Activity Linear Model (CALM), which uses activity-induced line shape changes in the spectra without requiring timing information. Using the CALM framework, we performed a joint fit with RVs and transits using exofastv2 and find Mp = $6.3_{-1.4}^{+1.4}$  $\, M_{\hbox{$\oplus $}}$ and Rp = $2.33^{+0.17}_{-0.15}$  $\, R_{\hbox{$\oplus $}}$, which places K2-167 b at the upper edge of the radius valley. We also find hints of a secondary companion at a ∼22 d period, but confirmation requires additional RVs. Although characterizing lower mass planets like K2-167 b is often impeded by stellar variability, these systems especially help probe the formation physics (i.e. photoevaporation, core-powered mass-loss) of the radius valley. In the future, CALM or similar techniques could be widely applied to FGK-type stars, help characterize a population of exoplanets surrounding the radius valley, and further our understanding of their formation.

     
    more » « less
  3. null (Ed.)
    Context. Young open clusters (ages of less than 200 Myr) have been observed to exhibit several peculiarities in their chemical compositions. These anomalies include a slightly sub-solar iron content, super-solar abundances of some atomic species (e.g. ionised chromium), and atypical enhancements of [Ba/Fe], with values up to ~0.7 dex. Regarding the behaviour of the other s -process elements like yttrium, zirconium, lanthanum, and cerium, there is general disagreement in the literature: some authors claim that they follow the same trend as barium, while others find solar abundances at all ages. Aims. In this work we expand upon our previous analysis of a sample of five young open clusters (IC 2391, IC 2602, IC 4665, NGC 2516, and NGC 2547) and one star-forming region (NGC 2264), with the aim of determining abundances of different neutron-capture elements, mainly Cu  I , Sr  I , Sr  II , Y  II , Zr  II , Ba  II , La  II , and Ce  II . For NGC 2264 and NGC 2547 we present the measurements of these elements for the first time. Methods. We analysed high-resolution, high signal-to-noise spectra of 23 solar-type stars observed within the Gaia -ESO survey. After a careful selection, we derived abundances of isolated and clean lines via spectral synthesis computations and in a strictly differential way with respect to the Sun. Results. We find that our clusters have solar [Cu/Fe] within the uncertainties, while we confirm that [Ba/Fe] is super-solar, with values ranging from +0.22 to +0.64 dex. Our analysis also points to a mild enhancement of Y, with [Y/Fe] ratios covering values between 0 and +0.3 dex. For the other s -process elements we find that [X/Fe] ratios are solar at all ages. Conclusions. It is not possible to reconcile the anomalous behaviour of Ba and Y at young ages with standard stellar yields and Galactic chemical evolution model predictions. We explore different possible scenarios related to the behaviour of spectral lines, from the dependence on the different ionisation stages and the sensitivity to the presence of magnetic fields (through the Landé factor) to the first ionisation potential effect. We also investigate the possibility that they may arise from alterations of the structure of the stellar photosphere due to the increased levels of stellar activity that affect the spectral line formation, and consequently the derived abundances. These effects seem to be stronger in stars at ages of less than ~ 100 Myr. However, we are still unable to explain these enhancements, and the Ba puzzle remains unsolved. With the present study we suggest that other elements, for example Sr, Zr, La, and Ce, might be more reliable tracer of the s -process at young ages, and we strongly encourage further critical observations. 
    more » « less
  4. We present the discovery of an Earth-mass planet (Mbsini= 1.26 ± 0.21M) on a 15.6 d orbit of a relatively nearby (d ~9.6 pc) and low-mass (0.167 ± 0.011M) M5.0 V star, Wolf 1069. Sitting at a separation of 0.0672 ± 0.0014 au away from the host star puts Wolf 1069 b in the habitable zone (HZ), receiving an incident flux ofS= 0.652 ± 0.029S. The planetary signal was detected using telluric-corrected radial-velocity (RV) data from the CARMENES spectrograph, amounting to a total of 262 spectroscopic observations covering almost four years. There are additional long-period signals in the RVs, one of which we attribute to the stellar rotation period. This is possible thanks to our photometric analysis including new, well-sampled monitoring campaigns undergone with the OSN and TJO facilities that supplement archival photometry (i.e., from MEarth and SuperWASP), and this yielded an updated rotational period range ofProt= 150–170 d, with a likely value at 169.3−3.6+3.7. The stellar activity indicators provided by the CARMENES spectra likewise demonstrate evidence for the slow rotation period, though not as accurately due to possible factors such as signal aliasing or spot evolution. Our detectability limits indicate that additional planets more massive than one Earth mass with orbital periods of less than 10 days can be ruled out, suggesting that perhaps Wolf 1069 b had a violent formation history. This planet is also the sixth closest Earth-mass planet situated in the conservative HZ, after Proxima Centauri b, GJ 1061 d, Teegarden’s Star c, and GJ 1002 b and c. Despite not transiting, Wolf 1069 b is nonetheless a very promising target for future three-dimensional climate models to investigate various habitability cases as well as for sub-m s−1RV campaigns to search for potential inner sub-Earth-mass planets in order to test planet formation theories.

     
    more » « less
  5. Abstract

    We present a radial velocity (RV) analysis of TOI-1136, a bright Transiting Exoplanet Survey Satellite (TESS) system with six confirmed transiting planets, and a seventh single-transiting planet candidate. All planets in the system are amenable to transmission spectroscopy, making TOI-1136 one of the best targets for intra-system comparison of exoplanet atmospheres. TOI-1136 is young (∼700 Myr), and the system exhibits transit timing variations (TTVs). The youth of the system contributes to high stellar variability on the order of 50 m s−1, much larger than the likely RV amplitude of any of the transiting exoplanets. Utilizing 359 High Resolution Echelle Spectrometer and Automated Planet Finder RVs collected as part of the TESS-Keck Survey, and 51 High-Accuracy Radial velocity Planetary Searcher North RVs, we experiment with a joint TTV-RV fit. With seven possible transiting planets, TTVs, more than 400 RVs, and a stellar activity model, we posit that we may be presenting the most complex mass recovery of an exoplanet system in the literature to date. By combining TTVs and RVs, we minimized Gaussian process overfitting and retrieved new masses for this system: (mb−g=3.500.7+0.8,6.321.3+1.1,8.351.6+1.8,6.071.01+1.09,9.73.7+3.9,5.63.2+4.1M). We are unable to significantly detect the mass of the seventh planet candidate in the RVs, but we are able to loosely constrain a possible orbital period near 80 days. Future TESS observations might confirm the existence of a seventh planet in the system, better constrain the masses and orbital properties of the known exoplanets, and generally shine light on this scientifically interesting system.

     
    more » « less