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&

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

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

Improving the open cluster census. I. Comparison of clustering algorithms applied to Gaia DR2 data

The census of open clusters in the Milky Way is in a never-before seen state of flux. Recent works have reported hundreds of new open clusters thanks to the incredible astrometric quality of the Gaia satellite, but other works have also reported that many open clusters discovered in the pre Gaia era may be associations. We aim to conduct a comparison of clustering algorithms used to detect open clusters, attempting to statistically quantify their strengths and weaknesses by deriving the sensitivity, specificity, and precision of each as well as their true positive rate against a larger sample. We selected DBSCAN, HDBSCAN, and Gaussian mixture models for further study, owing to their speed and appropriateness for use with Gaia data. We developed a preprocessing pipeline for Gaia data and developed the algorithms further for the specific application to open clusters. We derived detection rates for all 1385 open clusters in the fields in our study as well as more detailed performance statistics for 100 of these open clusters. DBSCAN was sensitive to 50% to 62% of the true positive open clusters in our sample, with generally very good specificity and precision. HDBSCAN traded precision for a higher sensitivity of up to 82%, especially across different distances and scales of open clusters. Gaussian mixture models were slow and only sensitive to 33% of open clusters in our sample, which tended to be larger objects. Additionally, we report on 41 new open cluster candidates detected by HDBSCAN, three of which are closer than 500 pc. When used with additional post-processing to mitigate its false positives, we have found that HDBSCAN is the most sensitive and effective algorithm for recovering open clusters in Gaia data. Our results suggest that many more new and already reported open clusters have yet to be detected in Gaia data.Comment: 28 pages, 13 figures, and 8 tables. Accepted in A&A. Supporting data is available on request until archiving at the CDS is complete

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

Improving the open cluster census. II. An all-sky cluster catalogue with Gaia DR3

Data from the Gaia satellite are revolutionising our understanding of the Milky Way. With every new data release, there is a need to update the census of open clusters. We aim to conduct a blind, all-sky search for open clusters using 729 million sources from Gaia DR3 down to magnitude $G\sim20$, creating a homogeneous catalogue of clusters including many new objects. We used the Hierarchical Density-Based Spatial Clustering of Applications with Noise (HDBSCAN) algorithm to recover clusters. We validated our clusters using a statistical density test and a Bayesian convolutional neural network for colour-magnitude diagram classification. We inferred basic astrometric parameters, ages, extinctions, and distances for the clusters in the catalogue. We recovered 7167 clusters, 2387 of which are candidate new objects and 4782 of which crossmatch to objects in the literature, including 134 globular clusters. A more stringent cut of our catalogue contains 4105 highly reliable clusters, 739 of which are new. Owing to the scope of our methodology, we are able to tentatively suggest that many of the clusters we are unable to detect may not be real, including 1152 clusters from the Milky Way Star Cluster (MWSC) catalogue that should have been detectable in Gaia data. Our cluster membership lists include many new members and often include tidal tails. Our catalogue's distribution traces the galactic warp, the spiral arm structure, and the dust distribution of the Milky Way. While much of the content of our catalogue contains bound open and globular clusters, as many as a few thousand of our clusters are more compatible with unbound moving groups, which we will classify in an upcoming work. We have conducted the largest search for open clusters to date, producing a single homogeneous star cluster catalogue which we make available with this paper.Comment: 31 pages, 17 figures. Accepted in A&A. Updated version that solves crossmatching issues to dwarf galaxies. Before archiving at the CDS completes, supporting data are available at https://drive.google.com/drive/folders/1-9R9132g9FGq6xfcjYI5qMWDFYnMttUG?usp=share_lin

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

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&

Precise radial velocities of giant stars VI. A possible 2:1 resonant planet pair around the K giant star η\eta Cet

We report the discovery of a new planetary system around the K giant $\eta$ Cet (HIP 5364, HD 6805) based on 118 high-precision optical radial velocities taken at Lick Observatory since July 2000. Since October 2011 an additional nine near-infrared Doppler measurements have been taken using the ESO CRIRES spectrograph (VLT, UT1). The visible data set shows two clear periodicities. Although we cannot completely rule out that the shorter period is due to rotational modulation of stellar features, the infrared data show the same variations as in the optical, which strongly supports that the variations are caused by two planets. Assuming the mass of $\eta$ Cet to be 1.7 $M_\odot$, the best edge-on coplanar dynamical fit to the data is consistent with two massive planets ($m_b\sin i$ = 2.6 $\pm$ 0.2 $M_{\mathrm{Jup}}$, $m_c\sin i$ = 3.3 $\pm$ 0.2 $M_{\mathrm{Jup}}$), with periods of $P_b$ = 407 $\pm$ 3 days and $P_c$ = 740 $\pm$ 5 days and eccentricities of $e_b$ = 0.12 $\pm$ 0.05 and $e_c$ = 0.08 $\pm$ 0.03. We tested a wide variety of edge-on coplanar and inclined planetary configurations for stability, which agree with the derived radial velocities. We find that in certain coplanar orbital configurations with moderate $e_b$ eccentricity, the planets can be effectively trapped in an anti-aligned 2:1 mean motion resonance. A much larger non-resonant stable region exists in low-eccentricity parameter space, although it appears to be much farther from the best fit than the 2:1 resonant region. In all other cases, the system is categorized as unstable or chaotic. Another conclusion from the coplanar inclined dynamical test is that the planets can be at most a factor of $\sim$ 1.4 more massive than their suggested minimum masses. This stability constraint on the inclination excludes the possibility of two brown dwarfs, and strongly favors a planetary system.Comment: 15 pages, 11 figures, accepted for publication in A&A on June 20, 201