Issue |
A&A
Volume 625, May 2019
|
|
---|---|---|
Article Number | A22 | |
Number of page(s) | 15 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201834028 | |
Published online | 01 May 2019 |
Precise radial velocities of giant stars
XII. Evidence against the proposed planet Aldebaran b⋆
1
Landessternwarte, Zentrum für Astronomie der Universität Heidelberg, Königstuhl 12, 69117 Heidelberg, Germany
e-mail: katja.reichert@uni-heidelberg.de
2
Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstraße 12-14, 69120 Heidelberg, Germany
3
Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
Received:
5
August
2018
Accepted:
1
March
2019
Context. Radial-velocity variations of the K giant star Aldebaran (α Tau) were first reported in the early 1990s. After subsequent analyses, the radial-velocity variability with a period of ∼629 d has recently been interpreted as caused by a planet of several Jovian masses.
Aims. We want to further investigate the hypothesis of an extrasolar planet around Aldebaran.
Methods. We combine 165 new radial-velocity measurements from Lick Observatory with seven already published data sets comprising 373 radial-velocity measurements. We perform statistical analyses and investigate whether a Keplerian model properly fits the radial velocities. We also perform a dynamical stability analysis for a possible two-planet solution. Furthermore, the possibility of oscillatory convective modes as cause for the observed radial-velocity variability is discussed.
Results. As best Keplerian fit to the combined radial-velocity data we obtain an orbit for the hypothetical planet with a smaller period (P = 607 d) and a larger eccentricity (e = 0.33 ± 0.04) than the previously proposed one. However, the residual scatter around that fit is still large, with a standard deviation of 117 ms−1. In 2006/2007, the statistical power of the ∼620 d period showed a temporary but significant decrease. Plotting the growth of power in reverse chronological order reveals that a period around 620 d is clearly present in the newest data but not in the data taken before ∼2006. Furthermore, an apparent phase shift between radial-velocity data and orbital solution is observable at certain times. A two-planet Keplerian fit matches the data considerably better than a single-planet solution, but poses severe dynamical stability issues.
Conclusions. The radial-velocity data from Lick Observatory do not further support but in fact weaken the hypothesis of a substellar companion around Aldebaran. Oscillatory convective modes might be a plausible alternative explanation of the observed radial-velocity variations.
Key words: stars: individual: α Tau / planets and satellites: detection / techniques: radial velocities / instrumentation: spectrographs
© ESO 2019
1. Introduction
The number of detected extrasolar planets has grown continuously since the early 1990s but has especially increased since 2014 due to the very successful Kepler space mission. Most of the detected extrasolar planets orbit main-sequence (MS) stars. Only a small fraction of them have been found around giant stars; the first one (around ι Draconis) was discovered in 2002 by Frink et al. (2002). The number of detected planetary companions around giant stars has almost doubled within recent years. Meanwhile, 112 substellar companions1 around 102 G and K giant stars have been found.
Due to the effects of the evolution of the host stars into giants on the orbits of their companions, such systems are of special interest as they allow us to get an idea about the dramatic destiny of our own solar system. Furthermore, evolved G or K giant stars are suitable targets for radial-velocity (RV) planet search programs due to their numerous absorption lines, compared to their F or A progenitor stars (Quirrenbach et al. 2011). In addition, detecting planets around their fast-rotating progenitor stars on the MS using spectroscopic methods is usually very challenging, whereas the further evolved giant stars enable the search for extrasolar planets around stars in a higher stellar mass range.
Hatzes & Cochran (1993) were the first to report low-amplitude variations with a period of 643 days in the RV data of the K giant star Aldebaran (HIP 21421, HD 29139), also known as α Tau. Even two years before the first detection of an extrasolar planet around a solar-type star, they already assumed that one explanation for these variations could be a planetary companion with a minimum mass of 11.4 Jupiter masses (assuming a stellar mass of 2.5 M⊙), since periods in this range are much too long to be caused by radial pulsations. However, they could not exclude that the variations are caused by rotational modulation of surface features or non-radial pulsations.
Furthermore, they found these variations to be consistent in amplitude and phase for at least ten years by comparison with previous measurements, as part of their McDonald Observatory Planetary Search (MOPS) program, and measurements of Walker et al. (1989). In addition, there were no hints for a period of ∼643 d in a subsequent analysis of the spectral line shapes by Hatzes & Cochran (1998). However, they found a period of 50 days in the line bisector measurements and assumed it to be due to stellar oscillations. Astrometric measurements of Gatewood (2008) could not refute the planet hypothesis of Hatzes & Cochran (1998), though Gatewood (2008) notes that the observed motion is “marginal at best and smaller than suggested by Hatzes & Cochran (1998)”. Gatewood (2008) also mentioned that the possible companion is likely to have a mass of 3 to 4 Jovian masses and not as large as the mass of 11 MJup which Hatzes & Cochran (1998) had derived.
Hatzes et al. (2015; hereafter H15) combined data taken between 1980 and 2013 at the Landessternwarte Tautenburg (TLS), the McDonald Observatory (McD), the Canada-France-Hawaii Telescope (CFHT), the Dominion Astrophysical Observatory (DAO), and at the Bohyunsan Optical Astronomy Observatory (BOAO), to re-investigate the nature of the long-period and long-lived RV variations in Aldebaran. They used three different data sets from McD, the original measurements using the 2.1 m telescope (hereafter “McD-2.1m” data set), measurements with the coude spectrograph in the so-called cs11 focus (hereafter “McD-CS11” data set), and measurements using the Tull Spectrograph (hereafter “McD-Tull” data set).
They were able to confirm the presence of a period in the range of ∼629 d in the RV data but did not find similar periods in any activity indicator. Therefore, they concluded that these variations are not caused by stellar activity and assumed them to be due to a substellar companion. Their orbital solution yielded a period P = 628.96 ± 0.90 d, an eccentricity e = 0.10 ± 0.05, and a RV amplitude K = 142.1 ± 7.2 ms−1 as parameters for the potential planetary companion.
Although the RV data published by H15 already extend over a time span of almost 34 years we add available data from Lick Observatory to those published by H15, as the Lick data were partially taken during a time with sparse sampling in the H15 data sets, to re-examine the planet hypothesis.
We note that Farr et al. (2018) have recently presented 302 new SONG RVs for Aldebaran obtained over a time span of about 120 days. They model their new RVs together with those from H15 with an orbiting planet plus a Gaussian process term due to acoustic oscillations. We show in Appendix A that their modelling of the acoustic oscillations (with a dominating period around 6 days) does not help with the interpretation of the long-term RV variation (many hundreds of days). Their model neither refutes nor confirms the planet hypothesis, although it does help with the asteroseismic determination of stellar parameters.
2. Stellar parameters of Aldebaran
Aldebaran is a K5 III giant star (Gray et al. 2006) and is located at a distance of about 20 pc (van Leeuwen 2007) in the constellation of Taurus. With an average apparent magnitude of mV = 0.87 (ESA 1997) it is one of the brightest stars in the sky.
We determined the stellar parameters of Aldebaran using evolutionary tracks by Bressan et al. (2012) and a Bayesian inference method based on the approach by Jørgensen & Lindegren (2005) with some improvements; see Stock et al. (2018) for details. Our method is capable of providing a probability for each of the two degenerate post-MS evolutionary stages in the HRD. Additionally, we provide asymmetric 1σ confidence intervals, as they describe the non-symmetrical probability density functions more properly. Our determined parameters are the mode values (maximum) of the probability density functions. The mode is less biased than the mean, especially if the distribution is bi-modal or truncated (see also Jørgensen & Lindegren 2005). As input parameters we used the metallicity [Fe/H] = −0.36 ± 0.10 as given by Hekker & Meléndez (2007), the HIPPARCOS parallax of 48.94 ± 0.77 mas (van Leeuwen 2007), and B, V photometry provided by the HIPPARCOS catalog (ESA 1997).
Aldebaran is not included in the first or second Gaia data release. Gaia DR2 however provides a parallax of 47.34 ± 0.11 mas (Gaia Collaboration 2018) for a star in the close vicinity of Aldebaran, α Tau B. It is a red dwarf with a mean magnitude of 11.9 mag in the G band (Gaia Collaboration 2018). Its angular distance to Aldebaran of only about 650 AU as well as its proper motion, which is comparable to that of Aldebaran, suggest that the red dwarf is gravitationally bound to Aldebaran. It would thus make sense to adopt the Gaia parallax of the companion as the overall parallax of the system. However, since the parallax of Aldebaran was determined already very precisely in van Leeuwen (2007), we used this value for determining its stellar parameters.
To avoid any biases which can occur due to converting parallax measurements to absolute magnitudes, we used the astrometry-based luminosity (Arenou & Luri 1999), a quantity which is linear in the parallax.
We determine that most probably (P = 97.8%) Aldebaran is a very late RGB star with a current mass of M⊙, which includes a Reimer’s mass loss of M* = 0.02 M⊙ calculated using η = 0.2 as suggested by Miglio et al. (2012). For the radius we obtain
M⊙ and for the age
Gyr. Table 1 also lists the stellar parameter for the unlikely case that Aldebaran is already on the HB.
Stellar parameters of Aldebaran.
Our determined mass is smaller than that of H15, who determined M* = 1.13 ± 0.11 M⊙ as well as a radius of R* = 36.68 ± 2.46 R⊙ and an age of τ = 6.6 ± 2.4 Gyr. The different stellar parameters are caused by the lower metallicity estimate of [Fe/H] = −0.36 ± 0.1 by Hekker & Meléndez (2007) compared to the estimate by H15 of [Fe/H] = −0.27 ± 0.05. Mozurkewich et al. (2003) interferometrically measured a diameter of 21.10 ± 0.21 mas for the limb-darkened case, corresponding to a radius of 46.3 ± 0.9 R⊙, while Richichi & Roccatagliata (2005) determine 20.58 ± 0.03 mas, corresponding to 45.2 ± 0.7 R⊙. Our value for the radius of Aldebaran is thus in excellent agreement with the interferometrically determined values.
3. The radial-velocity data from Lick Observatory
According to H15, a companion is the most likely explanation for the ∼629 d period in the RV measurements. However, it is still unclear whether giant stars show phenomena such as long-term non-radial pulsation patterns; hence, the planet hypothesis should be treated with caution. If a substellar companion is present, its period should be recovered in independent RV measurements. Therefore, we analysed available RV data from Lick Observatory and searched for RV variations in the range of ∼629 d. Furthermore, we added those data to the published ones of H15 and conducted the same analyses for the combined data sets as well as for the Lick data separately.
3.1. The radial-velocity survey at Lick Observatory
We monitored Aldebaran between 2000 and 2011 using the Hamilton Spectrograph at the 60 cm CAT Telescope at Lick Observatory as part of our RV survey of evolved stars (see, e.g., Frink et al. 2001). The survey started in 1999 and comprises a sample of 373 G and K giants. Almost all stars are of luminosity class III and brighter than V = 6 mag. For the survey, the 0.6 m Coudé Auxiliary Telescope (CAT) and the Hamilton Échelle Spectrograph were used, together with an iodine cell providing a stable reference. In these data, several extrasolar planets were already discovered (e.g., Frink et al. 2002; Reffert et al. 2006). In addition, these data provided valuable insights into the properties of (giant) extrasolar planets around giant stars in a statistical sense (e.g., Reffert et al. 2015).
3.2. The radial-velocity data of Aldebaran
A total of 165 RV measurements of α Tau were taken at Lick Observatory between JD = 2451781 (August 24, 2000) and JD = 2455896 (November 30, 2011). The RVs varied between −249.1 ms−1 and 325.5 ms−1 (mean RV subtracted), with a median uncertainty for the measurements of ∼3.8 ms−1. The observed dispersion around the mean of all data points was 129.3 ms−1. The first 69 measurements in particular are of great importance due to the small number of observations in the data of H15 between early 2000 and late 2002. Unfortunately, even with the complementary Lick data, there is still a large gap of more than six years between August 1993 and January 2000 in which no RV measurements are available.
The aim of combining the data sets of H15 with the Lick data is to obtain a long time baseline with as many measurements as possible in order to extract the long-period RV variation with confidence. The seven data sets of H15 contain 373 RV measurements in total, taken in the time ranges 1980–1993 and 2000–2013. By adding the Lick data we obtain a total number of 538 RV measurements.
4. Analysis of the radial-velocity data
4.1. Lomb-Scargle analysis
Since our RV data are unevenly sampled, we performed a statistical analysis corresponding to the least-squares spectral analysis of Lomb (1976) and Scargle (1982). Lomb-Scargle (L-S) periodograms showing all trial periods with their statistical power for the Lick data, the seven H15 data sets, and for the combined data sets (Lick + H15), respectively, are shown in Fig. 1.
![]() |
Fig. 1. Top: L-S periodogram for the Lick data. The most significant peak is at a period of 587.84 d with a FAP of 1.6 × 10−17. Middle: L-S periodogram for the H15 data. The most significant period is 629.08 d and has a FAP of less than 10−41. Bottom: L-S periodogram for the Lick data combined with the data of H15. The most significant peak is found at a period of 625.02 d with a FAP of less than 10−41. In all three panels, the significance thresholds represented by the vertical dashed lines correspond to FAPs of 0.1%, 1%, and 5%, respectively. |
The most significant peak in the Lick data has a period of 587.8 d, and a false-alarm probability (FAP) of 1.6 × 10−17. The period therefore deviates considerably from the most significant period of 629.1 d in the middle panel that represents the H15 data. Further periods around 766 d, 14 d, 998 d, and 263 d are also apparent in the top panel of Fig. 1. The left part of Table 2 lists all significant periods in the Lick data which have a FAP below 0.1%.
Significant periods and their FAPs in the L-S periodogram of the Lick data (left) and Lick + H15 data (right).
By reanalysing the H15 data we find additional periods which could not be found in the Lick data. Apart from the most prominent period of 629.1 d, which coincides with the period of the orbital solution for a single substellar companion presented in H15, further RV variations with periods around 687 d, 231 d, 304 d, and 873 d occur in the H15 data which are not present in the Lick data.
By adding the Lick data to the H15 data we find that now the most significant peak represents a period of 625.02 d. The L-S periodogram for the combined data is shown in the bottom panel of Fig. 1, while the exact values of the periods and their corresponding FAPs are listed in the right part of Table 2. The 766 d period from the Lick data has slightly increased and is now located at 775 d. Additional peaks are visible, for example at 687 d and 873 d, which are lacking in the Lick data and can solely be found in the H15 data. We again have a long-period RV variation near 1030 d.
4.2. Orbital solution for the Lick data
Deriving an orbital solution for a single companion for the Lick data yields a period similar to the most significant one from the L-S periodogram in the top panel of Fig. 1. However, since the classical L-S periodogram does not take into account eccentric orbits, the values for the periods do not coincide precisely. The fitted solution, which is plotted in the top-left panel of Fig. 2, has a period of P = 593.80 ± 2.75 d and deviates more than 12σ from the value of H15 of PH15 = 628.96 ± 0.90 d. Taking a stellar jitter of ∼90 ms−1 into account, the reduced χ2 scales down to almost unity (without jitter, ). With an eccentricity of e = 0.39 ± 0.06, the orbit of the potential planet is much more eccentric than in the orbital solution of H15 (eH15 = 0.10 ± 0.05) and deviates almost 4σ from this value.
![]() |
Fig. 2. Top left: RV measurements for Aldebaran from Lick Observatory (165 in total). The orbital solution for a substellar companion with a period of ∼594 d is represented by the black solid line. Bottom left: residual RV variations for the Lick RV measurements of Aldebaran after removing the orbital solution (black dashed line) for the Lick data of Table 3. Top right: RV measurements for Aldebaran from Lick Observatory. The orbital solution of H15 for a substellar companion with a period of 628.96 ± 0.90 d is overplotted and is represented by the black solid line. Bottom right: residual RV variations for the Lick RV measurements of Aldebaran after removing the orbital solution (black dashed line) of H15. |
In Table 3, we compare the orbital parameters of the Lick solution (second column) with those of the solution of H15 (last column).
Comparison of the parameters of our orbital solution for the Lick data, and for the combined data, respectively, and those of H15 for the substellar companion to Aldebaran.
Removing the orbital solution reveals a large scatter of the residual RVs around the fit over the whole time range, which can be seen in the bottom-left panel of Fig. 2. The maximum absolute residual RV with a value of 213.9 ms−1 is reached at JD = 2451931 (January 21, 2001). The rms value of the residual RVs is about 83 ms−1 and therefore of the same order as the expected stellar jitter which is around 91 ms−1, as estimated from its B − V colour index (see Trifonov et al. 2014).
The top-right panel of Fig. 2 shows the Lick RV data with the orbital solution of H15 overplotted. The orbital solution of H15 has a reduced of 1052.99 with respect to the Lick data, while our solution has
(in both cases without accounting for jitter).
H15 pointed out that there are time spans where their orbital solution does not fit the data points very well due to an apparent phase shift. This can be seen in Fig. 2 of H15, for example between JD = 2452600 (≈November 2002) and JD = 2453100 (≈April 2004). Such an apparent phase shift around late 2003 and early 2004 between RV measurements and orbital solution is also observable for the Lick data (for our Keplerian fit as well as for the overplotted H15 fit).
This indicates that the signal during that time period is not consistent with a long-lived, coherent one. We can therefore conclude that the Lick data do not further support the hypothesis of a single substellar companion, as they cannot be described very well by a single-planet orbital solution and in particular not by the solution of H15.
4.3. Orbital solution for the combined data
Figure 3 shows the RVs of the Lick data together with the seven H15 data sets as a function of time. There are still no data available between late 1993 and early 2000. The dispersion around the mean of all data points is 142.74 ms−1, while the median measurement error is 4.2 ms−1. The best fitting orbital solution has period P = 607.22 ± 1.17 d, eccentricity e = 0.33 ± 0.04, RV amplitude K = 156.98 ms−1, and a reduced χ2 of 519.55; a stellar jitter of about 160 m s−1 would be required for a reduced of the order of 1, not too far from the expected value. The value of H15 of 628.96 ± 0.90 days deviates more than 15σ from the value determined for the combined data. The discrepancy of the eccentricities (now e = 0.33, and thus much higher than the H15 value of eH15 = 0.10) is also not negligible.
![]() |
Fig. 3. RV measurements for Aldebaran from Lick Observatory and all the data sets of H15. The orbital solution for a substellar companion with a period of 607.22 d is represented by the black solid line in the top panel. Lower panel: residuals with respect to this one-planet fit. |
Removing the orbital solution again reveals high residual RVs with a maximum absolute value of 435.3 ms−1 at JD = 2448201 (November 5, 1990), as shown in the bottom panel of Fig. 3. The rms value for the residual RVs is 117 ms−1.
Checking the RV data for alias frequencies is crucial as they often interfere with the originally observed frequencies. For further explanation, see for example Dawson & Fabrycky (2010). We find that the ∼620 d period is not caused by a linear combination of other frequencies, but that it is responsible for the 230 day peak in the L-S periodogram (f230 d ≈ f1 yr + f620 d). Removing the orbital solution therefore causes the 230.66 d period to disappear (see Fig. 4). Very significant peaks at 775 d and 687 d are still present. A period around 873 d is also still observable in Fig. 4 above the significance level, although its power decreased relative to the highest peaks. It might be the corresponding alias to the 230 d period, since f873 d ≈ f1 yr − f620 d. The same can be observed for the Lick data (left part of Table 2), where the 225 d and 988 d periods could be 1 yr aliases of the dominant 588 d period (f225 d ≈ f1 yr + f588 d and f998 d ≈ f1 yr − f588 d).
![]() |
Fig. 4. Lomb-Scargle periodogram for the Lick data combined with the data of H15 after removing the 607.22 d period of the orbital solution. The alias period 230.66 d has disappeared as well. The significance thresholds represented by the vertical dashed lines correspond to FAPs of 0.1%, 1%, and 5%, respectively. |
4.4. Growth of the Lomb-Scargle power
An exoplanet orbiting a star should cause a strictly periodic RV signal whose statistical L-S power increases linearly with the number of measurements. For simulated data sampled in the same manner as the actual 538 RV measurements, this is shown by the blue circles in both panels of Fig. 5. The simulated data have the same measurement errors as the real data, and a stellar jitter of 91 ms−1 has been added.
![]() |
Fig. 5. Lomb-Scargle power of the 625 d period as a function of the number of data points (N). The red squared data points represent the real data while the blue circular data points represent simulated data. Top: growth of the Lomb-Scargle power from 1980–2013. Data points were added in steps of 20 in chronological order. A decrease in statistical power is observable for the real data between N = 300 and N = 320. Bottom: growth of the Lomb-Scargle power from 2013–1980. Data points were added in steps of 20 in reversed chronological order. Before N = 240 (which represent the newest RV data) the ∼620 d period is clearly visible as the growth in power shows the same linear slope as the simulated data. The older RV data (taken before 2006) however cause a decrease in the slope and a drop in power at N = 340. |
The actual measured data from H15 and the Lick survey, added in chronological order and represented by the red squares in the top panel of Fig. 5, have a shallower slope, and at some point show a non-linear behaviour; the L-S power temporarily decreases between the 300th and the 320th data point, corresponding to the time range between November 2006 and February 2007. In Fig. 5 in H15, such a decrease in power is also observable in the same time span (we note that H15 erroneously assign it to the time range between 2002 and early 2004, where they observed an apparent phase shift between the RV measurements and their orbital solution). The systematic decrease in L-S power during this time span weakens the hypothesis of a substellar companion.
Adding the data in reverse chronological order (see bottom panel of Fig. 5) reveals that the newest data (∼2006–2013), which represent in this panel N = 1 to N = 240, show a coherent RV variation in the range of 620 d; the L-S power linearly increases with the number of measurements. After N = 240 the slope decreases significantly, which is a hint that the 620 d signal is not as clearly present in the data taken before 2006 as it should be if it is caused by a substellar companion.
Since the generalized Lomb-Scargle periodogram (GLS) of Zechmeister & Kürster (2009) has a different normalization (it is normalized to unity) and can show a slightly different behaviour we compare the growth of the power in the classical L-S periodogram (from the top panel in Fig. 5) with the growth of power in the GLS. While for the combined data the L-S power remains almost unity for the simulated data, it shows multiple strong drops for the real data, which can be seen in Fig. 6. The most significant period in the GLS even changes while adding further RV measurements, which is indicated by the different colours. After the first 60 data points, the period in the range of ∼620 d arises and is most significant. After 140–160 measurements, two further periods around 550 d and 760 d become dominant. The latter one replaces the ∼620 d peak as most significant period from N = 180 to N = 240, while the ∼550 d RV signal dominates the periodogram between N = 320 and N = 500. Although the ∼620 d period completely vanishes between N = 400 and N = 500 (most likely due to its proximity to the ∼550 d period), it again becomes the most significant period after 500 measurements.
![]() |
Fig. 6. Growth of the L-S power in the GLS as a function of the number of data points (N). Data points were added in steps of 20 in chronological order. The red squared data points represent the ∼620 d period, while the green triangles represent the ∼750 d period and the yellow triangles the ∼550 d period. The blue circular data points represent simulated data. See the text for more details. |
4.5. Orbital solution with two planetary companions
As we still have two very significant periods (775 d and 687 d) left in Fig. 4, it is worth also investigating solutions with multiple planets. The best Keplerian two-planet solution we find has orbital periods of 614.10 d and 772.83 d, where the period of the potential outer planet corresponds to the most significant period of Fig. 4. This two-planet solution for all eight data sets is shown in the upper panel of Fig. 7. Especially for the old RV data taken between 1980 and 1993, the scatter around the two-planet fit is still very large, as can be seen in the lower panel of Fig. 7. The orbital parameters of both potential planets as well as the reduced χ2 and the rms scatter of the residual RVs around the two-planet fit are listed in Table 4.
![]() |
Fig. 7. RV measurements for Aldebaran from Lick Observatory and all data sets of H15. The orbital solution for two substellar companions with periods of 614.10 d and 772.83 d is represented by the black solid line in the top panel. Lower panel: residuals with respect to this two-planet fit. |
Planetary parameters of the orbital solution for two substellar companions to Aldebaran for the Lick + H15 data sets.
However, as both planets would have several Jovian masses and quite similar semi-major axes it is unlikely that this corresponds to a stable configuration. A dynamical fit to the data was not possible as the fit did not converge, or the solution was not stable over the timescale of the observations. In order to find possible dynamically stable solutions we create 11 million samples of orbital parameters for both planets, 10 million samples based on Gaussian distributions around the 3σ parameter space of each of our Keplerian best-fit parameters, namely period, mean anomaly, eccentricity, longitude of periastron and RV semi-amplitude. An additional 1 million samples were uniformly sampled around the 10σ parameter space of each of our Keplerian best-fit parameters in order to find possible dynamically stable islands which in case of significant dynamical interaction could lie far away from the best Keplerian fit parameters. The errors of our Keplerian best-fit parameters for the hypothetical two-planet solution that were used for the Monte-Carlo (MC) sampling are determined by bootstrapping and are slightly larger than the errors based on the covariance matrix shown in Table 4. The bootstrap uncertainties for all solutions presented in this paper are given in Table A.1.
We use co-planar configurations with i = 90° and keep the mass of the host star fixed at 0.91 M⊙. We test all 11 million samples for long-term dynamical stability using the Wisdom-Holman integrator (also known as MVS; Wisdom & Holman 1991), which was modified to work with Jacobi input elements. We choose 1 Myr as the maximum integration time as well as an integration time step of 1 day. The fit is regarded unstable if at any time the separation between the star and one of the companions exceeds 5 AU (planet ejection) or falls below 0.001 AU (star-planet collision, since this limit is well below the stellar radius).
We find only a very small number of dynamically stable solutions around our best Keplerian fit. Out of the 10 million Gaussian sampled combinations of orbital parameters investigated only 2016 result in stable configurations according to our definition, which corresponds to a fraction of only 0.02%. From the 1 million uniformly sampled combinations of orbital parameters around the 10σ parameter space of the best Keplerian fit we find 455 stable configurations (0.05%). These numbers would further decrease if we allow inclination angles less than 90°.
Despite the small fraction of stable configurations, these represent an interesting possibility for a stable two-planet system. To see if these configurations could explain the observed RVs we reconstruct the radial stellar motion induced by these synthesised stable configurations over the past ∼35 yr and we test it against the data. For this purpose we adopted a dynamical model (Tan et al. 2013; Trifonov et al. 2014) fed with fixed initial orbital parameters from the stable samples and adjustable RV offsets for each dataset. By studying the quality and the overall dynamics of these models, we find the following:
-
(i)
Twelve out of the 2016 Gaussian sampled stable solutions are in a 4:3 MMR, but the quality of these models is poor as they show rather large residuals. The rms of the residuals of the best 4:3 MMR model is 127.72 ms−1, which is larger than the rms of the simpler one-planet model for the combined data, which has 117.37 ms−1.
-
(ii)
The rest of the Gaussian sampled stable solutions were found in 1:1 MMR configurations, which induce an RV curve practically indistinguishable from a one-planet model (e.g. see Laughlin & Chambers 2002). Therefore, the quality of these considerably more complex 1:1 MMR two-planet models is comparable, or in most cases even worse than the one-planet model, which once again speaks against the two-planet hypothesis.
-
(iii)
Investigating the dynamical models of the 455 stable solutions, which were uniformly sampled around the 10σ parameter space of the best Keplerian fit, reveals that these fits are not consistent with the observed RV data. For the vast majority of stable solutions, the eccentricity of the potential second planet (P = 772.83 d) is higher than e = 0.4, resulting in RV fits that describe the observed RV variations poorly and show clear systematics in the residuals. The rms of the best-fitting dynamical model out of the 455 stable solutions is 116.9 ms−1, which is no significant improvement compared to the simpler one-planet model.
Further, almost all of the stable solutions have crossing orbits; it is rather unclear whether and how such a system could form. Given the tiny fraction of stable solutions in our MC samples, the presence of crossing orbits in almost all stable solutions, the presence of systematics in the residuals of the RV data, and the poor quality of these stable configurations, we conclude that the two-planet model is not a viable solution to explain the RV data of Aldebaran.
5. Oscillatory convective modes?
Long secondary periods, which are of the order of 400–1500 days and are found in about 25–30% of AGB stars (Wood et al. 2004), have been known for a long time, but not theoretically understood. Saio et al. (2015) suggested that non-adiabatic g− modes present in the deep convective envelopes of luminous giant stars, so-called oscillatory convective modes, could be responsible for the observed periods. They have presented calculations which show that, depending on the adopted value for the mixing length parameter, oscillatory convective modes would be present in the same region of period-luminosity space where the stars with long secondary periods are found.
In order to test whether oscillatory convective modes could possibly explain the observed RV variations in Aldebaran, we compare its location in period-luminosity space with those of other stars known to exhibit long secondary periods, similar to the diagram presented by Saio et al. (2015), in Fig. 8. The small dots represent stars with photometric periods determined from the OGLE survey of the Large Magellanic Cloud (LMC; Soszyński et al. 2009). The stars have been identified in 2MASS and dereddened, so that their extinction corrected colours best match expected intrinsic colours of giant stars as given by Bessell & Brett (1988). The scaling factors for the absorption in the various photometric bands were taken from Schlegel et al. (1998). We then used the infrared flux method (Blackwell & Lynas-Gray 1994) to obtain luminosities.
We only show in Fig. 8 those stars for which the 2MASS photometry had errors smaller than 0.05 mag in K band, and for which the absorption could be relatively well determined (good match to expected intrinsic colours after correction for absorption, and formal error of the fitted absorption less than 0.5 mag in V band). These criteria eliminated many of the initial 91 644 stars in the OGLE sample with 2MASS photometry, leaving 55 832 stars. Of those, 32 601 have photometric amplitudes smaller than 0.015 mag and are shown as green points in Fig. 8. We choose an upper limit on the photometric amplitude, since the stars in our Lick survey are known not to be photometrically variable based on HIPPARCOS photometry. The subset of black points are those 3833 stars which have been identified by Soszyński et al. (2009) to fall on sequence D, that is, those stars with long secondary periods. Whether or not the Lick survey stars are photometrically variable with periods similar to these sequence-D stars and small photometric amplitudes (mmag level) is not known, since the cadence of the HIPPARCOS photometry would not have allowed this kind of photometric periodicity to be detected.
![]() |
Fig. 8. Period-luminosity diagram to illustrate the position of Aldebaran and γ Dra (star symbols) with respect to published exoplanets from our Lick survey (red circles, labeled with HIP number) and candidate exoplanets (blue circles). The small dots are stars with periods identified by the OGLE survey in the LMC; the subset of black dots are those stars identified as showing long secondary periods in the OGLE photometry. The thick solid lines show model calculations by Saio et al. (2015) for the periods of oscillatory convective modes for stars of various masses, and for a mixing length parameter of 1.2. Aldebaran and γ Dra’s periods are rather close to the predicted periods of oscillatory convective modes, especially if one were to extrapolate the models for 1.0 and 1.3 M⊙. |
We overplotted the computed periods for stars of 1.0, 1.3, and 2.0 M⊙ and various luminosities from Saio et al. (2015; solid lines), for a mixing length parameter of 1.2. These periods coincide in general with the order of magnitude of the observed periods and stellar luminosities in the OGLE sample. A mixing length parameter value of 1.9 would provide a poorer fit to the sequence-D stars (black), but in turn a better fit to the radial modes of AGB stars on sequences A, B, C, and C′ (green points). Since the models by Saio et al. (2015) for the calculation of oscillatory convective modes do not take into account turbulent pressure and possible overshooting, the precise shape and location of the models could change. In particular, Saio et al. (2015) note that a better theory of convection could extend the range in which oscillatory convective modes are found.
Furthermore, we note that the theoretical calculations of the oscillatory convective mode periods have been performed for stellar models with slightly sub-solar metallicity (Z = 0.008, corresponding to [Fe/H] ≈ −0.3 dex). For the comparison with data from the LMC, one should take into account that the average metallicity is somewhat smaller in the LMC than in the solar neighbourhood; Cole et al. (2000) for instance determined an average metallicity of around −0.6 dex for a sample of field giants in the LMC. On the other hand, the median metallicity of stars in the Lick sample is −0.09 dex, while the stars orbited by giant planets (and shown in Fig. 8) have preferentially higher metallicities. Thus, we caution that different underlying metallicities could affect the comparison between data and theory in Fig. 8 to some extent.
We also indicate in Fig. 8 the locations of Aldebaran and γ Dra (star symbols), as well as published exoplanets (red circles) from our Lick survey. Blue circles indicate stars with more or less periodic RV patterns with uncertain interpretation at the moment; the RV variations could either be caused by an orbiting planet or by non-radial pulsations. There is a clear separation in luminosity visible between those categories: published planet discoveries from our survey are typically found around stars with log(L/L⊙) smaller than 2.5 (L < 300 L⊙), while the more doubtful cases have log(L/L⊙) larger than 2.5. It would be tempting now to attribute this to the onset of oscillatory convective modes for the higher-luminosity stars, but it could also be the result of the well-known fact that higher-luminosity stars show larger RV jitter. This is very well understood and is due to solar-like (i.e. radial) pulsations, with periods much smaller (typically hours do days) than the long secondary periods discussed here. The larger RV jitter makes it harder to distinguish true exoplanets with Keplerian RV curves from less regular RV curves which show periods, amplitude, and/or phase changes, such as Aldebaran. In any case, a much larger number of RV measurements taken over a longer time period is necessary for an unambiguous planet detection around evolved stars of high luminosity.
Hatzes et al. (2018) already noted that the period of γ Dra is very close to the period predicted for oscillatory convective modes by Saio et al. (2015), and the same is true for Aldebaran, which has very similar stellar parameters and evolutionary stage to γ Dra. Although we presently have no way to unambiguously confirm or reject oscillatory convective modes in Aldebaran, Fig. 8 cautions us to be very careful with exoplanet claims around stars with luminosities larger than about 300 L⊙.
6. Discussion and conclusion
The number of confirmed detections of extrasolar planets around giant stars is small compared to the total number of confirmed extrasolar planets. Therefore, it is of special interest to have a clean sample as a basis for proper statistics of planets around giant stars. Here we have further investigated the evidence for a Jupiter-mass planet around the K giant Aldebaran, which had been suggested by Hatzes et al. (2015).
We conclude that the RV data from Lick Observatory do not support the existence of this planet. By fitting an orbital solution to the Lick data combined with those of Hatzes et al. (2015) we obtain large residual RVs up to 435.3 ms−1. Even if we take into account the fact that Aldebaran is an evolved K giant with a stellar jitter in the range of 100 ms−1, this large systematic scatter around the orbital solution weakens the planet hypothesis.
Plotting the power of the most significant period in the Lomb-Scargle periodogram while continuously adding further measurements in chronological order reveals a strong decrease in statistical power after 300 measurements, contrary to expectations for a strictly periodic signal. Using the generalized L-S periodogram instead yields several instances where the power declines, and further reveals that the claimed planet period is not continuously present; it changes considerably in significance depending on the number of measurements.
A system with two planets formally fits the RV data better than a single planet. However, both planets would have at least several Jupiter masses and similar semi-major axes (around 1.5 AU), leading to large gravitational interactions. We perform a dynamical analysis and search for stable solutions in the vicinity of the best Keplerian fit, integrating for up to 1 Myr. We find that only 0.02% of the Gaussian sampled and only 0.05% of the uniformly sampled configurations are stable over this time range, and many of these formally stable solutions involve crossing orbits. We thus conclude that a two-planet configuration cannot explain the RV observations of Aldebaran.
The reason for the observable RV variations remains unknown, but an intrinsic effect related to the star itself appears more likely than an extrinsic one. In general, long-period variations in K giants could be caused by non-radial pulsations, spots (rotational modulation of surface features), or substellar companions. Since Aldebaran is photometrically stable at the level of 10 mmag (ESA 1997), large stellar spots that could in principle explain the large-amplitude RV variations cannot be present. Radial pulsations in K giants take place at timescales that are much too short to explain the long periods present in the periodograms; based on the scaling relations from Kjeldsen & Bedding (1995) the radial pulsation period for Aldebaran should be about 6.6 days.
An apparent phase shift between the RV data and the orbital solution, for example around 2004, as described in Hatzes et al. (2015), is also observable in the Lick data as well as in the combined RV data sets. Hatzes et al. (2018) recently described a similar phenomenon for the K5 III star γ Dra, where RV measurements taken between 2003 and 2011 suggest the presence of a planetary companion. New measurements however taken between 2011 and 2017 show that the RV variations abruptly cease around 2011 and return phase-shifted after three years with the same amplitude and period. The authors argue that this could possibly be related to dipole oscillatory convective modes, as proposed by Saio et al. (2015). If this is indeed the case for γ Dra, it could also be considered a likely explanation for Aldebaran, as the two stars are remarkably similar (both have spectral type K5 III, Teff between 3900 and 4000 K, log(g [cm s−2]) between 1.2 and 1.7 and log(L/L⊙) between 2.6 and 2.7, and both stars are nearing the tip of the RGB and will start burning helium soon). Since the predicted period-luminosity relations for dipole oscillatory convection modes roughly match the observations of γ Dra, the same is true for Aldebaran (although the mass of Aldebaran is slightly smaller than that of γ Dra).
We cannot completely rule out the existence of one or multiple substellar companions around the K giant star Aldebaran. However, it seems that neither a single- nor a two-planet solution is a sufficient explanation for the observed RV variations. One could however imagine that stellar variability is superimposed on the RV signal of one planet, but it would still be difficult to explain observed phenomena such as the apparent phase shift around 2004 between the RV data and the Keplerian fit. Nevertheless, the most likely cause seems to be a type of stellar oscillation in highly evolved K giant stars that we have not yet discovered. Further investigations on the nature of giant stars, especially of the evolved K5 III-type giants, are therefore needed before conclusions can be made on the true cause of the RV variations in giant stars like Aldebaran.
Acknowledgments
We kindly thank Hideyuki Saio for making his computed periods of oscillatory convective modes available to us. We also want to thank Paul Heeren for his technical assistance with some parts of the dynamical analysis. This work makes use of Systemic2 (Meschiari et al. 2009). K.R. and S.R. acknowledge the support of the DFG priority program SPP 1992 “Exploring the Diversity of Extrasolar Planets” (SP 345/20-1 and RE 2694/5-1, respectively). S.S. and S.R. further acknowledge support by the DFG Research Unit FOR2544 “Blue Planets around Red Stars” (RE 2694/4-1).
References
- Arenou, F., & Luri, X. 1999, in Harmonizing Cosmic Distance Scales in a Post-HIPPARCOS Era, eds. D. Egret, & A. Heck, ASP Conf. Ser., 167, 13 [NASA ADS] [Google Scholar]
- Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 [NASA ADS] [CrossRef] [Google Scholar]
- Blackwell, D. E., & Lynas-Gray, A. E. 1994, A&A, 282, 899 [NASA ADS] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Cole, A. A., Smecker-Hane, T. A., & Gallagher, III, J. S. 2000, AJ, 120, 1808 [NASA ADS] [CrossRef] [Google Scholar]
- Dawson, R. I., & Fabrycky, D. C. 2010, ApJ, 722, 937 [NASA ADS] [CrossRef] [Google Scholar]
- ESA 1997, in The HIPPARCOS and TYCHO Catalogues, ESA SP, 1200 [Google Scholar]
- Farr, W. M., Pope, B. J. S., Davies, G. R., et al. 2018, ApJ, 865, L20 [NASA ADS] [CrossRef] [Google Scholar]
- Frink, S., Quirrenbach, A., Fischer, D., Röser, S., & Schilbach, E. 2001, PASP, 113, 173 [NASA ADS] [CrossRef] [Google Scholar]
- Frink, S., Mitchell, D. S., Quirrenbach, A., et al. 2002, ApJ, 576, 478 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gatewood, G. 2008, AJ, 136, 452 [NASA ADS] [CrossRef] [Google Scholar]
- Gray, R. O., Corbally, C. J., Garrison, R. F., et al. 2006, ApJ, 132, 161 [Google Scholar]
- Hatzes, A. P., & Cochran, W. D. 1993, ApJ, 413, 339 [NASA ADS] [CrossRef] [Google Scholar]
- Hatzes, A. P., & Cochran, W. D. 1998, MNRAS, 293, 469 [NASA ADS] [CrossRef] [Google Scholar]
- Hatzes, A. P., Cochran, W. D., Endl, M., et al. 2015, A&A, 580, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hatzes, A. P., Endl, M., Cochran, W. D., et al. 2018, AJ, 155, 120 [NASA ADS] [CrossRef] [Google Scholar]
- Hekker, S., & Meléndez, J. 2007, A&A, 475, 1003 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
- Jørgensen, B. R., & Lindegren, L. 2005, A&A, 436, 127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kjeldsen, H., & Bedding, T. R. 1995, A&A, 293, 87 [NASA ADS] [Google Scholar]
- Laughlin, G., & Chambers, J. E. 2002, AJ, 124, 592 [NASA ADS] [CrossRef] [Google Scholar]
- Lomb, N. R. 1976, Ap&SS, 39, 447 [NASA ADS] [CrossRef] [Google Scholar]
- Meschiari, S., Wolf, A. S., Rivera, E., et al. 2009, PASP, 121, 1016 [NASA ADS] [CrossRef] [Google Scholar]
- Miglio, A., Brogaard, K., Stello, D., et al. 2012, MNRAS, 419, 2077 [NASA ADS] [CrossRef] [Google Scholar]
- Mozurkewich, D., Armstrong, J. T., Hindsley, R. B., et al. 2003, AJ, 126, 2502 [NASA ADS] [CrossRef] [Google Scholar]
- Quirrenbach, A., Reffert, S., & Bergmann, C. 2011, in AIP Conf. Ser., eds. S. Schuh, H. Drechsel, & U. Heber, 1331, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Reffert, S., Quirrenbach, A., Mitchell, D. S., et al. 2006, ApJ, 652, 661 [NASA ADS] [CrossRef] [Google Scholar]
- Reffert, S., Bergmann, C., Quirrenbach, A., Trifonov, T., & Künstler, A. 2015, A&A, 574, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Richichi, A., & Roccatagliata, V. 2005, A&A, 433, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saio, H., Wood, P. R., Takayama, M., & Ita, Y. 2015, MNRAS, 452, 3863 [NASA ADS] [CrossRef] [Google Scholar]
- Scargle, J. D. 1982, ApJ, 263, 835 [NASA ADS] [CrossRef] [Google Scholar]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
- Soszyński, I., Udalski, A., Szymański, M. K., et al. 2009, Acta Astron., 59, 239 [Google Scholar]
- Stock, S., Reffert, S., & Quirrenbach, A. 2018, A&A, 616, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tan, X., Payne, M. J., Lee, M. H., et al. 2013, ApJ, 777, 101 [NASA ADS] [CrossRef] [Google Scholar]
- Trifonov, T., Reffert, S., Tan, X., Lee, M. H., & Quirrenbach, A. 2014, A&A, 568, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Walker, G. A. H., Yang, S., Campbell, B., & Irwin, A. W. 1989, ApJ, 343, L21 [NASA ADS] [CrossRef] [Google Scholar]
- Wisdom, J., & Holman, M. 1991, AJ, 102, 1528 [NASA ADS] [CrossRef] [Google Scholar]
- Wood, P. R., Olivier, E. A., & Kawaler, S. D. 2004, ApJ, 604, 800 [NASA ADS] [CrossRef] [Google Scholar]
- Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Comments on Farr et al. (2018)
Farr et al. (2018) fitted a planet to all the available RV data sets of H15 plus their own high-cadence SONG data, while simultaneously applying a Gaussian Process-based Continuous Auto-Regressive Moving Average (CARMA) method to model the acoustic oscillations which happen on much shorter timescales. The planet period which they fitted is 629 days, while the stellar oscillations have their maximum frequency at 2.24 μHz, corresponding to a period of 5.1 days (Farr et al. 2018).
It is of course possible to model any signal with the CARMA formalism. However, the stochastic component must then contain power at the periods of the residuals to the Keplerian component. This is apparently why Farr et al. (2018) achieved a good fit to the data, but is inconsistent with their interpretation of the stochastic component with short-period oscillations as we will show in the following.
In Fig. A.1 we show the full H15 data set supplemented both with the SONG data from Farr et al. (2018) and our own Lick data presented in this paper, together with the Keplerian fit from Farr et al. (2018) (since the values for the periastron time and the longitude of the periastron were not given, we adjusted these parameters ourselves). The lower panel of Fig. A.1 shows the residuals with respect to the Keplerian fit, and Fig. A.2 shows the periodogram of these residuals.
The periodogram of the residuals is quite illustrative. It is clearly visible that there are several period regimes where significant power is still present in the residuals. As expected, there is power at various periods around 5 days, which corresponds to the acoustic oscillations which we did not model here. The power around periods of 5 days is mainly due to the high-cadence SONG data; these peaks in the L-S periodogram are absent in Fig. 4 where the SONG data have not been used.
The largest peaks however occur at periods which are of the same order as the fitted Keplerian period of 629 days. The two strongest peaks in Fig. A.2 are seen at periods of 586 and 766 days, above and below the assumed Keplerian period. The CARMA modelling will not be able to reduce these peaks, since they are located in a completely different period range (around 5 days, as opposed to longer than 500 days, respectively). These additional long periods are also seen in Fig. 4, where the power at the acoustic periods is absent, so that the long additional peaks cannot be aliases resulting from combinations of similar shorter periods.
However, these peaks that are of the same order as the assumed orbital period are the ones that challenge the existence of the planet. They are responsible for the observed deviations from a strictly periodic Keplerian orbit, such as an apparent phase shift in 2004 or the decrease in the L-S power of the assumed orbital period in 2006/2007, which we discussed in the main body of the paper and led us to conclude that the planet claimed to orbit Aldebaran might not exist.
Thus, the modelling of the acoustic modulations does not help with the interpretation of the much longer periods in the RV data of Aldebaran, and it will also not enable higher sensitivity for planet detection around giant stars in general (in particular since most Lick stars are much less evolved than Aldebaran and thus have even shorter acoustic oscillation periods, while at the same time lacking high cadence SONG data).
![]() |
Fig. A.1. All available RV data sets for Aldebaran, i.e. all the H15 data sets plus the ones from Lick (in blue) and Song (in black), with the planet fit from Farr et al. (2018) overplotted in the upper panel. Lower panel: residuals to that fit. As one can see, the phase of the periodic planet signal does not match during the time frame from about 2002 to 2005, shortly after the long gap in the overall RV time series. This mismatch can also be seen in the residuals plot; during that same time frame from about 2002 to 2005 the residuals are larger than at other times. |
Bootstrap uncertainties for the one-planet solution for the Lick data and for the combined RV data of Lick and H15, as well as for the two-planet solution.
![]() |
Fig. A.2. Periodogram of the residuals to the planet fit from Farr et al. (2018), shown in Fig. A.1. One can clearly see that there is excess power at periods around 5–6 days as well as even higher and very significant power at periods of 550–800 days, close to the period of 629 days attributed to the putative planet. |
Appendix B: Radial-velocity data from Lick Observatory for Aldebaran
Radial-velocity measurements with error from Lick Observatory for Aldebaran.
All Tables
Significant periods and their FAPs in the L-S periodogram of the Lick data (left) and Lick + H15 data (right).
Comparison of the parameters of our orbital solution for the Lick data, and for the combined data, respectively, and those of H15 for the substellar companion to Aldebaran.
Planetary parameters of the orbital solution for two substellar companions to Aldebaran for the Lick + H15 data sets.
Bootstrap uncertainties for the one-planet solution for the Lick data and for the combined RV data of Lick and H15, as well as for the two-planet solution.
All Figures
![]() |
Fig. 1. Top: L-S periodogram for the Lick data. The most significant peak is at a period of 587.84 d with a FAP of 1.6 × 10−17. Middle: L-S periodogram for the H15 data. The most significant period is 629.08 d and has a FAP of less than 10−41. Bottom: L-S periodogram for the Lick data combined with the data of H15. The most significant peak is found at a period of 625.02 d with a FAP of less than 10−41. In all three panels, the significance thresholds represented by the vertical dashed lines correspond to FAPs of 0.1%, 1%, and 5%, respectively. |
In the text |
![]() |
Fig. 2. Top left: RV measurements for Aldebaran from Lick Observatory (165 in total). The orbital solution for a substellar companion with a period of ∼594 d is represented by the black solid line. Bottom left: residual RV variations for the Lick RV measurements of Aldebaran after removing the orbital solution (black dashed line) for the Lick data of Table 3. Top right: RV measurements for Aldebaran from Lick Observatory. The orbital solution of H15 for a substellar companion with a period of 628.96 ± 0.90 d is overplotted and is represented by the black solid line. Bottom right: residual RV variations for the Lick RV measurements of Aldebaran after removing the orbital solution (black dashed line) of H15. |
In the text |
![]() |
Fig. 3. RV measurements for Aldebaran from Lick Observatory and all the data sets of H15. The orbital solution for a substellar companion with a period of 607.22 d is represented by the black solid line in the top panel. Lower panel: residuals with respect to this one-planet fit. |
In the text |
![]() |
Fig. 4. Lomb-Scargle periodogram for the Lick data combined with the data of H15 after removing the 607.22 d period of the orbital solution. The alias period 230.66 d has disappeared as well. The significance thresholds represented by the vertical dashed lines correspond to FAPs of 0.1%, 1%, and 5%, respectively. |
In the text |
![]() |
Fig. 5. Lomb-Scargle power of the 625 d period as a function of the number of data points (N). The red squared data points represent the real data while the blue circular data points represent simulated data. Top: growth of the Lomb-Scargle power from 1980–2013. Data points were added in steps of 20 in chronological order. A decrease in statistical power is observable for the real data between N = 300 and N = 320. Bottom: growth of the Lomb-Scargle power from 2013–1980. Data points were added in steps of 20 in reversed chronological order. Before N = 240 (which represent the newest RV data) the ∼620 d period is clearly visible as the growth in power shows the same linear slope as the simulated data. The older RV data (taken before 2006) however cause a decrease in the slope and a drop in power at N = 340. |
In the text |
![]() |
Fig. 6. Growth of the L-S power in the GLS as a function of the number of data points (N). Data points were added in steps of 20 in chronological order. The red squared data points represent the ∼620 d period, while the green triangles represent the ∼750 d period and the yellow triangles the ∼550 d period. The blue circular data points represent simulated data. See the text for more details. |
In the text |
![]() |
Fig. 7. RV measurements for Aldebaran from Lick Observatory and all data sets of H15. The orbital solution for two substellar companions with periods of 614.10 d and 772.83 d is represented by the black solid line in the top panel. Lower panel: residuals with respect to this two-planet fit. |
In the text |
![]() |
Fig. 8. Period-luminosity diagram to illustrate the position of Aldebaran and γ Dra (star symbols) with respect to published exoplanets from our Lick survey (red circles, labeled with HIP number) and candidate exoplanets (blue circles). The small dots are stars with periods identified by the OGLE survey in the LMC; the subset of black dots are those stars identified as showing long secondary periods in the OGLE photometry. The thick solid lines show model calculations by Saio et al. (2015) for the periods of oscillatory convective modes for stars of various masses, and for a mixing length parameter of 1.2. Aldebaran and γ Dra’s periods are rather close to the predicted periods of oscillatory convective modes, especially if one were to extrapolate the models for 1.0 and 1.3 M⊙. |
In the text |
![]() |
Fig. A.1. All available RV data sets for Aldebaran, i.e. all the H15 data sets plus the ones from Lick (in blue) and Song (in black), with the planet fit from Farr et al. (2018) overplotted in the upper panel. Lower panel: residuals to that fit. As one can see, the phase of the periodic planet signal does not match during the time frame from about 2002 to 2005, shortly after the long gap in the overall RV time series. This mismatch can also be seen in the residuals plot; during that same time frame from about 2002 to 2005 the residuals are larger than at other times. |
In the text |
![]() |
Fig. A.2. Periodogram of the residuals to the planet fit from Farr et al. (2018), shown in Fig. A.1. One can clearly see that there is excess power at periods around 5–6 days as well as even higher and very significant power at periods of 550–800 days, close to the period of 629 days attributed to the putative planet. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.