Advanced
PLANETARY COMPANION IN K GIANT σ PERSEI*
PLANETARY COMPANION IN K GIANT σ PERSEI*
Journal of The Korean Astronomical Society. 2014. Mar, 47(2): 69-76
Copyright © 2014, null
  • Received : February 19, 2014
  • Accepted : March 04, 2014
  • Published : March 27, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
BYEONG-CHEOL LEE
Astronomy and Space Science Major, University of Science and Technology, Gajeong-ro Yuseong-gu, Daejeon 305-333, Korea
INWOO HAN
Korea Astronomy and Space Science Institute 776, Daedeokdae-ro, Yuseong-gu, Daejeon 305-348, KoreaE-mail:bclee@kasi.re.kr,iwhan@kasi.re.kr,tlotv@kasi.re.kr,kmkim@kasi.re.kr
MYEONG-GU PARK
Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Daegu 702-701, KoreaE-mail:mgp@knu.ac.kr
DAVID E. MKRTICHIAN
Crimean Astrophysical Observatory, Taras Shevchenko National University of Kyiv, 98409, Nauchny, Crimea, UkraineE-mail:davidmkrt@gmail.com
GWANGHUI JEONG
Astronomy and Space Science Major, University of Science and Technology, Gajeong-ro Yuseong-gu, Daejeon 305-333, Korea
KANG-MIN KIM
Korea Astronomy and Space Science Institute 776, Daedeokdae-ro, Yuseong-gu, Daejeon 305-348, KoreaE-mail:bclee@kasi.re.kr,iwhan@kasi.re.kr,tlotv@kasi.re.kr,kmkim@kasi.re.kr
GENNADY VALYAVIN
Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnii Arkhyz, 369167, RussiaE-mail:gvalyavin@sao.ru

Abstract
We report the detection of an exoplanet candidate in orbit around σ Persei from a radial velocity (RV) survey. The system exhibits periodic RV variations of 579.8 ± 2.4 days. The purpose of the survey is to search for low-amplitude and long-period RV variations in giants and examine the origin of the variations using the fiber-fed Bohyunsan Observatory Echelle Spectrograph installed at the 1.8-m telescope of Bohyunsan Optical Astronomy Observatory in Korea. We present high-accuracy RV measurements of σ Per made from December 2003 to January 2014. We argue that the RV variations are not related to the surface inhomogeneities but instead a Keplerian motion of the planetary companion is the most likely explanation. Assuming a stellar mass of 2.25 ± 0.5 M , we obtain a minimum planetary companion mass of 6.5 ± 1.0 M Jup , with an orbital semi-major axis of 1.8 ± 0.1 AU, and an eccentricity of 0.3 ± 0.1 around σ Per.
Keywords
1. INTRODUCTION
Ever since the discovery of the first exoplanet around a main-sequence (MS) star by Mayor & Queloz (1995) using the radial velocity (RV) method, approximately 1000 exoplanets have been found so far. It has been demonstrated that the RV method is more efficient for late spectral type stars having an ample number of spectral lines. Planet searches generally focus on late F, G, and K low-mass dwarfs. Stellar and hence planet masses in such systems can be precisely determined. The RV search for exoplanets around early class B and A-type stars is difficult due to the high-rotational velocities of these stars and the low accuracy of RV measurements.
For giant stars, on the other hand, RV shifts can be measured with high precision due to their sharp spectral lines. Lack of knowledge about planet formation makes RV surveys of giant stars an important endeavor. Compared to MS stars, however, giants may show more complex RV variations due to the influence of various surface processes on the line profiles, such as chromospheric activities, spots, large convection cells, and stellar pulsations. It is thus difficult to identify the origin of RV variations in giant star systems. Also, it is difficult to determine accurate masses of the host giant stars, because giant stars of different masses are located in roughly the same region of the H-R Diagram. Finally, there appears to exist a selection bias in that only massive planets with relatively long orbital periods have been detected around giant stars, whereas surveys of MS stars tend to find low-mass planets with short orbital periods.
Frink et al. (2002) discovered the first planetary companion around the K giant star ι Dra (K2). Since then, several companions around K giant stars have been detected through the RV method. Recognizing the importance of exoplanet formation and evolution, several exoplanet survey groups began exploring evolved K giants (Frink et al. 2002; Setiawan et al. 2003; Hatzes et al. 2005; D¨olinger et al. 2007; Han et al. 2010).
We initiated a precise RV survey using the 1.8-m telescope at Bohyunsan Optical Astronomy Observatory (BOAO) in 2003 to search for exoplanets and to study the oscillations for 6 F giants, 55 K giants, and 10 M giants. In the course of the program we have reported the discovery of six new exoplanets (Han et al. 2010; Lee et al. 2011, 2012b, 2012c, 2013), an exoplanet candidate (Lee et al. 2012a), and confirmed the exoplanet around Pollux (Han et al. 2008). The K giant σ Per is one of the 55 K giants in our sample for which we have obtained precise RV measurements in the past 11 years at BOAO. Here we report the detection of low-amplitude and long-period RV variations for σ Per, possibly caused by a planetary companion. In Section 2, we describe the observations and data reduction. In Section 3, the stellar characteristics of the host star are derived. The RV variation measurements and possible origins are presented in Section 4. Finally, in Section 5, we discuss the main results from our study.
2. OBSERVATIONS AND REDUCTION
The observations were carried out as a part of the RV variation study of 55 K giants. We used the fiberfed, high-resolution ( R = 90000) Bohyunsan Observatory Echelle Spectrograph (BOES; Kim et al. 2007) installed at the 1.8-m telescope of BOAO, Korea. The BOES spectra covered the wavelength coverage from 3500 to 10500 Å. To provide precise RV measurements, an iodine absorption (I 2 ) cell was used with a wavelength region of 4900−6000 Å. The characteristic signal-to-noise (S/N) for the I 2 region was about 200 at typical exposure times ranging from 5 to 15 minutes.
We recorded 71 spectra for σ Per from December 2003 to January 2014 (52 nights in total). The basic reduction of spectra was performed using the IRAF software package and the DECH (Galazutdinov 1992) code. The precise RV measurements related to the I 2 analysis were undertaken using the RVI2CELL (Han et al. 2007) code, which is based on a method by Butler et al. (1996) and Valenti et al. (1995). However, for the modeling of the instrument profile, we used the matrix formula described by Endl et al. (2000). We solved the matrix equation using singular value decomposition instead of the maximum entropy method adopted by Endl et al. (2000).
The long-term stability of the BOES was demonstrated by observing the RV standard star τ Ceti, which was constant with an rms scatter of 6.8 m s −1 (Lee et al. 2013). The RV measurements for σ Per are listed in Table 1 and shown in Fig. 1 .
PPT Slide
Lager Image
— RV measurements (top) and rms scatter of the residuals (bottom) for 𝜎 Per taken at BOAO from December 2003 to January 2014. The solid line is the orbital solution with a period of 579.8 days and an eccentricity of 0.3.
3. STELLAR CHARACTERISTICS
Evolved stars are known to exhibit low-amplitude and long-period RV variations produced by rotation and surface activity such as spots, plages, or filaments. The investigation of the stellar characteristics thus is important in discerning the nature of any RV variations.
- 3.1 Fundamental Parameters
We assumed the main photometric parameters for σ Per (= HD 21552 = HR 1052 = HIP 16335) from the HIPPARCOS catalog (ESA 1997). σ Per is a cool evolved star of type K3 III and has an apparent magnitude V = 4.36. The astrometric parallax (π) was taken from the improved HIPPARCOS value by van Leeuwen (2007), 9.07 ± 0.26 mas. Other stellar parameters are taken from Luck & Heiter (2007): T eff = 4367 K, [Fe/H] = – 0.22 ± 0.23, and log g = 1.8 ± 0.1. We also determined atmospheric parameters directly from our spectra by measuring 274 equivalent widths (EW) of Fe I and Fe II lines. We estimated T eff = 4201 ± 22 K, [Fe/H] = – 0.21 ± 0.06, log g = 1.77 ± 0.09, and vmicro = 1.5 ± 0.1 using the program TGVIT (Takeda et al. 2005). Our estimated values were in agreement with Luck & Heiter (2007). Table 2 summarizes the basic stellar parameters of σ Per.
- 3.2 Stellar Mass
Due to disagreement in the mass of the host star as measured by different authors and methods, the planet mass is subject to some uncertainty. In order to derive stellar mass, one needs stellar parameters, such as the effective temperature. Luck & Heiter (2007) applied two basic methods of mass determination for nearby giants including σ Per: photometric and spectroscopic. The photometric mass was derived using the isochrones of Bertelli et al. (1994), interpolated using the absolute V magnitude and photometric temperature. The authors then derived the luminosity and surface gravity using standard relations. The spectroscopic mass was determined using the direct spectroscopic value for the surface gravity and the luminosity using standard relations. Both mass estimates are shown in Table 2 . The mass derived from the spectroscopic data is larger than that derived from the photometry. The difference is due to the larger surface gravity measured by the spectroscopic analysis, which demands larger mass. The use of different temperatures (i.e., photometric vs. spectroscopic) has little effect on the physical determination of the mass. We considered results from both the photometric and spectroscopic approaches.
RV measurements for 𝜎 Per from December 2003 to January 2014 using the BOES.
PPT Slide
Lager Image
RV measurements for 𝜎 Per from December 2003 to January 2014 using the BOES.
Stellar parameters for 𝜎 Per assumed in the present paper
PPT Slide
Lager Image
References. 1 HIPPARCOS (ESA 1997); 2a Luck & Heiter (2007) estimated by spectroscopic method; 2b Luck & Heiter (2007) estimated by photometric method; 3 van Leeuwen (2007); 4 This work; 5 Pasinetti-Fracassini et al. (2001); 6 McDonald et al. (2012); 7 de Medeiros & Mayor (1999); 8 Hekker & Mel´endez (2007)
- 3.3 Rotational Velocity
In evolved stars, low-amplitude and long-period RV variations may arise from the rotation (Lee et al. 2008, 2012a). Thus, the rotational period is very important in distinguishing between RV variations and rotational modulation. In order to estimate the stellar rotational velocity, we used a line-broadening model by Takeda et al. (2008). The observed stellar spectrum was fitted by convolving the intrinsic spectrum model and a total macrobroadening function. The intrinsic spectrum model can be computed if a model atmosphere characterized by the microturbulence, v micro , and an elemental abundance are given. The total macrobroadening function is a convolution of three component functions: the instrumental broadening dispersion ( v ip ), the rotational broadening velocity ( v rot ), and the macroturbulence ( v macro ). These component functions broaden spectral lines without altering their EW. Distinguishing between the broadening components is difficult for spectra of slowly rotating late type stars because they have similar intrinsic line profiles.
For the determination of line broadening, we used the automatic spectrum-fitting technique (Takeda 1995) for the spectrum within the wavelength range 6080− 6089 Å. This technique employs seven free parameters that specify the best fitting solution: the abundances of six elements (Si, Ti, V, Fe, Co, and Ni) and the total macrobroadening parameter ( v M ) which is given as
PPT Slide
Lager Image
where vip can be described as
PPT Slide
Lager Image
With ( c/R ) = 2 km s −1 in the case of R = 90000.
We used the Gray (1989) assumption that the macroturbulence depends only on the surface gravity. We thus use
PPT Slide
Lager Image
We know the value of v M from the observed spectrum. So the rotational velocity can be evaluated from Eq. 1.
In this work, we measured v M to be 3.9 km s −1 and estimated the rotational velocity equal to 1.3 km s −1 . In the literature, the v rot sin( i ) values are taken from various sources obtained by different methods. De Medeiros & Mayor (1999) measured rotational velocities using the cross-correlation technique on spectra obtained from the CORAVEL spectrometer with an uncertainty of about 1.0 km s −1 . Hekker & Meléndez (2007) determined rotational velocities from FWHM measurements using the Coudé Auxiliary Telescope (CAT) at Lick observatory. Lastly, Massarotti et al. (2008) used the cross-correlation of the observed spectra against templates drawn from a library of synthetic spectra calculated by Kurucz (1992) for different stellar atmospheres, even though σ Per was not calculated in the study.
As stated above, estimated rotational velocities depend critically on the adopted templates and methods used for the calibration of the line broadening. Thus it would be inappropriate to choose any one value as the “correct” one. Using these determinations, we can define a range of v rot sin( i ) measurements of 1.0−2.21 km s −1 . Based on the rotational velocities and the stellar radius, we derived the range of upper limits for the rotational period of
Prot = 2πR/[vrot sin(i)] = 641−1288 days.
4. RADIAL VELOCITY VARIATIONS AND ORIGIN
Evolved stars exhibit pulsations as well as surface activity, which result in low-amplitude RV variability on different time scales. Generally, short-term RV variations have been known to be the result of stellar pulsations (Hatzes & Cochran 1998), whereas long-term (hundreds of days) variations with a low-amplitude may be caused by three kinds of phenomena: planetary companions, stellar oscillations, or rotational modulations. To establish the origin of the RV variabilities for σ Per, we examined 1) the orbital fit, 2) the stellar chromospheric activity, 3) the HIPPARCOS photometry, and 4) the spectral line bisectors.
- 4.1 Orbital Solutions
Period analysis was applied to assess the significance of periodic trends or variations of the time series. The Lomb-Scargle periodogram (Lomb 1976; Scargle 1982), which is a useful tool for investigating long-period variations for unequally spaced data, was used to calculate the RV time series. Our RV measurements for σ Per are listed in Table 1 and shown in Fig. 1 (top panel). The Lomb-Scargle periodogram for σ Per shows a significant power peak at f 1 =0.00172 c d −1 and no significant peak in the residuals after subtracting the major period ( Fig. 2 ). The best fit Keplerian orbit is found with to have a period P = 579.8 ± 2.4 days, a semi-amplitude K = 96.0 ± 6.5 m s −1 , and an eccentricity e = 0.3 ± 0.1. Assuming a stellar mass of 2.25 ± 0.5 M (Luck & Heiter 2007), which includes the range of mass (1.82 M ) estimated by the photometric method, we derived a minimum mass of a planetary companion m sin( i ) = 6.5 ± 1.0 M Jup at a distance a = 1.8 AU from the host. The RV phase diagram for the orbit is shown in Fig. 3 . All the orbital elements are listed in Table 3 .
PPT Slide
Lager Image
— The Lomb-Scargle periodogram of the RV mea- surements for 𝜎 Per. The periogram shows a significant power peak at a frequency of 0.00172 c d−1 corresponding to a period of 579.8 days (top) and no significant power peak at the residual (bottom).
PPT Slide
Lager Image
— Phase diagram for 𝜎 Per. The solid line is phased with the orbital solution that fits the data with a period of 579.8 days and an eccentricity of 0.3.
Orbital parameters for 𝜎 Per b.
PPT Slide
Lager Image
Orbital parameters for 𝜎 Per b.
The statistical significance of this signal was calculated using the bootstrap randomization technique (Kürster et al. 1999). The measured RV values were randomly shuffled keeping the observed times fixed. The Lomb-Scargle periodogram was then computed for 200000 “random” datasets. The fraction of periodograms having power higher than the amplitude peak in the range of 0 < f < 0.02 c d −1 represents the false-alarm probability (FAP) that noise would create the detected signal. A small FAP is therefore an indication of a significant peak. We found that among 200000 trials the highest peak of the randomly generated periodograms never exceeded that in the original data, thus the FAP is < 10 −5 . We could not find any significant power peaks in the residual periodogram (FAP of 4 × 10 −2 ).
- 4.2 Chromospheric Activities
The EW variations of Ca II H & K, Hα, and Ca II 8662 Å lines are frequently used a chromospheric activity indicators because they are sensitive to the stellar atmospheric activity. Such activity could have a decisive effect on the RV variations. While the Hα is sensitive to the atmospheric stellar activity (Kürster et al. 2003) and useful to measure the variations, it is not easy to estimate the Hα EW due to line blending in stellar spectra and telluric lines if any. Another activity indicator, the Ca II 8662 Å line, is also not suitable because significant fringing and saturation of our CCD spectra is seen at wavelengths longer than 7500 Å.
The emission of the Ca II H & K comes from the chromosphere and the core of the spectral line often exhibits central reversal in the presence of chromospheric activity (Saar & Donahue 1997). It is common in cool stars and is intimately connected to the convective envelope and magnetic activity. Unfortunately, the spectra are not clear enough to resolve the emission feature in the Ca II K line core for σ Per. Instead, we use the relatively stronger Ca II H region shown in Fig. 4. There are weak central emission lines with weak reversal in the core of Ca II H.We used the averaged Ca II H line profiles at the upper, lower, and near zero part of the RV curve to identify any systematic difference related to the RV variations. For σ Per, these correspond to the positive deviation on JD-2454038.348295 (RV = 123.6 m s −1 ), the negative deviation on JD- 2455454.324212 (RV = − 91.8 m s −1 ), and the zero deviation on JD-2455932.979094 (RV = 18.3 m s −1 ) of the RV curve. Core emissions and reversals in the Ca II H central region are very weak and there are no systematic differences in the line profiles at these three extremal parts of the RV curves. This means that σ Per, at the moment of observations, exhibited at most a modest chromospheric activity.
- 4.3 Photometric Variations
We analyzed the HIPPARCOS photometry data to search for possible brightness variations which might be caused by the rotational modulation of cool stellar spots. The available photometry database comprises 77 HIPPARCOS measurements for σ Per from January 1990 to February 1993. They maintained photometric stability down to rms scatters of 0.006 magnitude, which correspond to 0.14% variations over the time span of the observations. Fig. 5 shows no significant photometric variations.
PPT Slide
Lager Image
— The Ca II H spectral region for 𝜎 Per at different phases. The Ca II H core features exhibit weak emission at the line center. Line profiles do not show any significantly different features related to the RV variations.
PPT Slide
Lager Image
— The HIPPARCOS photometric variations from January 1990 to February 1993 for 𝜎 Per.
- 4.4 Line Bisector
The RV variations due to planetary companions should not produce any changes in the spectral line shape whereas the surface inhomogeneities do. Stellar rotational modulations of surface inhomogeneities can create variable asymmetries in the spectral line profiles. Thus, the analysis of the shapes of the spectral lines, in particular the analysis of the line bisectors, may identify the origin of observed RV variations and prove or disprove the existence of planetary companions (e.g., Queloz et al. 2001).
The RV differences of the central values at high and low flux levels of the line profiles are defined as a bisector velocity span (BVS). To calculate the BVS, following Hatzes et al. (2005) and Lee et al. (2013), we selected the unblended strong line of Ni I 6643.6Å that is located beyond the I 2 absorption region. We estimated the BVS of the profile between two different flux levels at 0.8 and 0.4 of the central depth as the span points, which avoided the spectral core and wings where the differences in the bisector measurements are large and noisy. The BVS variations for σ Per as a function of RV are shown in Fig. 6 . There is no obvious correlation between the RV and BVS measurements.
PPT Slide
Lager Image
— Relation between the RV and BVS variations for the Ni I 6643.6 Å spectral line. The solid line marks the slope of the spectral line.
5. DISCUSSION
Our analysis of the RV measurements for σ Per revealed a 579 day periodic variation that persisted between 2003 and 2014 for almost six and half cycles. In evolved stars, low-amplitude and long-term RV variations, such as in hundreds of days, may be caused by stellar pulsations, rotational modulations of inhomogeneous surface features, or planetary companions. Of these, RV variations due to pulsations are unlikely because the period of the fundamental radial mode pulsation is of several days and two orders of magnitude longer compared to observed RV variations. The Ca II H line was used to monitor chromospheric activities, but does not show any significant variations. The HIPPARCOS photometry does not show any detectable variations nor the BVSs show any correlation with the RV measurements. All these null results from the analysis of the chromospheric activities or surface inhomogeneities argue against the rotational modulation of surface activities as the origin of observed RV variations and in favor of the planetary companion hypothesis.
Our precise RV measurements for σ Per have established the presence of variations with a 579.8 ± 2.4 day period. Due to the uncertainty in the mass of the host star, the planetary companion mass can range from 5.5−7.5 M Jup with an orbital semi-major axis of 1.8 ± 0.1 AU, and an eccentricity of 0.3 ± 0.1. The planetary companion to σ Per has similar properties to other companions around K giant stars. They are mostly massive planets with the average mass of ~ 5 M Jup and an orbital radius of ~ 1.5 AU. The low metallicity of − 0.2 for σ Per is also consistent with the average value estimated for other planet-hosting stars ([Fe/H] ~ − 0.12).
The residual RV variations after removing the 579.8- day signal are 33.7 m s −1 . This value of scatter is significantly greater than the RV precision for the RV standard star τ Ceti (6.8 m s −1 ) and the typical internal error of individual RV accuracies for K giants from our survey of ~ 7.8 m s −1 . The significant residual appears to be caused by systematics in RV measurement. Otherwise, the variations, whether periodic or not, are typical in evolved K giants (Setiawan et al. 2003; Hatzes et al. 2005; Dölinger et al. 2007; de Medeiros et al. 2009; Han et al. 2010; Lee et al. 2012c) and can be possibly attributed to pulsational and convection activity. Furthermore, residual RV variations tend to increase toward later spectral types.
Acknowledgements
BCL acknowledges partial support by the KASI (Korea Astronomy and Space Science Institute) grant 2014-1-400-06. MGP was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education, Science and Technology (2012R1A1A4A010 13596). DEM acknowledges his work as part of the research activity of the National Astronomical Research Institute of Thailand (NARIT), which is supported by the Ministry of Science and Technology of Thailand. This research made use of the SIMBAD database, operated at the CDS, Strasbourg, France. We thank the developers of the Bohyunsan Observatory Echelle Spectrograph (BOES) and all staff of the Bohyunsan Optical Astronomy Observatory (BOAO).
References
Bertelli G. , Bressan A. , Chiosi C. , Fagotto F. , Nasi E. 1994 Theoretical Isochrones from Models with New Radiative Opacities A&AS 106 275 -
Butler R. P. , Marcy G. W. , Williams E. , McCarthy C. , Dosanjh P. , Vogt S. S. 1996 Attaining Doppler Precision of 3 ms-1 PASP 108 500 -    DOI : 10.1086/133755
de Medeiros J. R. , Mayor M. 1999 A Catalog of Rotational and Radial Velocities for Evolved Stars A&AS 139 433 -    DOI : 10.1051/aas:1999401
de Medeiros J. R. , Setiawan J. , Hatzes A. P. , Pasquini L. , Girardi L. , Udry S. , Döllinger M. P. , da Silva L. 2009 A Planet around the Evolved Intermediate-Mass Star HD 110014 A&A 504 617 -    DOI : 10.1051/0004-6361/200911658
Döllinger M. P. , Hatzes A. P. , Pasquini L. , Guenther E. W. , Hartmann M. , Girardi L. , Esposito M. 2007 Discovery of a Planet around the K Giant Star 4 Ursae Majoris A&A 472 649 -    DOI : 10.1051/0004-6361:20066987
Endl M. , Kürster M. , Els S. 2000 The Planet Search Program at the ESO Coudé Echelle Spectrometer. I. Data Modeling Technique and Radial Velocity Precision Tests A&A 362 585 -
Frink S. , Mitchell D. S. , Quirrenbach A. , Fischer D. A. , Marcy G.W. , Butler R. P. 2002 Discovery of a Substellar Companion to the K2 III Giant ι Draconis ApJ 576 478 -    DOI : 10.1086/341629
Galazutdinov G. A. 1992 Special Astrophysical Observatory Preprint 92 SAO Nizhnij Arkhyz
Gray D. F. 1989 The Rotational Break for G Giants ApJ 347 1021 -    DOI : 10.1086/168192
Han I. , Kim K.-M. , Lee B.-C. , Valyavin G. 2007 Development of RVI2CELL -A Precise Radial Velocity Estimation Program with BOES Data PKAS 22 75 -
Han I. , Lee B.-C. , Kim K.-M. , Mkrtichian D. E. 2008 Confirmation of the Exoplanet around 𝛽 GEM from the RV Observations Using BOES JKAS 41 59 -
Han I. , Lee B.-C. , Kim K. M. , Mkrtichian D. E. , Hatzes A. P. , Valyavin G. 2010 Detection of a Planetary Companion around the Giant Star γ1 Leonis A&A 509 A24 -    DOI : 10.1051/0004-6361/200912536
Hatzes A. P. , Cochran W. D. 1998 Stellar Oscillations in K Giant Stars ASPC 154 311 -
Hatzes A. P. , Guenther E. W. , Endl M. , Cochran W. D. , Döllinger M. P. , Bedalov A. 2005 A Giant Planet around the Massive Giant Star HD 13189 A&A 437 743 -    DOI : 10.1051/0004-6361:20052850
Hekker S. , Meléndez J. 2007 Precise Radial Velocities of Giant Stars. III. Spectroscopic Stellar Parameters A&A 475 1003 -    DOI : 10.1051/0004-6361:20078233
Kim K.-M. , Han I. , Valyavin G. G. , Plachinda S. , Jang J. G. , Jang B.-H. , Seong H. C. , Lee B.-C. , Kang D.-I. , Park B.-G. , Yoon T. S. , Vogt S. S. 2007 The BOES Spectropolarimeter for Zeeman Measurements of Stellar Magnetic Fields PASP 119 1052 -    DOI : 10.1086/521959
Kürster M. , Hatzes A. P. , Cochran W. D. 1999 Precise Radial Velocities of Proxima Centauri. Strong Constraints on a Substellar Companion A&A 344 L5 -
Kürster M. , Endl M. , Rouesnel F. 2003 The Low-Level Radial Velocity Variability in Barnard’s Star (= GJ 699). Secular Acceleration, Indications for Convective Redshift, and Planet Mass Limits A&A 403 1077 -    DOI : 10.1051/0004-6361:20030396
Kurucz R. L. 1992 Model Atmospheres for Population Synthesis, The Stellar Populations of Galaxies IAUS 149 225 -
Lee B.-C. , Mkrtichian D. E. , Han I. , Park M.-G. , Kim K.-M 2008 Precise Radial Velocities of Polaris: Detection of Amplitude Growth AJ 135 2240 -    DOI : 10.1088/0004-6256/135/6/2240
Lee B.-C. , Mkrtichian D. E. , Han I. , Kim K.-M. , Park M.-G. 2011 A Likely Exoplanet Orbiting the Oscillating K-Giant α Arietis A&A 529 A134 -    DOI : 10.1051/0004-6361/201016293
Lee B.-C. , Han I. , Park M.-G. , Kim K.-M. , Mkrtichian D. E. 2012 Detection of the 128-Day Radial Velocity Variations in the Supergiant α Persei. Rotational Modulations, Pulsations, or a Planet? A&A 543 A37 -    DOI : 10.1051/0004-6361/201118539
Lee B.-C. , Han I. , Park M.-G. , Mkrtichian D. E. , Kim K.-M 2012 A Planetary Companion around the K Giant 𝜖 Corona Borealis A&A 546 A5 - 5    DOI : 10.1051/0004-6361/201219347
Lee B.-C. , Mkrtichian D. E. , Han I. , Park M.-G. , Kim K.-M 2012 Detection of an Exoplanet around the Evolved K Giant HD 66141 A&A 548 A118 -    DOI : 10.1051/0004-6361/201118014
Lee B.-C. , Han I. , Park M.-G 2013 Planetary Companions Orbiting M Giants HD 208527 and HD 220074 A&A 549 A2 -    DOI : 10.1051/0004-6361/201220301
Lomb N. R. 1976 Least-Squares Frequency Analysis of Unequally Spaced Data Ap&SS 39 447 -    DOI : 10.1007/BF00648343
Luck R. E. , Heiter U. 2007 Giants in the Local Region AJ 133 2464 -    DOI : 10.1086/513194
Massarotti A. , Latham D. W. , Stefanik R. P. , Fogel J. 2008 Rotational and Radial Velocities for a Sample of 761 HIPPARCOS Giants and the Role of Binarity AJ 135 209 -    DOI : 10.1088/0004-6256/135/1/209
Mayor M. , Queloz D. 1995 A Jupiter-Mass Xompanion to a Aolar-Rype Star Nature 378 355 -    DOI : 10.1038/378355a0
McDonald I. , Zijlstra A. A. , Boyer M. L. 2012 Fundamental Parameters and Infrared Excesses of Hipparcos Stars MNRAS 427 343 -    DOI : 10.1111/j.1365-2966.2012.21873.x
Pasinetti-Fracassini L. E. , Pastori L. , Covino S. , Pozzi A. 2001 Catalogue of Apparent Diameters and Absolute Radii of Stars (CADARS) - Third Edition - Comments and Statistics A&A 367 521 -    DOI : 10.1051/0004-6361:20000451
Queloz D. , Henry G. W. , Sivan J. P. 2001 No Planet for HD 166435 A&A 379 279 -    DOI : 10.1051/0004-6361:20011308
Saar S. H. , Donahue R. A. 1997 Activity-Related Radial Velocity Variation in Cool Stars ApJ 485 319 -    DOI : 10.1086/304392
Scargle J. D. 1982 Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data ApJ 263 835 -    DOI : 10.1086/160554
Setiawan J. , Pasquini L. , da Silva L. , von der Lühe O. , Hatzes A. 2003 Precise Radial Velocity Measurements of G and K Giants. First Results A&A 397 1151 -    DOI : 10.1051/0004-6361:20021559
Takeda Y. , Ohkubo M. , Sato B. , Kambe E. , Sadakane K. 2005 Spectroscopic Study on the Atmospheric Parameters of Nearby F-K Dwarfs and Subgiants PASJ 57 27 -
Takeda Y. , Sato B. , Murata D. 2008 Stellar Parameters and Elemental Abundances of Late-G Giants PASJ 60 781 -
Valenti J. A. , Butler R. P. , Marcy G. W. 1995 Determining Spectrometer Instrumental Profiles Using FTS Reference Spectra PASP 107 966 -    DOI : 10.1086/133645
van Leeuwen F. 2007 Validation of the New Hipparcos Reduction A&A 474 653 -    DOI : 10.1051/0004-6361:20078357