Precise radial velocities of giant stars. X. Bayesian stellar parameters and evolutionary stages for 372 giant stars from the Lick planet search

The determination of accurate stellar parameters of giant stars is essential for our understanding of such stars in general and as exoplanet host stars in particular. Precise stellar masses are vital for determining the lower mass limit of potential substellar companions with the radial velocity method. Our goal is to determine stellar parameters, including mass, radius, age, surface gravity, effective temperature and luminosity, for the sample of giants observed by the Lick planet search. Furthermore, we want to derive the probability of these stars being on the horizontal branch (HB) or red giant branch (RGB), respectively. We compare spectroscopic, photometric and astrometric observables to grids of stellar evolutionary models using Bayesian inference. We provide tables of stellar parameters, probabilities for the current post-main sequence evolutionary stage, and probability density functions for 372 giants from the Lick planet search. We find that $81\%$ of the stars in our sample are more probably on the HB. In particular, this is the case for 15 of the 16 planet host stars in the sample. We tested the reliability of our methodology by comparing our stellar parameters to literature values and find very good agreement. Furthermore, we created a small test sample of 26 giants with available asteroseismic masses and evolutionary stages and compared these to our estimates. The mean difference of the stellar masses for the 24 stars with the same evolutionary stages by both methods is only $\langle\Delta M\rangle=0.01\pm0.20\;\mathrm{M_\odot}$. We do not find any evidence for large systematic differences between our results and estimates of stellar parameters based on other methods. In particular we find no significant systematic offset between stellar masses provided by asteroseismology to our Bayesian estimates based on evolutionary models.Comment: 15 pages, 7 figures, accepted for publication in A&

Precise radial velocities of giant stars VIII. Testing for the presence of planets with CRIRES Infrared Radial Velocities

We have been monitoring 373 very bright (V < 6 mag) G and K giants with high precision optical Doppler spectroscopy for more than a decade at Lick Observatory. Our goal was to discover planetary companions around those stars and to better understand planet formation and evolution around intermediate-mass stars. However, in principle, long-term, g-mode nonradial stellar pulsations or rotating stellar features, such as spots, could effectively mimic a planetary signal in the radial velocity data. Our goal is to compare optical and infrared radial velocities for those stars with periodic radial velocity patterns and to test for consistency of their fitted radial velocity semiamplitudes. Thereby, we distinguish processes intrinsic to the star from orbiting companions as reason for the radial velocity periodicity observed in the optical. Stellar spectra with high spectral resolution have been taken in the H-band with the CRIRES near-infrared spectrograph at ESO's VLT for 20 stars of our Lick survey. Radial velocities are derived using many deep and stable telluric CO2 lines for precise wavelength calibration. We find that the optical and near-infrared radial velocities of the giant stars in our sample are consistent. We present detailed results for eight stars in our sample previously reported to have planets or brown dwarf companions. All eight stars passed the infrared test. We conclude that the planet hypothesis provides the best explanation for the periodic radial velocity patterns observed for these giant stars.Comment: 14 pages, 6 figures, 3 tables, accepted by Astronomy & Astrophysic

Precise radial velocities of giant stars. XI. Two brown dwarfs in 6:1 mean motion resonance around the K giant star ν\nu Ophiuchi

We present radial-velocity (RV) measurements for the K giant $\nu$ Oph (= HIP88048, HD163917, HR6698), which reveal two brown dwarf companions with a period ratio close to 6:1. For our orbital analysis we use 150 precise RV measurements taken at Lick Observatory between 2000 and 2011, and we combine them with RV data for this star available in the literature. Using a stellar mass of $M = 2.7\,M_\odot$ for $\nu$ Oph and applying a self-consistent N-body model we estimate the minimum dynamical companion masses to be $m_1\sin i \approx 22.2\,M_{\mathrm{Jup}}$ and $m_2\sin i \approx 24.7\,M_{\mathrm{Jup}}$, with orbital periods $P_1 \approx 530$ d and $P_2 \approx 3185$ d. We study a large set of potential orbital configurations for this system, employing a bootstrap analysis and a systematic $\chi_{\nu}^2$ grid-search coupled with our dynamical fitting model, and we examine their long-term stability. We find that the system is indeed locked in a 6:1 mean motion resonance (MMR), with $\Delta \omega$ and all six resonance angles $\theta_{1}, \ldots, \theta_{6}$ librating around 0$^\circ$. We also test a large set of coplanar inclined configurations, and we find that the system will remain in a stable resonance for most of these configurations. The $\nu$ Oph system is important for probing planetary formation and evolution scenarios. It seems very likely that the two brown dwarf companions of $\nu$ Oph formed like planets in a circumstellar disk around the star and have been trapped in a MMR by smooth migration capture.Comment: 17 pages, 9 figures. New version with corrected number in title. No other change

Disentangling 2:1 resonant radial velocity orbits from eccentric ones and a case study for HD 27894

In radial velocity observations, a pair of extrasolar planets near a 2:1 orbital resonance can be misinterpreted as a single eccentric planet, if data are sparse and measurement precision insufficient to distinguish between these models. We determine the fraction of alleged single-planet RV detected systems for which a 2:1 resonant pair of planets is also a viable model and address the question of how the models can be disentangled. By simulation we quantified the mismatch arising from applying the wrong model. Model alternatives are illustrated using the supposed single-planet system HD 27894 for which we also study the dynamical stability of near-2:1 resonant solutions. From the data scatter around the fitted single-planet Keplerians, we find that for $74\%$ of the $254$ putative single-planet systems, a 2:1 resonant pair cannot be excluded as a viable model, since the error due to the wrong model is smaller than the scatter. For $187$ stars $\chi ^2$-probabilities can be used to reject the Keplerian models with a confidence of $95\%$ for $54\%$ of the stars and with $99.9\%$ for $39\%$ of the stars. For HD 27894 a considerable fit improvement is obtained when adding a low-mass planet near half the orbital period of the known Jovian planet. Dynamical analysis demonstrates that this system is stable when both planets are initially placed on circular orbits. For fully Keplerian orbits a stable system is only obtained if the eccentricity of the inner planet is constrained to $<0.3$. A large part of the allegedly RV detected single-planet systems should be scrutinized in order to determine the fraction of systems containing near-2:1 resonant pairs of planets. Knowing the abundance of such systems will allow us to revise the eccentricity distribution for extrasolar planets and provide direct constraints for planetary system formation.Comment: 12 pages, 8 figures, one of them composed by two files, accepted by A&A, citations may appear in a non-standard way (double brackets) due to reformatting needs. Abstract slightly adjuste

Dynamical analysis of the circumprimary planet in the eccentric binary system HD59686

We present a detailed orbital and stability analysis of the HD~59686 binary-star planet system. HD~59686 is a single-lined moderately close ($a_{B} = 13.6\,$AU) eccentric ($e_{B} = 0.73$) binary, where the primary is an evolved K giant with mass $M = 1.9 M_{\odot}$ and the secondary is a star with a minimum mass of $m_{B} = 0.53 M_{\odot}$. Additionally, on the basis of precise radial velocity (RV) data a Jovian planet with a minimum mass of $m_p = 7 M_{\mathrm{Jup}}$, orbiting the primary on a nearly circular S-type orbit with $e_p = 0.05$ and $a_p = 1.09\,$AU, has recently been announced. We investigate large sets of orbital fits consistent with HD 59686's radial velocity data by applying bootstrap and systematic grid-search techniques coupled with self-consistent dynamical fitting. We perform long-term dynamical integrations of these fits to constrain the permitted orbital configurations. We find that if the binary and the planet in this system have prograde and aligned coplanar orbits, there are narrow regions of stable orbital solutions locked in a secular apsidal alignment with the angle between the periapses, $\Delta \omega$, librating about $0^\circ$. We also test a large number of mutually inclined dynamical models in an attempt to constrain the three-dimensional orbital architecture. We find that for nearly coplanar and retrograde orbits with mutual inclination $145^\circ \lesssim \Delta i \leq 180^\circ$, the system is fully stable for a large range of orbital solutions.Comment: 17 pages, 11 figures, accepted for publication by A

A comprehensive examination of the Eps Eri system -- Verification of a 4 micron narrow-band high-contrast imaging approach for planet searches

Due to its proximity, youth, and solar-like characteristics with a spectral type of K2V, Eps Eri is one of the most extensively studied systems in an extrasolar planet context. Based on radial velocity, astrometry, and studies of the structure of its circumstellar debris disk, at least two planetary companion candidates to Eps Eri have been inferred in the literature (Eps Eri b, Eps Eri c). Some of these methods also hint at additional companions residing in the system. Here we present a new adaptive optics assisted high-contrast imaging approach that takes advantage of the favourable planet spectral energy distribution at 4 microns, using narrow-band angular differential imaging to provide an improved contrast at small and intermediate separations from the star. We use this method to search for planets at orbits intermediate between Eps Eri b (3.4 AU) and Eps Eri c (40 AU). The method is described in detail, and important issues related to the detectability of planets such as the age of Eps Eri and constraints from indirect measurements are discussed. The non-detection of companion candidates provides stringent upper limits for the masses of additional planets. Using a combination of the existing dynamic and imaging data, we exclude the presence of any planetary companion more massive than 3 Mjup anywhere in the Eps Eri system. Specifically, with regards to the possible residual linear radial velocity trend, we find that it is unlikely to correspond to a real physical companion if the system is as young as 200 Myr, whereas if it is as old as 800 Myr, there is an allowed semi-major axis range between about 8.5 and 25 AU.Comment: 11 pages, 8 figures, A&A accepte

Precise Radial Velocities of Giant Stars VII. Occurrence Rate of Giant Extrasolar Planets as a Function of Mass and Metallicity

(abridged) We have obtained precise radial velocities for a sample of 373 G and K type giants at Lick Observatory regularly over more than 12 years. Planets have been identified around 15 giant stars; an additional 20 giant stars host planet candidates. We investigate the occurrence rate of substellar companions around giant stars as a function of stellar mass and metallicity. We probe the stellar mass range from about 1 to beyond 3 M_Sun, which is not being explored by main-sequence samples. We fit the giant planet occurrence rate as a function of stellar mass and metallicity with a Gaussian and an exponential distribution, respectively. We find strong evidence for a planet-metallicity correlation among the secure planet hosts of our giant star sample, in agreement with the one for main-sequence stars. However, the planet-metallicity correlation is absent for our sample of planet candidates, raising the suspicion that a good fraction of them might indeed not be planets. Consistent with the results obtained by Johnson for subgiants, the giant planet occurrence rate increases in the stellar mass interval from 1 to 1.9 M_Sun. However, there is a maximum at a stellar mass of 1.9 +0.1/-0.5 M_Sun, and the occurrence rate drops rapidly for masses larger than 2.5-3.0 M_Sun. We do not find any planets around stars more massive than 2.7 M_Sun, although there are 113 stars with masses between 2.7 and 5 M_Sun in our sample (corresponding to a giant planet occurrence rate < 1.6% at 68.3% confidence in that stellar mass bin). We also show that this result is not a selection effect related to the planet detectability being a function of the stellar mass. We conclude that giant planet formation or inward migration is suppressed around higher mass stars, possibly because of faster disk depletion coupled with a longer migration timescale.Comment: 13 pages plus long table appendix, accepted by A&

The young binary HD 102077: Orbit, spectral type, kinematics, and moving group membership

The K-type binary star HD 102077 was proposed as a candidate member of the TW Hydrae Association (TWA) which is a young (5-15 Myr) moving group in close proximity (~50 pc) to the solar system. The aim of this work is to verify this hypothesis by different means. We first combine diffraction-limited observations from the ESO NTT 3.5m telescope in SDSS-i' and -z' passbands and ESO 3.6m telescope in H-band with literature data to obtain a new, amended orbit fit, estimate the spectral types of both components, and reanalyse the Hipparcos parallax and proper motion taking the orbital motion into account. Moreover, we use two high-resolution spectra of HD 102077 obtained with the fibre-fed optical echelle spectrograph FEROS at the MPG/ESO 2.2m telescope to determine the radial velocity and the lithium equivalent width of the system. The trajectory of HD 102077 is well constrained and we derive a total system mass of $2.6 \pm 0.8\,$ M$_{\odot}$ and a semi-major axis of $14.9 \pm 1.6\,$AU. From the i'-z' colours we infer an integrated spectral type of K2V, and individual spectral types of K0 +/- 1 and K5 +/- 1. The radial velocity corrected for the orbital motion of the system is $17.6 \pm 2\,$km/s. Even though the parallax determination from the Hipparcos data is not influenced by the orbital motion, the proper motion changes to $\mu_\alpha*\cos(\delta) = -137.84 \pm 1.26\,$ mas/yr and $\mu_\delta = -33.53 \pm 1.45 \,$mas/yr. With the resultant space motion, the probability of HD 102077 being a member of TWA is less than 1%. Furthermore, the lithium equivalent width of $200 \pm 4\,$m\AA $\,$ is consistent with an age between 30 Myr and 120 Myr and thus older than the predicted age of TWA. In conclusion, HD 102077's age, galactic space motion, and position do not fit TWA or any other young moving group

On the Transit Potential of the Planet Orbiting iota Draconis

Most of the known transiting exoplanets are in short-period orbits, largely due to the bias inherent in detecting planets through the transit technique. However, the eccentricity distribution of the known radial velocity planets results in many of those planets having a non-negligible transit probability. One such case is the massive planet orbiting the giant star iota Draconis, a situation where both the orientation of the planet's eccentric orbit and the size of the host star inflate the transit probability to a much higher value than for a typical hot Jupiter. Here we present a revised fit of the radial velocity data with new measurements and a photometric analysis of the stellar variability. We provide a revised transit probability, an improved transit ephemeris, and discuss the prospects for observing a transit of this planet from both the ground and space.Comment: 6 pages, 7 figures, accepted for publication in ApJ. Radial velocities will be made available in the on-line version and through the NASA Star and Exoplanet Database (NStED). Minor corrections from ApJ proof have been applie