We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

Recent years have seen increasing interest in the characterization of sub-Neptune sized planets because of their prevalence in the Galaxy, contrasted with their absence in our solar system. HD 97658 is one of the brightest stars hosting a planet of this kind, and we present the transmission spectrum of this planet by combining four HST transits, twelve Spitzer/IRAC transits, and eight MOST transits of this system. Our transmission spectrum has higher signal to noise ratio than that from previous works, and the result suggests that the slight increase in transit depth from wavelength 1.1 to 1.7 microns reported in previous works on the transmission spectrum of this planet is likely systematic. Nonetheless, our atmospheric modeling results are not conclusive as no model provides an excellent match to our data. Nonetheless we find that atmospheres with high C/O ratios (C/O >~ 0.8) and metallicities of >~ 100x solar metallicity are favored. We combine the mid-transit times from all the new Spitzer and MOST observations and obtain an updated orbital period of P=9.489295 +/- 0.000005 d, with a best-fit transit time center at T_0 = 2456361.80690 +/- 0.00038 (BJD). No transit timing variations are found in this system. We also present new measurements of the stellar rotation period (34 +/- 2 d) and stellar activity cycle (9.6 yr) of the host star HD 97658. Finally, we calculate and rank the Transmission Spectroscopy Metric of all confirmed planets cooler than 1000 K and with sizes between 1 and 4 R_Earth. We find that at least a third of small planets cooler than 1000 K can be well characterized using JWST, and of those, HD 97658b is ranked fifth, meaning it remains a high-priority target for atmospheric characterization.

The orbital eccentricities of directly imaged exoplanets and brown dwarf companions provide clues about their formation and dynamical histories. We combine new high-contrast imaging observations of substellar companions obtained primarily with Keck/NIRC2 together with astrometry from the literature to test for differences in the population-level eccentricity distributions of 27 long-period giant planets and brown dwarf companions between 5-100 AU using hierarchical Bayesian modeling. Orbit fits are performed in a uniform manner for companions with short orbital arcs; this typically results in broad constraints for individual eccentricity distributions, but together as an ensemble these systems provide valuable insight into their collective underlying orbital patterns. The shape of the eccentricity distribution function for our full sample of substellar companions is approximately flat from e=0-1. When subdivided by companion mass and mass ratio, the underlying distributions for giant planets and brown dwarfs show significant differences. Low mass ratio companions preferentially have low eccentricities, similar to the orbital properties of warm Jupiters found with radial velocities and transits. We interpret this as evidence for in situ formation on largely undisturbed orbits within massive, extended disks. Brown dwarf companions exhibit a broad peak at e $\approx$ 0.6-0.9 with evidence for a dependence on orbital period. This closely resembles the orbital properties and period-eccentricity trends of wide (1-200 AU) stellar binaries, suggesting that brown dwarfs in this separation range predominantly form in a similar fashion. We also report evidence that the "eccentricity dichotomy" observed at small separations extends to planets on wide orbits: the mean eccentricity for the multi-planet system HR 8799 is lower than for systems with single planets. (Abridged)

Discoveries of exoplanets using the radial velocity method are progressively reaching out to increasingly longer orbital periods as the duration of surveys continues to climb. The improving sensitivity to potential Jupiter analogs is revealing a diversity of orbital architectures that are substantially different from that found in our solar system. An excellent example of this is the recent discovery of HR 5183b; a giant planet on a highly eccentric ($e = 0.84$) $\sim$75~year orbit. The presence of such giant planet orbits are intrinsically interesting from the perspective of the dynamical history of planetary systems, and also for examining the implications of on-going dynamical stability and habitability of these systems. In this work, we examine the latter, providing results of dynamical simulations that explore the stable regions that the eccentric orbit of the HR 5183 giant planet allows to exist within the Habitable Zone of the host star. Our results show that, despite the incredible perturbing influence of the giant planet, there remain a narrow range of locations within the Habitable Zone where terrestrial planets may reside in long-term stable orbits. We discuss the effects of the giant planet on the potential habitability of a stable terrestrial planet, including the modulation of terrestrial planet eccentricities and the periodically spectacular view of the giant planet from the terrestrial planet location.

orbitize! is an open-source, object-oriented software package for fitting the orbits of directly-imaged objects. It packages the Orbits for the Impatient (OFTI) algorithm and a parallel-tempered Markov Chain Monte Carlo (MCMC) algorithm into a consistent and intuitive Python API. orbitize! makes it easy to run standard astrometric orbit fits; in less than 10 lines of code, users can read in data, perform one fit using OFTI and another using MCMC, and make two publication-ready figures. Extensive pedagogical tutorials, intended to be navigable by both orbit-fitting novices and seasoned experts, are available on our documentation page. We have designed the orbitize! API to be flexible and easy to use/modify for unique cases. orbitize! was designed by members of the exoplanet imaging community to be a central repository for algorithms, techniques, and know-how developed by this community. We intend for it to continue to expand and change as the field progresses and new techniques are developed, and call for community involvement in this process. Complete and up-to-date documentation is available at orbitize.info.

Based on two decades of radial velocity (RV) observations using Keck/HIRES and McDonald/Tull, and more recent observations using the Automated Planet Finder, we found that the nearby star HR 5183 (HD 120066) hosts a 3$M_J$ minimum mass planet with an orbital period of $74^{+43}_{-22}$ years. The orbit is highly eccentric (e$\simeq$0.84), shuttling the planet from within the orbit of Jupiter to beyond the orbit of Neptune. Our careful survey design enabled high cadence observations before, during, and after the planet's periastron passage, yielding precise orbital parameter constraints. We searched for stellar or planetary companions that could have excited the planet's eccentricity, but found no candidates, potentially implying that the perturber was ejected from the system. We did identify a bound stellar companion more than 15,000 au from the primary, but reasoned that it is currently too widely separated to have an appreciable effect on HR 5183 b. Because HR 5183 b's wide orbit takes it more than 30 au (1") from its star, we also explored the potential of complimentary studies with direct imaging or stellar astrometry. We found that a Gaia detection is very likely, and that imaging at 10 $\mu$m is a promising avenue. This discovery highlights the value of long-baseline RV surveys for discovering and characterizing long-period, eccentric Jovian planets. This population may offer important insights into the dynamical evolution of planetary systems containing multiple massive planets.

Short-period (P<1 hour) white dwarf binaries will be the most numerous sources for the space-based gravitational wave detector LISA. Based on thousands of resolved systems, we will be able to constrain binary evolution and provide a new map of the Milky Way and its surroundings. Here we predict the main properties of populations of different types of detached white dwarf binaries detected by LISA. For the first time, we combine a high-resolution cosmological simulation of a Milky Way-mass galaxy (from the FIRE project) with a binary population synthesis model for low and intermediate mass stars. Our model therefore provides a cosmologically realistic star formation and metallicity history for the galaxy and naturally produces its different components such as the thin and thick disk, the bulge, the stellar halo, and satellite galaxies and streams. With the simulation, we show how different galactic components contribute differently to the gravitational wave signal, due to their typical age and distance distributions. We find that the dominant LISA sources will be He-He systems and He-CO systems with important contributions from the thick disk and bulge but also a few systems in the stellar halo. The resulting sky map of the sources is different from previous models, with important consequences for the searches for electromagnetic counterparts and data analysis. We also emphasize that much of the science-enabling information regarding white dwarf binaries, such as the chirp mass and the sky localisation, becomes increasingly rich with long observations, including an extended mission up to 8 years.

We report on the discovery of a white dwarf companion to the nearby late G dwarf star, HD 159062. The companion is detected in 14 years of precise radial velocity (RV) data, and in high-resolution imaging observations. RVs of HD 159062 from 2003-2018 reveal an acceleration of $-13.3\pm0.12\ \rm{m s}^{-1}$, indicating that it hosts a companion with a long-period orbit. Subsequent imaging observations with the ShaneAO system on the Lick Observatory 3-meter Shane telescope, the PHARO AO system on the Palomar Observatory 5-meter telescope, and the NIRC2 AO system at the Keck II 10-meter telescope reveal a faint companion 2.7'' from the primary star. We performed relative photometry, finding $\Delta J = 10.09 \pm 0.38$ magnitudes, $\Delta K_{S} = 10.06 \pm 0.22$ magnitudes, and $\Delta L' = 9.67\pm0.08$ magnitudes for the companion from these observations. Analysis of the radial velocities, astrometry, and photometry reveals that the combined data set can only be reconciled for the scenario where HD 159062 B is a white dwarf. A full Bayesian analysis of the RV and imaging data to obtain the cooling age, mass, and orbital parameters of the white dwarf indicates that the companion is an old $M_{B} = 0.65^{+0.12}_{-0.04} M_{\odot}$ white dwarf with an orbital period of $P = 250^{+130}_{-76}$ years, and a cooling age of $\tau = 8.2^{+0.3}_{-0.5}$ Gyr.

K2-291 (EPIC 247418783) is a solar-type star with a radius of R_star = 0.899 $\pm$ 0.034 R_sun and mass of M_star=0.934 $\pm$ 0.038 M_sun. From K2 C13 data, we found one super-Earth planet (R_p = 1.589+0.095-0.072 R_Earth) transiting this star on a short period orbit (P = 2.225177 +6.6e-5 -6.8e-5 days). We followed this system up with adaptive-optic imaging and spectroscopy to derive stellar parameters, search for stellar companions, and determine a planet mass. From our 75 radial velocity measurements using HIRES on Keck I and HARPS-N on Telescopio Nazionale Galileo, we constrained the mass of EPIC 247418783b to M_p = 6.49 $\pm$ 1.16 M_Earth. We found it necessary to model correlated stellar activity radial velocity signals with a Gaussian process in order to more accurately model the effect of stellar noise on our data; the addition of the Gaussian process also improved the precision of this mass measurement. With a bulk density of 8.84+2.50-2.03 g cm-3, the planet is consistent with an Earth-like rock/iron composition and no substantial gaseous envelope. Such an envelope, if it existed in the past, was likely eroded away by photo-evaporation during the first billion years of the star's lifetime.

We report improved masses, radii, and densities for four planets in two bright M-dwarf systems, K2-3 and GJ3470, derived from a combination of new radial velocity and transit observations. Supplementing K2 photometry with follow-up Spitzer transit observations refined the transit ephemerides of K2-3 b, c, and d by over a factor of 10. We analyze ground-based photometry from the Evryscope and Fairborn Observatory to determine the characteristic stellar activity timescales for our Gaussian Process fit, including the stellar rotation period and activity region decay timescale. The stellar rotation signals for both stars are evident in the radial velocity data and are included in our fit using a Gaussian process trained on the photometry. We find the masses of K2-3 b, K2-3 c and GJ3470 b to be 6.48$^{+0.99}_{-0.93}$, 2.14$^{+1.08}_{-1.04}$, and 12.58$^{+1.31}_{-1.28}$ M$_\oplus$ respectively. K2-3 d was not significantly detected and has a 3-$\sigma$ upper limit of 2.80 M$_\oplus$. These two systems are training cases for future TESS systems; due to the low planet densities ($\rho$ $<$ 3.7 g cm$^{-3}$) and bright host stars (K $<$ 9 mag), they are among the best candidates for transmission spectroscopy in order to characterize the atmospheric compositions of small planets.

We present SCExAO/CHARIS high-contrast imaging/$JHK$ integral field spectroscopy of $\kappa$ And b, a directly-imaged low-mass companion orbiting a nearby B9V star. We detect $\kappa$ And b at a high signal-to-noise and extract high precision spectrophotometry using a new forward-modeling algorithm for (A-)LOCI complementary to KLIP-FM developed by Pueyo (2016). $\kappa$ And b's spectrum best resembles that of a low-gravity L0--L1 dwarf (L0--L1$\gamma$). Its spectrum and luminosity are very well matched by 2MASSJ0141-4633 and several other 12.5--15 $M_{\rm J}$ free floating members of the 40 $Myr$-old Tuc-Hor Association, consistent with a system age derived from recent interferometric results for the primary, a companion mass at/near the deuterium-burning limit (13$^{+12}_{-2}$ M$_{\rm J}$), and a companion-to-primary mass ratio characteristic of other directly-imaged planets ($q$ $\sim$ 0.005$^{+0.005}_{-0.001}$). We did not unambiguously identify additional, more closely-orbiting companions brighter and more massive than $\kappa$ And b down to $\rho$ $\sim$ 0.3" (15 au). SCExAO/CHARIS and complementary Keck/NIRC2 astrometric points reveal clockwise orbital motion. Modeling points towards a likely eccentric orbit: a subset of acceptable orbits include those that are aligned with the star's rotation axis. However, $\kappa$ And b's semimajor axis is plausibly larger than 75 au and in a region where disk instability could form massive companions. Deeper $\kappa$ And high-contrast imaging and low-resolution spectroscopy from extreme AO systems like SCExAO/CHARIS and higher resolution spectroscopy from Keck/OSIRIS or, later, IRIS on the Thirty Meter Telescope could help clarify $\kappa$ And b's chemistry and whether its spectrum provides an insight into its formation environment.

RadVel is an open source Python package for modeling Keplerian orbits in radial velocity (RV) time series. RadVel provides a convenient framework to fit RVs using maximum a posteriori optimization and to compute robust confidence intervals by sampling the posterior probability density via Markov Chain Monte Carlo (MCMC). RadVel allows users to float or fix parameters, impose priors, and perform Bayesian model comparison. We have implemented realtime MCMC convergence tests to ensure adequate sampling of the posterior. RadVel can output a number of publication-quality plots and tables. Users may interface with RadVel through a convenient command-line interface or directly from Python. The code is object-oriented and thus naturally extensible. We encourage contributions from the community. Documentation is available at http://radvel.readthedocs.io.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

Our Keck/NIRC2 imaging survey searches for stellar companions around 144 systems with radial velocity (RV) detected giant planets to determine whether stellar binaries influence the planets' orbital parameters. This survey, the largest of its kind to date, finds eight confirmed binary systems and three confirmed triple systems. These include three new multi-stellar systems (HD 30856, HD 86081, and HD 207832) and three multi-stellar systems with newly confirmed common proper motion (HD 43691, HD 116029, and HD 164509). We combine these systems with seven RV planet-hosting multi-stellar systems from the literature in order to test for differences in the properties of planets with semimajor axes ranging between 0.1-5 au in single vs multi-stellar systems. We find no evidence that the presence or absence of stellar companions alters the distribution of planet properties in these systems. Although the observed stellar companions might influence the orbits of more distant planetary companions in these systems, our RV observations currently provide only weak constraints on the masses and orbital properties of planets beyond 5 au. In order to aid future efforts to characterize long period RV companions in these systems, we publish our contrast curves for all 144 targets. Using four years of astrometry for six hierarchical triple star systems hosting giant planets, we fit the orbits of the stellar companions in order to characterize the orbital architecture in these systems. We find that the orbital plane of the secondary and tertiary companions are inconsistent with an edge-on orbit in four out of six cases.

We describe a Bayesian rejection sampling algorithm designed to efficiently compute posterior distributions of orbital elements for data covering short fractions of long-period exoplanet orbits. Our implementation of this method, Orbits for the Impatient (OFTI), converges up to several orders of magnitude faster than two implementations of MCMC in this regime. We illustrate the efficiency of our approach by showing that OFTI calculates accurate posteriors for all existing astrometry of the exoplanet 51 Eri b up to 100 times faster than a Metropolis-Hastings MCMC. We demonstrate the accuracy of OFTI by comparing our results for several orbiting systems with those of various MCMC implementations, finding the output posteriors to be identical within shot noise. We also describe how our algorithm was used to successfully predict the location of 51 Eri b six months in the future based on less than three months of astrometry. Finally, we apply OFTI to ten long-period exoplanets and brown dwarfs, all but one of which have been monitored over less than 3% of their orbits, producing fits to their orbits from astrometric records in the literature.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

We present a new nearby young moving group (NYMG) kinematic membership analysis code, LocAting Constituent mEmbers In Nearby Groups (LACEwING), a new Catalog of Suspected Nearby Young Stars, a new list of bona fide members of moving groups, and a kinematic traceback code. LACEwING is a convergence-style algorithm with carefully vetted membership statistics based on a large numerical simulation of the Solar Neighborhood. Given spatial and kinematic information on stars, LACEwING calculates membership probabilities in 13 NYMGs and three open clusters within 100 pc. In addition to describing the inputs, methods, and products of the code, we provide comparisons of LACEwING to other popular kinematic moving group membership identification codes. As a proof of concept, we use LACEwING to reconsider the membership of 930 stellar systems in the Solar Neighborhood (within 100 pc) that have reported measurable lithium equivalent widths. We quantify the evidence in support of a population of young stars not attached to any NYMGs, which is a possible sign of new as-yet-undiscovered groups or of a field population of young stars.

We have conducted an angular differential imaging survey with NIRC2 at Keck in search of close-in substellar companions to a sample of seven systems with confirmed planetary-mass companions (PMCs) on wide orbits (>50 AU). These wide-separation PMCs pose significant challenges to all three possible formation mechanisms: core accretion plus scattering, disk instability, and turbulent fragmentation. We explore the possibility that these companions formed closer in and were scattered out to their present-day locations by searching for other massive bodies at smaller separations. The typical sensitivity for this survey is \Delta K ~ 12.5 at 1". We identify eight candidate companions, whose masses would reach as low as one Jupiter mass if gravitationally bound. From our multi-epoch astrometry we determine that seven of these are conclusively background objects, while the eighth near DH Tau is ambiguous and requires additional monitoring. We rule out the presence of >7 M$_{\rm Jup}$ bodies in these systems down to 15 - 50 AU that could be responsible for scattering. This result combined with the totality of evidence suggests that dynamical scattering is unlikely to have produced this population of PMCs. We detect orbital motion from the companions ROXs 42B b and ROXs 12 b, and from this determine 95% upper limits on the companions' eccentricities of 0.58 and 0.83 respectively. Finally, we find that the 95% upper limit on the occurrence rate of additional planets with masses between 5 - 15 M$_{\rm Jup}$ outside of 40 AU in systems with PMCs is 54%.

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

We present basic observational data and association membership analysis for 45 young and active low-mass stellar systems from the ongoing RECONS photometry and astrometry program at the Cerro Tololo Inter-American Observatory. Most of these systems have saturated X-ray emission (log(Lx/Lbol) > -3.5) based on X-ray fluxes from the ROSAT All-Sky Survey, and many are significantly more luminous than main-sequence stars of comparable color. We present parallaxes and proper motions, Johnson-Kron-Cousins VRI photometry, and multiplicity observations from the CTIOPI program on the CTIO 0.9m telescope. To this we add low-resolution optical spectroscopy and line measurements from the CTIO 1.5m telescope, and interferometric binary measurements from the Hubble Space Telescope Fine Guidance Sensors. We also incorporate data from published sources: JHKs photometry from the 2MASS point source catalog; X-ray data from the ROSAT All-Sky Survey; and radial velocities from literature sources. Within the sample of 45 systems, we identify 21 candidate low-mass pre-main-sequence members of nearby associations, including members of beta Pictoris, TW Hydrae, Argus, AB Doradus, two ambiguous 30 Myr old systems, and one object that may be a member of the Ursa Major moving group. Of the 21 candidate young systems, 14 are newly identified as a result of this work, and six of those are within 25 parsecs of the Sun.