The most abundant stars in the Galaxy, M dwarfs, are very commonly hosts to diverse systems of low-mass planets. Their abundancy implies that the general occurrence rate of planets is dominated by their occurrence rate around such M dwarfs. In this article, we combine the M dwarf surveys conducted with the HIRES/Keck, PFS/Magellan, HARPS/ESO, and UVES/VLT instruments supported with data from several other instruments. We analyse the radial velocities of an approximately volume- and brightness-limited sample of 426 nearby M dwarfs in order to search for Doppler signals of cadidate planets. In addition, we analyse spectroscopic activity indicators and ASAS photometry to rule out radial velocity signals corresponding to stellar activity as Doppler signals of planets. We calculate estimates for the occurrence rate of planets around the sample stars and study the properties of this occurrence rate as a function of stellar properties. Our analyses reveal a total of 118 candidate planets orbiting nearby M dwarfs. Based on our results accounting for selection effects and sample detection threshold, we estimate that M dwarfs have on average at least 2.39$^{+4.58}_{-1.36}$ planets per star orbiting them. Accounting for the different sensitivities of radial velocity surveys and Kepler transit photometry implies that there are at least 3.0 planets per star orbiting M dwarfs. We also present evidence for a population of cool mini-Neptunes and Neptunes with indications that they are found an order of magnitude more frequently orbiting the least massive M dwarfs in our sample.

We present a new precision radial velocity dataset that reveals a multi-planet system orbiting the G0V star HD 34445. Our 18-year span consists of 333 precision radial velocity observations, 56 of which were previously published, and 277 which are new data from Keck Observatory, Magellan at Las Campanas Observatory, and the Automated Planet Finder at Lick Observatory. These data indicate the presence of six planet candidates in Keplerian motion about the host star with periods of 1057, 215, 118, 49, 677, and 5700 days, and minimum masses of 0.63, 0.17, 0.1, 0.05, 0.12 and 0.38 Jupiter masses respectively. The HD 34445 planetary system, with its high degree of multiplicity, its long orbital periods, and its induced stellar radial velocity half-amplitudes in the range $2 \,{\rm m\, s^{-1}} \lesssim K \lesssim 5\,{\rm m\, s^{-1}}$ is fundamentally unlike either our own solar system (in which only Jupiter and Saturn induce significant reflex velocities for the Sun), or the Kepler multiple-transiting systems (which tend to have much more compact orbital configurations)

Upper Scorpius is a subgroup of the nearest OB association, Scorpius--Centaurus. Its young age makes it an important association to study star and planet formation. We present parallaxes to 52 low mass stars in Upper Scorpius, 28 of which have full kinematics. We measure ages of the individual stars by combining our measured parallaxes with pre-main sequence evolutionary tracks. We find there is a significant difference in the ages of stars with and without circumstellar disks. The stars without disks have a mean age of 4.9+/-0.8 Myr and those with disks have an older mean age of 8.2+/-0.9 Myr. This somewhat counterintuitive result suggests that evolutionary effects in young stars can dominate their apparent ages. We also attempt to use the 28 stars with full kinematics (i.e.\ proper motion, radial velocity, and parallax) to trace the stars back in time to their original birthplace to obtain a trackback age. We find, as expected given large measurement uncertainties on available radial velocity measurements, that measurement uncertainties alone cause the group to diverge after a few Myr.

Transit photometry of the M8V dwarf star TRAPPIST-1 (2MASS J23062928-0502285) has revealed the presence of at least seven planets with masses and radii similar to that of Earth orbiting at distances that might allow liquid water to be present on their surfaces. We have been following TRAPPIST-1 since 2011 with the CAPSCam astrometric camera on the 2.5-m du Pont telescope at the Las Campanas Observatory in Chile. In 2016 we noted that TRAPPIST-1 lies slightly farther away than previously thought, at 12.49 pc, rather than 12.1 pc. Here we examine fifteen epochs of CAPSCam observations of TRAPPIST-1, spanning the five years from 2011 to 2016, and obtain a revised trigonometric distance of $12.56 \pm 0.12$ pc. The astrometric data analysis pipeline shows no evidence for a long-period astrometric wobble of TRAPPIST-1. After proper motion and parallax are removed, residuals at the level of $\pm 1.3$ millarcsec (mas) remain. The amplitude of these residuals constrains the masses of any long-period gas giant planets in the TRAPPIST-1 system: no planet more massive than $\sim 4.6 M_{Jup}$ orbits with a 1 yr period, and no planet more massive than $\sim 1.6 M_{Jup}$ orbits with a 5 yr period. Further refinement of the CAPSCam data analysis pipeline, combined with continued CAPSCam observations, should either detect any long-period planets, or put an even tighter constraint on these mass upper limits.

We describe a 20-year survey carried out by the Lick-Carnegie Exoplanet Survey Team (LCES), using precision radial velocities from HIRES on the Keck-I telescope to find and characterize extrasolar planetary systems orbiting nearby F, G, K, and M dwarf stars. We provide here 60,949 precision radial velocities for 1,624 stars contained in that survey. We tabulate a list of 357 significant periodic signals that are of constant period and phase, and not coincident in period and/or phase with stellar activity indices. These signals are thus strongly suggestive of barycentric reflex motion of the star induced by one or more candidate exoplanets in Keplerian motion about the host star. Of these signals, 225 have already been published as planet claims, 60 are classified as significant unpublished planet candidates that await photometric follow-up to rule out activity-related causes, and 54 are also unpublished, but are classified as "significant" signals that require confirmation by additional data before rising to classification as planet candidates. Of particular interest is our detection of a candidate planet with a minimum mass of 3.9 Earth masses and an orbital period of 9.9 days orbiting Lalande 21185, the fourth-closest main sequence star to the Sun. For each of our exoplanetary candidate signals, we provide the period and semi-amplitude of the Keplerian orbital fit, and a likelihood ratio estimate of its statistical significance. We also tabulate 18 Keplerian-like signals that we classify as likely arising from stellar activity.

The TW Hya Association (TWA) is a nearby stellar association with an age of $\sim$ 5-10 Myr. This is an important age for studying the late stages of star and planet formation. We measure the parallaxes of 14 candidate members of TWA. That brings to 38 the total number of individual stars with fully measured kinematics, i.e. proper motion, radial velocity, and parallax, to describe their motions through the Galaxy. We analyze these kinematics to search for convergence to a smaller volume in the past, but we find the association is never much more compact than it is at present. We show that it is difficult to measure traceback ages for associations such as TWA that have expected velocity dispersions of 1-2 km s$^{-1}$ with typical measurement uncertainties. We also use our stellar distances and pre-main-sequence evolutionary tracks to find the average age of the association of 7.9 $\pm$ 1.0 Myr. Additionally, our parallax measurement of TWA 32 indicates it should be considered a bona fide member of TWA. Two new candidate members have high membership probabilities, and we assign them TWA numbers: TWA 45 for 2MASS J11592786-4510192 and TWA 46 for 2MASS J12354615-4115531.

We report trigonometric parallaxes for 134 low mass stars and brown dwarfs, of which 38 have no previously published measurement and 79 more have improved uncertainties. Our survey targeted nearby targets, so 119 are closer than 30 pc. Of the 38 stars with new parallaxes, 14 are within 20 pc and seven are likely brown dwarfs (spectral types later than L0). These parallaxes are useful for studies of kinematics, multiplicity, and spectrophotometric calibration. Two objects with new parallaxes are confirmed as young stars with membership in nearby young moving groups: LP 870-65 in AB Doradus and G 161-71 in Argus. We also report the first parallax for the planet-hosting star GJ 3470; this allows us to refine the density of its Neptune-mass planet. One T-dwarf, 2MASS J12590470-4336243, previously thought to lie within 4 pc, is found to be at 7.8 pc, and the M-type star 2MASS J01392170-3936088 joins the ranks of nearby stars as it is found to be within 10 pc. Five stars that are over-luminous and/or too red for their spectral types are identified and deserve further study as possible young stars.

We present a set of 109 new, high-precision Keck/HIRES radial velocity (RV) observations for the solar-type star HD 32963. Our dataset reveals a candidate planetary signal with a period of 6.49 $\pm$ 0.07 years and a corresponding minimum mass of 0.7 $\pm$ 0.03 Jupiter masses. Given Jupiter's crucial role in shaping the evolution of the early Solar System, we emphasize the importance of long-term radial velocity surveys. Finally, using our complete set of Keck radial velocities and correcting for the relative detectability of synthetic planetary candidates orbiting each of the 1,122 stars in our sample, we estimate the frequency of Jupiter analogs across our survey at approximately 3%.

We report initial performance results emerging from 600 hours of observations with the Automated Planet Finder (APF) telescope and Levy Spectrometer located at UCO/Lick Observatory. We have obtained multiple spectra of 80 G, K and M-type stars, which comprise 4,954 individual Doppler radial velocity (RV) measurements with a median internal uncertainty of 1.35 ms$^{-1}$. We find a strong, expected correlation between the number of photons accumulated in the 5000-6200$\AA$ iodine region of the spectrum, and the resulting internal uncertainty estimates. Additionally, we find an offset between the population of G and K stars and the M stars within the data set when comparing these parameters. As a consequence of their increased spectral line densities, M-type stars permit the same level of internal uncertainty with 2x fewer photons than G-type and K-type stars. When observing M stars, we show that the APF/Levy has essentially the same speed-on-sky as Keck/HIRES for precision RVs. In the interest of using the APF for long-duration RV surveys, we have designed and implemented a dynamic scheduling algorithm. We discuss the operation of the scheduler, which monitors ambient conditions and combines on-sky information with a database of survey targets to make intelligent, real-time targeting decisions.

We present new, high-precision Doppler radial velocity (RV) data sets for the nearby K3V star HD 219134. The data include 175 velocities obtained with the HIRES Spectrograph at the Keck I Telescope, and 101 velocities obtained with the Levy Spectrograph at the Automated Planet Finder Telescope (APF) at Lick Observatory. Our observations reveal six new planetary candidates, with orbital periods of P=3.1, 6.8, 22.8, 46.7, 94.2 and 2247 days, spanning masses of msini=3.8, 3.5, 8.9, 21.3, 10.8 and 108 M_earth respectively. Our analysis indicates that the outermost signal is unlikely to be an artifact induced by stellar activity. In addition, several years of precision photometry with the T10 0.8~m automatic photometric telescope (APT) at Fairborn Observatory demonstrated a lack of brightness variability to a limit of ~0.0002 mag, providing strong support for planetary-reflex motion as the source of the radial velocity variations. The HD 219134 system, with its bright (V=5.6) primary provides an excellent opportunity to obtain detailed orbital characterization (and potentially follow-up observations) of a planetary system that resembles many of the multiple-planet systems detected by Kepler, and which are expected to be detected by NASA's forthcoming TESS Mission and by ESA's forthcoming PLATO Mission.

Both astronomical observations of the interaction of Type II supernova remnants (SNR) with dense interstellar clouds as well as cosmochemical studies of the abundances of daughter products of short-lived radioisotopes (SLRIs) formed by supernova nucleosynthesis support the hypothesis that the Solar Systems SLRIs may have been derived from a supernova. This paper continues a series devoted to examining whether such a shock wave could have triggered the dynamical collapse of a dense, presolar cloud core and simultaneously injected sufficient abundances of SLRIs to explain the cosmochemical evidence. Here we examine the effects of shock waves striking clouds whose spin axes are oriented perpendicular, rather than parallel, to the direction of propagation of the shock front. The models start with 2.2 solar mass cloud cores and shock speeds of 20 or 40 km/sec. Central protostars and protoplanetary disks form in all models, though with disk spin axes aligned somewhat randomly. The disks derive most of their angular momentum not from the initial cloud rotation, but from the Rayleigh-Taylor fingers that also inject shock wave SLRIs. Injection efficiencies, fi, the fraction of the incident shock wave material injected into the collapsing cloud core, are 0.04 - 0.1 in these models, similar to when the rotation axis is parallel to the shock propagation direction. Evidently altering the rotation axis orientation has only a minor effect on the outcome, strengthening the case for this scenario as an explanation for the Solar Systems SLRIs.

We present the results of a large suite of three-dimensional (3D) models of the collapse of magnetic molecular cloud cores using the adaptive mesh refinement (AMR) code Enzo2.2 in the ideal magnetohydrodynamics (MHD) approximation. The cloud cores are initially either prolate or oblate, centrally condensed clouds with masses of 1.73 or 2.73 $M_\odot$, respectively. The radial density profiles are Gaussian, with central densities 20 times higher than boundary densities. A barotropic equation of state is used to represent the transition from low density, isothermal phases, to high density, optically thick phases. The initial magnetic field strength ranges from 6.3 to 100 $\mu$G, corresponding to clouds that are strongly to marginally supercritical, respectively, in terms of the mass to magnetic flux ratio. The magnetic field is initially uniform and aligned with the clouds' rotation axes, with initial ratios of rotational to gravitational energy ranging from $10^{-4}$ to 0.1. Two significantly different outcomes for collapse result: (1) formation of single protostars with spiral arms, and (2) fragmentation into multiple protostar systems. The transition between these two outcomes depends primarily on the initial magnetic field strength, with fragmentation occurring for mass to flux ratios greater than about 14 times the critical ratio for prolate clouds. Oblate clouds typically fragment into several times more clumps than prolate clouds. Multiple, rather than binary, system formation is the general rule in either case, suggesting that binary stars are primarily the result of the orbital dissolution of multiple protostar systems.

A key test of the supernova triggering and injection hypothesis for the origin of the solar system's short-lived radioisotopes is to reproduce the inferred initial abundances of these isotopes. We present here the most detailed models to date of the shock wave triggering and injection process, where shock waves with varied properties strike fully three dimensional, rotating, dense cloud cores. The models are calculated with the FLASH adaptive mesh hydrodynamics code. Three different outcomes can result: triggered collapse leading to fragmentation into a multiple protostar system; triggered collapse leading to a single protostar embedded in a protostellar disk; or failure to undergo dynamic collapse. Shock wave material is injected into the collapsing clouds through Rayleigh-Taylor fingers, resulting in initially inhomogeneous distributions in the protostars and protostellar disks. Cloud rotation about an axis aligned with the shock propagation direction does not increase the injection efficiency appreciably, as the shock parameters were chosen to be optimal for injection even in the absence of rotation. For a shock wave from a core-collapse supernova, the dilution factors for supernova material are in the range of $\sim 10^{-4}$ to $\sim 3 \times 10^{-4}$, in agreement with recent laboratory estimates of the required amount of dilution for $^{60}$Fe and $^{26}$Al. We conclude that a type II supernova remains as a promising candidate for synthesizing the solar system's short-lived radioisotopes shortly before their injection into the presolar cloud core by the supernova's remnant shock wave.

The Automated Planet Finder (APF) is a facility purpose-built for the discovery and characterization of extrasolar planets through high-cadence Doppler velocimetry of the reflex barycentric accelerations of their host stars. Located atop Mt. Hamilton, the APF facility consists of a 2.4-m telescope and its Levy spectrometer, an optical echelle spectrometer optimized for precision Doppler velocimetry. APF features a fixed format spectral range from 374 nm - 970 nm, and delivers a "Throughput" (resolution * slit width product) of 114,000 arc-seconds, with spectral resolutions up to 150,000. Overall system efficiency (fraction of photons incident on the primary mirror that are detected by the science CCD) on blaze at 560 nm in planet-hunting mode is 15%. First-light tests on the RV standard stars HD 185144 and HD 9407 demonstrate sub-meter per second precision (RMS per observation) held over a 3-month period. This paper reviews the basic features of the telescope, dome, and spectrometer, and gives a brief summary of first-light performance.

A variety of stellar sources have been proposed for the origin of the short-lived radioisotopes that existed at the time of the formation of the earliest Solar System solids, including Type II supernovae, AGB and super-AGB stars, and Wolf-Rayet star winds. Our previous adaptive mesh hydrodynamics models with the FLASH2.5 code have shown which combinations of shock wave parameters are able to simultaneously trigger the gravitational collapse of a target dense cloud core and inject significant amounts of shock wave gas and dust, showing that thin supernova shocks may be uniquely suited for the task. However, recent meteoritical studies have weakened the case for a direct supernova injection to the presolar cloud, motivating us to re-examine a wider range of shock wave and cloud core parameters, including rotation, in order to better estimate the injection efficiencies for a variety of stellar sources. We find that supernova shocks remain as the most promising stellar source, though planetary nebulae resulting from AGB star evolution cannot be conclusively ruled out. Wolf-Rayet star winds, however, are likely to lead to cloud core shredding, rather than to collapse. Injection efficiencies can be increased when the cloud is rotating about an axis aligned with the direction of the shock wave, by as much as a factor of $\sim 10$. The amount of gas and dust accreted from the post-shock wind can exceed that injected from the shock wave, with implications for the isotopic abundances expected for a supernova source.

Magnetic fields are important contributers to the dynamics of collapsing molecular cloud cores, and can have a major effect on whether collapse results in a single protostar or fragmentation into a binary or multiple protostar system. New models are presented of the collapse of magnetic cloud cores using the adaptive mesh refinement (AMR) code Enzo2.0. The code was used to calculate the ideal magnetohydrodynamics (MHD) of initially spherical, uniform density and rotation clouds with density perturbations, i.e., the Boss and Bodenheimer (1979) standard isothermal test case for three dimensional (3D) hydrodynamics (HD) codes. After first verifying that Enzo reproduces the binary fragmentation expected for the non-magnetic test case, a large set of models was computed with varied initial magnetic field strengths and directions with respect to the cloud core axis of rotation (parallel or perpendicular), density perturbation amplitudes, and equations of state. Three significantly different outcomes resulted: (1) contraction without sustained collapse, forming a denser cloud core, (2) collapse to form a single protostar with significant spiral arms, and (3) collapse and fragmentation into binary or multiple protostar systems, with multiple spiral arms. Comparisons are also made with previous MHD calculations of similar clouds with barotropic equations of state. These results for the collapse of initially uniform density spheres illustrate the central importance of both magnetic field direction and field strength for determining the outcome of dynamic protostellar collapse.

A supernova is a likely source of short-lived radioisotopes (SLRIs) that were present during the formation of the earliest solar system solids. A suitably thin and dense supernova shock wave may be capable of triggering the self-gravitational collapse of a molecular cloud core while simultaneously injecting SLRIs. Axisymmetric hydrodynamics models have shown that this injection occurs through a number of Rayleigh-Taylor (RT) rings. Here we use the FLASH adaptive mesh refinement (AMR) hydrodynamics code to calculate the first fully three dimensional (3D) models of the triggering and injection process. The axisymmetric RT rings become RT fingers in 3D. While ~ 100 RT fingers appear early in the 3D models, only a few RT fingers are likely to impact the densest portion of the collapsing cloud core. These few RT fingers must then be the source of any SLRI spatial heterogeneity in the solar nebula inferred from isotopic analyses of chondritic meteorites. The models show that SLRI injection efficiencies from a supernova several pc away fall at the lower end of the range estimated for matching SLRI abundances, perhaps putting them more into agreement with recent reassessments of the level of 60Fe present in the solar nebula.

The short-lived radioisotope $^{60}$Fe requires production in a core collapse supernova or AGB star immediately before its incorporation into the earliest solar system solids. Shock waves from a somewhat distant supernova, or a relatively nearby AGB star, have the right speeds to simultaneously trigger the collapse of a dense molecular cloud core and to inject shock wave material into the resulting protostar. A new set of FLASH2.5 adaptive mesh refinement hydrodynamical models shows that the injection efficiency depends sensitively on the assumed shock thickness and density. Supernova shock waves appear to be thin enough to inject the amount of shock wave material necessary to match the short-lived radioisotope abundances measured for primitive meteorites. Planetary nebula shock waves from AGB stars, however, appear to be too thick to achieve the required injection efficiencies. These models imply that a supernova pulled the trigger that led to the formation of our solar system.

The discovery of decay products of a short-lived radioisotope (SLRI) in the Allende meteorite led to the hypothesis that a supernova shock wave transported freshly synthesized SLRI to the presolar dense cloud core, triggered its self-gravitational collapse, and injected the SLRI into the core. Previous multidimensional numerical calculations of the shock-cloud collision process showed that this hypothesis is plausible when the shock wave and dense cloud core are assumed to remain isothermal at ~10 K, but not when compressional heating to ~1000 K is assumed. Our two-dimensional models (Boss et al. 2008) with the FLASH2.5 adaptive mesh refinement (AMR) hydrodynamics code have shown that a 20 km/sec shock front can simultaneously trigger collapse of a 1 solar mass core and inject shock wave material, provided that cooling by molecular species such as H2O, CO, and H2 is included. Here we present the results for similar calculations with shock speeds ranging from 1 km/sec to 100 km/sec. We find that shock speeds in the range from 5 km/sec to 70 km/sec are able to trigger the collapse of a 2.2 solar mass cloud while simultaneously injecting shock wave material: lower speed shocks do not achieve injection, while higher speed shocks do not trigger sustained collapse. The calculations continue to support the shock-wave trigger hypothesis for the formation of the solar system, though the injection efficiencies in the present models are lower than desired.

Cosmochemical evidence for the existence of short-lived radioisotopes (SLRI) such as $^{26}$Al and $^{60}$Fe at the time of the formation of primitive meteorites requires that these isotopes were synthesized in a massive star and then incorporated into chondrites within $\sim 10^6$ yr. A supernova shock wave has long been hypothesized to have transported the SLRI to the presolar dense cloud core, triggered cloud collapse, and injected the isotopes. Previous numerical calculations have shown that this scenario is plausible when the shock wave and dense cloud core are assumed to be isothermal at $\sim 10$ K, but not when compressional heating to $\sim 1000$ K is assumed. We show here for the first time that when calculated with the FLASH2.5 adaptive mesh refinement (AMR) hydrodynamics code, a 20 km/sec shock wave can indeed trigger the collapse of a 1 $M_\odot$ cloud while simultaneously injecting shock wave isotopes into the collapsing cloud, provided that cooling by molecular species such as H$_2$O, CO$_2$, and H$_2$ is included. These calculations imply that the supernova trigger hypothesis is the most likely mechanism for delivering the SLRI present during the formation of the solar system.