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.

We present a time series of the highest resolution spectra yet published for the nearby Type Ia supernova (SN) 2014J in M82. They were obtained at 11 epochs over 33 days around peak brightness with the Levy Spectrograph (resolution R~110,000) on the 2.4m Automated Planet Finder telescope at Lick Observatory. We identify multiple Na I D and K I absorption features, as well as absorption by Ca I H & K and several of the more common diffuse interstellar bands (DIBs). We see no evolution in any component of Na I D, Ca I, or in the DIBs, but do establish the dissipation/weakening of the two most blueshifted components of K I. We present several potential physical explanations, finding the most plausible to be photoionization of circumstellar material, and discuss the implications of our results with respect to the progenitor scenario of SN 2014J.

Exoplanets of a few Earth masses can be now detected around nearby low-mass stars using Doppler spectroscopy. In this paper, we investigate the radial velocity variations of Kapteyn's star, which is both a sub-dwarf M-star and the nearest halo object to the Sun. The observations comprise archival and new HARPS, HIRES and PFS Doppler measurements. Two Doppler signals are detected at periods of 48 and 120 days using likelihood periodograms and a Bayesian analysis of the data. Using the same techniques, the activity indicies and archival ASAS-3 photometry show evidence for low-level activity periodicities of the order of several hundred days. However, there are no significant correlations with the radial velocity variations on the same time-scales. The inclusion of planetary Keplerian signals in the model results in levels of correlated and excess white noise that are remarkably low compared to younger G, K and M dwarfs. We conclude that Kapteyn's star is most probably orbited by two super-Earth mass planets, one of which is orbiting in its circumstellar habitable zone, becoming the oldest potentially habitable planet known to date. The presence and long-term survival of a planetary system seems a remarkable feat given the peculiar origin and kinematic history of Kapteyn's star. The detection of super-Earth mass planets around halo stars provides important insights into planet-formation processes in the early days of the Milky Way.

Precision radial velocities from the Automated Planet Finder and Keck/HIRES reveal an M*sin(i) =18 +/- 2 Earth mass planet orbiting the nearby M3V star GJ 687. This planet has an orbital period, P = 38.14 days, and a low orbital eccentricity. Our Stromgren b and y photometry of the host star suggests a stellar rotation signature with a period of P = 60 days. The star is somewhat chromospherically active, with a spot filling factor estimated to be several percent. The rotationally{induced 60-day signal, however, is well-separated from the period of the radial velocity variations, instilling confidence in the interpretation of a Keplerian origin for the observed velocity variations. Although GJ 687b produces relatively little specific interest in connection with its individual properties, a compelling case can be argued that it is worthy of remark as an eminently typical, yet at a distance of 4.52 pc, a very nearby representative of the galactic planetary census. The detection of GJ 687b indicates that the APF telescope is well suited to the discovery of low-mass planets orbiting low-mass stars in the as-yet relatively un-surveyed region of the sky near the north celestial pole.

We present precision radial velocity (RV) data sets from Keck-HIRES and from Lick Observatory's new Automated Planet Finder Telescope and Levy Spectrometer on Mt. Hamilton that reveal a multiple-planet system orbiting the nearby, slightly evolved, K-type star HD 141399. Our 91 observations over 10.5 years suggest the presence of four planets with orbital periods of 94.35, 202.08, 1070.35, and 3717.35 days and minimum masses of 0.46, 1.36, 1.22, and 0.69 Jupiter masses respectively. The orbital eccentricities of the three inner planets are small, and the phase curves are well sampled. The inner two planets lie just outside the 2:1 resonance, suggesting that the system may have experienced dissipative evolution during the protoplanetary disk phase. The fourth companion is a Jupiter-like planet with a Jupiter-like orbital period. Its orbital eccentricity is consistent with zero, but more data will be required for an accurate eccentricity determination.

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.

The K2.5 dwarf HD 40307 has been reported to host three super-Earths. The system lacks massive planets and is therefore a potential candidate for having additional low-mass planetary companions. We re-derive Doppler measurements from public HARPS spectra of HD 40307 to confirm the significance of the reported signals using independent data analysis methods. We also investigate these measurements for additional low-amplitude signals. We used Bayesian analysis of our radial velocities to estimate the probability densities of different model parameters. We also estimated the relative probabilities of models with differing numbers of Keplerian signals and verified their significance using periodogram analyses. We investigated the relation of the detected signals with the chromospheric emission of the star. As previously reported for other objects, we found that radial velocity signals correlated with the S-index are strongly wavelength dependent. We identify two additional clear signals with periods of 34 and 51 days, both corresponding to planet candidates with minimum masses a few times that of the Earth. An additional sixth candidate is initially found at a period of 320 days. However, this signal correlates strongly with the chromospheric emission from the star and is also strongly wavelength dependent. When analysing the red half of the spectra only, the five putative planetary signals are recovered together with a very significant periodicity at about 200 days. This signal has a similar amplitude as the other new signals reported in the current work and corresponds to a planet candidate with M sin i = 7 Me (HD 40307 g). ...

Keck/HIRES precision radial velocities of HD 207832 indicate the presence of two Jovian-type planetary companions in Keplerian orbits around this G star. The planets have minimum masses of 0.56 and 0.73 Jupiter-masses with orbital periods of ~162 and ~1156 days, and eccentricities of 0.13 and 0.27, respectively. Stromgren b and y photometry reveals a clear stellar rotation signature of the host star with a period of 17.8 days, well separated from the period of the radial velocity variations, reinforcing their Keplerian origin. The values of the semimajor axes of the planets suggest that these objects have migrated from the region of giant planet formation to closer orbits. In order to examine the possibility of the existence of additional (small) planets in the system, we studied the orbital stability of hypothetical terrestrial-sized objects in the region between the two planets and interior to the orbit of the inner body. Results indicated that stable orbits exist only in a small region interior to planet b. However, the current observational data offer no evidence for the existence of additional objects in this system.

We report the detection of two new planets from the Anglo-Australian Planet Search. These planets orbit two stars each previously known to host one planet. The new planet orbiting HD 142 has a period of 6005\pm427 days, and a minimum mass of 5.3M_Jup. HD142c is thus a new Jupiter analog: a gas-giant planet with a long period and low eccentricity (e = 0.21 \pm 0.07). The second planet in the HD 159868 system has a period of 352.3\pm1.3 days, and m sin i=0.73\pm0.05 M_Jup. In both of these systems, including the additional planets in the fitting process significantly reduced the eccentricity of the original planet. These systems are thus examples of how multiple-planet systems can masquerade as moderately eccentric single-planet systems.

We re-analyze 4 years of HARPS spectra of the nearby M1.5 dwarf GJ 667C available through the ESO public archive. The new radial velocity (RV) measurements were obtained using a new data analysis technique that derives the Doppler measurement and other instrumental effects using a least-squares approach. Combining these new 143 measurements with 41 additional RVs from the Magellan/PFS and Keck/HIRES spectrometers, reveals 3 additional signals beyond the previously reported 7.2-day candidate, with periods of 28 days, 75 days, and a secular trend consistent with the presence of a gas giant (Period sim 10 years). The 28-day signal implies a planet candidate with a minimum mass of 4.5 Mearth orbiting well within the canonical definition of the star's liquid water habitable zone, this is, the region around the star at which an Earth-like planet could sustain liquid water on its surface. Still, the ultimate water supporting capability of this candidate depends on properties that are unknown such as its albedo, atmospheric composition and interior dynamics. The 75-day signal is less certain, being significantly affected by aliasing interactions among a potential 91-day signal, and the likely rotation period of the star at 105 days detected in two activity indices. GJ 667C is the common proper motion companion to the GJ 667AB binary, which is metal poor compared to the Sun. The presence of a super-Earth in the habitable zone of a metal poor M dwarf in a triple star system, supports the evidence that such worlds should be ubiquitous in the Galaxy.

We have obtained precision astrometry of the planet hosting M dwarf GJ 317 in the framework of the Carnegie Astrometric Planet Search project. The new astrometric measurements give a distance determination of 15.3 pc, 65% further than previous estimates. The resulting absolute magnitudes suggest it is metal rich and more massive than previously assumed. This result strengthens the correlation between high metallicity and the presence of gas giants around low mass stars. At 15.3 pc, the minimal astrometric amplitude for planet candidate GJ 317b is 0.3 milliarcseconds (edge-on orbit), just below our astrometric sensitivity. However, given the relatively large number of observations and good astrometric precision, a Bayesian Monte Carlo Markov Chain analysis indicates that the mass of planet b has to be smaller than twice the minimum mass with a 99% confidence level, with a most likely value of 2.5 Mjup. Additional RV measurements obtained with Keck by the Lick-Carnegie Planet search program confirm the presence of an additional very long period planet candidate, with a period of 20 years or more. Even though such an object will imprint a large astrometric wobble on the star, its curvature is yet not evident in the astrometry. Given high metallicity, and the trend indicating that multiple systems are rich in low mass companions, this system is likely to host additional low mass planets in its habitable zone that can be readily detected with state-of-the-art optical and near infrared RV measurements.

We present new precise HIRES radial velocity (RV) data sets of five nearby stars obtained at Keck Observatory. HD 31253, HD 218566, HD 177830, HD 99492 and HD 74156 are host stars of spectral classes F through K and show radial velocity variations consistent with new or additional planetary companions in Keplerian motion. The orbital parameters of the candidate planets in the five planetary systems span minimum masses of M sin i = 27.43 M_{earth} to M sin i = 8.28 M_{jup}, periods of 17.05 to 4696.95 days and eccentricities ranging from circular to extremely eccentric (e ~ 0.63). The 5th star, HD 74156, was known to have both a 52-day and a 2500-day planet, and was claimed to also harbor a 3rd planet at 336d, in apparent support of the "Packed Planetary System" hypothesis. Our greatly expanded data set for HD 74156 provides strong confirmation of both the 52-day and 2500-d planets, but strongly contradicts the existence of a 336-day planet, and offers no significant evidence for any other planets in the system.

We present 11 years of HIRES precision radial velocities (RV) of the nearby M3V star Gliese 581, combining our data set of 122 precision RVs with an existing published 4.3-year set of 119 HARPS precision RVs. The velocity set now indicates 6 companions in Keplerian motion around this star. Differential photometry indicates a likely stellar rotation period of ~94 days and reveals no significant periodic variability at any of the Keplerian periods, supporting planetary orbital motion as the cause of all the radial velocity variations. The combined data set strongly confirms the 5.37-day, 12.9-day, 3.15-day, and 67-day planets previously announced by Bonfils et al. (2005), Udry et al. (2007), and Mayor et al (2009). The observations also indicate a 5th planet in the system, GJ 581f, a minimum-mass 7.0 M_Earth planet orbiting in a 0.758 AU orbit of period 433 days and a 6th planet, GJ 581g, a minimum-mass 3.1 M_Earth planet orbiting at 0.146 AU with a period of 36.6 days. The estimated equilibrium temperature of GJ 581g is 228 K, placing it squarely in the middle of the habitable zone of the star and offering a very compelling case for a potentially habitable planet around a very nearby star. That a system harboring a potentially habitable planet has been found this nearby, and this soon in the relatively early history of precision RV surveys, indicates that eta_Earth, the fraction of stars with potentially habitable planets, is likely to be substantial. This detection, coupled with statistics of the incompleteness of present-day precision RV surveys for volume-limited samples of stars in the immediate solar neighborhood suggests that eta_Earth could well be on the order of a few tens of percent. If the local stellar neighborhood is a representative sample of the galaxy as a whole, our Milky Way could be teeming with potentially habitable planets.

(Abreviated) Continued radial velocity monitoring of the nearby M4V red dwarf star GJ~876 with Keck/HIRES has revealed the presence of a Uranus-mass fourth planetary companion in the system. The new planet has a mean period of $P_e=126.6$ days (over the 12.6-year baseline of the radial velocity observations), and a minimum mass of $m_e\sin{i_e}=12.9\pm 1.7\,M_{\oplus}$. Self-consistent, N-body fits to the radial velocity data set show that the four-planet system has an invariable plane with an inclination relative to the plane of the sky of $i=59.5^{\circ}$. The fit is not significantly improved by the introduction of a mutual inclination between the planets ``b'' and ``c,'' but the new data do confirm a non-zero eccentricity, $e_d=0.207\pm0.055$ for the innermost planet, ``d.'' In our best-fit coplanar model, the mass of the new component is $m_e=14.6\pm1.7\,M_{\oplus}$. Our best-fitting model places the new planet in a 3-body resonance with the previously known giant planets (which have mean periods of $P_c=30.4$ and $P_b=61.1$ days). The critical argument, $\varphi_{\rm Laplace}=\lambda_c-3\lambda_b+2\lambda_e$, for the Laplace resonance librates with an amplitude of $\Delta\varphi_{\rm Laplace}=40\pm13^{\circ}$ about $\varphi_{\rm Laplace}=0^{\circ}$. Numerical integration indicates that the four-planet system is stable for at least a billion years (at least for the coplanar cases). This resonant configuration of three giant planets orbiting an M-dwarf primary differs from the well-known Laplace configuration of the three inner Galilean satellites of Jupiter, which are executing very small librations about $\varphi_{\rm Laplace}=180^{\circ}$, and which never experience triple conjunctions. The GJ~876 system, by contrast, comes close to a triple conjunction between the outer three planets once per every orbit of the outer planet, ``e.''

Precision radial velocities from Keck/HIRES reveal a Saturn-mass planet orbiting the nearby M4V star HIP 57050. The planet has a minimum mass of 0.3 Jupiter-mass, an orbital period of 41.4 days, and an orbital eccentricity of 0.31. V-band photometry reveals a clear stellar rotation signature of the host star with a period of 98 days, well separated from the period of the radial velocity variations and reinforcing a Keplerian origin for the observed velocity variations. The orbital period of this planet corresponds to an orbit in the habitable zone of HIP 57050, with an expected planetary temperature of approximately 230 K. The star has a metallicity of [Fe/H] = 0.32+/-0.06 dex, of order twice solar and among the highest metallicity stars in the immediate solar neighborhood. This newly discovered planet provides further support that the well-known planet-metallicity correlation for F, G, and K stars also extends down into the M-dwarf regime. The a priori geometric probability for transits of this planet is only about 1%. However, the expected eclipse depth is ~7%, considerably larger than that yet observed for any transiting planet. Though long on the odds, such a transit is worth pursuing as it would allow for high quality studies of the atmosphere via transmission spectroscopy with HST. At the expected planetary effective temperature, the atmosphere may contain water clouds.

We present precision radial velocity (RV) data that reveal a multiple exoplanet system orbiting the bright nearby G5V star 61 Virginis. Our 4.6 years of combined Keck/HIRES and Anglo-Australian Telescope precision RVs indicate the hitherto unknown presence of at least three planets orbiting this well-studied star. These planets are all on low-eccentricity orbits with periods of 4.2, 38.0, and 124.0 days, and projected masses (M sini) of 5.1, 18.2, and 24.0 M_Earth, respectively. Test integrations of systems consistent with the RV data suggest that the configuration is dynamically stable. Depending on the effectiveness of tidal dissipation within the inner planet, the inner two planets may have evolved into an eccentricity fixed-point configuration in which the apsidal lines of all three planets corotate. This conjecture can be tested with additional observations. We present a 16-year time series of photometric observations of 61 Virginis, which comprise 1194 individual measurements, and indicate that it has excellent photometric stability. No significant photometric variations at the periods of the proposed planets have been detected. This new system is the first known example of a G-type Sun-like star hosting a Super-Earth mass planet. It joins HD 75732 (55 Cnc), HD 69830, GJ 581, HD 40307, and GJ 876 as a growing group of exoplanet systems that have multiple planets orbiting with periods less than an Earth-year. The ubiquity of such systems portends that space-based transit-search missions such as KEPLER and COROT will find many multi-transiting systems.

The Doppler wobble induced by the extra-solar planet HD 134987b was first detected by data from the Keck Telescope nearly a decade ago, and was subsequently confirmed by data from the Anglo-Australian Telescope. However, as more data have been acquired for this star over the years since, the quality of a single Keplerian fit to that data has been getting steadily worse. The best fit single Keplerian to the 138 Keck and AAT observations now in hand has an root-mean-square (RMS) scatter of 6.6 m/s. This is significantly in excess of both the instrumental precision achieved by both the Keck and Anglo-Australian Planet Searches for stars of this magnitude, and of the jitter expected for a star with the properties of HD134987. However, a double Keplerian (i.e. dual planet) fit delivers a significantly reduced RMS of 3.3 m/s. The best-fit double planet solution has minimum planet masses of 1.59 and 0.82Mjup, orbital periods of 258 and 5000d, and eccentricities of 0.23 and 0.12 respectively. We find evidence that activity-induced jitter is a significant factor in our fits and do not find evidence for asteroseismological p-modes. We also present seven years of photometry at a typical precision of 0.003mag with the T8 0.8m automatic photometric telescope at Fairborn observatory. These observations do not detect photometric variability and support the inference that the detected radial-velocity periods are due to planetary mass companions rather than due to photospheric spots and plages.

We present precision radial velocity data that reveal a Super-Earth mass planet and two probable additional planets orbiting the bright nearby G0V star HD 1461. Our 12.8 years of Keck HIRES precision radial velocities indicate the presence of a 7.4M_Earth planet on a 5.77-day orbit. The data also suggest, but cannot yet confirm, the presence of outer planets on low-eccentricity orbits with periods of 446.1 and 5017 days, and projected masses (M sin i) of 27.9 and 87.1M_Earth, respectively. Test integrations of systems consistent with the radial velocity data suggest that the configuration is dynamically stable. We present a 12.2-year time series of photometric observations of HD 1461, which comprise 799 individual measurements, and indicate that it has excellent long-term photometric stability. However, there are small amplitude variations with periods comparable to those of the suspected 2nd and 3rd signals in the radial velocities near 5000 and 446 days, thus casting some suspicion on those periodicities as Keplerian signals. If the 5.77-day companion has a Neptune-like composition, then its expected transit depth is of order ~0.5 millimags. The geometric a priori probability of transits is ~8%. Phase-folding of the ground-based photometry shows no indication that transits of the 5.77-day companion are occurring, but high-precision follow-up of HD 1461 during upcoming transit phase windows will be required to definitively rule out or confirm transits. This new system joins a growing list of solar-type stars in the immediate galactic neighborhood that are accompanied by at least one Neptune- (or lower) mass planets having orbital periods of 50 days or less.

We present the systemic Console, a new all-in-one, general-purpose software package for the analysis and combined multiparameter fitting of Doppler radial velocity (RV) and transit timing observations. We give an overview of the computational algorithms implemented in the Console, and describe the tools offered for streamlining the characterization of planetary systems. We illustrate the capabilities of the package by analyzing an updated radial velocity data set for the HD128311 planetary system. HD128311 harbors a pair of planets that appear to be participating in a 2:1 mean motion resonance. We show that the dynamical configuration cannot be fully determined from the current data. We find that if a planetary system like HD128311 is found to undergo transits, then self-consistent Newtonian fits to combined radial velocity data and a small number of timing measurements of transit midpoints can provide an immediate and vastly improved characterization of the planet's dynamical state.

We investigate the possibility that the large orbital eccentricity of the transiting Neptune-mass planet Gliese 436b is maintained in the face of tidal dissipation by a second planet in the system. We find that the currently observed configuration can be understood if Gliese 436b and a putative companion have evolved to a quasi-stationary fixed point in which the planets' orbital apses are co-linear and in which secular variations in the orbital eccentricities of the two planets have been almost entirely damped out. We adopt an octopole-order secular theory based on a Legendre expansion in the semi-major axis ratio to delineate well-defined regions of (P_c, M_c, e_c) space that can be occupied by a perturbing companion. We incorporate the evolutionary effect of tidal dissipation into our secular model of the system, and solve the resulting initial value problems for a large sample of the allowed configurations. We then polish the stationary configurations derived from secular theory with full numerical integrations. We present our results in the form of candidate companion planets to Gliese 436b. For these candidates, radial velocity half-amplitudes, K_c, are of order 3 m/s, and the maximum amplitude of orbit-to-orbit transit timing variations are of order Delta t=1 s to Delta t=5s. For the particular example case of a perturber with orbital period, P_c=40 d, mass, M_c=8.5 M_Earth, and eccentricity, e_c=0.58, we confirm our semi-analytic calculations with a full numerical 3-body integration of the orbital decay that includes tidal damping and spin evolution. Additionally, we discuss the possibility of many-perturber stationary configurations, utilizing modified Laplace-Lagrange secular theory.

Precision Doppler measurements from an intensive 48 night "Rocky Planet Search" observing campaign on the Anglo-Australian Telescope (AAT) have revealed the presence of a low-mass exoplanet orbiting the G1 dwarf HD16417. Subsequent Doppler observations with the AAT, as well as independent observations obtained by the Keck Planet Search, have confirmed this initial detection and refine the orbital parameters to period 17.24+/-0.01 d, eccentricity 0.20+/-0.09, orbital semi-major axis 0.14+/-0.01 AU and minimum planet mass 22.1+/-2.0 Mearth. HD 16417 raises the number of published exoplanets with minimum masses of less than 25 Mearth to eighteen. Interestingly, the distribution of detected sub-25 Mearth planets over the spectral types G, K and M is almost uniform. The detection of HD 16417b by an intensive observing campaign clearly demonstrates the need for extended and contiguous observing campaigns when aiming to detect low-amplitude Doppler planets in short period orbits. Perhaps most critically it demonstrates that the search for low-mass Doppler planets will eventually require these traditional "bright-time" projects to extend throughout dark lunations.

We simulate the formation of planetary systems around Alpha Centauri B. The N-body accretionary evolution of a 1/r disk populated with 400-900 lunar-mass protoplanets is followed for 200 Myr. All simulations lead to the formation of multiple-planet systems with at least one planet in the 1-2 MEarth mass range at 0.5-1.5 AU. We examine the detectability of our simulated planetary systems by generating synthetic radial velocity observations including noise based on the radial velocity residuals to the recently published three planet fit to the nearby K0V star HD 69830. Using these synthetic observations, we find that we can reliably detect a 1.8 MEarth planet in the habitable zone of Alpha Centauri B after only three years of high cadence observations. We also find that the planet is detectable even if the radial velocity precision is 3 m/s, as long as the noise spectrum is white. Our results show that the greatest uncertainty in our ability to detect rocky planets in the Alpha Centauri system is the unknown magnitude of ultra-low frequency stellar noise.

We report the results of a globally coordinated photometric campaign to search for transits by the P ~ 30 d and P ~ 60 d outer planets of the 3-planet system orbiting the nearby M-dwarf Gl 876. These two planets experience strong mutual perturbations, which necessitate use of a dynamical (four-body) model to compute transit ephemerides for the system. Our photometric data have been collected from published archival sources, as well as from our photometric campaigns that were targeted to specific transit predictions. Our analysis indicates that transits by planet "c" (P ~ 30 d) do not currently occur, in concordance with the best-fit i = 50 degree co-planar configuration obtained by dynamical fits to the most recent radial velocity data for the system. Transits by planet "b" (P ~ 60 d) are not entirely ruled out by our observations, but our data indicate that it is very unlikely that they occur. Our experience with the Gl 876 system suggests that a distributed ground-based network of small telescopes can be used to search for transits of very low mass M-stars by terrestrial-sized planets.

High precision, high cadence radial velocity monitoring over the past 8 years at the W. M. Keck Observatory reveals evidence for a third planet orbiting the nearby (4.69 pc) dM4 star GJ 876. The residuals of three-body Newtonian fits, which include GJ 876 and Jupiter mass companions b and c, show significant power at a periodicity of 1.9379 days. Self-consistently fitting the radial velocity data with a model that includes an additional body with this period significantly improves the quality of the fit. These four-body (three-planet) Newtonian fits find that the minimum mass of companion ``d'' is m sin i = 5.89 +- 0.54 Earth masses and that its orbital period is 1.93776 (+- 7x10^-5) days. Assuming coplanar orbits, an inclination of the GJ 876 planetary system to the plane of the sky of ~ 50 degrees gives the best fit. This inclination yields a mass for companion d of m = 7.53 +- 0.70 Earth masses, making it by far the lowest mass companion yet found around a main sequence star other than our Sun. Precise photometric observations at Fairborn Observatory confirm low-level brightness variability in GJ 876 and provide the first explicit determination of the star's 96.7-day rotation period. Even higher precision short-term photometric measurements obtained at Las Campanas imply that planet d does not transit GJ 876.

We present results of extensive numerical studies of the stability of non-interacting particles in the planetary systems of stars Upsilon Andromedae, GJ 876, 47 UMa, and 55 Cancri. We discuss the possibility of the existence of islands of stability and/or instability at different regions in these multi-body systems, and their probable correspondence to certain mean-motion resonances. The results of our study may be applied to questions concerning the stability of terrestrial planets in these systems and also the trapping of particles in resonances with the planets.

In this paper, we present preliminary results on the stability of massless particles in two and three-planet systems. The results of our study may be used to address questions concerning the stability of terrestrial planets in these systems and also the trapping of particles in resonances with the planets. The possibility of the existence of islands of stability and/or instability at different regions in multi-body systems and their probable correspondence to certain resonances are also discussed.