New stellar encounters discovered in the second Gaia data release
AAstronomy & Astrophysics manuscript no. stellar_encounters_gdr2 c (cid:13)
ESO 201826 June 2018
New stellar encounters discovered in the second Gaia data release
C.A.L. Bailer-Jones, J. Rybizki, R. Andrae, M. Fouesneau
Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, GermanySubmitted 19 May 2018. Revised 17 June 2018. Accepted 19 June 2018.
ABSTRACT
Passing stars may play an important role in the evolution of our solar system. We search for close stellar encounters to the Sun amongall 7.2 million stars in Gaia DR2 that have six-dimensional phase space data. We characterize encounters by integrating their orbitsthrough a Galactic potential and propagating the correlated uncertainties via a Monte Carlo resampling. After filtering to removespurious data, we find 694 stars that have median (over uncertainties) closest encounter distances within 5 pc, all occurring within15 Myr from now. 26 of these have at least a 50% chance of coming closer than 1 pc (and 7 within 0.5 pc), all but one of which arenewly discovered here. We confirm some and refute several other previously-identified encounters, confirming suspicions about theirdata. The closest encounter in the sample is Gl 710, which has a 95% probability of coming closer than 0.08 pc (17 000 AU). Takingmass estimates obtained from Gaia astrometry and multiband photometry for essentially all encounters, we find that Gl 710 also hasthe largest impulse on the Oort cloud. Using a Galaxy model, we compute the completeness of the Gaia DR2 encountering sample asa function of perihelion time and distance. Only 15% of encounters within 5 pc occurring within ± / or cool stars. Accounting for the incompleteness, we infer the present rate ofencounters within 1 pc to be 19 . ± . ff ect is hard to quantify. Key words.
Oort cloud – methods: analytical, statistical – solar neighbourhood – stars: kinematics and dynamics – surveys: Gaia
1. Introduction
The first convincing evidence for relative stellar motion camefrom Edmund Halley in 1718. Since then – if not well be-fore – people have wondered how close other stars may cometo our own. Several studies over the past quarter century haveused proper motions and radial velocities to answer this ques-tion (Matthews 1994; Mülläri & Orlov 1996; García-Sánchezet al. 1999, 2001; Dybczy´nski 2006; Bobylev 2010a,b; Jiménez-Torres et al. 2011; Bailer-Jones 2015a; Dybczy´nski & Berski2015; Mamajek et al. 2015; Berski & Dybczy´nski 2016; Bobylev& Bajkova 2017; Bobylev 2017; Bailer-Jones 2018).Other than being interesting in their own right, close encoun-ters may have played a significant role in the evolution of oursolar system, in particular of the Oort cloud. This may also havehad implications for the development of life, since a strong per-turbation of the Oort cloud by an encountering star could pushcomets into the inner solar system. An ensuing collision with theEarth could be catastrophic enough to cause a mass extinction.Such a fate probably befell the dinosaurs 65 Myr ago. Studies ofstellar encounters and comet impacts on the Earth can also beused to learn about the general hazards for life on exoplanets.The recent publication of the second Gaia data release(Gaia DR2; Gaia Collaboration 2018) is a boon to encounterstudies. It contains six-dimensional (6D) kinematic data – po-sition, parallax, proper motion, and radial velocity – for 7.2 mil-lion stars. This is 22 times larger than our previous encounterstudy using TGAS in the first Gaia data release. Furthermore,TGAS contained no radial velocities, so these had to be obtainedby cross matching to external catalogues.Here we report on results of looking for close encountersin Gaia DR2. Our approach follows very closely that taken by Bailer-Jones (2015a) (paper 1) for Hipparcos and Bailer-Jones(2018) (paper 2) for TGAS. The main di ff erences here are: – We use only Gaia data. This avoids the complications ofhaving to obtain radial velocities from several di ff erent cata-logues, which resulted in some data heterogeneity and dupli-cate sources. – Using mass estimates from multiband photometry andGaia DR2 astrometry in Fouesneau et al. (in preparation), weestimate the momentum transfer from an encountering star toOort cloud comets for essentially every encounter. – We account for the incompleteness of our encounter sam-ple by sampling from a self-consistent spatial and kinematicGalaxy model. This is more realistic that the analytic modeldeveloped in paper 2.The main limitation in the number of encounter candidates inthis study is still the availability of radial velocities. WhileGaia DR2 contains five-parameter astrometry for 1.33 billionstars, only 7.2 million have published radial velocities. These arepredominantly brighter than G =
14 mag, and are also limited tothe approximate T e ff range of 3550 < T e ff / K < T e ff estimates in Gaia DR2 are described in Andrae et al.2018).In section 2 we describe how we select our sample and inferthe distribution of the encounter parameters for each candidate.We analyse the results in section 3, discussing new (and dubious)cases, and highlighting disagreement with earlier results. Issuesof spurious data and imperfect filtering we discuss in subsec-tions 2.2 and 3.2. In section 4 we introduce the new complete-ness model and use this to derive the completeness-corrected en-counter rate. We conclude in section 5 with a brief discussion. Article number, page 1 of 12 a r X i v : . [ a s t r o - ph . S R ] J un & A proofs: manuscript no. stellar_encounters_gdr2
2. Identification and characterization of closeencounters
We searched the Gaia archive for all stars which would approachwithin 10 pc of the Sun under the assumption that they move onunaccelerated paths relative to the solar system (the so-called“linear motion approximation”, LMA, defined in paper 1). Thearchive ADQL query for this is in appendix A. This yielded 3865encounter candidates, which we call the unfiltered sample . Inboth this selection and the orbital integration used later, we esti-mate distance using inverse parallax rather than doing a properinference (Bailer-Jones 2015b; Bailer-Jones et al. 2018). This isacceptable here, because for the filtered sample we define later,95% have fractional parallax uncertainties below 0.08 (99% be-low 0.14; the largest is 0.35), meaning the simple inversion is areasonably good approximation.The Gaia data are used as they are, other than we added0.029 mas to the parallaxes to accommodate the global paral-lax zeropoint Lindegren et al. (2018). (Neglecting this would re-move 83 sources from the unfiltered sample.) There is no evi-dence for a proper motion zeropoint o ff set (Arenou et al. 2018).While there may be a systematic di ff erence between the Gaia ra-dial velocities and other catalogues of up to 0.5 km s − , the ori-gin of this is unclear (Katz et al. 2018). This is of the order ofthe gravitational redshifts which are also not corrected for. Noneof the Gaia uncertainties have been adjusted to accommodatea possible over- or underestimate of their values (e.g. missingRMS systematics).To get more precise encounter parameters for the set of can-didates, we integrated the orbits of the unfiltered sample througha smooth Galactic potential forward and backwards in time.(It was shown in paper 1 that the deviation of an orbit dueto perturbations by individual stars can be neglected.) We usethe same procedure and model as described in paper 2. TheGalactic potential, described in detail in paper 1, is a three-component axisymmetric model. The bar and spiral arms are notincluded, partly because their properties are not well determined,but mostly because using the same potential as in our previousstudies eases comparison of results. As the orbit segments up toencounter are generally short compared to the scale lengths inthe model, the exact choice of potential will have only a smallimpact on the encounter parameters for most stars.In order to accommodate and propagate the uncertainties inthe data, we draw 2000 samples from the 6D covariant proba-bility density function (PDF) over the data – position, parallax,proper motion, and radial velocity – for each star and integratethe orbits of each of these “surrogates” through the potential.The distribution of the perihelion time, distance, and speed overthe surrogates for each star is used to characterize the encounter(in the next section). Comparisons of encounter parameters com-puted with the LMA, orbital integration of the nominal data, andorbital integration of the surrogates were shown in papers 1 and2. Due to the nonlinear transformation from astrometric mea-surements to perihelion parameters, neglecting the full PDF canlead to erroneous results (and not just erroneous uncertainties),in particular for stars with long travel times. There are many reasons why astrometric solutions in Gaia DR2may be “wrong” for some stars, in the sense that the reporteduncertainties may not be representative of the true uncertainties. The main reasons are: neglect of accelerated motions (i.e. unseencompanions); cross-matching errors leading to the inclusion ofobservations of other sources (spurious data); a poor correctionof the so-called “DOF bug” (see appendix A of Lindegren et al.2018). These can lead to erroneous estimates of the quantities ortheir uncertainties.Various metrics on the astrometric solution are reported inGaia DR2 to help identify good solutions. Lindegren et al. (2018)discuss some of these and give an example of a set of cutswhich may be used to define a conservative sample, i.e. an agres-sive removal of poor solutions (e.g. in their Figure C.2). Thisis not appropriate for our work, however, because we are look-ing to determine the encounter rate, not just find the most re-liable encounters. While quantiles on the various metrics areeasily measured, there is no good model for the expected dis-tribution of these for only non-spurious results. Concepts like“the reduced χ should be about one” are simplistic at best (andstatistically questionable), and also don’t tell us whether devi-ations from this are “wrong” or just have mildly deviant as-trometry or slightly underestimated uncertainties. Even a highlysignificant astrometric excess noise of a few mas may be oflittle consequence if the parallax and proper motion are large.It is therefore di ffi cult to make a reliable cut. The main qual-ity metric of interest here is the “unit weight error”, definedin appendix A of Lindegren et al. (2018) as u = ( χ /ν ) / ,where χ is the metric astrometric_chi2_al and ν is the de-grees of freedom, equal to astrometric_n_good_obs_al − u correlates quite strongly with the astrometric_excess_noise and reasonably well, but lesstightly, with astrometric_excess_noise_sig . There is nocorrelation between u and visibility_periods_used . Fig-ure C.2 of Lindegren et al. (2018) plots u against G . We see u increasing for brighter sources, albeit it with a lot of scat-ter. The full range of u in our unfiltered sample is 0.63 to122, and the brighest star in our unfiltered sample (Gaia DR25698015743040715264 = rho Puppis) has G = .
67 mag.This corresponds to u =
35 for the selection in Figure C.2,so we only retain sources with u <
35. This moderatelyliberal cut reduces the sample size to 3465. Rho Puppis it-self has u = astrometric_excess_noise = astrometric_excess_noise_sig = visibility_periods_used is at least 8. Accord-ing to Arenou et al. (2018), this should help to remove the mostspurious proper motions. This cuts the sample down to 3379.Large uncertainties in the data, provided they are represen-tative of the true uncertainty, are not a problem per se becausethey are accommodated by the resampling we used to map thePDF of the perihelion parameters. We will also accommodatethis PDF when we compute the completeness and encounter ratein section 4. We therefore do not filter out sources due to largeuncertainties.One should be very careful about filtering out extreme val-ues of the data just because they are extreme. Gaia DR2 is knownto include sources with implausibly large parallaxes. For exam-ple, there are 21 sources with (cid:36) > G >
19 mag none of them appear in any of our se-lections. Other, less extreme, spurious values are impossible toidentify without using additional information. Spuriously largeproper motions are less of a problem, because close encounter
Article number, page 2 of 12ailer-Jones et al.: Close encounters to the Sun in Gaia-DR2 . . . . . . perihelion time (t ph ) / Myr pe r i he li on d i s t an c e ( d ph ) / p c Fig. 1.
The distribution of the perihelion parameters for the 2000surrogates used to characterize the encounter of Gl 710 = Gaia DR24270814637616488064. stars generally have relatively small proper motions (unless theyare currently very close to encounter). In contrast, stars whichcome close tend to be those that currently have large radial ve-locities (compared to their transverse velocities). The radial ve-locity processing for Gaia DR2 may have produced spuriouslylarge radial velocities (the pipeline can produce values up to ± − ).Katz et al. (2018) report they were hard to ver-ify due to the absence of observations of standards with radialvelocities larger than 550 kms. They further state that the preci-sions are lower for stars with | v r | >
175 km s − , but not that theradial velocities themselves are problematic. Moreover, they vi-sually inspected all results with | v r | >
500 km s − and removedsuspicious results. At this point in the filtering we have onlyseven stars with | v r | >
500 km s − , six of which were deter-mined by just two focal plane transits ( rv_nb_transits = | σ (v r ) | <
20 km s − . It is tempting to filter out stars with a smallnumber of transits, but this also removes valid measurements,such as that for the closest known encounter, Gl 710, which hastwo transits and a radial velocity consistent with non-Gaia mea-surements. We do not filter on rv_nb_transits .We define the remaning set of 3379 stars as the filteredsample . It comprises all stars that approach within 10 pc ofthe Sun according to the LMA (but not necessarily the or-bital integration) and which satisfy the filters u <
35 and visibility_periods_used ≥
8. When it comes to consider-ing the completeness-corrected encounter rate (section 4.3), weshall further limit this sample to stars with G < .
3. Close encounters found in Gaia DR2
We first look at the overall results, then discuss individual starsfound approaching within 1 pc, and finally those close encoun-ters from papers 1 and 2 not found in the present study.
The perihelion for each star is described by the distribution ofthe 2000 surrogates in the perihelion time, t ph , distance, d ph , and Table 1.
The number of stars in the filtered sample found by the orbitintegration to have d medph < d maxph (for any t medph ). Stars with potentiallyproblematic data have not been excluded. d maxph / pc No. stars ∞ . .
25 3 −20 −10 0 10 20 perihelion time (t phmed ) / Myr pe r i he li on d i s t an c e ( d ph m ed ) / p c llll ll ll ll ll lll ll l lllll l lll ll lll lll lll ll lllll lll llll ll l lll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll lllll lll ll ll l ll ll l lll l lll l ll l l llll llll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll llll lll ll ll lll l ll l lllll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll lll llll ll lllll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lllll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lllll llll l ll ll llllll l ll lll l ll lll lll l lll ll l ll ll lll l llll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll llll ll llll ll lll lll l ll llllll lllll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llllll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l lll llll ll ll l lll lllll ll ll ll lll ll lll l ll lll l llll ll l ll ll ll ll ll ll llll l l l lll lll llll l lll lll l ll ll l l ll llllll ll ll ll llllll ll llll llll llll l llll l l llllll l l lll l ll lll l ll ll llllll lll ll lll ll llll l ll l lll l lll ll ll llll llll ll lllll lllll lll lll ll l lllll llll lll l ll llll lll ll lll l l lllll ll ll ll ll ll l ll l lll lll l lll lll lll ll ll lll ll lll lll ll l lll l lll l lll ll ll lll ll lll lll lllll lll l ll l ll lll l ll ll llll lll l lll ll ll l lll l ll lll ll l ll lll ll llll ll l llll llll lll lll ll l ll ll l ll llll l lll lll ll l ll ll lll llll llll l lll l llll l ll ll ll l lll l ll lll l llll llll ll ll l lll ll lll ll ll ll lll ll llll l ll ll llll l lll llll ll l lll l l ll ll ll ll l lll lll l lll lll ll llll l ll lll l l ll ll l ll ll l lll l ll l ll lllll ll lll llll l ll ll ll ll ll ll l lll l lll ll l llll lll l ll lll ll lllll l llll l l lll ll ll l ll lll lll l ll lll ll lll ll ll lll l ll llll lll l lll llll l ll ll l ll ll ll ll l l lll l l ll lll l l lll llllll ll ll l llll lllll lll lll lll ll lll lll lll lll lll ll ll lll l ll ll llll llll ll ll lll ll llll lll l l llll ll ll lllll l lll lllll lll ll llll l ll ll ll lll ll lll l llll ll lll l l lll lllll ll lll lll ll lll ll lll l lll ll llll ll lllll llllll l ll l l ll llllll lllll ll llll l l ll l ll lll l lll ll llll ll llll ll lll llll llllll l ll l l ll ll llllll llll llllll ll l l lll ll l lll lll ll l ll l ll l l l llllllll lll llll ll ll l ll l ll l ll ll llll l l lll ll llllll ll lll l l lll ll llll lll l ll l ll lll lll ll lllll ll ll llll ll ll ll ll l llll ll ll lll ll ll llll ll llll ll lll lll ll lll llll llll lll ll lll ll llll l ll l lll ll l llllll lll ll ll ll lll ll lllll ll ll lll ll l ll ll l lll l llll llllll l lll ll l llll l ll l ll l l lll l lll ll l ll ll llllll lll llll llll l ll ll l ll lll llll ll l ll llll l l ll lll lll ll ll lll ll lll l ll ll ll ll lll lll l lllll ll ll lll ll ll ll lll l llll l ll ll ll l ll ll llll lll ll l ll lll l lll llllll ll l lll l lll llll lll ll ll lll ll lll lll l ll l lll llll lll l llll l ll lll ll lll l ll lll ll ll ll ll l lll lllllll ll ll l ll ll lll lll ll ll llll ll llll lll l lll ll lllll ll ll l llll ll ll ll ll ll lll llll l ll lll ll llll lllllll ll lll llll ll ll ll llll ll ll l lll l l l llll lll ll lll ll ll ll ll ll lll l lllll lll l ll ll lll l ll ll lll lll ll l llll ll lll ll llll llll l lllll lll ll ll l ll lll l lll ll l llll ll l ll lll llll l ll lll ll lllll llll ll llll l ll ll l l lll ll lll l lll ll ll lll l l llll ll lll l lllll lll ll l l lll ll ll lll l llll ll lll ll ll lll lll l l l lllll ll lll ll l ll llll ll l ll lll l lllll ll lll l ll ll ll llll l ll l ll llll l lll lll ll ll l l lll l llllllll llll lll ll ll lll ll llll l ll ll ll ll ll ll ll ll ll ll ll llll lll ll llll l lll lll ll llll l lll lll ll lll lll llllll lll lll ll lll ll ll llllll ll lllll lll ll ll l ll ll lll ll l ll ll l l ll l ll llll lll llll ll l ll lll ll lll ll lllll l ll ll lll lll l ll lll ll lll l lll ll lll ll lll lll ll l lll ll ll l llll ll llll lll llll l ll lll lll llll lll ll l ll l ll llll l ll lll lll l lll ll l ll l lll ll ll llll ll l lll ll ll l lll ll l ll l l ll ll llll lll l ll lll llll llll llllll l lll l llll ll l ll lll l lll l ll lll lllll ll lll l l lll l lll l ll l l ll ll l l l llll l l lll l ll lllll ll l ll l ll l llll l ll l l ll lll ll lll l l ll lll ll l ll ll llll ll ll ll ll lll ll l lllll l lll ll lll lll lll ll lllll lll llll ll l lll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll lllll lll ll ll l ll ll l lll l lll l ll l l llll llll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll llll lll ll ll lll l ll l lllll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll lll llll ll lllll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lllll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lllll llll l ll ll llllll l ll lll l ll lll lll l lll ll l ll ll lll l llll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll llll ll llll ll lll lll l ll llllll lllll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llllll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l lll llll ll ll l lll lllll ll ll ll lll ll lll l ll lll l llll ll l ll ll ll ll ll ll llll l l l lll lll llll l lll lll l ll ll l l ll llllll ll ll ll llllll ll llll llll llll l llll l l llllll l l lll l ll lll l ll ll llllll lll ll lll ll llll l ll l lll l lll ll ll llll llll ll lllll lllll lll lll ll l lllll llll lll l ll llll lll ll lll l l lllll ll ll ll ll ll l ll l lll lll l lll lll lll ll ll lll ll lll lll ll l lll l lll l lll ll ll lll ll lll lll lllll lll l ll l ll lll l ll ll llll lll l lll ll ll l lll l ll lll ll l ll lll ll llll ll l llll llll lll lll ll l ll ll l ll llll l lll lll ll l ll ll lll llll llll l lll l llll l ll ll ll l lll l ll lll l llll llll ll ll l lll ll lll ll ll ll lll ll llll l ll ll llll l lll llll ll l lll l l ll ll ll ll l lll lll l lll lll ll llll l ll lll l l ll ll l ll ll l lll l ll l ll lllll ll lll llll l ll ll ll ll ll ll l lll l lll ll l llll lll l ll lll ll lllll l llll l l lll ll ll l ll lll lll l ll lll ll lll ll ll lll l ll llll lll l lll llll l ll ll l ll ll ll ll l l lll l l ll lll l l lll llllll ll ll l llll lllll lll lll lll ll lll lll lll lll lll ll ll lll l ll ll llll llll ll ll lll ll llll lll l l llll ll ll lllll l lll lllll lll ll llll l ll ll ll lll ll lll l llll ll lll l l lll lllll ll lll lll ll lll ll lll l lll ll llll ll lllll llllll l ll l l ll llllll lllll ll llll l l ll l ll lll l lll ll llll ll llll ll lll llll llllll l ll l l ll ll llllll llll llllll ll l l lll ll l lll lll ll l ll l ll l l l llllllll lll llll ll ll l ll l ll l ll ll llll l l lll ll llllll ll lll l l lll ll llll lll l ll l ll lll lll ll lllll ll ll llll ll ll ll ll l llll ll ll lll ll ll llll ll llll ll lll lll ll lll llll llll lll ll lll ll llll l ll l lll ll l llllll lll ll ll ll lll ll lllll ll ll lll ll l ll ll l lll l llll llllll l lll ll l llll l ll l ll l l lll l lll ll l ll ll llllll lll llll llll l ll ll l ll lll llll ll l ll llll l l ll lll lll ll ll lll ll lll l ll ll ll ll lll lll l lllll ll ll lll ll ll ll lll l llll l ll ll ll l ll ll llll lll ll l ll lll l lll llllll ll l lll l lll llll lll ll ll lll ll lll lll l ll l lll llll lll l llll l ll lll ll lll l ll lll ll ll ll ll l lll lllllll ll ll l ll ll lll lll ll ll llll ll llll lll l lll ll lllll ll ll l llll ll ll ll ll ll lll llll l ll lll ll llll lllllll ll lll llll ll ll ll llll ll ll l lll l l l llll lll ll lll ll ll ll ll ll lll l lllll lll l ll ll lll l ll ll lll lll ll l llll ll lll ll llll llll l lllll lll ll ll l ll lll l lll ll l llll ll l ll lll llll l ll lll ll lllll llll ll llll l ll ll l l lll ll lll l lll ll ll lll l l llll ll lll l lllll lll ll l l lll ll ll lll l llll ll lll ll ll lll lll l l l lllll ll lll ll l ll llll ll l ll lll l lllll ll lll l ll ll ll llll l ll l ll llll l lll lll ll ll l l lll l llllllll llll lll ll ll lll ll llll l ll ll ll ll ll ll ll ll ll ll ll llll lll ll llll l lll lll ll llll l lll lll ll lll lll llllll lll lll ll lll ll ll llllll ll lllll lll ll ll l ll ll lll ll l ll ll l l ll l ll llll lll llll ll l ll lll ll lll ll lllll l ll ll lll lll l ll lll ll lll l lll ll lll ll lll lll ll l lll ll ll l llll ll llll lll llll l ll lll lll llll lll ll l ll l ll llll l ll lll lll l lll ll l ll l lll ll ll llll ll l lll ll ll l lll ll l ll l l ll ll llll lll l ll lll llll llll llllll l lll l llll ll l ll lll l lll l ll lll lllll ll lll l l lll l lll l ll l l ll ll l l l llll l l lll l ll lllll ll l ll l ll l llll l ll l l ll lll ll lll l l ll lll ll l ll ll llll ll ll ll ll lll ll l lllll l lll ll lll lll lll ll lllll lll llll ll l lll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll lllll lll ll ll l ll ll l lll l lll l ll l l llll llll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll llll lll ll ll lll l ll l lllll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll lll llll ll lllll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lllll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lllll llll l ll ll llllll l ll lll l ll lll lll l lll ll l ll ll lll l llll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll llll ll llll ll lll lll l ll llllll lllll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llllll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l lll llll ll ll l lll lllll ll ll ll lll ll lll l ll lll l llll ll l ll ll ll ll ll ll llll l l l lll lll llll l lll lll l ll ll l l ll llllll ll ll ll llllll ll llll llll llll l llll l l llllll l l lll l ll lll l ll ll llllll lll ll lll ll llll l ll l lll l lll ll ll llll llll ll lllll lllll lll lll ll l lllll llll lll l ll llll lll ll lll l l lllll ll ll ll ll ll l ll l lll lll l lll lll lll ll ll lll ll lll lll ll l lll l lll l lll ll ll lll ll lll lll lllll lll l ll l ll lll l ll ll llll lll l lll ll ll l lll l ll lll ll l ll lll ll llll ll l llll llll lll lll ll l ll ll l ll llll l lll lll ll l ll ll lll llll llll l lll l llll l ll ll ll l lll l ll lll l llll llll ll ll l lll ll lll ll ll ll lll ll llll l ll ll llll l lll llll ll l lll l l ll ll ll ll l lll lll l lll lll ll llll l ll lll l l ll ll l ll ll l lll l ll l ll lllll ll lll llll l ll ll ll ll ll ll l lll l lll ll l llll lll l ll lll ll lllll l llll l l lll ll ll l ll lll lll l ll lll ll lll ll ll lll l ll llll lll l lll llll l ll ll l ll ll ll ll l l lll l l ll lll l l lll llllll ll ll l llll lllll lll lll lll ll lll lll lll lll lll ll ll lll l ll ll llll llll ll ll lll ll llll lll l l llll ll ll lllll l lll lllll lll ll llll l ll ll ll lll ll lll l llll ll lll l l lll lllll ll lll lll ll lll ll lll l lll ll llll ll lllll llllll l ll l l ll llllll lllll ll llll l l ll l ll lll l lll ll llll ll llll ll lll llll llllll l ll l l ll ll llllll llll llllll ll l l lll ll l lll lll ll l ll l ll l l l llllllll lll llll ll ll l ll l ll l ll ll llll l l lll ll llllll ll lll l l lll ll llll lll l ll l ll lll lll ll lllll ll ll llll ll ll ll ll l llll ll ll lll ll ll llll ll llll ll lll lll ll lll llll llll lll ll lll ll llll l ll l lll ll l llllll lll ll ll ll lll ll lllll ll ll lll ll l ll ll l lll l llll llllll l lll ll l llll l ll l ll l l lll l lll ll l ll ll llllll lll llll llll l ll ll l ll lll llll ll l ll llll l l ll lll lll ll ll lll ll lll l ll ll ll ll lll lll l lllll ll ll lll ll ll ll lll l llll l ll ll ll l ll ll llll lll ll l ll lll l lll llllll ll l lll l lll llll lll ll ll lll ll lll lll l ll l lll llll lll l llll l ll lll ll lll l ll lll ll ll ll ll l lll lllllll ll ll l ll ll lll lll ll ll llll ll llll lll l lll ll lllll ll ll l llll ll ll ll ll ll lll llll l ll lll ll llll lllllll ll lll llll ll ll ll llll ll ll l lll l l l llll lll ll lll ll ll ll ll ll lll l lllll lll l ll ll lll l ll ll lll lll ll l llll ll lll ll llll llll l lllll lll ll ll l ll lll l lll ll l llll ll l ll lll llll l ll lll ll lllll llll ll llll l ll ll l l lll ll lll l lll ll ll lll l l llll ll lll l lllll lll ll l l lll ll ll lll l llll ll lll ll ll lll lll l l l lllll ll lll ll l ll llll ll l ll lll l lllll ll lll l ll ll ll llll l ll l ll llll l lll lll ll ll l l lll l llllllll llll lll ll ll lll ll llll l ll ll ll ll ll ll ll ll ll ll ll llll lll ll llll l lll lll ll llll l lll lll ll lll lll llllll lll lll ll lll ll ll llllll ll lllll lll ll ll l ll ll lll ll l ll ll l l ll l ll llll lll llll ll l ll lll ll lll ll lllll l ll ll lll lll l ll lll ll lll l lll ll lll ll lll lll ll l lll ll ll l llll ll llll lll llll l ll lll lll llll lll ll l ll l ll llll l ll lll lll l lll ll l ll l lll ll ll llll ll l lll ll ll l lll ll l ll l l ll ll llll lll l ll lll llll llll llllll l lll l llll ll l ll lll l lll l ll lll lllll ll lll l l lll l lll l ll l l ll ll l l l llll l l lll l ll lllll ll l ll l ll l llll l ll l l ll lll ll lll l l ll lll ll l ll ll Fig. 2.
Perihelion times and distances computed by orbit integration forthe filtered sample (278 points lie outside the plotting range). Open cir-cles show the median of the perihelion time and distance distributions.The error bars show the limits of the 5% and 95% percentiles. speed v ph . The distribution for one particular star is shown inFigure 1. As in papers 1 and 2 we summarize these using themedian, and characterize their uncertainty using the 5th and 95thpercentiles (which together form a 90% confidence interval, CI).The number of stars coming within various perihelion distancesis shown in Table 1 (see also Figure 10 later). Summary perihel-ion data for those stars with d medph < d medph <
10 pc had they been subject to the orbital inte-gration, but were never selected because they had d ph >
10 pcfrom the LMA. This latter omission will in principle lead to anunderestimate of the derived encounter rate. However, as shownin both papers 1 and 2, this mostly a ff ects stars that will en-counter further in the past / future and / or nearer to the edge ofthe 10 pc distance limit. As stars at such large distances canhardly be considered encountering, we will only be interested in Article number, page 3 of 12 & A proofs: manuscript no. stellar_encounters_gdr2
Table 2.
Perihelion parameters for all stars with a median perihelion distance (median of the samples; d medph ) below 1 pc, sorted by this value. Thefirst column is the Gaia DR2 source ID. Columns 2, 5, and 8 are t medph , d medph , and v medph respectively. The columns labelled 5% and 95% are thebounds of corresponding confidence intervals. Columns 11–16 list the parallax ( (cid:36) , plus the 0.029 mas zeropoint o ff set), total proper motion ( µ ),and radial velocity (v r ) along with their 1-sigma uncertainties. Column 17 is the estimated mass of the star from Fouesneau et al. (in preparation);NA indicates it is missing. The formal 1-sigma uncertainties in the masses are a few percent (systematics are likely to be higher). Those encounterresults we consider bogus are marked with the dagger symbol in the final column. These are discussed in the text. Other may also be dubious:quality metrics from Gaia DR2 can be found in Table 3. The online table at CDS includes all 3379 stars in the filtered sample and reports somecolumns to a higher numerical precision. t ph / kyr d ph / pc v ph / km s − (cid:36) σ ( (cid:36) ) µ σ ( µ ) v r σ (v r ) M med 5% 95% med 5% 95% med 5% 95% mas mas yr − km s − M (cid:12) † † † † † −15 −10 −5 0 5 10 15 perihelion time (t phmed ) / Myr pe r i he li on d i s t an c e ( d ph m ed ) / p c llll ll ll ll ll lll ll l lllll l lll ll ll l lll lll ll lllll lll llll ll l l ll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll ll lll lll ll ll l ll ll l lll l lll l ll l l ll ll l lll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll ll ll lll ll ll lll l ll l ll lll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll ll l llll ll ll lll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lll ll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lll ll llll l ll ll lll l ll l ll lll l ll lll lll l lll ll l ll ll lll l ll ll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll ll ll ll llll ll lll lll l ll llllll l l lll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llll ll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l ll llll ll ll ll ll lll ll l lllll l lll ll ll l lll lll ll lllll lll llll ll l l ll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll ll lll lll ll ll l ll ll l lll l lll l ll l l ll ll l lll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll ll ll lll ll ll lll l ll l ll lll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll ll l llll ll ll lll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lll ll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lll ll llll l ll ll lll l ll l ll lll l ll lll lll l lll ll l ll ll lll l ll ll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll ll ll ll llll ll lll lll l ll llllll l l lll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llll ll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l ll llll ll ll ll ll lll ll l lllll l lll ll ll l lll lll ll lllll lll llll ll l l ll lll l llll lll l l ll ll ll ll lllll l llll lll lll l ll ll ll ll lll lll ll ll l ll ll l lll l lll l ll l l ll ll l lll l lll l l ll ll lll ll l l lll l lll lll ll ll l llll lll l llll ll ll lll ll ll lll l ll l ll lll ll lll l lll ll ll llll ll lll llll l llll l lll l ll ll ll l ll lll ll l llll ll ll lll ll l llll lll ll ll llll ll ll ll l l l lll lllll ll ll ll lll ll ll lll ll llll ll ll ll l lll l ll llll lll l lll llll l ll l lll ll l ll ll lllll ll l l lll llll l lll ll llll ll lll ll l l ll lll ll llll l ll ll lll l ll l ll lll l ll lll lll l lll ll l ll ll lll l ll ll l ll l l ll l llll ll ll ll l l ll l llll ll llll ll l ll ll l ll ll lll ll l lll ll ll ll ll ll llll ll lll lll l ll llllll l l lll l llll ll ll ll ll lll l lll ll ll l lll lll ll l ll lll ll lll lll lllll l lll ll l ll lll ll ll l ll l ll llll ll llll ll lll lll ll ll ll lll ll ll lll lll l ll lllll ll lll ll ll ll l ll ll lll l lll l ll Fig. 3.
As Figure 2, but just showing the 694 stars with d medph < encounters within 5 pc from now on. Moreover, when we later compute the (completeness-corrected) encounter rate, we willlimit the sample to a narrower time window.The encounters with d medph < | t ph | . This is primarily a consequence of the magnitude limit inthe sample (95% brighter than G = / future generally correspond tostars that are currently more distant, and so more likely to be be-low the limiting magnitude. The e ff ective time limit of this studyis 5–10 Myr. Comparing this with Figure 3 of paper 2, we seethat while Gaia DR2 has found many more encounters than ourTGAS study (by a factor of about seven within 5 pc), Gaia DR2does not allow us to probe much further into the past / future, be-cause of the similar magnitude limits on the samples.We also see in Figure 3 a slight reduction in the density ofencounters very close to the present time. (The e ff ect is not quiteas strong as first appears, as it is partly an illusion produced bythe shorter error bars.) This was also seen with in paper 2, wherewe argued this was due to two things: missing bright stars in theGaia catalogue; and the ever smaller volume available for en-counters to occur at arbitrarily near times. Both of these applyfor Gaia DR2, although we argue later that this is also a conse-quence of the limited T e ff range for stars with radial velocities inGaia DR2. Article number, page 4 of 12ailer-Jones et al.: Close encounters to the Sun in Gaia-DR2
Table 3.
Additional data from Gaia DR2 for the close encounters listed in Table 2. They are all taken directly from the catalogue, except for theastrometric “unit weight error” u which is calculated as sqrt[astrometric_chi2_al/(astrometric_n_good_obs_al-5)] . The online tableat CDS includes all 3379 stars in the filtered sample. Gaia DR2 source ID G BP-RP u No. astrometric astrometric No. l b mag mag visibility excess excess RVS deg degperiods noise noise sig transits4270814637616488064 9.06 1.70 1.22 10 0.00 0.00 2 27 6955098506408767360 12.41 0.76 16.19 10 2.49 1661.28 5 176 105571232118090082816 11.79 1.50 1.25 15 0.00 0.00 10 249 -252946037094755244800 12.34 1.49 27.41 10 3.12 4225.87 22 228 -74071528700531704704 12.44 0.78 23.52 10 3.92 3138.25 3 7 -11510911618569239040 8.88 0.77 1.57 17 0.00 0.00 6 126 0154460050601558656 15.37 NA 6.47 10 1.85 359.52 3 174 -116608946489396474752 12.28 1.44 0.89 10 0.00 0.00 6 23 -613376241909848155520 12.52 0.78 28.53 11 5.46 8824.25 8 190 51791617849154434688 11.00 1.09 1.49 11 0.00 0.00 9 70 -184265426029901799552 12.20 2.07 18.88 10 3.44 2533.02 3 32 05261593808165974784 12.69 2.02 1.39 18 0.00 0.00 14 285 -275896469620419457536 13.55 1.98 1.06 11 0.03 0.23 5 313 64252068750338781824 12.10 0.92 20.42 9 2.53 1453.10 5 26 -31949388868571283200 13.12 NA 1.34 11 0.13 5.81 2 86 -131802650932953918976 12.69 1.02 0.99 14 0.00 0.00 11 52 -113105694081553243008 12.31 1.24 25.05 9 4.78 6025.61 6 215 -25231593594752514304 12.03 1.87 1.22 16 0.00 0.00 2 293 -94472507190884080000 12.90 0.87 13.00 9 2.17 1206.88 5 29 153996137902634436480 11.74 0.87 32.65 9 5.07 7026.13 4 212 643260079227925564160 11.73 2.13 1.98 13 0.12 8.12 7 188 -335700273723303646464 11.96 0.84 32.64 15 5.26 8991.52 8 242 75551538941421122304 13.10 1.71 1.23 16 0.06 1.31 11 257 -212924378502398307840 12.62 1.25 1.10 14 0.00 0.00 14 232 -166724929671747826816 11.97 1.31 13.55 9 1.87 1105.66 2 352 -113972130276695660288 9.88 2.18 1.58 8 0.00 0.00 3 227 655163343815632946432 12.92 1.25 1.21 10 0.00 0.00 5 194 -502926732831673735168 9.56 0.72 1.41 14 0.00 0.00 9 230 -122929487348818749824 11.21 0.78 1.54 12 0.00 0.00 4 233 -5939821616976287104 9.91 1.43 1.50 8 0.00 0.00 2 181 183458393840965496960 12.07 0.63 11.71 9 3.07 2046.73 2 172 9 −15 −10 −5 0 5 10 15 perihelion time (t ph ) / Myr pe r i he li on d i s t an c e ( d ph ) / p c Fig. 4.
As Figure 3, but now showing the individual surrogates from theorbital integrations (but just plotting 100 per star rather than all 2000 perstar). Surrogates from stars with median perihelion parameters outsidethe plotting range (and so not shown in Figure 3) are shown here.
The uncertainties in t ph and d ph are correlated. This can beseen in Figure 4, where we plot the individual surrogates insteadof the summary median and axis-parallel error bars. Each staris generally represented as an ellipsoidal shape pointing roughlytowards t ph = d ph = perihelion speed (v phmed ) / km s - pe r i he li on d i s t an c e ( d ph m ed ) / p c l l llll ll lll ll l lll llll ll lll l ll ll llll l ll l lll lll l lll l ll ll lll l ll l ll l ll l l llll l ll lll ll ll ll ll l lll ll ll ll ll l ll l ll lll llll lll llll lll l l ll lll l ll ll ll ll l ll ll l l llll l l l lll l ll ll l lll ll llll ll ll lll llll l l lll l ll l ll ll l ll lll ll l ll ll l lllll ll l ll ll ll ll ll ll ll ll ll ll ll ll l llll l llll ll lll ll ll lllll lllllll ll l lll lll l ll ll l lll ll ll l l ll lll llll lll ll ll ll ll ll lll lll ll lll ll ll lll ll lll l ll ll ll ll lll ll l ll llll lll ll l ll l lll lll lll l ll l lllll lll l l llll lll ll lll ll ll ll lll l lll lll lllll l ll lll l lll l l ll l ll lllll ll lll ll l llll llll l ll ll llll lll l l l llll l lll ll ll l lll lll ll l l lll l lll lll l lll l ll ll ll ll lll l l ll ll lll l ll ll l ll ll ll ll lllll lll lll llll ll ll ll ll ll ll lll l ll l ll lll ll l lll lll ll ll lll ll l ll lll l lll lll lllll l llll ll llll ll l llll l ll l l lll lll l ll lll lll l lll lll lll l l l ll ll l ll ll ll l ll ll lll l llll ll lll ll l lll llll ll lll l ll ll llll l ll l lll lll l lll l ll ll lll l ll l ll l ll l l llll l ll lll ll ll ll ll l lll ll ll ll ll l ll l ll lll llll lll llll lll l l ll lll l ll ll ll ll l ll ll l l llll l l l lll l ll ll l lll ll llll ll ll lll llll l l lll l ll l ll ll l ll lll ll l ll ll l lllll ll l ll ll ll ll ll ll ll ll ll ll ll ll l llll l llll ll lll ll ll lllll lllllll ll l lll lll l ll ll l lll ll ll l l ll lll llll lll ll ll ll ll ll lll lll ll lll ll ll lll ll lll l ll ll ll ll lll ll l ll llll lll ll l ll l lll lll lll l ll l lllll lll l l llll lll ll lll ll ll ll lll l lll lll lllll l ll lll l lll l l ll l ll lllll ll lll ll l llll llll l ll ll llll lll l l l llll l lll ll ll l lll lll ll l l lll l lll lll l lll l ll ll ll ll lll l l ll ll lll l ll ll l ll ll ll ll lllll lll lll llll ll ll ll ll ll ll lll l ll l ll lll ll l lll lll ll ll lll ll l ll lll l lll lll lllll l llll ll llll ll l llll l ll l l lll lll l ll lll lll l lll lll lll l l l ll ll l ll ll ll l ll ll lll l llll ll lll ll l lll llll ll lll l ll ll llll l ll l lll lll l lll l ll ll lll l ll l ll l ll l l llll l ll lll ll ll ll ll l lll ll ll ll ll l ll l ll lll llll lll llll lll l l ll lll l ll ll ll ll l ll ll l l llll l l l lll l ll ll l lll ll llll ll ll lll llll l l lll l ll l ll ll l ll lll ll l ll ll l lllll ll l ll ll ll ll ll ll ll ll ll ll ll ll l llll l llll ll lll ll ll lllll lllllll ll l lll lll l ll ll l lll ll ll l l ll lll llll lll ll ll ll ll ll lll lll ll lll ll ll lll ll lll l ll ll ll ll lll ll l ll llll lll ll l ll l lll lll lll l ll l lllll lll l l llll lll ll lll ll ll ll lll l lll lll lllll l ll lll l lll l l ll l ll lllll ll lll ll l llll llll l ll ll llll lll l l l llll l lll ll ll l lll lll ll l l lll l lll lll l lll l ll ll ll ll lll l l ll ll lll l ll ll l ll ll ll ll lllll lll lll llll ll ll ll ll ll ll lll l ll l ll lll ll l lll lll ll ll lll ll l ll lll l lll lll lllll l llll ll llll ll l llll l ll l l lll lll l ll lll lll l lll lll lll l l l ll ll l ll ll ll l ll ll ll Fig. 5.
Median perihelion velocities from the orbit integrations for thoseencounters shown in Figure 3 The velocity axis is a logarithmic scale.
The perihelion speeds for encounters with d medph < medph <
100 km s − . The very fastencounters are almost entirely due to stars with large radial ve-locities.The masses of the encountering stars are listed in Table 2.The 5th, 50th, and 95th quantiles of the mass distribution overthe entire filtered sample are 0.53, 0.91, and 1.64 M (cid:12) respec- Article number, page 5 of 12 & A proofs: manuscript no. stellar_encounters_gdr2 −15 −10 −5 0 5 10 15 perihelion time (t phmed ) / Myr pe r i he li on d i s t an c e ( d ph m ed ) / p c l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l Fig. 6.
As Figure 3, but now plotting each star as a circle, the area ofwhich is proportional to M / (v medph d medph ). −15 −10 −5 0 5 10 15 perihelion time (t phmed ) / Myr pe r i he li on d i s t an c e ( d ph m ed ) / p c l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l Fig. 7.
As Figure 3, but now plotting each star as a circle, the area ofwhich is proportional to M / (v medph ( d medph ) ). tively. The lack of massive stars in the sample is a consequenceof the T e ff filtering on radial velocities published in Gaia DR2.Close encounters are interesting not least because they per-turb the Oort cloud. The degree of perturbation, or rather the im-pulse which comets can receive, depends not only on the distanceof the encountering star, but also on its speed and the mass. Ac-cording to the simple impulse approximation (Öpik 1932; Oort1950; Rickman 1976; Dybczynski 1994) the impulse transfer isgiven by M v − d − α ph where α is 1 for very close encounters (onthe order of the comet–Sun separation), and 2 otherwise. Theimpulse of the encountering stars is visualized in Figures 6 andFigures 7 for α equal to 1 and 2 respectively. Regardless of whichimpulse approximation we use, it is the closest encounters whichhave the greatest impact. Figure 8 shows the colour–absolute magnitude diagram (CMD)for the close encounters on the assumption of zero interstellarextinction. The upper panel shows the 26 non-bogus encounterswith d medph < BP−RP [mag] G + l og v + [ m ag ] l l l ll l lllll ll l llll llll lll l lll llll ll l lll ll lll l ll ll l lll ll l ll ll ll ll l l ll ll l ll l ll lll l ll lllll ll ll ll l ll l ll llll lll ll ll ll ll ll l ll ll ll llll l l lll lll ll l lll ll ll l ll ll ll lll ll llll lll lll l ll l ll l ll l l ll ll l l ll lll ll ll l ll ll ll ll lll l l ll ll lll l l ll l lll lll ll lll l ll ll lll lll ll llll ll lll ll ll ll l l lll ll llll l ll ll ll lll ll lll l lll l llll lll lll ll lll ll l lll lll ll ll l l l lll l l l l lll ll ll lll ll l lll ll llll ll ll l ll ll ll ll lll lll ll ll lll ll l llll lll lll ll llll lllll llll l l ll l l llll ll ll ll l ll l ll ll lll l l llll ll ll l l lll l lll l ll lllll l l lll l ll ll l ll l l ll lll l ll ll ll ll l ll l l ll lll l l l ll l lll llll llll ll l lll l l ll ll ll l l ll ll ll lll l ll ll l llllll l l ll ll ll llll l ll l ll ll llll l ll l ll l ll l lll ll l lll l lll l lll l ll ll l ll lll ll lll ll l l lll lll ll l ll l l lll ll l lll l ll l lll ll ll lllll lll ll l ll llll ll ll ll l l lll l l l l lllll ll ll ll l lll lll ll l l l llllll l ll llll ll ll ll ll ll l ll ll l l lll ll ll ll lll lll l ll lll l ll ll l lll l l lll l l ll l ll lll ll l ll lll lll ll llll l ll l l lll l ll ll lll l ll l lll l ll ll ll ll l ll l lll l ll llll ll ll l l l ll l ll lll lllll l ll l lll ll ll ll ll l llll ll ll llll l lll ll ll ll lll ll l lll l lll l l ll l llll lll l ll ll lll l ll ll l lll lll l ll l lll ll ll lll l lll ll l ll l ll l ll lll ll ll llll lll l ll llll l l llll l ll ll ll l ll lll ll lll lll ll ll lll ll l ll lll ll lll ll l ll l llll lll ll ll ll lll llll l ll lll l lll ll l ll l ll ll ll ll ll l lll ll lll l ll lll llll l ll l ll lll ll lll ll ll ll l ll l l ll l lll l ll lllll l l l l llll ll lll l l ll l ll l ll l ll ll ll ll l lll lll l ll ll ll lll l ll ll ll lll lll l ll ll l l ll l ll l ll ll llll l lll l ll l ll ll l llll lll ll l ll l l ll l l ll l l lll lll l l l l ll lll llll ll l ll ll l ll l ll lll ll ll ll lll ll lll l lll ll ll l ll ll l ll l lll l l ll ll ll ll l ll lll lll l ll l ll l ll ll l ll ll l ll ll l lll ll ll ll ll ll ll ll lll l lll lll ll llll llll l l ll ll ll ll ll lll l l ll l ll ll ll ll lll lll lll llll ll l lll ll llll lll lll l lll lll l l ll ll lll ll lll l ll ll l ll ll ll lll l l l l ll lll ll l lll ll lll l l lll lll llll ll l l ll l ll l lll lll ll l l lll lll l ll ll ll lll ll l lll ll llll l l l lll ll l ll l l l ll ll ll ll lll l ll l ll ll l ll ll ll l llll l llll lll ll lll l l l ll l l llll ll l lll l l ll ll l ll lll l lllll ll llll lll lll ll l ll llll ll lll l lll l l ll ll ll l l ll lll l lll ll l l llll ll lll ll ll l lll llll l lll l l l l ll llll lll l l lll lll ll l l ll ll l l ll ll lllll ll llll ll ll lll l llll l ll ll l ll llll ll l llll l lll llll lll ll l l ll l l lll lll lll ll l l lll ll ll ll ll lll ll ll l llll llll ll lll ll l ll l lll ll ll l ll l lll lll l l ll ll lll l llll l llll ll lll lll lll ll ll l l ll lll l l ll llll lll lll l lll l lllll u BP−RP [mag] G + l og v + [ m ag ] l l l ll l lllll ll l llll llll lll l lll llll ll l lll ll lll l ll ll l lll ll l ll ll ll ll l l ll ll l ll l ll lll l ll lllll ll ll ll l ll l ll llll lll ll ll ll ll ll l ll ll ll llll l l lll lll ll l lll ll ll l ll ll ll lll ll llll lll lll l ll l ll l ll l l ll ll l l ll lll ll ll l ll ll ll ll lll l l ll ll lll l l ll l lll lll ll lll l ll ll lll lll ll llll ll lll ll ll ll l l lll ll llll l ll ll ll lll ll lll l lll l llll lll lll ll lll ll l lll lll ll ll l l l lll l l l l lll ll ll lll ll l lll ll llll ll ll l ll ll ll ll lll lll ll ll lll ll l llll lll lll ll llll lllll llll l l ll l l llll ll ll ll l ll l ll ll lll l l llll ll ll l l lll l lll l ll lllll l l lll l ll ll l ll l l ll lll l ll ll ll ll l ll l l ll lll l l l ll l lll llll llll ll l lll l l ll ll ll l l ll ll ll lll l ll ll l llllll l l ll ll ll llll l ll l ll ll llll l ll l ll l ll l lll ll l lll l lll l lll l ll ll l ll lll ll lll ll l l lll lll ll l ll l l lll ll l lll l ll l lll ll ll lllll lll ll l ll llll ll ll ll l l lll l l l l lllll ll ll ll l lll lll ll l l l llllll l ll llll ll ll ll ll ll l ll ll l l lll ll ll ll lll lll l ll lll l ll ll l lll l l lll l l ll l ll lll ll l ll lll lll ll llll l ll l l lll l ll ll lll l ll l lll l ll ll ll ll l ll l lll l ll llll ll ll l l l ll l ll lll lllll l ll l lll ll ll ll ll l llll ll ll llll l lll ll ll ll lll ll l lll l lll l l ll l llll lll l ll ll lll l ll ll l lll lll l ll l lll ll ll lll l lll ll l ll l ll l ll lll ll ll llll lll l ll llll l l llll l ll ll ll l ll lll ll lll lll ll ll lll ll l ll lll ll lll ll l ll l llll lll ll ll ll lll llll l ll lll l lll ll l ll l ll ll ll ll ll l lll ll lll l ll lll llll l ll l ll lll ll lll ll ll ll l ll l l ll l lll l ll lllll l l l l llll ll lll l l ll l ll l ll l ll ll ll ll l lll lll l ll ll ll lll l ll ll ll lll lll l ll ll l l ll l ll l ll ll llll l lll l ll l ll ll l llll lll ll l ll l l ll l l ll l l lll lll l l l l ll lll llll ll l ll ll l ll l ll lll ll ll ll lll ll lll l lll ll ll l ll ll l ll l lll l l ll ll ll ll l ll lll lll l ll l ll l ll ll l ll ll l ll ll l lll ll ll ll ll ll ll ll lll l lll lll ll llll llll l l ll ll ll ll ll lll l l ll l ll ll ll ll lll lll lll llll ll l lll ll llll lll lll l lll lll l l ll ll lll ll lll l ll ll l ll ll ll lll l l l l ll lll ll l lll ll lll l l lll lll llll ll l l ll l ll l lll lll ll l l lll lll l ll ll ll lll ll l lll ll llll l l l lll ll l ll l l l ll ll ll ll lll l ll l ll ll l ll ll ll l llll l llll lll ll lll l l l ll l l llll ll l lll l l ll ll l ll lll l lllll ll llll lll lll ll l ll llll ll lll l lll l l ll ll ll l l ll lll l lll ll l l llll ll lll ll ll l lll llll l lll l l l l ll llll lll l l lll lll ll l l ll ll l l ll ll lllll ll llll ll ll lll l llll l ll ll l l ll llll ll ll l lll llll ll lll llllll ll llll ll ll ll ll lll ll l ll l l lll lll ll l ll lllll ll l llll l ll llll ll lll l llll lllll ll l ll l lllll l lll ll ll l ll llll llll ll l lll l ll ll ll ll l ll ll lll ll ll lll l l lll ll lll l l lllll l l ll lllllll l llll ll l lll l ll lll lll l ll l ll l llll ll llllll ll llll lll llll ll l ll l lll ll ll l lllll l ll l lll lll ll llll l llll lll l llll lll ll ll l l llll lll lll ll lll ll ll lllll lll l lll l l llll l ll l llll ll lll ll lll lll ll llll l l ll l ll l ll l ll lllll lll l llll ll l lll ll ll lll lll ll ll lll ll ll lll ll llll lll lll llll l llll ll lll ll l lll l l lll l ll llll lll ll l l ll ll l lllll l l llll ll ll lll ll ll lll l ll ll llll ll lll ll ll l ll ll l llllll lllll llll lll ll l l ll lll lll ll l ll ll l lll ll lllll lll ll ll ll l l llll lll lll ll lll l lll l lll l l l lllll ll l l ll l lll ll llll ll ll ll l l llllll llll ll ll lllll lll ll lll ll ll llllll ll l l l l ll l l ll l ll ll l ll ll lll l ll l ll lll lll lll lll l ll ll lll lll ll l l lll ll ll ll ll ll ll lll lll lll l ll ll ll lllll llll l ll l l lll ll ll lll ll l ll llll ll ll l llll ll lllll l ll ll l lll ll ll l ll llll l ll ll lll lll lllll ll ll lll l l ll ll llll ll llll ll ll ll lll l lll l ll ll lll l l ll l llll lll ll l ll lll lll l ll lll ll ll lll lll l lllll ll lll lllll l lllll l ll ll l l l lll ll l ll lll l ll ll l ll ll ll ll ll l ll llll ll ll l lllll ll lll ll l lllll ll l ll l ll lll lll ll llll ll ll ll ll lll ll l ll lll ll llll ll ll ll ll lll ll ll ll lllll l l lll ll l lllll llll ll lll l ll ll ll lll l ll lll ll ll ll ll ll lll ll ll ll ll lll lll ll ll ll ll ll lllll l ll l ll lll l ll l l ll llll ll ll lll l l llll l lll l lll l ll lllll lll ll ll ll lll l lll l llll lll ll ll lll l lll ll ll ll ll lll ll l lll l lll lll l ll l ll l l ll lll l ll ll l ll ll lll l lll lll llll l l ll lll l lll ll ll l ll l lll ll l l l lll lll ll ll ll ll llll ll llll l llll l llllll ll l ll lll l ll l lll ll lll ll lll l ll ll llll ll ll ll lll ll lll llll llll ll l ll ll llll llll l ll l l ll lll l llll ll ll ll lll ll l ll ll ll lll l l l ll ll ll ll ll l lll lll ll ll llll l lll ll l l ll lll ll lll l lll l lll lll ll lllll l lll lll l ll l lll lll l llll l ll ll lll llll l lll lll ll llllll l lll l ll llll lll l l lll l llll lll lll l lll lll l ll llll l ll ll l ll l ll lll ll l lll ll lll l lll lllll ll ll ll l ll ll lllll lll llllll ll lll lll l l llll ll lll lllll llll lll l ll ll l l lll lll lll lll ll llll lll lll ll lll lll lll lll l lll lll llll ll l lll l lll lll l lll lll ll ll l ll l ll ll ll l l ll l llll lllllll ll ll lll l ll l llll ll ll ll ll l ll l l lll llllll l ll l ll lll ll ll ll l ll lll ll lll ll lll ll ll ll l lll l ll l llll ll llll lll llllll lll lll ll llll l llll l lll lll l llll l ll ll ll lll l lll lll ll l lll l ll l l ll llll l l lll llll lll ll l ll l lll l llll l lll l l llll ll lll ll l l lll ll ll l ll ll ll ll l lll ll ll l ll ll lll lll lll l lll lll ll l lll lll lll l llll ll l lll ll ll ll ll ll ll lll l l lll ll l ll l ll ll l lll ll ll llll ll ll ll ll l lll ll lll ll ll l ll ll l lll l llll ll llll ll l ll lll lllll ll ll llll lll ll l ll lll ll lll ll lll ll lll l ll lll l lll ll ll llll ll lll ll l l llllll ll ll ll l l ll l l lll ll lll l l ll ll llll l lllll ll lll l l lll ll lllll ll ll lll l ll lll ll ll llll ll l ll ll l lll l ll ll ll lll l ll ll lll ll l ll lll l l ll l ll ll lll l llll ll ll ll lll l ll l ll ll ll lll l lll l lll ll l l lll ll l lll ll ll ll llll ll ll l l lll l l l lll l l lll l ll l ll ll lll lll ll l lllll l ll lll l ll llll ll ll lll ll l ll ll ll l ll ll lll ll llll lllll l lll l ll l ll lll llll ll l llll l l ll lll ll lll l ll ll l lll ll ll lll l ll ll ll lll ll ll lll lll ll llll l ll ll ll l l ll lllll ll l l ll ll ll ll lll l lll l llll lll ll lll lll lll ll l ll llll l llll l ll l l ll lll ll l lll ll ll ll l l ll lll ll l lll ll lll ll ll l lll l l lll ll l ll ll l lll ll l lllll ll llll ll l lll lll ll ll ll l lllll lll llll ll llll l ll lll ll llll l ll lll ll lll l l lll llll ll ll lll ll l l lll l lll l ll lll lll lll ll l lll l l llll l lll ll l llll ll l lll ll lll ll llll llll lll lll ll ll llll lll ll ll l l ll ll llll l ll lll ll lll lllll l lll l lll l ll lll lll llll ll ll ll ll lll lll l l llll ll lll ll ll lll lll l ll ll ll lll llll l ll ll l llll lll ll llll lll l ll ll ll l ll l llll ll lll ll lll llll ll l lll ll ll ll l ll llll lll l ll ll ll lll lll l ll ll ll llll ll ll l llll ll ll ll ll l llll ll ll ll ll l l ll ll ll llll lll ll lll l lll ll lll ll l ll ll ll l ll ll lllll l l lll lll l lll llll lll lll l l l lllllll l ll lll lllll ll ll ll ll lll ll ll ll ll lll l ll lll lll ll l ll ll lll lll ll l ll l llll lll ll ll ll l lll ll llll lll l ll l l l l ll l l lll lll lll ll l l lll ll ll ll ll lll ll ll l llll llll ll lll ll l ll l lll ll ll l ll l lll lll l l ll ll lll l llll l llll ll lll lll lll ll ll l l ll lll l l ll llll lll lll l lll l lllll u Fig. 8.
Colour–magnitude diagram for stellar encounters coloured ac-cording to the unit weight error u in Table 3. The vertical axis equals M G assuming extinction is zero. For orientation, the grey lines are unred-dened solar metallicity PARSEC isochrones for 1 Gyr (solid) and 10 Gyr(dashed) from Marigo et al. (2017), the grey points are white dwarfswithin 20 pc of the Sun identified by Hollands et al. (2018) (many lieoutside the range of the plot), and the small blue dots show all sources inGaia DR2 with (cid:36) >
50 mas (plotted on the bottom layer, so most are ob-scured by coloured and grey points, especially in the lower panel). Theupper panel shows all reliable (non-bogus) encounters with d medph < u on higher layers so as to betterlocate larger values of u . The colour bar spans the full range of valuesin both panels. error u >
10 (the others have u <
2) and lie some way below andleft of the main sequence (MS), yet well above the white dwarf(WD) sequence. The dimensionless quantity u measures howwell the five-parameter solution describes the astrometric data(see section 2.2). If we assumed that all the encounters were sin-gle main sequence stars, then one interpretation of this diagramis that the parallaxes of these 11 sources are wrong; specificallythat they are all overestimated by one or two orders of magni-tude (i.e. the sources should be much further away). This seemsunlikely, although it is the case that a search for encounters willpreferentially include sources with spuriously large parallaxesrather than sources with spuriously small parallaxes (because theformer are more likely to have apparent close encounters). As wesee a smaller proportion of sources between the MS and WDsamong nearby stars (blue points in Figure 8), it at first seems un- Article number, page 6 of 12ailer-Jones et al.: Close encounters to the Sun in Gaia-DR2 likely that we would get such a large proportion of such sourcesamong the closest encounters (upper panel). The magnitude ofthis selection e ff ect is impossible to estimate without knowingwhich parallaxes are actually spurious. On the other hand, if welook instead at all encounters (lower panel), then we find that just14% of all encounters within the plotted range lie below the MS,whereas this figure is 24% for the nearby stars. In other words,there is in fact a decreased tendency for encounters to be belowthe MS compared to nearby stars. This comparison is not ideal,however, because some of these nearby stars are WDs – essen-tially absent from the encountering sample due to the magnitudelimit – and the encounters are drawn from a much larger volumeWe noted in section 2.2 that a large value of u does notmean the parallax is wrong. Indeed, a relatively poor fit of thefive-parameter solution is in principle expressed by larger val-ues of the astrometric uncertainties (and these are quite corre-lated with u ). Here, the 11 sources in the upper panel of Fig-ure 8 have σ ( (cid:36) ) between 0.5 and 1.1 mas, whereas the other 15have σ ( (cid:36) ) < .
06 mas. We have already accommodated the un-certainties in the identification of the close encounters. Further-more, although the parallax uncertainties are larger, they pro-duce uncertainties in the absolute magnitude of no more than0.08 mag, which is smaller than the size of the points plot-ted in Figure 8. Contamination of the astrometric solution bynearby (sub-arcsecond) sources is possible, and although mostof the sources are near to the Galactic plane (as are unques-tioned sources), they are all bright, so less likely to be a ff ected bycrowding. Visual inspection of non-Gaia images of these targetsshows no obvious contaminants. They do generally have fewervisibility periods in the solutions, however: 9, 10, or 11, as op-posed to 8–18 for the sources lying near the main sequence.Given their location in the CMD, it is possible that someor all of these 11 sources are subdwarfs or MS-WD binaries.The latter, of which many have been discovered (e.g. Rebassa-Mansergas et al. 2016), can in principle lie in much of the regionbetween the MS and WD sequence. Such binarity could poten-tially explain elevated values of u , although only if the periodswere short enough and the amplitudes large enough to be de-tectable over Gaia DR2’s 22-month baseline. Even then the par-allax should be reasonably accurate, because the di ff erent astro-metric displacements at di ff erent times should not create a largebias (Lennart Lindegren, private communication). Furthermore,four of the 11 sources have much larger radial velocity uncertain-ties than the other encounters in Table 2 ( >
14 km s − as opposedto < − ). As this uncertainty is computed from the stan-dard deviation of measurements at di ff erent epochs (equation 1of Katz et al. 2018), this could be an indication of binarity.In conclusion, the correlation of the CMD location with u issuspicious, and would be consistent with a selection e ff ect thatpreferentially includes spuriously large parallaxes in a search forstellar encounters. On the other hand, we do not have specificevidence to claim that these stars have spurious solutions or un-accounted for errors, and it is quite possible that some or all areMS-WD binaries. Spectroscopy may help identify the nature ofthese sources, and astrometry over a longer timebase (e.g. fromGaia DR3) should help to determine if the astrometric solutionsare good. We find 31 stars with d medph < u , astrometric_excess_noise and / or astrometric_excess_noise_sig , as can be seen in Table 3.Only one of the encounters has been discovered by the previousstudies listed in the introduction. 21 do not have Hipparcos orTycho IDs, so could not have been discovered in papers 1 or 2.Of the remaining nine stars, seven have obvious Tycho matchesbut were not found in paper 2 for reasons discussed below. Thelast two stars are problematic and are also discussed below.The closest approaching star, and the only one to havebeen found already, is Gaia DR2 4270814637616488064, bet-ter known as Gl 710 (Tyc 5102-100-1, Hip 89825). This is aK7 dwarf known from many previous studies to be a very closeencounter. In paper 1 we found a median encounter distance of0.267 pc (90% CI 0.101–0.444 pc). A much smaller – and moreprecise – proper motion measured by TGAS gave a closer en-counter, with a median distance of 0.076 pc (90% CI 0.048–0.104 pc) (paper 2 and Berski & Dybczy´nski 2016). Gaia DR2has slightly decreased this distance – 0.0676 pc (13 900 AU) –and has narrowed the uncertainty – 90% CI 0.0519–0.0842 pc(10 700–17 400 AU) – although this is statistically consistentwith the TGAS result. This is well within the Oort cloud. Due tothe slightly more negative radial velocity published in Gaia DR2than used in the earlier studies, the encounter time is slightlyearlier (but no more precisely determined). We note, however,that the radial velocity in Gaia DR2 comes from just two Gaiafocal plane transits, the minimum for a value to be reported inthe catalogue. Despite its low mass, Gl 710 imparts the biggestimpulse according to either impact approximation (Figures 6and 7), in part because of its low velocity. Using the same data,de la Fuente Marcos & de la Fuente Marcos (2018) get a slightlycloser mean (and median) distance of 0.052 pc (with a symmetricstandard deviation of 0.010 pc). Their estimate is close to whatone gets for a constant gravitational potential, namely 0.0555 pc(computed by propagating the set of surrogates using the LMAand taking the mean). This suggests they used a very di ff erentpotential from ours. Note that ignoring the parallax zeropointand / or ignoring the astrometric correlations changes the estimateby less than 0.0001 pc, so these are not the cause of the di ff er-ence.The second closest approach in Table 2, Gaia DR2955098506408767360, is also the one with the second largestimpulse. It is one of the most massive of the closest en-counters, with a perihelion in the relatively recent past. Fourfainter 2MASS sources lie within 8”, but there is no goodreason to think these have interfered with the Gaia solution.The most recent reliable encounter in the table is Gaia DR23376241909848155520, which has a 50% chance of havingpassed within 0.5 pc, around 450 kyr ago.Among the 21 sources without Tycho matches, most havematches to 2MASS (and a couple of others to other surveys).A few of these are in very crowded regions (e.g at low Galac-tic latitude) and have close companions, but in all but one casethese companions are much fainter, and there is no specific ev-idence to suggest a problem with the Gaia measurements. Theone exception is Gaia DR2 5700273723303646464, which has acompanion 4” away that is just 1.6 mag fainter in the J-band (thebest match to this in Gaia DR2 is 0.7 mag fainter in the G band,and has a completely di ff erent parallax and proper motion). Thisis a borderline case. Its measurements are not suspicious, butsome of the quality metrics make it questionable. We decide toflag this as bogus.There are seven encounters with obvious Tycho matches thatwere not found as encounters in paper 2, either because theylacked a radial velocity, or because their astrometry has changed Article number, page 7 of 12 & A proofs: manuscript no. stellar_encounters_gdr2 substantially from TGAS to Gaia DR2. These are as follows(with notes on the problematic cases). – Gaia DR2 510911618569239040 = TYC 4034-1077-1 . – Gaia DR2 154460050601558656 = TYC 1839-310-1. Thisis rather faint, G = .
37 mag, and has a rather high radialvelocity based on just 3 transits. It also has no BP-RP colourand a large astrometric_excess_noise_sig of 360, sowe consider this encounter to be bogus. – Gaia DR2 1791617849154434688 = TYC 1662-1962-1. – Gaia DR2 1949388868571283200 = TYC 2730-1701-1. Thelarge radial velocity is based on just two transits, so we con-sider this to be bogus. – Gaia DR2 3972130276695660288 = TYC 1439-2125-1 = GJ 3649. – Gaia DR2 2926732831673735168 = TYC 5960-2077-1. – Gaia DR2 2929487348818749824 = TYC 5972-2542-1.The final two sources have ambiguous matches. Simbad listsfour objects within 5” of Gaia DR2 5231593594752514304 ( G = .
03 mag), although in reality these four are probably just twoor three unique sources. Gaia DR2 itself identifies one of these asGaia DR2 5231593599052768896, which with G = .
12 mag isalmost certainly Hip 53534 with V = .
21 mag. Whichever otherstar our close encounter is, its RVS spectrum is surely contami-nated by this much brighter star just a few arcseconds away, sug-gesting that the whoppingly large radial velocity of −
715 km s − measured by RVS is in fact spurious. We therefore consider thisencounter to be bogus.Finally, Gaia DR2 939821616976287104 would appear tomatch with TYC 2450-1618-1. However the 2MASS imageclearly shows this to be a double star, and this surely has contam-inated the RVS spectrum. In that case we do not trust the highradial velocity reported (568 km s − ) and conservatively considerthis to be a bogus encounter too. In paper 2 we found two stars with d medph < d medph = .
103 pc (90% CI 0.041–0.205 pc). It was deemed du-bious due to the inconsistency of its magnitude and Hipparcosparallax with its apparent spectral type. This suspicion is con-firmed with Gaia DR2, which has a much smaller parallax of1 . ± .
027 mas compared to 146 . ± .
81 in Hipparcos-2(the proper motions and radial velocities are consistent within1 σ ). This now places the closest approach at 420 pc. It was con-jectured in paper 1 that the cause for the incorrect parallax is abinary companion at a very problematic separation for the Hip-parcos detectors. Hip 63721. Its best match, Gaia DR21459521872495435904, has a parallax of 5 . ± .
04 mas,compared to 216 . ± .
53 mas in Hipparcos. This Hipparcossolution was doubted in paper 1 on the grounds of inconsistencywith spectral type and it being a double star.Hip 91012. The parallax and proper motion of its best match,Gaia DR2 4258121978471892864, are consistent with Hippar-cos, but the Gaia DR2 radial velocity is now much smaller, at − . ± .
61 km s − , compared to the RAVE value of − . ± . > − ) nearby (113 mas) star – Gaia DR2 andHipparcos agree – but the Gaia DR2 radial velocity is muchsmaller: 21 . ± .
15 km s − compared to the implausibly largevalue of − . ± . . ± . − versus 0 . ± . − ),meaning it no longer comes near to the Sun.Hip 55606. The proper motion of its best match, Gaia DR23590767623539652864, is consistent with Hipparcos (parallaxless so), but the Gaia DR2 radial velocity is now much smaller, at − . ± .
70 km s − , compared to the implausibly large value of − . ± . σ of the Hipparcosastrometry (which was rather imprecise, with a proper motionconsistent with zero). Using the (large) RAVE radial velocityof 368 km s − from paper 1 with the Gaia DR2 astrometry, theLMA gives a closest approach of 2.1 pc 660 kyr in the past.Hip 103738 (gamma Microscopii). This G6 giant was thepotentially most massive encounter coming within 1 pc in pa-per 1, found to have d medph = .
83 pc (90% CI 0.35–1.34 pc)based on the Hipparcos-2 proper motion of 1 . ± .
35 mas yr − .(A second entry in the catalogue based on the larger Tycho-2proper motion put it at 3.73 pc; 90% CI 2.28–5.22.) This staris Gaia DR2 6781898461559620480 with G = .
37 mag. Theproper motion in Gaia DR2 is ten times (17 σ ) larger, 17 . ± . − , with no suggestion from the quality metrics thatthis is a poor solution. (Both the parallaxes and radial veloci-ties agree with 1 σ .) Either there is a significant, unaccounted foracceleration, or the proper motion in Hipparcos-2 or Gaia DR2(or both) is wrong. Using the Gaia DR2 data we find the en-counter to be at t medph = − − − d medph = . α Centauri A). This is not in Gaia DR2 on ac-count of its brightness.Hip 71681 ( α Centauri B). This is not in Gaia DR2 on ac-count of its brightness.Hip 70890 (Proxima Centauri). This is Gaia DR25853498713160606720. It has no radial velocity in Gaia DR2,but the parallax and proper motion agree with Hipparcos towithin 2 σ . Given their high precision, the perihelion parametersusing the new astrometry together with the old radial velocityare consistent with the result in paper 1 (now more precise).Hip 3829 (van Maanen’s star). This is Gaia DR22552928187080872832. It has no radial velocity in Gaia DR2,but the parallax and proper motion agree with Hipparcos Article number, page 8 of 12ailer-Jones et al.: Close encounters to the Sun in Gaia-DR2 to better than 1 σ . Given their high precision, the perihelionparameters using the new astrometry with the old radial velocityare consistent with the result in paper 1 (now more precise).Hip 42525. Two Gaia DR2 sources with the same par-allax, proper motion, and radial velocity as each other(within their uncertainties) match this. They are Gaia DR2913394034663258752 and Gaia DR2 913394034663259008.The proper motion and radial velocity agree with the data usedin paper 1, but the Gaia DR2 parallax of 6 . ± . . ± . . ± .
4. Completeness correction and the encounter rate
Gaia DR2 does not identify all encounters, primarily because ofits faint-end magnitude limit for radial velocities. To computean encounter rate over a specified time and distance window wemust correct for this. This was done in paper 2 using a simpleanalytic model. Although insightful and easy to work with, thatmodel made severe assumptions of isotropic spatial and homoge-neous velocity distributions. Here we take the same conceptualapproach to building a completeness map, but replace the ana-lytic model with a more realistic Galaxy simulation, in which thespatial and kinematic distributions are also self-consistent. Wefirst simulate the positions and velocities of all stars in the nearbyGalaxy (“mock full Galaxy”), and trace the orbits backwards andforwards in time (with LMA) to determine the distribution ofencounters in perihelion coordinates. Call this F mod ( t ph , d ph ), thenumber of encounters per unit perihelion time and distance. Wethen repeat this including only the stars which would have full6D kinematic information in Gaia DR2 (“mock Gaia-observedGalaxy”), then trace orbits to give F exp ( t ph , d ph ). The ratio ofthese two distributions gives the completeness map, C ( t ph , d ph ),which can be interpreted as the fraction of encounters at a spe-cific perihelion time and distance that will actually be observed(in our sample). To generate our mock full Galaxy we use Galaxia (Sharma et al.2011) to sample stars following a Besancon-like (Robin et al.2003) model distribution of the Galaxy, in a similar way to howwe generated the mock Gaia DR2 catalogue described in Rybizkiet al. (2018). The main di ff erence is that we now apply no mag-nitude cut, but instead a distance cut, to keep the computationsmanageable. We use 3 kpc, as only stars with radial velocitiesabove 200 km s − would reach the Sun within 15 Myr. (Our fi-nal completeness-corrected encounter rate only uses the modelout to 5 Myr, so only mock stars beyond 3 kpc with radial ve-locities above 600 km s − will be neglected.) We checked on asubsample of the data that increasing this limit to 10 kpc madenegligible di ff erence to the distribution of encounters. We thencomputed the perihelion parameters of all stars using the LMA.Although the LMA is not an accurate model for long paths inthe Galaxy, the inaccuracy introduced by this for the distribu-tion as a whole will be small. The impact is further diminished −15 −10 −5 0 5 10 15 perihelion time (t ph ) / Myr pe r i he li on d i s t an c e ( d ph ) / p c Fig. 9.
Perihelion times and distances for star samples from the mockfull Galaxy (no Gaia selection). Each star is represented by a singlesample. given that we are ultimately interested only in the ratio of thetwo encounter distributions from these two models.Figure 9 shows the positions in perihelion parameter spaceof individual stars from the sampled mock full Galaxy. The den-sity of encounters shows essentially no dependence on t ph overthe range shown. This is di ff erent from what was found from themodel in paper 2 which showed a drop o ff to large | t ph | . Thiswas in part due to the spatial distribution adopted (and its incon-sistency with the velocity distribution adopted). In the presentmodel we also do not see the drop in density towards t ph =
0. Theencounter density in Figure 9 does show a strong dependency on d ph . Closer inspection shows that the number of encounters perunit distance varies linearly with d ph , i.e. the number of encoun-ters within d ph varies as d . This is what the simple model in pa-per 2 predicts (see section 4.2 of that paper, which also explainswhy it does not vary as d , as one may initially expect). To con-struct the completeness map we therefore replace the 2D distri-bution in Figure 9 with a 1D distribution F exp ( t ph , d ph ) = a d ph ,where the constant a is fit from the simulated data. The distribu-tion of the real encounters also shows this linear variation out toseveral pc, as shown in Figure 10.To produce the mock Gaia-observed Galaxy we query themock catalogue of Rybizki et al. (2018) using the followingADQL query SELECT parallax, pmra, pmdec, radial_velocityFROM gdr2mock.mainWHERE phot_g_mean_mag <= 12.5AND teff_val > 3550 AND teff_val < 6900
The T e ff limit simulates the limit on the radial velocities pub-lished in Gaia DR2 (Katz et al. 2018), and the magnitude limit isthe selection we will apply to the observed encounter sample tocompute the completeness-corrected encounter rate below. Wedo not specifically account for the incompleteness at the brightend in Gaia DR2 because it is poorly defined. There are so fewbright stars in the model that this is anyway a minor source oferror. This above query delivers 4.4 million stars. (As a compar-ison, the number of stars in Gaia DR2 with complete 6D kine-matic data brighter than G = . Article number, page 9 of 12 & A proofs: manuscript no. stellar_encounters_gdr2 perihelion distance (d ph ) / pc nu m be r o f en c oun t e r s pe r p c l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l ll l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l Fig. 10.
The variation of the number of observed encounters per unitperihelion distance as a function of perihelion distance (for the filteredsample for all t ph ). The open circles count the encountering stars dis-cretely, using d medph as the distance estimate for each star. The filled cir-cles count each surrogate for each encountering star separately. Thegreen line is the distribution we expect under the assumption that thedistribution is linear. It has a gradient of 50.0 pc − . (The number of en-counters within some distance is the integral of this, i.e. ∝ d .) Thedata do not follow this near to the limit (at 10 pc) because not all starsfound in the LMA-based selection actually come within 10 pc after do-ing the orbital integration. The error bars show the Poisson noise com-puted from the theoretical relation (only attached to the open points toavoid crowding). in the latter. In other words, measurement errors in the real Gaiaobservations can lead to stars being preferentially scattered intoor out of some part of the perihelion parameter space. We ac-commodate this in our model by replacing each star in the mockGaia-observed Galaxy with a set of 100 noise-perturbed sam-ples generated using a simple error model. For this we use a 4D-Gaussian PDF (the Gaia-uncertainties on RA and Declinationare negligible), with the following standard deviations, obtainedfrom inspection of the Gaia DR2 uncertainties shown in Linde-gren et al. (2018): σ ( (cid:36) ) = .
068 mas; σ ( µ α ∗ ) = .
059 mas yr − ; σ ( µ δ ) = .
041 mas yr − ; σ (v r ) = . − . We neglect corre-lations as there is no reliable model for these. Just as was donewith the surrogates for the real data, the orbits of each of thesamples are traced (but here using LMA rather than a potential).Figure 11 shows the resulting positions in perihelion param-eter space of the samples from the mock Gaia-observed Galaxy.In principle this is a model for the observations shown in Fig-ure 4. They are not the same due to shot noise, imperfect mod-elling of the (complex) Gaia DR2 selection function, and in par-ticular the fact that Galaxia is not an exact model of our Galaxy.A particular di ff erence is that the mock catalogue shows moreencounters at larger perihelion times. There are several possiblecauses for this. One possibility is di ff ering velocity distributions:if stars in the mock catalogue were generally slower, then themost distant visible stars would arrive at larger perihelion times.An identical reproduction of the observational distribution is notnecessary, however, because the completeness is dependent pri-marily on the change from the mock full Galaxy to the mockGaia-observed Galaxy.As an aside, we can use the mock Gaia-observed simula-tion to try to explain features in the observed perihelion distri-bution (e.g. Figure 3). In particular, the gap at small | t ph | might −15 −10 −5 0 5 10 15 perihelion time (t ph ) / Myr pe r i he li on d i s t an c e ( d ph ) / p c Fig. 11.
Perihelion times and distances for star samples from the mockGaia-observed Galaxy. Each star is represented by 100 samples. This isjust a window on a much larger set of encounters, so includes samples ofstars which are predominantly outside the range shown (e.g. the medianperihelion position is outside). −15 −10 −5 0 5 10 15 perihelion time (t ph ) / Myr pe r i he li on d i s t an c e ( d ph ) / p c Fig. 12.
The completeness map, C ( t ph , d ph ), which can be interpretedas the probability of observing an encounter in the selected Gaia DR2sample as a function of its perihelion time and distance. “Selected” heremeans for G < . be explained by the T e ff cut. This cut removes cool M dwarfs,which are also intrinsically faint. To be observable now, theymust therefore be relatively near, in which case they will allencounter relatively soon (past or future). M dwarfs are highlyabundant, so would dominate the encounters at small | t ph | . Yetthey have been removed by the T e ff cut. A similar argument inprinciple also applies to the brightest stars – many of which arebright because they are close by, and thus near to encounter now– which are also missing in Gaia DR2, but these are fewer, sotheir net contribution is smaller.We can now build the completeness map. We bin the dis-tribution in Figure 11, and divide by the number of samples,to get a binned distribution of the (fractional) number of en-countering stars per bin. We then divide this by the expected(continuous) distribution from the mock full Galaxy, which wasjust F exp ( t ph , d ph ) = a d ph , to give our completeness map. Thisis shown in Figure 12. The map was computed and fitted overthe range −
15 to +
15 Myr and 0 to 10 pc to improve the fit to
Article number, page 10 of 12ailer-Jones et al.: Close encounters to the Sun in Gaia-DR2 F exp ( t ph , d ph ) (by reducing the shot noise), but is only shownover a smaller d ph range (and will only be used over a nar-rower t ph range too). The completeness ranges between nearlyzero and 0.48, with an average value (over the bins) of 0.09 for | t ph | <
10 Myr and 0.14 for | t ph | < ff ered from edge e ff ects (which could be miti-gated by expanding the grid, but at the expense of sensitivity tolow completeness values at larger times). If a star is observed to be encountering at position t ph , d ph , thenthe completeness-corrected (fractional) number of encounterscorresponding to this is simply 1 / C ( t ph , d ph ). When each staris instead represented by a set of N sur surrogates, the complete-nesses at which are { C i } , the completeness-corrected number ofencounters is n cor = N sur (cid:88) i C i . (1)This same expression applies when we have a set of stars repre-sented by surrogates over some range (“window”) of t ph and d ph with corresponding completenesses { C i } . N sur is still the numberof surrogates per stars. n cor is then the completeness-correctednumber of encounters in the window.Random errors in n cor come from two main sources: (i) Pois-son noise from the finite number of encounters observed; (ii)noise in the completeness map. The first source is easily ac-commodated. The fractional number of encounters observed in awindow, n enc , is just the sum of the fraction of all surrogates perstar which lie in that window. The signal-to-noise ratio in this is √ n enc , so the standard deviation in n cor due to this alone wouldbe n cor / √ n enc . This does not include any noise due arising fromthe original sample selection (e.g. from changes in the filtering).The second source can be treated with a first order propagationof errors in equation 1 δ n cor = N sur (cid:88) i C i δ C i C i . (2)In practice it is di ffi cult to determine δ C i . We have experimentedwith varying the approach to constructing the completeness map.We approximate the uncertainties arising thereby as a constantfractional uncertainty of f c = δ C i / C i = .
1. In that case equation2 can be written δ n cor = f c n cor . Combining this with the Poissonterm for (i) we may approximate the total random uncertainty in n cor as σ ( n cor ) = n cor (cid:32) n enc + f c (cid:33) / . (3)To compute n enc we use just the filtered encounter samplewith G < . | t ph | < d ph < n enc = uncorrected rate of 46 encounters per Myrwithin 5 pc. This compares to 639 stars with t medph and d medph in thisrange (i.e. neglecting the uncertainties gives an uncorrected rateof 64 stars per Myr within 5 pc.) Applying the completeness cor-rection as outlined above yields n cor = ±
542 encounters,which corresponds to 491 ±
54 encounters per Myr within 5 pc.We could calculate the corrected rates for smaller upper limits on d ph , but such results are sensitive to the use of fewer bins in thecompleteness map and are more a ff ected by the scattering of thesurrogates. We instead scale the value found for 5 pc using theexpectation that the number of encounters within some distancegrows quadratically with distance. (Figure 10 shows that out to5 pc the scaling is as expected, una ff ected by the drop-o ff near to10 pc.) This gives encounter rates of 78 . ± . . ± . ±
59 perMyr (which scales to 21 . ± . . ± . | t ph | < . d ph < n enc = ±
44 per Myr within 5 pc. This is 1.5 σ timessmaller than that obtained with the larger time window. This maysuggest some time variability in the encounter rate, although thisis hard to distinguish given the di ffi culty of propagating all theuncertainties. It is also clear that equation 1 is rather sensitive tosmall values of the completeness. Over the window | t ph | < d ph < C i < .
01. It is an unfortunate and unavoidable fact that it isthe low completeness regions which contribute the largest un-certainty to any attempt to correct the encounter rate.
5. Conclusions
We have identified the closest stellar encounters to the Sun fromamong the 7.2 million stars in Gaia DR2 that have 6D phasespace information. Encounters were found by integrating theirorbits in a smooth gravitational potential. The correlated uncer-tainties were accounted for by a Monte Carlo resampling of the6D likelihood distribution of the data and integrating a swarmof surrogate particles. The resulting distributions over perihel-ion time, distance, and speed are generally asymmetric, and aresummarized by their 5th, 50th, and 95th percentiles.We find 31, 8, and 3 stars which have come – or which willcome – within 1 pc, 0.5 pc, and 0.25 pc of the Sun, respectively.These numbers drop to 26, 7, and 3 when we remove likely in-correct results (“bogus encounters”) following a subjective anal-ysis (including visual inspection of images). Quality metrics inthe Gaia catalogue are not calibrated and are hard to use to iden-tify sources with reliable data. Thus some of the stars in our en-counter list are sure to be bogus, and we are sure to have missedsome others for the same reason. In particular, a number of en-counters with unexpected positions in the CMD have large val-ues of the astrometric unit weight error. These are not neces-sarily poor astrometric solutions. At least some could be mainsequence–white dwarf binaries.The closest encounter found is Gl 710, long known to bea close encounter, now found to come slightly closer and withslightly better determined perihelion parameters. Most of the
Article number, page 11 of 12 & A proofs: manuscript no. stellar_encounters_gdr2 other encounters found are discovered here, including 25 within1 pc. Using newly available masses for 98% of the sample com-puted by Fouesneau et al. (in preparation) from Gaia astrometryand multiband photometry, we compute the impulse transfer tothe Oort cloud using the impulse approximation. For both the1 / d ph or 1 / d dependencies, Gl 710 induces the largest impulse.Berski & Dybczy´nski (2016) studied the impact of this star onthe Oort cloud in some detail using the encounter parametersfrom the first Gaia data release. They found that it would injecta large flux of Oort cloud comets toward the inner solar sys-tem. Given that the encounter parameters are not significantlychanged in Gaia DR2, this conclusion still holds. It remains tobe studied what the cumulative e ff ect is of the many more en-counters found in our study.The main factors limiting the accuracy of our resulting peri-helion parameters are: for most stars, the accuracy of the radialvelocities; for distant stars, the accuracy of our potential model;and for some stars, the neglect of possible binarity.Our sample is not complete. The main limitation is theavailability of radial velocities in Gaia DR2. 99% of the stars(whether in the unfiltered or filtered sample) are brighter than G = . T e ff range 3550–6900 K, thereby lim-iting the number of (numerous) late-type stars as well as mas-sive stars. Both T e ff and G should be extended in subsequentGaia data releases. Correcting for this incompleteness using amodel based on the Galaxia simulation, we infer the encounterrate averaged over the past / future 5 Myr to be 491 ±
54 Myr − within 5 pc. When scaled to encounters within 1 pc, this is19 . ± . − . We caution, however, that the accuracy of thisrate is limited by the completeness model assumptions and fit-ting the resulting completeness map, as well as the distributionof the actual encounters. It may also be overestimated if there isa large fraction of sources with spuriously large parallax valuesthat are not accounted for by their formal uncertainties.The other source of incompleteness is missing bright starsin Gaia DR2. Bright stars saturate even with the shortest CCDgate on Gaia, making them harder to calibrate. Astrometry forthese will only be provided in later data releases. Although fewin number (and therefore acceptably neglected by our complete-ness correction), some bright stars are currently nearby and / ormassive, so encounter parameters more accurate than those ob-tained with Hipparcos could ultimately reveal some importantencounters. The case of gamma Microscopii – no longer a closeencounter in Gaia DR2 – has already been discussed. Anothercase is the A2 dwarf zeta Leporis, computed in paper 1 to have d medph = .
30 pc 850 kyr ago. It is in Gaia DR2 but without a radialvelocity. Using the old radial velocity we find its perihelion pa-rameters to be rather similar ( t medph = −
860 kyr, d medph = .
43 pc).Sirius, Altair, and Algol are not in Gaia DR2.There are no doubt many more close – and probably closer –encounters to be discovered in future Gaia data releases.
Acknowledgements.
This work is based on data from the European SpaceAgency (ESA) mission Gaia ( ), pro-cessed by the Gaia Data Processing and Analysis Consortium (DPAC, ). Funding for theDPAC has been provided by national institutions, in particular the institutionsparticipating in the Gaia Multilateral Agreement. We thank Lennart Lindegrenfor discussions about the Gaia astrometric solution and Ted von Hippel for dis-cussions about white dwarfs. This work was funded in part by the DLR (Germanspace agency) via grant 50 QG 1403. We are grateful for the availability of theSimbad object database and Aladin sky atlas, both developed and provided byCDS, Strasbourg.
References
Andrae, R., Fouesneau, M., Creevey, O., et al. 2018, ArXiv e-prints,arXiv:1804.09374Arenou, F., Luri, X., Babusiaux, C., et al. 2018, ArXiv e-prints,arXiv:1804.09375Bailer-Jones, C. A. L. 2015a, A&A, 575, A35Bailer-Jones, C. A. L. 2015b, PASP, 127, 994Bailer-Jones, C. A. L. 2018, A&A, 609Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R.2018, ArXiv e-prints, arXiv:1804.10121Berski, F. & Dybczy´nski, P. A. 2016, A&A, 595, L10Bobylev, V. V. 2010a, Astronomy Letters, 36, 220Bobylev, V. V. 2010b, Astronomy Letters, 36, 816Bobylev, V. V. 2017, Astronomy Reports, 61, 883Bobylev, V. V. & Bajkova, A. T. 2017, Astronomy Letters, 43, 559de la Fuente Marcos, R. & de la Fuente Marcos, C. 2018, ArXiv e-prints,arXiv:1805.02644Dybczynski, P. A. 1994, Celestial Mechanics and Dynamical Astronomy, 58, 139Dybczy´nski, P. A. 2006, A&A, 449, 1233Dybczy´nski, P. A. & Berski, F. 2015, MNRAS, 449, 2459Gaia Collaboration. 2018, ArXiv e-prints, arXiv:1804.09365García-Sánchez, J., Preston, R. A., Jones, D. L., et al. 1999, AJ, 117, 1042García-Sánchez, J., Weissman, P. R., Preston, R. A., et al. 2001, A&A, 379, 634Hollands, M. A., Tremblay, P. E., Gaensicke, B. T., Gentile-Fusillo, N. P., &Toonen, S. 2018, ArXiv e-prints, arXiv:1805.12590Jiménez-Torres, J. J., Pichardo, B., Lake, G., & Throop, H. 2011, MNRAS, 418,1272Katz, D., Sartoretti, P., Cropper, M., et al. 2018, ArXiv e-prints,arXiv:1804.09372Lindegren, L., Hernandez, J., Bombrun, A., et al. 2018, ArXiv e-prints,arXiv:1804.09366Mamajek, E. E., Barenfeld, S. A., Ivanov, V. D., et al. 2015, ApJ, 800, L17Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77Matthews, R. A. J. 1994, QJRAS, 35, 1Mülläri, A. A. & Orlov, V. V. 1996, Earth Moon and Planets, 72, 19Oort, J. H. 1950, Bull. Astron. Inst. Netherlands, 11, 91Rebassa-Mansergas, A., Ren, J. J., Parsons, S. G., et al. 2016, MNRAS, 458,3808Rickman, H. 1976, Bulletin of the Astronomical Institutes of Czechoslovakia,27, 92Robin, A. C., Reylé, C., Derrière, S., & Picaud, S. 2003, A&A, 409, 523Rybizki, J., Demleitner, M., Fouesneau, M., et al. 2018, ArXiv e-prints,arXiv:1804.01427Sharma, S., Bland-Hawthorn, J., Johnston, K. V., & Binney, J. 2011, ApJ, 730, 3van Leeuwen, F., ed. 2007, Astrophysics and Space Science Library, Vol. 350,Hipparcos, the New Reduction of the Raw DataÖpik, E. 1932, Proceedings of the American Academy of Arts and Sciences, 67,169
Appendix A: Gaia archive query
Below is the ADQL query used to select stars from Gaia DR2which have perihelion distances less than 10 pc according to thelinear motion approximation (specified by equation 4 in paper 1).