The nearby (6.5 pc) star HD 219134 was recently shown by Motalebi et al. (2015) and Vogt et al. (2015) to host several planets, the innermost of which is transiting. We present twenty-seven years of radial velocity observations of this star from the McDonald Observatory Planet Search program, and nineteen years of stellar activity data. We detect a long-period activity cycle measured in the Ca II $S_{HK}$ index, with a period of $4230 \pm 100$ days (11.7 years), very similar to the 11-year Solar activity cycle. Although the period of the Saturn-mass planet HD 219134 h is close to half that of the activity cycle, we argue that it is not an artifact due to stellar activity. We also find a significant periodicity in the $S_{HK}$ data due to stellar rotation with a period of 22.8 days. This is identical to the period of planet f identified by Vogt et al. (2015), suggesting that this radial velocity signal might be caused by rotational modulation of stellar activity rather than a planet. Analysis of our radial velocities allows us to detect the long-period planet HD 219134 h and the transiting super-Earth HD 219134 b. Finally, we use our long time baseline to constrain the presence of longer-period planets in the system, excluding to $1\sigma$ objects with $M\sin i>0.36 M_J$ at 12 years (corresponding to the orbital period of Jupiter) and $M\sin i>0.72 M_J$ at a period of 16.4 years (assuming a circular orbit for an outer companion).

We report the detection of two new long-period giant planets orbiting the stars HD 95872 and HD 162004 (psi1 Draconis B) by the McDonald Observatory planet search. The planet HD 95872b has a minimum mass of 4.6 M_Jup and an orbital semi-major axis of 5.2 AU. The giant planet psi1 Dra Bb has a minimum mass of 1.5 M_Jup and an orbital semi-major axis of 4.4 AU. Both of these planets qualify as Jupiter analogs. These results are based on over one and a half decades of precise radial velocity measurements collected by our program using the McDonald Observatory Tull Coude spectrograph at the 2.7 m Harlan J. Smith telescope. In the case of psi1 Draconis B we also detect a long-term non-linear trend in our data that indicates the presence of an additional giant planet, similar to the Jupiter-Saturn pair. The primary of the binary star system, psi1 Dra A, exhibits a very large amplitude radial velocity variation due to another stellar companion. We detect this additional member using speckle imaging. We also report two cases - HD 10086 and HD 102870 (beta Virginis) - of significant radial velocity variation consistent with the presence of a planet, but that are probably caused by stellar activity, rather than reflexive Keplerian motion. These two cases stress the importance of monitoring the magnetic activity level of a target star, as long-term activity cycles can mimic the presence of a Jupiter-analog planet.

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 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.

We carried out photometric observations for HBC 722 in SDSS r, i and z bands from 2011 April to 2013 May with a Camera for Quasars in Early uNiverse attached to the 2.1m Otto Struve telescope at McDonald Observatory. The post-outburst phenomena were classified into five phases according to not only brightness but also color variations, which might be caused by physical changes in the emitting regions of optical and near-infrared bands. A series of spectral energy distribution (SED) is presented to support color variations and track the time evolution of SED in optical/near-infrared bands after the outburst. Given two years of data, possible periodicities of r, i and z bands were checked. We found out three families of signals around ~6, ~10 and ~1 days in three bands, which is broadly consistent with Green et al. (2013). We also examined short term variability (intra-day and day scales) to search for evidences of flickering by using the micro-variability method. We found clear signs of day scale variability and weak indications of intra-day scale fluctuations, which implies that the flickering event occurs in HBC 722 after outburst.

Recent simulations have shown that the formation of planets in circumbinary configurations (such as those recently discovered by Kepler) is dramatically hindered at the planetesimal accretion stage. The combined action of the binary and the protoplanetary disk acts to raise impact velocities between km-sized planetesimals beyond their destruction threshold, halting planet formation within at least 10 AU from the binary. It has been proposed that a primordial population of "large" planetesimals (100 km or more in size), as produced by turbulent concentration mechanisms, would be able to bypass this bottleneck; however, it is not clear whether these processes are viable in the highly perturbed circumbinary environments. We perform two-dimensional hydrodynamical and N-body simulations to show that km-sized planetesimals and collisional debris can drift and be trapped in a belt close to the central binary. Within this belt, planetesimals could initially grow by accreting debris, ultimately becoming "indestructible" seeds that can accrete other planetesimals in-situ despite the large impact speeds. We find that large, indestructible planetesimals can be formed close to the central binary within $10^5$ years, therefore showing that even a primordial population of "small" planetesimals can feasibly form a planet.

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.

We present the detection of day-timescale periodic variability in the r-band lightcurve of newly outbursting FU Orionis-type object HBC 722, taken from > 42 nights of observation with the CQUEAN instrument on the McDonald Observatory 2.1m telescope. The optical/near-IR lightcurve of HBC 722 shows a complex array of periodic variability, clustering around 5.8 day (0.044 mag amplitude) and 1.28 day (0.016 mag amplitude) periods, after removal of overall baseline variation. We attribute the unusual number of comparable strength signals to a phenomenon related to the temporary increase in accretion rate associated with FUors. We consider semi-random "flickering", magnetic braking/field compression and rotational asymmetries in the disk instability region as potential sources of variability. Assuming the 5.8 day period is due to stellar rotation and the 1.28 day period is indicative of Keplerian rotation at the inner radius of the accretion disk (at 2 R(star)), we derive a B-field strength of 2.2-2.7 kG, slightly larger than typical T Tauri stars. If instead the 5.8 day signal is from a disk asymmetry, the instability region has an outer radius of 5.4 R(star), consistent with models of FUor disks. Further exploration of the time domain in this complicated source and related objects will be key to understanding accretion processes.

The existence of planets born in environments highly perturbed by a stellar companion represents a major challenge to the paradigm of planet formation. In numerical simulations, the presence of a close binary companion stirs up the relative velocity between planetesimals, which is fundamental in determining the balance between accretion and erosion. However, the recent discovery of circumbinary planets by Kepler establishes that planet formation in binary systems is clearly viable. We perform N-body simulations of planetesimals embedded in a protoplanetary disk, where planetesimal phasing is frustrated by the presence of stochastic torques, modeling the expected perturbations of turbulence driven by the magnetorotational instability (MRI). We examine perturbation amplitudes relevant to dead zones in the midplane (conducive to planet formation in single stars), and find that planetesimal accretion can be inhibited even in the outer disk (4-10 AU) far from the central binary, a location previously thought to be a plausible starting point for the formation of circumbinary planets.

The recently discovered circumbinary planets (Kepler-16 b, Kepler-34 b, Kepler-35 b) represent the first direct evidence of the viability of planet formation in circumbinary orbits. We report on the results of N-body simulations investigating planetesimal accretion in the Kepler-16 b system, focusing on the range of impact velocities under the influence of both stars' gravitational perturbation and friction from a putative protoplanetary disk. Our results show that planet formation might be effectively inhibited for a large range in semi-major axis (1.75 < a_P < 4 AU), suggesting that the planetary core must have either migrated from outside 4 AU, or formed in situ very close to its current location.

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.

(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.''

Transit timing variations - deviations from strict periodicity between successive passages of a transiting planet - can be used to probe the structure and dynamics of multiple-planet systems. In this paper, we examine prospects for numerically solving the so-called inverse problem, the determination of the orbital elements of a perturbing body from the transit timing variations it induces. We assume that the planetary systems under examination have a limited number of Doppler velocity measurements, and show that a more extensive radial velocity characterization with precision comparable to the semiamplitude of the perturber may remove degeneracies in the solution. We examine several configurations of interest, including (1) a prototypical non-resonant system, modeled after HD40307 b and c, which contains multiple super-Earth mass planets, (2) a hypothetical system containing a transiting giant planet with a terrestrial-mass companion trapped in low-order mean motion resonance, and (3) the HAT-P-13 system, in which forced precession by an outer perturbing body that is well characterized by Doppler radial velocity measurements can give insight into the interior structure of a perturbing planet, and for which the determination of mutual inclination between the transiting planet and its perturber is a key issue.

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.

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 (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.

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.

In this letter, we briefly describe the evolution of a variety of self-gravitating protoplanetary disk models that contain annular grooves (e.g. gaps) in their surface density. These grooves are inspired by the density gaps that are presumed to open in response to the formation of a giant planet. Our work provides an extension of the previously studied groove modes that are known in the context of stellar disks. The emergence of spiral gravitational instabilities (GI) is predicted via a generalized eigenvalue code that performs a linear analysis, and confirmed with hydrodynamical simulations. We find the presence of a groove drives a fast-growing two-armed mode in moderately massive disks, and extends the importance of self-gravitating instabilities down to lower disk masses than for which they would otherwise occur. We discuss the potential importance of this instability in the context of planet formation, e.g. the modification of the torques driving Type II migration.