The next generation of exoplanet-hunting spectrographs should deliver up to an order of magnitude improvement in radial velocity precision over the standard 1 m/s state of the art. This advance is critical for enabling the detection of Earth-mass planets around Sun-like stars. New calibration techniques such as laser frequency combs and stabilized etalons ensure that the instrumental stability is well characterized. However, additional sources of error include stellar noise, undetected short-period planets, and telluric contamination. To understand and ultimately mitigate error sources, the contributing terms in the error budget must be isolated to the greatest extent possible. Here, we introduce a new high cadence radial velocity program, the EXPRES 100 Earths program, which aims to identify rocky planets around bright, nearby G and K dwarfs. We also present a benchmark case: the 62-d orbit of a Saturn-mass planet orbiting the chromospherically quiet star, HD 3651. The combination of high eccentricity (0.6) and a moderately long orbital period, ensures significant dynamical clearing of any inner planets. Our Keplerian model for this planetary orbit has a residual RMS of 58 cm/s over a $\sim 6$ month time baseline. By eliminating significant contributors to the radial velocity error budget, HD 3651 serves as a standard for evaluating the long term precision of extreme precision radial velocity (EPRV) programs.

`Oumuamua (I1 2017) was the first macroscopic ($l\sim100\,{\rm m}$) body observed to traverse the inner solar system on an unbound hyperbolic orbit. Its light curve displayed strong periodic variation, and it showed no hint of a coma or emission from molecular outgassing. Astrometric measurements indicate that 'Oumuamua experienced non-gravitational acceleration on its outbound trajectory, but energy balance arguments indicate this acceleration is inconsistent with a water ice sublimation-driven jet of the type exhibited by solar system comets. We show that all of `Oumaumua's observed properties can be explained if it contained a significant fraction of molecular hydrogen (H$_{2}$) ice. H$_{2}$ sublimation at a rate proportional to the incident solar flux generates a surface-covering jet that reproduces the observed acceleration. Mass wasting from sublimation leads to monotonic increase in the body axis ratio, explaining `Oumuamua's shape. Back-tracing `Oumuamua's trajectory through the Solar System permits calculation of its mass and aspect ratio prior to encountering the Sun. We show that H$_{2}$-rich bodies plausibly form in the coldest dense cores of Giant Molecular Clouds, where number densities are of order $n\sim10^5$, and temperatures approach the $T=3\,{\rm K}$ background. Post-formation exposure to galactic cosmic rays implies a $\tau \sim 100$ Myr age, explaining the kinematics of `Oumuamua's inbound trajectory.

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

The discovery of 'Oumuamua (1I/2017 U1), the first interstellar interloper, suggests an abundance of free-floating small bodies whose ejection into galactic space cannot be explained by the current population of confirmed exoplanets. Shortly after 'Oumuamua's discovery, observational results from the DSHARP survey illustrated the near-ubiquity of ring/gap substructures within protoplanetary disks, strongly suggesting the existence of a vast population of as-yet undetected wide-separation planets that are capable of efficiently ejecting debris from their environments. These planets have $a \gtrsim 5$ au and masses of order Neptune's or larger, and they may accompany $\sim$50% of newly formed stars (Zhang et al. 2018). We combine the DSHARP results with statistical constraints from current time-domain surveys to quantify the population of detectable icy planetesimals ejected by disk-embedded giant planets through gravity assists. Assessment of the expected statistical distribution of interstellar objects is critical to accurately plan for and interpret future detections. We show that the number density of interstellar objects implied by 'Oumuamua is consistent with 'Oumuamua itself having originated as an icy planetesimal ejected from a DSHARP-type system via gravity assists, with the caveat that 'Oumuamua's lack of observed outgassing remains in strong tension with a cometary origin. Under this interpretation, 'Oumuamua's detection points towards a large number of long-period giant planets in extrasolar systems, supporting the hypothesis that the observed gaps in protoplanetary disks are carved by planets. In the case that 'Oumuamua is an ejected cometary planetesimal, we conclude that LSST should detect up to a few interstellar objects per year of 'Oumuamua's size or larger and over 100 yr$^{-1}$ for objects with $r > 1\,{\rm m}$.

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.

Recent work suggests that many short-period extrasolar planets may have spin obliquities that are significantly tilted with respect to their orbital planes. These large obliquities are a natural outcome of ``secular spin-orbit resonance'', a configuration in which the planetary spin precession frequency matches the frequency of orbit nodal regression, or a Fourier component thereof. While exoplanet spin obliquities have not yet been measured directly, they may be detectable indirectly through their signatures in various observations, such as photometric measurements across the full phase of a planet's orbit. In this work, we employ a thermal radiative model to explore how large polar tilts affect full-phase light curves, and we discuss the range of unique signatures that are expected to result. We show that the well-studied short-period planets HD 149026 b, WASP-12 b, and CoRoT-2 b all exhibit phase curve features that may arise from being in high-obliquity states. We also constrain the parameters and assess the detectability of hypothetical perturbing planets that could maintain the planets in these states. Among the three planets considered, CoRoT-2 b has the tightest constraints on its proposed obliquity ($45.8^{\circ} \pm 1.4^{\circ}$) and axial orientation. For HD 149026 b, we find no significant evidence for a non-zero obliquity, and the phase curve of WASP-12 b is too complicated by strong tidal distortions for a conclusive assessment.

We discuss the feasibility of and present initial designs and approximate cost estimates for a large ($N\sim2000$) network of small photometric telescopes that is purpose-built to monitor $V \lesssim 15$ Gaia Mission program stars for occultations by minor solar system bodies. The implementation of this network would permit measurement of the solar system's tidal gravity field to high precision, thereby revealing the existence of distant trans-Neptunian objects such as the proposed "Planet Nine." As a detailed example of the network capabilities, we investigate how occultations by Jovian Trojans can be monitored to track the accumulation of gravitational perturbations, thereby constraining the presence of undetected massive solar system bodies. We also show that the tidal influence of Planet Nine can be discerned from that of smaller, nearer objects in the Kuiper belt. Moreover, ephemerides for all small solar system bodies observed in occultation could be significantly improved using this network, thereby improving spacecraft navigation and refining Solar System modeling. Finally, occultation monitoring would generate direct measurements of size distributions for asteroid populations, permitting a better understanding of their origins.

We show that the $P\sim8\,{\rm h}$ photometric period and the astrometrically measured $A_{\rm ng}\sim2.5\times10^{-4}\,{\rm cm\,s^{-2}}$ non-gravitational acceleration (at $r\sim1.4\,{\rm AU}$) of the interstellar object 1I/2017 (`Oumuamua) can be explained by a nozzle-like venting of volatiles whose activity migrated to track the sub-solar location on the object's surface. Adopting the assumption that `Oumuamua was an elongated $a\times b \times c$ ellipsoid, this model produces a pendulum-like rotation of the body and implies a long semi-axis $a\sim 5A_{\rm ng}P^2/4\pi^2 \sim 260\,{\rm m}$. This scale agrees with the independent estimates of `Oumuamua's size that stem from its measured brightness, assuming an albedo of $p\sim0.1$, appropriate to ices that have undergone long-duration exposure to the interstellar cosmic ray flux. Using ray-tracing, we generate light curves for ellipsoidal bodies that are subject to both physically consistent sub-solar torques and to the time-varying geometry of the Sun-Earth-`Oumuamua configuration. Our synthetic light curves display variations from chaotic tumbling and changing cross-sectional illumination that are consistent with the observations, while avoiding significant secular changes in the photometric periodicity. If our model is correct, `Oumuamua experienced mass loss that wasted $\sim 10\%$ of its total mass during the $\sim 100\,{\rm d}$ span of its encounter with the inner Solar System and had an icy composition with a very low $[{\rm C}/{\rm O}]\lesssim 0.003$. Our interpretation of `Oumuamua's behavior is consistent with the hypothesis that it was ejected from either the outer regions of a planetesimal disk after an encounter with an embedded $M_{\rm p} \sim M_{\rm Nep}$ planet or from an exo-Oort cloud.

NASA's Kepler mission revealed that $\sim 30\%$ of Solar-type stars harbor planets with sizes between that of Earth and Neptune on nearly circular and co-planar orbits with periods less than 100 days. Such short-period compact systems are rarely found with planet pairs in mean-motion resonances (MMRs) -- configurations in which the planetary orbital periods exhibit a simple integer ratio -- but there is a significant overabundance of planet pairs lying just wide of the first-order resonances. Previous work suggests that tides raised on the planets by the host star may be responsible for forcing systems into these configurations by draining orbital energy to heat. Such tides, however, are insufficient unless there exists a substantial and as-yet unidentified source of extra dissipation. Here we show that this cryptic heat source may be linked to "obliquity tides" generated when a large axial tilt (obliquity) is maintained by secular resonance-driven spin-orbit coupling. We present evidence that typical compact, nearly-coplanar systems frequently experience this mechanism, and we highlight additional features in the planetary orbital period and radius distributions that may be its signatures. Extrasolar planets that maintain large obliquities will exhibit infrared light curve features that are detectable with forthcoming space missions. The observed period ratio distribution can be explained if typical tidal quality factors for super-Earths and sub-Neptunes are similar to those of Uranus and Neptune.

Recent analyses have revealed a mystery. The orbital period of the highly inflated hot Jupiter, WASP-12b, is decreasing rapidly. The rate of inspiral, however, is too fast to be explained by either eccentricity tides or equilibrium stellar tides. While dynamical stellar tides are possible, they require a subgiant structure for the star, whereas stellar models point toward a main sequence host. Here, we show that these hitherto irreconcilable observations might be explained by planetary obliquity tides if planet b's spin vector is trapped in a high-obliquity state maintained by a secular spin-orbit resonance with an unseen exterior perturbing planet. We derive constraints on the obliquity ($\epsilon\gtrsim50^{\circ}$), reduced tidal quality factor ($Q^{\prime}\sim10^{6}-10^{7}$), and perturbing planet parameters ($M_{2}\sim10-20M_{\oplus}$, $a_2\lesssim0.04\,{\rm AU}$) required to generate the observed orbital decay. Direct N-body simulations that include tidal and spin dynamics reinforce the plausibility of the scenario. Furthermore, we show that the resonance could have been captured when planet b's obliquity was small, making the proposed sequence of events easy to explain. The hypothetical perturbing planet is within the limits of current radial velocity constraints on the system yet is also detectable. If it exists, it could provide evidence in favor of the in situ formation hypothesis for hot Jupiters.

Models of the Sun's long-term evolution suggest that its luminosity was substantially reduced 2-4 billion years ago, which is inconsistent with substantial evidence for warm and wet conditions in the geological records of both ancient Earth and Mars. Typical solutions to this so-called "faint young Sun paradox" consider changes in the atmospheric composition of Earth and Mars, and while attractive, geological verification of these ideas is generally lacking-particularly for Mars. One possible underexplored solution to the faint young Sun paradox is that the Sun has simply lost a few percent of its mass during its lifetime. If correct, this would slow, or potentially even offset the increase in luminosity expected from a constant-mass model. However, this hypothesis is challenging to test. Here, we propose a novel observational proxy of the Sun's ancient mass that may be readily measured from accumulation patterns in sedimentary rocks on Earth and Mars. We show that the orbital parameters of the Solar system planets undergo quasi-cyclic oscillations at a frequency, given by secular mode g_2-g_5, that scales approximately linearly with the Sun's mass. Thus by examining the cadence of sediment accumulation in ancient basins, it is possible distinguish between the cases of a constant mass Sun and a more massive ancient Sun to a precision of greater than about 1 per cent. This approach provides an avenue toward verification, or of falsification, of the massive early Sun hypothesis.

We report on the confirmation of a transiting giant planet around the relatively hot (Teff = 6801 $\pm$ 56 K) star HD2685, whose transit signal was detected in Sector 1 data of the TESS mission. We confirmed the planetary nature of the transit signal by using Doppler velocimetric measurements with CHIRON, CORALIE and FEROS, as well as photometric data with CHAT and LCOGT. From the photometry and radial velocities joint analysis, we derived the following parameters for HD2685 $b$: $P$=4.12692$\pm$0.00004 days, M$_P$=1.18 $\pm$ 0.09 $M_J$ and $R_P$=1.44 $\pm$ 0.01 $R_J$. This system is a typical example of an inflated transiting Hot-Jupiter in a circular orbit. Given the host star apparent visual magnitude ($V$ = 9.6 mag), this is one of the brightest known stars hosting a transiting Hot-Jupiter, and a good example of the upcoming systems that will be detected by TESS during the two-year primary mission. This is also an excellent target for future ground and space based atmospheric characterization as well as a good candidate for measuring the projected spin-orbit misalignment angle via the Rossiter-McLaughlin effect.

We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known.

We report the photometry of six transits of the hot Jupiter HAT-P-29b obtained from 2013 October to 2015 January. We analyze the new light curves, in combination with the published photometric, and Doppler velocimetric, and spectroscopic measurements, finding an updated orbital ephemeris for the HAT-P-29 system, $T_{\rm C}[0]= 2456170.5494(15)\,[\rm{BJD_{TDB}}]$ and $P=5.723390(13) \,{\rm days}$. It is $17.63\,{\rm s}$ ($4.0\,\sigma$) longer than the previously published value, amounting to errors exceeding $2.5\,\rm{hrs}$ at the time of writing (on UTC 2018 June 1). The measured transit mid-times for HAT-P-29b show no compelling evidence of timing anomalies from a linear model, which rules out the presence of a perturbers with masses greater than $0.6$, $0.7$, $0.5$, and $0.4\,{\rm M_\oplus}$ near the $1:2$, $2:3$, $3:2$, and $2:1$ resonances with HAT-P-29b, respectively.

Although the formation and the properties of hot Jupiters (with orbital periods $P<10\,$d) have attracted a great deal of attention, the origins of warm Jupiters ($10<P<100\,$d) are less well-studied. Using a transit timing analysis, we present the orbital parameters of five planetary systems containing warm Jupiters, Kepler 30, Kepler 117, Kepler 302, Kepler 487 and Kepler 418. Three of them, Kepler-30 c($M_p=549.4{\pm{5.6}}M_{\oplus}$), Kepler-117 c($M_p=702{\pm{63}}M_{\oplus}$) and Kepler 302 c($M_p=933\pm 527M_{\oplus}$), are confirmed to be real warm Jupiters based on their mass. Insights drawn from the radius-temperature relationship lead to the inference that hot Jupiters and warm Jupiters can be roughly separated by ${T_{\rm eff,c}} =1123.7\pm3.3$ K. Also, ${T_{\rm eff,c}}$ provides a good separation for Jupiters with companion fraction consistent with zero($T_{\rm eff}>T_{\rm eff,c}$) and those with companion fraction significantly different from zero ($T_{\rm eff}<T_{\rm eff,c}$).

We present a uniform assessment of existing near-infrared Spitzer Space Telescope observations of planet-bearing stars. Using a simple four-parameter blackbody thermal model, we analyze stars for which photometry in at least one of Spitzer's IRAC bands has been obtained over either the entirety or a significant fraction of the planetary orbit. Systems in this category comprise ten well-studied systems with Hot Jupiters on circular or near-circular orbits (HAT-P-7, HD 149026, HD 189733, HD 209458, WASP-12, WASP-14, WASP-18, WASP-19, WASP-33, and WASP-43), as well as three stars harboring planets on significantly eccentric orbits (GJ 436, HAT-P-2, and HD 80606). We find that our simple model, in almost all cases, accurately reproduces the minimum and maximum planetary emission, as well as the phase offsets of these extrema with respect to transits/secondary eclipses. For one notable exception, WASP-12 b, adding an additional parameter to account for its tidal distortion is not sufficient to reproduce its photometric features. Full-orbit photometry is available in multiple wavelengths for 10 planets. We find that the returned parameter values for independent fits to each band are largely in agreement. However, disagreements in night-side temperature suggest distinct atmospheric layers, each with their own characteristic minimum temperature. In addition, a diversity in albedos suggests variation in opacity of the photospheres. While previous works have pointed out trends in photometric features based on system properties, we cannot conclusively identify analogous trends for physical model parameters. To make the connection between full-phase data and physical models more robust, a higher signal-to-noise must come from both increased resolution and a careful treatment of instrumental systematics.

We present eight new light curves of the transiting extra-solar planet HAT-P-25b obtained from 2013 to 2016 with three telescopes at two observatories. We use the new light curves, along with recent literature material, to estimate the physical and orbital parameters of the transiting planet. Specifically, we determine the mid-transit times (T$_{C}$) and update the linear ephemeris, T$_{C[0]}$=2456418.80996$\pm$0.00025 [$\mathrm{BJD}_\mathrm{TDB}$] and P=3.65281572$\pm$0.00000095 days. We carry out a search for transit timing variations (TTVs), and find no significant TTV signal at the $\Delta T=$80 s-level, placing a limit on the possible strength of planet-planet interactions ($\mathrm{TTV_{G}}$). In the course of our analysis, we calculate the upper mass-limits of the potential nearby perturbers. Near the 1:2, 2:1, and 3:1 resonances with HAT-P-25b, perturbers with masses greater than 0.5, 0.3, and 0.5 $\mathrm{M_{\oplus}}$ respectively, can be excluded. Furthermore, based on the analysis of TTVs caused by light travel time effect (LTTE) we also eliminate the possibility that a long-period perturber exists with $M_{\rm p}> 3000 \,\mathrm{M_{J}}$ within $a=11.2\,{\rm AU}$ of the parent star.

A rapid accumulation of observations and interpretation have followed in the wake of 1I `Oumuamua's passage through the inner Solar System. We briefly outline the consequences that this first detection of an interstellar asteroid implies for the planet-forming process, and we assess the near-term prospects for detecting and observing (both remotely and in situ) future Solar System visitors of this type. Drawing on detailed heat-transfer calculations that take both `Oumuamua's unusual shape and its chaotic tumbling into account, we affirm that the lack of a detectable coma in deep images of the object very likely arises from the presence of a radiation-modified coating of high molecular weight material (rather than a refractory bulk composition). Assuming that `Oumuamua is a typical representative of a larger population with a kinematic distribution similar to Population I stars in the local galactic neighborhood, we calculate expected arrival rates, impact parameters and velocities of similar objects and assess their prospects for detection using operational and forthcoming facilities. Using `Oumuamua as a proof-of-concept, we assess the prospects for missions that intercept interstellar objects (ISOs) using conventional chemical propulsion. Using a "launch on detection" paradigm, we estimate wait times of order 10 years between favorable mission opportunities with the detection capabilities of the Large-Scale Synoptic Survey Telescope (LSST), a figure that will be refined as the population of interstellar asteroids becomes observationally better constrained.

Gliese 876 harbors one of the most dynamically rich and well-studied exoplanetary systems. The nearby M4V dwarf hosts four known planets, the outer three of which are trapped in a Laplace mean-motion resonance. A thorough characterization of the complex resonant perturbations exhibited by the orbiting planets, and the chaotic dynamics therein, is key to a complete picture of the system's formation and evolutionary history. Here we present a reanalysis of the system using six years of new radial velocity (RV) data from four instruments. This new data augments and more than doubles the size of the decades-long collection of existing velocity measurements. We provide updated estimates of the system parameters by employing a computationally efficient Wisdom-Holman N-body symplectic integrator, coupled with a Gaussian Process (GP) regression model to account for correlated stellar noise. Experiments with synthetic RV data show that the dynamical characterization of the system can differ depending on whether a white noise or correlated noise model is adopted. Despite there being a region of stability for an additional planet in the resonant chain, we find no evidence for one. Our new parameter estimates place the system even deeper into resonance than previously thought and suggest that the system might be in a low energy, quasi-regular double apsidal corotation resonance. This result and others will be used in a subsequent study on the primordial migration processes responsible for the formation of the resonant chain.

The Kepler-9 system harbors three known transiting planets. The system holds significant interest for several reasons. First, the outer two planets exhibit a period ratio that is close to a 2:1 orbital commensurability, with attendant dynamical consequences. Second, both planets lie in the planetary mass "desert" that is generally associated with the rapid gas agglomeration phase of the core accretion process. Third, there exist attractive prospects for accurately measuring both the sky-projected stellar spin-orbit angles as well as the mutual orbital inclination between the planets in the system. Following the original \textit{Kepler} detection announcement in 2010, the initially reported orbital ephemerides for Kepler-9~b and c have degraded significantly, due to the limited time base-line of observations on which the discovery of the system rested. Here, we report new ground-based photometric observations and extensive dynamical modeling of the system. These efforts allow us to photometrically recover the transit of Kepler-9~b, and thereby greatly improve the predictions for upcoming transit mid-times. Accurate ephemerides of this system are important in order to confidently schedule follow-up observations of this system, for both in-transit Doppler measurements as well as for atmospheric transmission spectra taken during transit.

We present spectroscopic measurements of the Rossiter-McLaughlin effect for the planet b of Kepler-9 multi-transiting planet system. The resulting sky-projected spin-orbit angle is $\lambda=-13^{\circ} \pm 16^{\circ}$, which favors an aligned system and strongly disfavors highly misaligned, polar, and retrograde orbits. Including Kepler-9, there are now a total of 4 Rossiter-McLaughlin effect measurements for multiplanet systems, all of which are consistent with spin-orbit alignment.

The arrival of the robustly hyperbolic asteroid A/2017 U1 has potentially interesting ramifications for the planet-formation process. Although extrapolations from a sample size of one are necessarily uncertain, order-of-magnitude estimates suggest that the Galaxy contains a substantial mass in similar bodies. We argue that despite its lack of Coma, A/2017 U1 likely contained a significant mass fraction of volatile components, and we argue that its presence can be used to infer a potentially large population of as-yet undetected Neptune-like extrasolar planets.

The widespread prevalence of close-in, nearly coplanar super-Earth- and sub-Neptune-sized planets in multiple-planet systems was one of the most surprising results from the Kepler mission. By studying a uniform sample of Kepler "multis" with mass measurements from transit timing variations (TTVs), we show that a given planetary system tends to harbor a characteristic type of planet. That is, planets in a system have both masses and radii that are far more similar than if the system were assembled randomly from planets in the population. This finding has two important ramifications. First, the large intrinsic compositional scatter in the planet mass-radius relation is dominated by system-to-system variance rather than intra-system variance. Second, if provided enough properties of the star and primordial protoplanetary disk, there may be a substantial degree of predictability in the outcome of the planet formation process. We show that stellar mass and metallicity account for of order $20\%$ of the variation in outcomes; the remainder is as-yet unknown.

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)

Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are generally amenable to analytic treatment only in idealized form. In this paper, we explore the evolution of compressible two-dimensional flows using an implicit dual-time hydrodynamical scheme that strictly conserves vorticity (if applied to simulate inviscid flows for which Kelvin's Circulation Theorem is applicable). The algorithm is based on the work of Lerat, Falissard & Side (2007), who proposed it in the context of terrestrial applications such as the blade-vortex interactions generated by helicopter rotors. We present several tests of Lerat et al.'s vorticity-preserving approach, which we have implemented to second-order accuracy, providing side-by-side comparisons with other algorithms that are frequently used in protostellar disk simulations. The comparison codes include one based on explicit, second-order van-Leer advection, one based on spectral methods, and another that implements a higher-order Godunov solver. Our results suggest that Lerat et al's algorithm will be useful for simulations of astrophysical environments in which vortices play a dynamical role, and where strong shocks are not expected.

The optical, full-phase photometric variations of a short-period planet provide a unique view of the planet's atmospheric composition and dynamics. The number of planets with optical phase curve detections, however, is currently too small to study them as an aggregate population, motivating an extension of the search to non-transiting planets. Here we present an algorithm for the detection of non-transiting, short-period giant planets in the Kepler field. The procedure uses the phase curves themselves as evidence for the planets' existence. We employ a supervised learning algorithm to recognize the salient time-dependent properties of synthetic phase curves; we then search for detections of signals that match these properties. After demonstrating the algorithm's capabilities, we classify 142,630 FGK Kepler stars without confirmed planets or KOIs and, for each one, assign a probability of a phase curve of a non-transiting planet being present. We identify 60 high-probability non-transiting hot Jupiter candidates. We also derive constraints on the candidates' albedos and offsets of the phase curve maxima. These targets are strong candidates for follow-up radial velocity confirmation and characterization. Once confirmed, the atmospheric information content in the phase curves may be studied in yet greater detail.

We present ten $R$-band photometric observations of eight different transits of the hot Jupiter HAT-P-33b, which has been targeted by our Transiting Exoplanet Monitoring Project (TEMP). The data were obtained by two telescopes at the Xinglong Station of National Astronomical Observatories of China (NAOC) from 2013 December through 2016 January, and exhibit photometric scatter of $1.6-3.0\,\rm{mmag}$. After jointly analyzing the previously published photometric data, radial-velocity (RV) measurements, and our new light curves, we revisit the system parameters and orbital ephemeris for the HAT-P-33b system. Our results are consistent with the published values except for the planet-to-star radius ratio ($R_{P}/R_{*}$), the ingress/egress duration ($\tau$) and the total duration ($T_{14}$), which together indicate a slightly shallower and shorter transit shape. Our results are based on more complete light curves, whereas the previously published work had only one complete transit light curve. No significant anomalies in Transit Timing Variations (TTVs) are found, and we place upper mass limits on potential perturbers, largely supplanting the loose constraints provided by the extant RV data. The TTV limits are stronger near mean-motion resonances, especially for the low-order commensurabilities. We can exclude the existence of a perturber with mass larger than 0.6, 0.3, 0.5, 0.5, and $0.3\,{\rm M_\oplus}$ near the 1:3, 1:2, 2:3, 3:2, and 2:1 resonances, respectively.

The newly detected TRAPPIST-1 system, with seven low-mass, roughly Earth-sized planets transiting a nearby ultra-cool dwarf, is one of the most important exoplanet discoveries to date. The short baseline of the available discovery observations, however, means that the planetary masses (obtained through measurement of transit timing variations of the planets of the system) are not yet well constrained. The masses reported in the discovery paper were derived using a combination of photometric timing measurements obtained from the ground and from the Spitzer spacecraft, and have uncertainties ranging from 30\% to nearly 100\%, with the mass of the outermost, $P=18.8\,{\rm d}$, planet h remaining unmeasured. Here, we present an analysis that supplements the timing measurements of the discovery paper with 73.6 days of photometry obtained by the K2 Mission. Our analysis refines the orbital parameters for all of the planets in the system. We substantially improve the upper bounds on eccentricity for inner six planets (finding $e<0.02$ for inner six known members of the system), and we derive masses of $0.79\pm 0.27 \,M_{\oplus}$, $1.63\pm 0.63\,M_{\oplus}$, $0.33\pm 0.15\,M_{\oplus}$, $0.24^{+0.56}_{-0.24}\,M_{\oplus}$, $0.36\pm 0.12\,M_{\oplus}$, $0.566\pm 0.038\,M_{\oplus}$, and $0.086\pm 0.084\,M_{\oplus}$ for planets b, c, d, e, f, g, and h, respectively.

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

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.

A number of authors have proposed that the statistically significant orbital alignment of the most distant Kuiper Belt Objects (KBOs) is evidence of an as-yet undetected planet in the outer solar system, now referred to colloquially a "Planet Nine". Dynamical simulations by Batygin & Brown (2016) have provided constraints on the range of the planet's possible orbits and sky locations. We extend these investigations by exploring the suggestion of Malhotra et al. (2016) that Planet Nine is in small integer ratio mean-motion resonances (MMRs) with several of the most distant KBOs. We show that the observed KBO semi-major axes present a set of commensurabilities with an unseen planet at $\sim 654 \ \mathrm{AU}$ ($P \sim 16,725 \ \mathrm{yr}$) that has a greater than $98\%$ chance of stemming from a sequence of MMRs rather than from a random distribution. We describe and implement a Monte-Carlo optimization scheme that drives billion-year dynamical integrations of the outer solar system to pinpoint the orbital properties of perturbers that are capable of maintaining the KBOs' apsidal alignment. This optimization exercise suggests that the unseen planet is most consistently represented with mass, $m \sim 6-12 M_{\oplus}$, semi-major axis, $a \sim 654 \ \mathrm{AU}$, eccentricity, $e\sim0.45$, inclination, $i \sim 30^{\circ}$, argument of periastron, $\omega \sim 150^{\circ}$, longitude of ascending node, $\Omega \sim 50^{\circ}$, and mean anomaly, $M \sim 180^{\circ}$. A range of sky locations relative to this fiducial ephemeris are possible. We find that the region $30^{\circ}\, \lesssim \mathrm{RA} \lesssim 50^{\circ}$, $-20^{\circ}\, \lesssim \mathrm{Dec} \lesssim 20^{\circ}$ is promising.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Hot Jupiters, giant extrasolar planets with orbital periods shorter than ~10 days, have long been thought to form at large radial distances, only to subsequently experience long-range inward migration. Here, we propose that in contrast with this picture, a substantial fraction of the hot Jupiter population formed in situ via the core accretion process. We show that under conditions appropriate to the inner regions of protoplanetary disks, rapid gas accretion can be initiated by Super-Earth type planets, comprising 10-20 Earth masses of refractory composition material. An in situ formation scenario leads to testable consequences, including the expectation that hot Jupiters should frequently be accompanied by additional low-mass planets with periods shorter than ~100 days. Our calculations further demonstrate that dynamical interactions during the early stages of planetary systems' lifetimes should increase the inclinations of such companions, rendering transits rare. High-precision radial velocity monitoring provides the best prospect for their detection.

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

Two decades after the discovery of 51 Peg b, the formation processes and atmospheres of short-period gas giants remain poorly understood. Observations of eccentric systems provide key insights on those topics as they can illuminate how a planet's atmosphere responds to changes in incident flux. We report here the analysis of multi-day multi-channel photometry of the eccentric (e~ 0.93) hot Jupiter HD 80606 b obtained with the Spitzer Space Telescope. The planet's extreme eccentricity combined with the long coverage and exquisite precision of new periastron-passage observations allow us to break the degeneracy between the radiative and dynamical timescales of HD 80606 b's atmosphere and constrain its global thermal response. Our analysis reveals that the atmospheric layers probed heat rapidly (~4 hr radiative timescale) from $\lt$500 to 1400 K as they absorb ~ 20% of the incoming stellar flux during the periastron passage, while the planet's rotation period is 93$\pm_{35}^{85}$ hr, which exceeds the predicted pseudo-synchronous period (40 hr).

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $\mu$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $\mu$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $\mu$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $\mu$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1\sigma$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.

We outline a photometric method for detecting the presence of a non-transiting short-period giant planet in a planetary system harboring one or more longer period transiting planets. Within a prospective system of the type that we consider, a hot Jupiter on an interior orbit inclined to the line-of-sight signals its presence through approximately sinusoidal full-phase photometric variations in the stellar light curve, correlated with astrometrically induced transit timing variations for exterior transiting planets. Systems containing a hot Jupiter along with a low-mass outer planet or planets on inclined orbits are a predicted hallmark of in situ accretion for hot Jupiters, and their presence can thus be used to test planetary formation theories. We outline the prospects for detecting non-transiting hot Jupiters using photometric data from typical \textit{Kepler} objects of interest (KOIs). As a demonstration of the technique, we perform a brief assessment of \textit{Kepler} candidates and identify a potential non-transiting hot Jupiter in the KOI-1822 system. Candidate non-transiting hot Jupiters can be readily confirmed with a small number of Doppler velocity observations, even for stars with $V\gtrsim14$.

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

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.

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.

We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 \sigma$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 \sigma$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We are still in the early days of exoplanet discovery. Astronomers are beginning to model the atmospheres and interiors of exoplanets and have developed a deeper understanding of processes of planet formation and evolution. However, we have yet to map out the full complexity of multi-planet architectures or to detect Earth analogues around nearby stars. Reaching these ambitious goals will require further improvements in instrumentation and new analysis tools. In this chapter, we provide an overview of five observational techniques that are currently employed in the detection of exoplanets: optical and IR Doppler measurements, transit photometry, direct imaging, microlensing, and astrometry. We provide a basic description of how each of these techniques works and discuss forefront developments that will result in new discoveries. We also highlight the observational limitations and synergies of each method and their connections to future space missions.

The statistics of extrasolar planetary systems indicate that the default mode of planet formation generates planets with orbital periods shorter than 100 days, and masses substantially exceeding that of the Earth. When viewed in this context, the Solar System is unusual. Here, we present simulations which show that a popular formation scenario for Jupiter and Saturn, in which Jupiter migrates inward from a > 5 AU to a ~ 1.5 AU before reversing direction, can explain the low overall mass of the Solar System's terrestrial planets, as well as the absence of planets with a < 0.4 AU. Jupiter's inward migration entrained s ~ 10-100 km planetesimals into low-order mean-motion resonances, shepherding and exciting their orbits. The resulting collisional cascade generated a planetesimal disk that, evolving under gas drag, would have driven any pre-existing short-period planets into the Sun. In this scenario, the Solar System's terrestrial planets formed from gas-starved mass-depleted debris that remained after the primary period of dynamical evolution.

Thousands of extrasolar planets have been discovered, and it is clear that the galactic planetary census draws on a diversity greatly exceeding that exhibited by the solar system's planets. We review significant landmarks in the chronology of extrasolar planet detection, and we give an overview of the varied observational techniques that are brought to bear. We then discuss the properties of the currently known distribution, using the mass-period diagram as a guide to delineating hot Jupiters, eccentric giant planets, and a third, highly populous, category that we term "ungiants", planets having masses less than 30 Earth masses and orbital periods less than 100 days. We then move to a discussion of the bulk compositions of the extrasolar planets. We discuss the long-standing problem of radius anomalies among giant planets, as well as issues posed by the unexpectedly large range in sizes observed for planets with masses somewhat greater than Earth's. We discuss the use of transit observations to probe the atmospheres of extrasolar planets; various measurements taken during primary transit, secondary eclipse, and through the full orbital period, can give clues to the atmospheric compositions, structures, and meteorologies. The extrasolar planet catalog, along with the details of our solar system and observations of star-forming regions and protoplanetary disks, provide a backdrop for a discussion of planet formation in which we review the elements of the favored pictures for how the terrestrial and giant planets were assembled. We conclude by listing several research questions that are relevant to the next ten years and beyond.

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$\mu$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $\mu$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $\mu$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

In contrast to photometric transits, whose peak signal occurs at mid-transit due to occultation of the brightest region of the disk, polarimetric transits provide a signal upon ingress and egress due to occultation of the polarized stellar limb. Limb polarization, the bright corollary to limb darkening, arises from the $90^\circ$ scattering angle and low optical depth experienced by photons at the limb. In addition to the ratio $R_{\rm p} / R_*$, the amplitude of a polarimetric transit is expected to be controlled by the strength and width of the stellar limb polarization profile, which depend on the scattering-to-total opacity ratio at the stellar limb. We present a short list of the systems providing the highest expected signal-to-noise ratio for detection of this effect, and we draw particular attention to HD 80606b. This planet is spin/orbit misaligned, has a three-hour ingress, and has a bright parent star, which make it an attractive target. We report on test observations of an HD 80606b ingress with the POLISH2 polarimeter at the Lick Observatory Shane 3-m telescope. We conclude that unmodeled telescope systematic effects prevented polarimetric detection of this event. We outline a roadmap for further refinements of exoplanet polarimetry, whose eventual success will require a further factor of ten reduction in systematic noise.

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.

Unlike hot Jupiters or other gas giants, super-Earths are expected to have a wide variety of compositions, ranging from terrestrial bodies like our own to more gaseous planets like Neptune. Observations of transiting systems, which allow us to directly measure planet masses and radii and constrain atmospheric properties, are key to understanding the compositional diversity of the planets in this mass range. Although Kepler has discovered hundreds of transiting super-Earth candidates over the past four years, the majority of these planets orbit stars that are too far away and too faint to allow for detailed atmospheric characterization and reliable mass estimates. Ground-based transit surveys focus on much brighter stars, but most lack the sensitivity to detect planets in this size range. One way to get around the difficulty of finding these smaller planets in transit is to start by choosing targets that are already known to contain super-Earth sized bodies detected using the radial velocity technique. Here we present results from a Spitzer program to observe six of the most favorable RV detected super-Earth systems, including HD 1461, HD 7924, HD 156668, HIP 57274, and GJ 876. We find no evidence for transits in any of their 4.5 micron flux light curves, and place limits on the allowed transit depths and corresponding planet radii that rule out even the most dense and iron-rich compositions for these objects. We also observed HD 97658, but the observation window was based on a possible ground-based transit detection (Henry et al. 2011) that was later ruled out; thus the window did not include the predicted time for the transit detection recently made by MOST (Dragomir et al. 2013).

We propose to use Kepler in 2-wheel mode to conduct a detailed search for Earth-sized planets orbiting ultra-cool stars and brown dwarfs (spectral types from M7 to L3). This population of objects presents several advantages for exoplanet surveys. First, ultra-cool stars and brown dwarfs are small and thus result in favorable planet-to-star area ratios. Second, because of their low effective temperature, the inner edge of their habitable zone is extremely close (2 to 3 days only). Third, our targets are bright at infrared wavelengths, which will enable detailed follow-up studies. Our program therefore represents a unique opportunity to find a transiting Earth-size exoplanet for which atmospheric features (including biosignatures) could be detected with near-to-come facilities such as JWST. Such exoplanet has not been discovered yet. Kepler in 2-wheel mode provides the required stability and photometric precision to make this survey successful. Our initial target sample includes 60 ultra-cool stars and brown dwarfs from which we expect to detect at least one transiting planet. We propose to monitor each source for 4 days, resulting in a total program duration of ~240 days.

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$\sigma$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

Exoplanetary science has reached a historic moment. The James Webb Space Telescope will be capable of probing the atmospheres of rocky planets, and perhaps even search for biologically produced gases. However this is contingent on identifying suitable targets before the end of the mission. A race therefore, is on, to find transiting planets with the most favorable properties, in time for the launch. Here, we describe a realistic opportunity to discover extremely favorable targets - rocky planets transiting nearby brown dwarfs - using the Spitzer Space Telescope as a survey instrument. Harnessing the continuous time coverage and the exquisite precision of Spitzer in a 5,400 hour campaign monitoring nearby brown dwarfs, we will detect a handful of planetary systems with planets as small as Mars. The survey we envision is a logical extension of the immense progress that has been realized in the field of exoplanets and a natural outcome of the exploration of the solar neighborhood to map where the nearest habitable rocky planets are located (as advocated by the 2010 Decadal Survey). Our program represents an essential step towards the atmospheric characterization of terrestrial planets and carries the compelling promise of studying the concept of habitability beyond Earth-like conditions. In addition, our photometric monitoring will provide invaluable observations of a large sample of nearby brown dwarfs situated close to the M/L transition. This is why, we also advocate an immediate public release of the survey data, to guarantee rapid progress on the planet search and provide a treasure trove of data for brown dwarf science.

The transit signature of exoplanets provides an avenue through which characterization of exoplanetary properties may be undertaken, such as studies of mean density, structure, and atmospheric composition. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is a program to expand the catalog of transiting planets around bright host stars by refining the orbits of known planets discovered with the radial velocity technique. Here we present results for the HD 38529 system. We determine fundamental properties of the host star through direct interferometric measurements of the radius and through spectroscopic analysis. We provide new radial velocity measurements that are used to improve the Keplerian solution for the two known planets, and we find no evidence for a previously postulated third planet. We also present 12 years of precision robotic photometry of HD 38529 that demonstrate the inner planet does not transit and the host star exhibits cyclic variations in seasonal mean brightness with a timescale of approximately six years.

Close-in super-Earths, with radii R = 2-5 R_Earth and orbital periods P < 100 days, orbit more than half, and perhaps nearly all Sun-like stars in the universe. We use this omnipresent population to construct the minimum-mass extrasolar nebula (MMEN), the circumstellar disk of solar-composition solids and gas from which such planets formed, if they formed near their current locations and did not migrate. In a series of back-of-the-envelope calculations, we demonstrate how in-situ formation in the MMEN is fast, efficient, and can reproduce many of the observed properties of close-in super-Earths, including their gas-to-rock fractions. Testable predictions are discussed.

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 \mu m bands of the Spitzer Space Telescope. The 3.6 and 4.5 \mu m data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 \mu m that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 \mu m bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

As part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS), we present new radial velocities and photometry of the HD 192263 system. Our analysis of the already available Keck-HIRES and CORALIE radial velocity measurements together with the five new Keck measurements we report in this paper results in improved orbital parameters for the system. We derive constraints on the size and phase location of the transit window for HD 192263b, a Jupiter-mass planet with a period of 24.3587 \pm 0.0022 days. We use 10 years of Automated Photoelectric Telescope (APT) photometry to analyze the stellar variability and search for planetary transits. We find continuing evidence of spot activity with periods near 23.4 days. The shape of the corresponding photometric variations changes over time, giving rise to not one but several Fourier peaks near this value. However, none of these frequencies coincides with the planet's orbital period and thus we find no evidence of star-planet interactions in the system. We attribute the ~23-day variability to stellar rotation. There are also indications of spot variations on longer (8 years) timescales. Finally, we use the photometric data to exclude transits for a planet with the predicted radius of 1.09 RJ, and as small as 0.79 RJ.

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6{\mu}m and 4.5{\mu}m for each target. XO-4b exhibits a stronger eclipse depth at 4.5{\mu}m than at 3.6{\mu}m, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6{\mu}m, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(\omega) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

The discovery of transiting planets around bright stars holds the potential to greatly enhance our understanding of planetary atmospheres. In this work we present the search for transits of HD 168443b, a massive planet orbiting the bright star HD 168443 V=6.92 with a period of 58.11 days. The high eccentricity of the planetary orbit e=0.53 significantly enhances the a-priori transit probability beyond that expected for a circular orbit, making HD 168443 a candidate for our ongoing Transit Ephemeris Refinement and Monitoring Survey (TERMS). Using additional radial velocities from Keck-HIRES, we refined the orbital parameters of this multi-planet system and derived a new transit ephemeris for HD 168443b. The reduced uncertainties in the transit window make a photometric transit search practicable. Photometric observations acquired during predicted transit windows were obtained on three nights. CTIO 1.0 m photometry acquired on 2010 September 7 had the required precision to detect a transit but fell just outside of our final transit window. Nightly photometry from the T8 0.8 m Automated Photometric Telescope (APT) at Fairborn Observatory, acquired over a span of 109 nights, demonstrates that HD 168443 is constant on a time scale of weeks. Higher-cadence photometry on 2011 April 28 and June 25 shows no evidence of a transit. We are able to rule out a non-grazing transit of HD 168443b.

The ongoing High Accuracy Radial velocity Planet Search (HARPS) has found that 30-50% of GK dwarfs in the solar neighborhood host planets with sub-Neptune masses in orbits of P < 50 days. At first glance, this overall occurrence rate seems inconsistent with the planet frequency measured during Q0-Q2 of the Kepler Mission, whose 1,235 detected planetary candidates imply that ~ 15% of main sequence dwarfs harbor short-period planets with R_pl < 4 R_Earth. A rigorous comparison between the two surveys is difficult, however, as they observe different stellar populations and measure different planetary properties. Here we report the results of a Monte Carlo study that can account for this discrepancy via plausible distributions of planetary compositions. We find that a population concurrently consisting of (1) dense silicate-iron planets and (2) low-density gas-dominated worlds provides a natural fit to the current data. In this scenario, the fraction of dense planets decreases with increasing mass, from f_rocky = 90% at M = 1 M_Earth to f_rocky = 10% at M = 17 M_Earth. Our best fit population has a total occurrence rate of 40% for 2 < P < 50 days and 1 < M < 17 M_Earth, and is characterized by simple power-law indices of the form N(M)dM ~ M^alpha dM and N(P)dP ~ P^beta dP with alpha = -1.0 and beta = 0.0. Our model population therefore contains four free parameters and is readily testable with future observations. Furthermore, our model's insistence that at least two distinct types of planets must exist in the survey data indicates that multiple formation mechanisms are at work to produce the population of planets commonly referred to as "super-Earths".

We find that Kepler exoplanet candidate (EC) host stars are preferentially metal-rich, including the low-mass stellar hosts of small-radius ECs. The last observation confirms a tentative hint that there is a correlation between the metallicity of low-mass stars and the presence of low-mass and small-radius exoplanets. In particular, we compare the J-H--g-r color-color distribution of Kepler EC host stars with a control sample of dwarf stars selected from the ~150,000 stars observed during Q1 and Q2 of the Kepler mission but with no detected planets. We find that at J-H = 0.30 characteristic of solar-type stars, the average g-r color of stars that host giant ECs is 4-sigma redder than the average color of the stars in the control sample. At the same time, the average g-r color of solar-type stars that host small-radius ECs is indistinguishable from the average color of the stars in the control sample. In addition, we find that at J-H = 0.62 indicative of late K dwarfs, the average g-r color of stars that host small-radius ECs is 4-sigma redder than the average color of the stars in the control sample. These offsets are unlikely to be caused by differential reddening, age differences between the two populations, or the presence of giant stars in the control sample. Stellar models suggest that the first color offset is due to a 0.2 dex enhancement in [Fe/H] of the giant EC host population at M_star = 1 M_Sun, while Sloan photometry of M 67 and NGC 6791 suggests that the second color offset is due to a similar [Fe/H] enhancement of the small-radius EC host population at M_star = 0.7 M_Sun. These correlations are a natural consequence of the core-accretion model of planet formation.

In this paper we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 micron bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 Feb 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7 sigma significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight, and that the planet's orbit is therefore likely to be misaligned. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. Our measured 8 micron secondary eclipse depths are consistent with a constant value, and we place a 1 sigma upper limit of 17% on changes in the planet's dayside flux in this band. Averaging over the eleven visits gives us an improved estimate of 0.0452% +/- 0.0027% for the secondary eclipse depth. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris, and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. [ABRIDGED]

HD 156846b is a Jovian planet in a highly eccentric orbit (e = 0.85) with a period of 359.55 days. The pericenter passage at a distance of 0.16 AU is nearly aligned to our line of sight, offering an enhanced transit probability of 5.4% and a potentially rich probe of the dynamics of a cool planetary atmosphere impulsively heated during close approach to a bright star (V = 6.5). We present new radial velocity (RV) and photometric measurements of this star as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS). The RV measurements from Keck-HIRES reduce the predicted transit time uncertainty to 20 minutes, an order of magnitude improvement over the ephemeris from the discovery paper. We photometrically monitored a predicted transit window under relatively poor photometric conditions, from which our non-detection does not rule out a transiting geometry. We also present photometry that demonstrates stability at the millimag level over its rotational timescale.

Long-term orbital evolution of multi-planet systems under tidal dissipation often converges to a stationary state, known as the tidal fixed point. The fixed point is characterized by a lack of oscillations in the eccentricities and apsidal alignment among the orbits. Quantitatively, the nature of the fixed point is dictated by mutual interactions among the planets as well as non-Keplerian effects. We show that if a roughly coplanar system hosts a hot, sub-Saturn mass planet, and is tidally relaxed, separation of planet-planet interactions and non-Keplerian effects in the equations of motion leads to a direct determination of the true masses of the planets. Consequently, a "snap-shot" observational determination of the orbital state resolves the sin(I) degeneracy, and opens up a direct avenue towards identification of the true lowest-mass exo-planets detected. We present an approximate, as well as a general, mathematical framework for computation of the line of sight inclination of secular systems, and apply our models illustratively to the 61 Vir system. We conclude by discussing the observability of planetary systems to which our method is applicable and we set our analysis into a broader context by presenting a current summary of the various possibilities for determining the physical properties of planets from observations of their orbital states.

We present a systematic evaluation of the agreement between the observed radii of 90 well-characterized transiting extrasolar giant planets and their corresponding model radii. Our model radii are drawn from previously published calculations of core-less giant planets that have attained their asymptotic radii, and which have been tabulated for a range of planet masses and equilibrium temperatures. (We report a two-dimensional polynomial fitting function that accurately represents the models). As expected, the model radii provide a statistically significant improvement over a null hypothesis that the sizes of giant planets are completely independent of mass and effective temperature. As is well known, however, fiducial models provide an insufficient explanation; the planetary radius anomalies are strongly correlated with planetary equilibrium temperature. We find that the radius anomalies have a best-fit dependence, ${\cal R}\propto T_{\rm eff}^{\alpha}$, with $\alpha=1.4\pm0.6$. Incorporating this relation into the model radii leads to substantially less scatter in the radius correlation. The extra temperature dependence represents an important constraint on theoretical models for Hot Jupiters. Using simple scaling arguments, we find support for the hypothesis of Batygin and Stevenson (2010) that this correlation can be attributed to a planetary heating mechanism that is mediated by magnetohydrodynamic coupling between the planetary magnetic field and near-surface flow that is accompanied by ohmic dissipation at adiabatic depth. Additionally, we find that the temperature dependence is likely too strong to admit kinetic heating as the primary source of anomalous energy generation within the majority of the observed transiting planets.

Cool M dwarfs outnumber sun-like G stars by ten to one in the solar neighborhood. Due to their proximity, small size, and low mass, M-dwarf stars are becoming attractive targets for exoplanet searches via almost all current search methods. But what planetary systems can form around M dwarfs? Following up on the Cool Stars~16 Splinter Session "Planet Formation Around M Dwarfs", we summarize here our knowledge of protoplanetary disks around cool stars, how they disperse, what planetary systems might form and can be detected with current and future instruments.

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 photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

Transiting exoplanets provide access to data to study the mass-radius relation and internal structure of extrasolar planets. Long-period transiting planets allow insight into planetary environments similar to the Solar System where, in contrast to hot Jupiters, planets are not constantly exposed to the intense radiation of their parent stars. Observations of secondary eclipses additionally permit studies of exoplanet temperatures and large-scale exo-atmospheric properties. We show how transit and eclipse probabilities are related to planet-star system geometries, particularly for long-period, eccentric orbits. The resulting target selection and observational strategies represent the principal ingredients of our photometric survey of known radial-velocity planets with the aim of detecting transit signatures (TERMS).

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

Transiting planet discoveries have largely been restricted to the short-period or low-periastron distance regimes due to the bias inherent in the geometric transit probability. Through the refinement of planetary orbital parameters, and hence reducing the size of transit windows, long-period planets become feasible targets for photometric follow-up. Here we describe the TERMS project that is monitoring these host stars at predicted transit times.

The search for a habitable extrasolar planet has long interested scientists, but only recently have the tools become available to search for such planets. In the past decades, the number of known extrasolar planets has ballooned into the hundreds, and with it the expectation that the discovery of the first Earth-like extrasolar planet is not far off. Here we develop a novel metric of habitability for discovered planets, and use this to arrive at a prediction for when the first habitable planet will be discovered. Using a bootstrap analysis of currently discovered exoplanets, we predict the discovery of the first Earth-like planet to be announced in the first half of 2011, with the likeliest date being early May 2011. Our predictions, using only the properties of previously discovered exoplanets, accord well with external estimates for the discovery of the first potentially habitable extrasolar planet, and highlights the the usefulness of predictive scientometric techniques to understand the pace of scientific discovery in many fields.

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

The location of M dwarfs in the V-K_s--M_Ks color-magnitude diagram (CMD) has been shown to correlate with metallicity. We demonstrate that previous empirical photometric calibrations of M dwarf metallicity exploiting this correlation systematically underestimate or overestimate metallicity at the extremes of their range. We improve upon previous calibrations in three ways. We use both a volume-limited and kinematically-matched sample of F and G dwarfs from the Geneva-Copehnagen Survey (GCS) to infer the mean metallicity of M dwarfs in the Solar Neighborhood, we use theoretical models of M dwarf interiors and atmospheres to determine the effect of metallicity on M dwarfs in the V-K_s--M_Ks CMD, and we base our final calibration purely on high-resolution spectroscopy of FGK primaries with M dwarf companions. As a result, we explain an order of magnitude more of the variance in the calibration sample than previous photometric calibrations. We non-parametrically quantify the significance of the observation that M dwarfs that host exoplanets are preferentially in a region of the V-K_s--M_Ks plane populated by metal-rich M dwarfs. We find that the probability p that planet-hosting M dwarfs are distributed across the V-K_s--M_Ks CMD in the same way as field M dwarfs is p = 0.06 +/- 0.008. Interestingly, the subsample of M dwarfs that host Neptune and sub-Neptune mass planets may also be preferentially located in the region of the V-K_s--M_Ks plane populated by high-metallicity M dwarfs. The probability of this occurrence by chance is p = 0.40 +/- 0.02, and this observation hints that low-mass planets may be more likely to be found around metal-rich M dwarfs. An increased rate of low-mass planet occurrence around metal-rich M dwarfs would be a natural consequence of the core-accretion model of planet formation. (abridged)

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.

This is a scientific strategy for the detection and characterization of extrasolar planets; that is, planets orbiting other stars. As such, it maps out over a 15-year horizon the techniques and capabilities required to detect and measure the properties of planets as small as Earth around stars as large as our own Sun. It shows how the technology pieces and their development fit together to achieve the primary goal of the strategy: if planets like Earth exist around stars within some tens of light years of our own Solar System, those planets will be found and their basic properties characterized. Essential to this strategy is not only the search for and examination of individual planets, but also a knowledge of the arrangement, or architecture, of planetary systems around as large a number of stars as possible; this is the second goal of the strategy. The final goal of the strategy is the study of disks around stars, important both to understand the implications of the variety of exoplanet systems for planet formation, and to determine how many nearby stars have environments around them clean enough of debris that planets may be sought and, if found, characterized.

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.

Transiting planet discoveries have yielded a plethora of information towards understanding the structure and atmospheres of extra-solar planets. These discoveries have been restricted to the short-period or low-periastron distance regimes due to the bias inherent in the geometric transit probability. Through the refinement of planetary orbital parameters, and hence reducing the size of transit windows, long-period planets become feasible targets for photometric follow-up. Here we describe the TERMS project which is monitoring these host stars at predicted transit times.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

Transiting planet discoveries have yielded a plethora of information regarding the internal structure and atmospheres of extra-solar planets. These discoveries have been restricted to the low-periastron distance regime due to the bias inherent in the geometric transit probability. Monitoring known radial velocity planets at predicted transit times is a proven method of detecting transits, and presents an avenue through which to explore the mass-radius relationship of exoplanets in new regions of period/periastron space. Here we describe transit window calculations for known radial velocity planets, techniques for refining their transit ephemerides, target selection criteria, and observational methods for obtaining maximum coverage of transit windows. These methods are currently being implemented by the Transit Ephemeris Refinement and Monitoring Survey (TERMS).

Doppler and transit surveys are finding extrasolar planets of ever smaller mass and radius, and are now sampling the domain of superEarths (1-3 Earth radii). Recent results from the Doppler surveys suggest that discovery of a transiting superEarth in the habitable zone of a lower main sequence star may be possible. We evaluate the prospects for an all-sky transit survey targeted to the brightest stars, that would find the most favorable cases for photometric and spectroscopic characterization using the James Webb Space Telescope (JWST). We use the proposed Transiting Exoplanet Survey Satellite (TESS) as representative of an all-sky survey. We couple the simulated TESS yield to a sensitivity model for the MIRI and NIRSpec instruments on JWST. We focus on the TESS planets with radii between Earth and Neptune. Our simulations consider secondary eclipse filter photometry using JWST/MIRI, comparing the 11- and 15-micron bands to measure CO2 absorption in superEarths, as well as JWST/NIRSpec spectroscopy of water absorption from 1.7-3.0 microns, and CO2 absorption at 4.3-microns. We project that TESS will discover about eight nearby habitable transiting superEarths. The principal sources of uncertainty in the prospects for JWST characterization of habitable superEarths are superEarth frequency and the nature of superEarth atmospheres. Based on our estimates of these uncertainties, we project that JWST will be able to measure the temperature, and identify molecular absorptions (water, CO2) in one to four nearby habitable TESS superEarths.

Tidal dissipation within a short-period transiting extrasolar planet perturbed by a companion object can drive orbital evolution of the system to a so-called tidal fixed point, in which the apsidal lines of the transiting planet and its perturber are aligned, and for which variations in the orbital eccentricities of both planet and perturber are damped out. Significant contributions to the apsidal precession rate are made by the secular planet-planet interaction, by general relativity, and by the gravitational quadropole fields created by the transiting planet's tidal and rotational distortions. The fixed-point orbital eccentricity of the inner planet is therefore a strong function of the planet's interior structure. We illustrate these ideas in the specific context of the recently discovered HAT-P-13 exo-planetary system, and show that one can already glean important insights into the physical properties of the inner transiting planet. We present structural models of the planet, which indicate that its observed radius can be maintained for a one-parameter sequence of models that properly vary core mass and tidal energy dissipation in the interior. We use an octopole-order secular theory of the orbital dynamics to derive the dependence of the inner planet's eccentricity, on its tidal Love number. We find that the currently measured eccentricity, implies 0.116 < k2_{b} < 0.425, 0 M_{Earth}<M_{core}<120 M_{Earth}$, and Q_{b} < 300,000. Improved measurement of the eccentricity will soon allow for far tighter limits to be placed on all three of these quantities, and will provide an unprecedented probe into the interior structure of an extrasolar planet.

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.

We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

We investigate an in-situ formation scenario for Earth-mass terrestrial planets in short-period, potentially habitable orbits around low-mass stars (M_star < 0.3 M_sun). We then investigate the feasibility of detecting these Earth-sized planets. Our simulations of terrestrial planet formation follow the growth of planetary embryos in an annular region around a fiducial M7 primary. Our simulations couple a semi-analytic model to a full N-body integration to follow the growth from ~3x10^21 g to the final planetary system configurations that generally consist of 3-5 planets with masses of order 0.1 - 1.0 M_earth in or near the habitable zone of the star. To obtain a concrete estimate of the detectability of the planets arising in our simulations, we present a detailed Monte-Carlo transit detection simulation. We find that detection of 1 R_earth planets around the local M-dwarfs is challenging for a 1m class ground-based photometric search, but that detection of planets of larger radius is a distinct possibility. The detection of Earth-sized planets is straightforward, however, with an all-sky survey by a low-cost satellite mission. Given a reduced correlated noise level of 0.45 mmag and an intermediate planetary ice-mass fraction of planets orbiting a target list drawn from the nearest late-type M dwarfs, a ground-based photometric search could detect, on average, 0.8 of these planets with an extended search. A space-based photometric search (similar to the TESS mission) should discover ~17 of these Earth-sized planets during it's two year survey, with an assumed occurrence fraction of 28%.

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.

AIMS : To improve the parameters of the HD 17156 system (peculiar due to the eccentric and long orbital period of its transiting planet) and constrain the presence of stellar companions. METHODS : Photometric data were acquired for 4 transits, and high precision radial velocity measurements were simultaneously acquired with SARG@TNG for one transit. The template spectra of HD 17156 was used to derive effective temperature, gravity, and metallicity. A fit of the photometric and spectroscopic data was performed to measure the stellar and planetary radii, and the spin-orbit alignment. Planet orbital elements and ephemeris were derived from the fit. Near infrared adaptive optic images was acquired with ADOPT@TNG. RESULTS: We have found that the star has a radius of R_S = 1.43+/-0.03 R_sun and the planet R_P =1.02+/-0.08 R_jup. The transit ephemeris is T_c = 2454\756.73134+/-0.00020+N*21.21663+/-0.00045 BJD. The analysis of the Rossiter-Mclaughlin effect shows that the system is spin orbit aligned with an angle lambda = 4.8 +/- 5.3 deg. The analysis of high resolution images has not revealed any stellar companion with projected separation between 150 and 1000 AU from HD 17156.

We present two new in situ core accretion simulations of Saturn with planet formation timescales of 3.37 Myr (model S0) and 3.48 Myr (model S1), consistent with observed protostellar disk lifetimes. In model S0, we assume rapid grain settling reduces opacity due to grains from full interstellar values (Podolak 2003). In model S1, we do not invoke grain settling, instead assigning full interstellar opacities to grains in the envelope. Surprisingly, the two models produce nearly identical formation timescales and core/atmosphere mass ratios. We therefore observe a new manifestation of core accretion theory: at large heliocentric distances, the solid core growth rate (limited by Keplerian orbital velocity) controls the planet formation timescale. We argue that this paradigm should apply to Uranus and Neptune as well.

We have developed a new two-dimensional hydrostatically-balanced isobaric hydrodynamic model for use in simulation of exoplanetary atmospheres. We apply this model to the infrared photosphere of the hot Jupiter HD 189733 b, for which an excellent 8-micron light curve has been obtained. For reasonable parameter choices, the results of our model are consistent with these observations. In our simulations, strongly turbulent supersonic flow develops, with wind speeds of approximately 5 km/s. This flow geometry causes chaotic variation of the temperature distribution, leading to observable variations in the light curve from one orbit to the next.

This paper considers the effects of turbulence on mean motion resonances in extrasolar planetary systems and predicts that systems rarely survive in a resonant configuration. A growing number of systems are reported to be in resonance, which is thought to arise from the planet migration process. If planets are brought together and moved inward through torques produced by circumstellar disks, then disk turbulence can act to prevent planets from staying in a resonant configuration. This paper studies this process through numerical simulations and via analytic model equations, where both approaches include stochastic forcing terms due to turbulence. We explore how the amplitude and forcing time intervals of the turbulence affect the maintenance of mean motion resonances. If turbulence is common in circumstellar disks during the epoch of planet migration, with the amplitudes indicated by current MHD simulations, then planetary systems that remain deep in mean motion resonance are predicted to be rare. More specifically, the fraction of resonant systems that survive over a typical disk lifetime of 1 Myr is of order 0.01. If mean motion resonances are found to be common, their existence would place tight constraints on the amplitude and duty cycle of turbulent fluctuations in circumstellar disks. These results can be combined by expressing the expected fraction of surviving resonant systems in the approximate form P_b = C / N_{orb}^{1/2}, where the dimensionless parameter C = 10 - 50 and where N_{orb} is the number of orbits for which turbulence is active.

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.

A long-term numerical integration of the classical Newtonian approximation to the planetary orbital motions of the full Solar System (sun + 8 planets), spanning 20 Gyr, was performed. The results showed no severe instability arising over this time interval. Subsequently, utilizing a bifurcation method described by Jacques Laskar, two numerical experiments were performed with the goal of determining dynamically allowed evolutions for the Solar System in which the planetary orbits become unstable. The experiments yielded one evolution in which Mercury falls onto the Sun at ~1.261Gyr from now, and another in which Mercury and Venus collide in ~862Myr. In the latter solution, as a result of Mercury's unstable behavior, Mars was ejected from the Solar System at ~822Myr. We have performed a number of numerical tests that confirm these results, and indicate that they are not numerical artifacts. Using synthetic secular perturbation theory, we find that Mercury is destabilized via an entrance into a linear secular resonance with Jupiter in which their corresponding eigenfrequencies experience extended periods of commensurability. The effects of general relativity on the dynamical stability are discussed. An application of the bifurcation method to the outer Solar System (Jupiter, Saturn, Uranus, and Neptune) showed no sign of instability during the course of 24Gyr of integrations, in keeping with an expected Uranian dynamical lifetime of 10^(18) years.

NGC 205, a close satellite of the M31 galaxy, is our nearest example of a dwarf elliptical galaxy. Photometric and kinematic observations suggest that NGC 205 is undergoing tidal distortion from its interaction with M31. Despite earlier attempts, the orbit and progenitor properties of NGC 205 are not well known. We perform an optimized search for these unknowns by combining a genetic algorithm with restricted N-body simulations of the interaction. This approach, coupled with photometric and kinematic observations as constraints, allows for an effective exploration of the parameter space. We represent NGC 205 as a static Hernquist potential with embedded massless test particles that serve as tracers of surface brightness. We explore 3 distinct, initially stable configurations of test particles: cold rotating disk, warm rotating disk, and hot, pressure-supported spheroid. Each model reproduces some, but not all, of the observed features of NGC 205, leading us to speculate that a rotating progenitor with substantial pressure support could match all of the observables. Furthermore, plausible combinations of mass and scale length for the pressure-supported spheroid progenitor model reproduce the observed velocity dispersion profile. For all 3 models, orbits that best match the observables place the satellite 11+/-9 kpc behind M31 moving at very large velocities: 300-500 km/s on primarily radial orbits. Given that the observed radial component is only 54 km/s, this implies a large tangential motion for NGC 205, moving from the NW to the SE. These results suggest NGC 205 is not associated with the stellar arc observed to the NE of NGC 205. Furthermore, NGC 205's velocity appears to be near or greater than its escape velocity, signifying that the satellite is likely on its first M31 passage.

We report high-cadence time-series photometry of the recently-discovered transiting exoplanet system HD 17156, spanning the time of transit on UT 2007 October 1, from three separate observatories. We present a joint analysis of our photometry, previously published radial velocity measurements, and times of transit center for 3 additional events. Adopting the spectroscopically-determined values and uncertainties for the stellar mass and radius, we estimate a planet radius of Rp = 1.01 +/- 0.09 RJup and an inclination of i = 86.5 +1.1 -0.7 degrees. We find a time of transit center of Tc = 2454374.8338 +/- 0.0020 HJD and an orbital period of P = 21.21691 +/- 0.00071 days, and note that the 4 transits reported to date show no sign of timing variations that would indicate the presence of a third body in the system. Our results do not preclude the existence of a secondary eclipse, but imply there is only a 9.2% chance for this to be present, and an even lower probability (6.9%) that the secondary eclipse would be a non-grazing event. Due to its eccentric orbit and long period, HD 17156b is a fascinating object for the study of the dynamics of exoplanet atmospheres. To aid such future studies, we present theoretical light curves for the variable infrared emission from the visible hemisphere of the planet throughout its orbit.

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 18 years of Doppler shift measurements of a nearby star, 55 Cancri, that exhibit strong evidence for five orbiting planets. The four previously reported planets are strongly confirmed here. A fifth planet is presented, with an apparent orbital period of 260 days, placing it 0.78 AU from the star in the large empty zone between two other planets. The velocity wobble amplitude of 4.9 \ms implies a minimum planet mass \msini = 45.7 \mearthe. The orbital eccentricity is consistent with a circular orbit, but modest eccentricity solutions give similar \chisq fits. All five planets reside in low eccentricity orbits, four having eccentricities under 0.1. The outermost planet orbits 5.8 AU from the star and has a minimum mass, \msini = 3.8 \mjupe, making it more massive than the inner four planets combined. Its orbital distance is the largest for an exoplanet with a well defined orbit. The innermost planet has a semi-major axis of only 0.038 AU and has a minimum mass, \msinie, of only 10.8 \mearthe, one of the lowest mass exoplanets known. The five known planets within 6 AU define a {\em minimum mass protoplanetary nebula} to compare with the classical minimum mass solar nebula. Numerical N-body simulations show this system of five planets to be dynamically stable and show that the planets with periods of 14.65 and 44.3 d are not in a mean-motion resonance. Millimagnitude photometry during 11 years reveals no brightness variations at any of the radial velocity periods, providing support for their interpretation as planetary.

We employ a two-dimensional grid-based hydrodynamic model to simulate upper atmospheric dynamics on extrasolar giant planets. Our model is well-suited to simulate the dynamics of the atmospheres of planets with high orbital eccentricity that are subject to widely-varying irradiation conditions. We identify six such planets, with eccentricities between $e=0.28$ and $e=0.93$ and semimajor axes ranging from $a=0.0508$ A.U. to $a=0.432$ A.U., as particularly interesting objects for study. For each of these planets, we determine the temperature profile and resulting infrared light curves in the 8-$\mu$m Spitzer bands. Especially notable are the results for HD 80606b, which has the largest eccentricity ($e=0.9321$) of any known planet, and HAT-P-2b, which transits its parent star, so that its physical properties are well-constrained. Despite the variety of orbital parameters, the atmospheric dynamics of these eccentric planets display a number of interesting common properties. In all cases, the atmospheric response is primarily driven by the intense irradiation at periastron. The resulting expansion of heated air produces high-velocity turbulent flow, including long-lived circumpolar vortices. Additionally, a superrotating acoustic front develops on some planets; the strength of this disturbance depends on both the eccentricity and the temperature gradient resulting from uneven heating. The specifics of the resulting infrared light curves depend strongly on the orbital geometry. We show, however, that the variations on HD 80606 b and HAT-P-2b should be readily detectable at 4.5 and 8 $\mu$m using the Spitzer Space Telescope. Indeed, these two objects present the most attractive observational targets of all known high-$e$ exoplanets.

We report the detection of transits by the 3.1 M_Jup companion to the V=8.17 G0V star HD 17156. The transit was observed by three independant observers on Sept. 9/10, 2007 (two in central Italy and one in the Canary Islands), who obtained detections at confidence levels of 3.0 sigma, 5.3 sigma, and 7.9 sigma, respectively. The observations were carried out under the auspices of the Transitsearch.org network, which organizes follow-up photometric transit searches of known planet-bearing stars during the time intervals when transits are expected to possibly occur. Analyses of the 7.9 sigma data set indicates a transit depth d=0.0062+/-0.0004, and a transit duration t=186+/-5 min. These values are consistent with the transit of a Jupiter-sized planet with an impact parameter b=a*cos(i)/R_star ~ 0.8. This planet occupies a unique regime among known transiting extrasolar planets, both as a result of its large orbital eccentricity (e=0.67) and long orbital period (P=21.2 d). The planet receives a 26-fold variation in insolation during the course of its orbit, which will make it a useful object for characterization of exoplanetary atmospheric dynamics.

We report the detection of five Jovian mass planets orbiting high metallicity stars. Four of these stars were first observed as part of the N2K program and exhibited low RMS velocity scatter after three consecutive observations. However, follow-up observations over the last three years now reveal the presence of longer period planets with orbital periods ranging from 21 days to a few years. HD 11506 is a G0V star with a planet of \msini = 4.74 \mjup in a 3.85 year orbit. HD 17156 is a G0V star with a 3.12 \mjup planet in a 21.2 day orbit. The eccentricity of this orbit is 0.67, one of the highest known for a planet with a relatively short period. The orbital period for this planet places it in a region of parameter space where relatively few planets have been detected. HD 125612 is a G3V star with a planet of \msini = 3.5 \mjup in a 1.4 year orbit. HD 170469 is a G5IV star with a planet of \msini = 0.67 \mjup in a 3.13 year orbit. HD 231701 is an F8V star with planet of 1.08 \mjup in a 142 day orbit. All of these stars have supersolar metallicity. Three of the five stars were observed photometrically but showed no evidence of brightness variability. A transit search conducted for HD 17156 was negative but covered only 25% of the search space and so is not conclusive.

We report the discovery of two new planets: a 1.94 M_Jup planet in a 1.8-year orbit of HD 5319, and a 2.51 M_Jup planet in a 1.1-year orbit of HD 75898. The measured eccentricities are 0.12 for HD 5319 b and 0.10 for HD 75898 b, and Markov Chain Monte Carlo simulations based on derived orbital parameters indicate that the radial velocities of both stars are consistent with circular planet orbits. With low eccentricity and 1 < a < 2 AU, our new planets have orbits similar to terrestrial planets in the solar system. The radial velocity residuals of both stars have significant trends, likely arising from substellar or low-mass stellar companions.

We report the results of infrared (8 micron) transit and secondary eclipse photometry of the hot Neptune exoplanet, GJ436b using Spitzer. The nearly photon-limited precision of these data allow us to measure an improved radius for the planet, and to detect the secondary eclipse. The transit (centered at HJD = 2454280.78149 +/-0.00016) shows the flat-bottomed shape typical of infrared transits, and it precisely defines the planet-to-star radius ratio (0.0839 +/-0.0005), independent of the stellar properties. However, we obtain the planetary radius, as well as the stellar mass and radius, by fitting to the transit curve simultaneously with an empirical mass-radius relation for M-dwarfs (M=R). We find Rs=Ms=0.47 +/-0.02 in solar units, and Rp=27,600 +/-1170 km (4.33 +/-0.18 Earth radii). This radius significantly exceeds the radius of a naked ocean planet, and requires a gasesous hydrogen-helium envelope. The secondary eclipse occurs at phase 0.587 +/-0.005, proving a significant orbital eccentricity (e=0.15 +/-0.012). The amplitude of the eclipse (5.7 +/-0.8e-4) indicates a brightness temperature for the planet of T=712 +/-36K. If this is indicative of the planet's physical temperature, it suggests the occurrence of tidal heating in the planet. An uncharacterized second planet likely provides ongoing gravitational perturbations, to maintain GJ436b's orbit eccentricity over long time scales.

We use a grid-based shallow water model to simulate the atmospheric dynamics of the transiting hot Jupiter HD 209458b. Under the usual assumption that the planet is in synchronous rotation with zero obliquity, a steady state is reached with a well-localized cold spot centered 76 degrees east of the antistellar point. This represents a departure from predictions made by previous simulations in the literature that used the shallow water formalism; we find that the disagreement is explained by the factor of 30 shorter radiative timescale used in our model. We also examine the case that the planet is in Cassini state 2, in which the expected obliquity is ~90 degrees. Under these circumstances, a periodic equilibrium is reached, with the temperature slightly leading the solar forcing. Using these temperature distributions, we calculate disk-integrated bolometric infrared light curves from the planet. The light curves for the two models are surprisingly similar, despite large differences in temperature patterns in the two cases. In the zero-obliquity case, the intensity at the minimum is 66% of the maximum intensity, with the minimum occuring 72 degrees ahead of transit. In the high-obliquity case, the minimum occurs 54 degrees ahead of transit, with an intensity of 58% of the maximum.

We report high-precision atmospheric parameters for 1907 stars in the N2K low-resolution spectroscopic survey, designed to identify metal-rich FGK dwarfs likely to harbor detectable planets. 284 of these stars are in the ideal temperature range for planet searches, T_eff <= 6000K, and have a 10% or greater probability of hosting planets based on their metallicities. The stars in the low-resolution spectroscopic survey should eventually yield >60 new planets, including 8-9 hot Jupiters. Short-period planets have already been discovered orbiting the survey targets HIP 14801 and HD 149143.

This paper considers general relativistic (GR) effects in currently observed extrasolar planetary systems. Although GR corrections are small, they can compete with secular interactions in these systems and thereby play an important role. Specifically, some of the observed multiple planet systems are close to secular resonance, where the dynamics is extremely sensitive to GR corrections, and these systems can be used as laboratories to test general relativity. For the three-planet solar system Upsilon Andromedae, secular interaction theory implies an 80% probability of finding the system with its observed orbital elements if GR is correct, compared with only a 2% probability in the absence of GR. In the future, tighter constraints can be obtained with increased temporal coverage.

The observational properties of brown dwarfs pose challenges to the theory of star formation. Because their mass is much smaller than the typical Jeans mass of interstellar clouds, brown dwarfs are most likely formed through secondary fragmentation processes, rather than through the direct collapse of a molecular cloud core. In order to prevent substantial post-formation mass accretion, young brown dwarfs must leave the high density formation regions in which they form. We propose here that brown dwarfs are formed in the optically thin outer regions of circumbinary disks. Through post-formation dynamical interaction with their host binary stars, young brown dwarfs are either scattered to large distance or removed, with modest speed, from their cradles.

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.

Using the most recent kinematic and radial velocity data in the literature, we calculate the binding energy of Proxima Centauri relative to the center of mass of the Alpha Centauri system. When we adopt the centroids of the observed data, we find that the three stars constitute a bound system, albeit with a semi-major axis that is on order the same size as Alpha Centauri AB's Hill radius in the galactic potential. We carry out a Monte Carlo simulation under the assumption that the errors in the observed quantities are uncorrelated. In this simulation, 44% of the trial systems are bound, and systems on the 1-3 sigma tail of the radial velocity distribution can have Proxima currently located near the apastron position of its orbit. Our analysis shows that a further, very significant improvement in the characterization of the system can be gained by obtaining a more accurate measurement of the radial velocity of Proxima Centauri.

This paper studies the long term evolution of planetary systems containing short-period planets, including the effects of tidal circularization, secular excitation of eccentricity by companion planets, and stellar damping. For planetary systems subject to all of these effects, analytic solutions (or approximations) are presented for the time evolution of the semi-major axes and eccentricities. Secular interactions enhance the inward migration and accretion of hot Jupiters, while general relativity tends to act in opposition by reducing the effectiveness of the secular perturbations. The analytic solutions presented herein allow us to understand these effects over a wide range of parameter space and to isolate the effects of general relativity in these planetary systems.

This paper studies the effects of dynamical interactions among the planets in observed extrasolar planetary systems, including hypothetical additional bodies, with a focus on secular perturbations. These interactions cause the eccentricities of the planets to explore a distribution of values over time scales that are long compared to observational time baselines, but short compared to the age of the systems. The same formalism determines the eccentricity forcing of hypothetical test bodies (terrestrial planets) in these systems and we find which systems allow for potentially habitable planets. Such planets would be driven to nonzero orbital eccentricity and we derive the distribution of stellar flux experienced by the planets over the course of their orbits. The general relativistic corrections to secular interaction theory are included in the analysis and such effects are important in systems with close planets ($\sim$4 day orbits). Some extrasolar planetary systems (e.g., Upsilon Andromedae) can be used as a test of general relativity, whereas in other systems, general relativity can be used to constrain the system parameters (e.g., $\sin i \gta 0.93$ for HD160691). For the case of hot Jupiters, we discuss how the absence of observed eccentricity implies the absence of companion planets.

We present a modification to the iodine cell Doppler technique that eliminates the need for an observed stellar template spectrum. For a given target star, we iterate toward a synthetic template spectrum beginning with an existing spectrum of a similar star. We then perturb the shape of this first-guess template to match the program observation of the target star taken through an iodine cell. The elimination of a separate template observation saves valuable telescope time, a feature that is ideally suited for the quick-look strategy employed by the ``Next 2000 Stars'' (N2K) planet search program. Tests using Keck/HIRES spectra indicate that synthetic templates yield a short-term precision of 3 m/s and a long-term, run-to-run precision of 5 m/s. We used this new Doppler technique to discover three new planets: a 1.5 Mjup planet in a 2.1375 d orbit around HD 86081; a 0.71 Mjup planet in circular, 26.73 d orbit around HD 224693; and a Saturn-mass planet in an 18.179 d orbit around HD 33283. The remarkably short period of HD 86081b bridges the gap between the extremely short-period planets detected in the OGLE survey and the 16 Doppler-detected hot jupiters (P < 15 d), which have an orbital period distribution that piles up at about three days. We have acquired photometric observations of two of the planetary host stars with the automated photometric telescopes at Fairborn Observatory. HD 86081 and HD 224693 both lack detectable brightness variability on their radial velocity periods, supporting planetary-reflex motion as the cause of the radial velocity variability. HD 86081 shows no evidence of planetary transits in spite of a 17.6% transit probability.

When extrasolar planets are observed to transit their parent stars, we are granted unprecedented access to their physical properties. It is only for transiting planets that we are permitted direct estimates of the planetary masses and radii, which provide the fundamental constraints on models of their physical structure. In particular, precise determination of the radius may indicate the presence (or absence) of a core of solid material, which in turn would speak to the canonical formation model of gas accretion onto a core of ice and rock embedded in a protoplanetary disk. Furthermore, the radii of planets in close proximity to their stars are affected by tidal effects and the intense stellar radiation. As a result, some of these "hot Jupiters" are significantly larger than Jupiter in radius. Precision follow-up studies of such objects (notably with the space-based platforms of the Hubble and Spitzer Space Telescopes) have enabled direct observation of their transmission spectra and emitted radiation. These data provide the first observational constraints on atmospheric models of these extrasolar gas giants, and permit a direct comparison with the gas giants of the Solar system. Despite significant observational challenges, numerous transit surveys and quick-look radial velocity surveys are active, and promise to deliver an ever-increasing number of these precious objects. The detection of transits of short-period Neptune-sized objects, whose existence was recently uncovered by the radial-velocity surveys, is eagerly anticipated. Ultra-precise photometry enabled by upcoming space missions offers the prospect of the first detection of an extrasolar Earth-like planet in the habitable zone of its parent star, just in time for Protostars and Planets VI.

We report the detection of a second exoplanet orbiting the G6V dwarf HD73526. This second planet has an orbital period of 377 d, putting it in a 2:1 resonance with the previously known exoplanet, the orbital period for which is updated to 188d. Dynamical modeling of the combined system allows solution for a self-consistent set of orbital elements for both components. HD73526 is the fourth exoplanetary system (of a total of 18 systems with 2 or more components currently known) to have components detected in 2:1 resonance. Finding such a large fraction of multiple planets (more than 20 per cent) in 2:1 resonance strongly suggests that orbital migration, halted by stabilisation in a trapping resonance, plays an important role in the evolution of exoplanets in multiple planet systems.

We present evidence that stars with planets exhibit statistically significant silicon and nickel enrichment over the general metal-rich population. We also present simulations which predict silicon enhancement of planet hosts within the context of the core-accretion hypothesis for giant planet formation. Because silicon and oxygen are both alpha elements, [Si/Fe] traces [O/Fe], so the silicon enhancement in planet hosts predicts that these stars are oxygen-rich as well. We present new numerical simulations of planet formation by core accretion that establish the timescale on which a Jovian planet reaches rapid gas accretion, t_rga, as a function of solid surface density sigma_solid: (t_rga / 1 Myr) = (sigma_solid / 25.0 g cm^{-2})^{-1.44}. This relation enables us to construct Monte Carlo simulations that predict the fraction of star-disk systems that form planets as a function of [Fe/H], [Si/Fe], disk mass, outer disk radius and disk lifetime. Our simulations reproduce both the known planet-metallicity correlation and the planet-silicon correlation reported in this paper. The simulations predict that 16% of Solar-type stars form Jupiter-mass planets, in agreement with 12% predicted from extrapolation of the observed planet frequency-semimajor axis distribution. Although a simple interpretation of core accretion predicts that the planet-silicon correlation should be much stronger than the planet-nickel correlation, we observe the same degree of silicon and nickel enhancement in planet hosts. If this result persists once more planets have been discovered, it might indicate a complexity in the chemistry of planet formation beyond the simple accumulation of solids in the core accretion theory.

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 have created a framework to facilitate the construction of specialized target lists for radial velocity surveys that are biased toward stars that (1) possess planets and (2) are easiest to observe with current detection techniques. We use a procedure that uniformly estimates fundamental stellar properties of Tycho 2 stars, with errors, using spline functions of broadband photometry and proper motion found in Hipparcos/Tycho 2 and 2MASS. We provide estimates of temperature and distance for 2.4 million Tycho 2 stars that lack trigonometric distances. For stars that appear to be FGK dwarfs according to estimated temperature and absolute magnitude, we also derive [Fe/H] and identify unresolved binary systems with mass ratios between 1.25 and 3. Our spline function models are trained on the unique Valenti & Fischer (2005) set, composed of 1000 dwarfs with precise stellar parameters estimated from HIRES spectroscopy. For FGK dwarfs with V photometric error less than 0.05 magnitudes, or V < 9, our temperature model gives a one-sigma error of +58.7/-65.9 K and our metallicity model gives a one-sigma error of +0.13/-0.14 dex. Our estimates of distance and spectral type enable us to isolate 354,822 Tycho 2 dwarfs, 321,996 of which are absent from Hipparcos, with giant and subgiant contamination at 2.6% and 7.2%, respectively. 2,500 of these FGK dwarfs are bright (V < 9.0) and metal-rich ([Fe/H] > 0.2). Our metallicity estimates have been used to identify targets for N2K (Fischer et al. 2005), a large-scale radial velocity search for Hot Jupiters, which has published the detection of 4 Hot Jupiters with one transit. The broadband filtering outlined here is the first screening tier for N2K; the second tier is a low-resolution spectroscopy program headed by S.E. Robinson (astro-ph/0510150).

We present empirical calibrations that estimate stellar metallicity, effective temperature and surface gravity as a function of Lick/IDS indices. These calibrations have been derived from a training set of 261 stars for which (1) high-precision measurements of [Fe/H], T_eff and log g have been made using spectral-synthesis analysis of HIRES spectra, and (2) Lick indices have also been measured. Our [Fe/H] calibration, which has precision 0.07 dex, has identified a number of bright (V < 9) metal-rich stars which are now being screened for hot Jupiter-type planets. Using the Yonsei-Yale stellar models, we show that the calibrations provide distance estimates accurate to 20% for nearby stars. This paper outlines the second tier of the screening of planet-search targets by the N2K Consortium, a project designed to identify the stars most likely to harbor extrasolar planets. Discoveries by the N2K Consortium include the transiting hot Saturn HD 149026 b (Sato et al. 2005, astro-ph/0507009) and HD 88133 b (Fischer et al. 2005). See Ammons et al. (2005, In Press) for a description of the first tier of N2K metallicity screening, calibrations using broadband photometry.

Doppler measurements from Subaru and Keck have revealed radial velocity variations in the V=8.15, G0IV star HD 149026 consistent with a Saturn-Mass planet in a 2.8766 day orbit. Photometric observations at Fairborn Observatory have detected three complete transit events with depths of 0.003 mag at the predicted times of conjunction. HD 149026 is now the second brightest star with a transiting extrasolar planet. The mass of the star, based on interpolation of stellar evolutionary models, is 1.3 +/- 0.1 solar masses; together with the Doppler amplitude, K=43.3 m s^-1, we derive a planet mass Msin(i)=0.36 Mjup, and orbital radius of 0.042 AU. HD 149026 is chromospherically inactive and metal-rich with spectroscopically derived [Fe/H]=+0.36, Teff=6147 K, log g=4.26 and vsin(i)=6.0 km s^-1. Based on Teff and the stellar luminosity of 2.72 Lsun, we derive a stellar radius of 1.45 Rsun. Modeling of the three photometric transits provides an orbital inclination of 85.3 +/- 1.0 degrees and (including the uncertainty in the stellar radius) a planet radius of 0.725 +/- 0.05 Rjup. Models for this planet mass and radius suggest the presence of a ~67 Mearth core composed of elements heavier than hydrogen and helium. This substantial planet core would be difficult to construct by gravitational instability.

Observations of TRES-1b transits made during the late summer and fall 2004 observing season reveal a statistically significant but low amplitude brightening event during egress.

We report precise Doppler measurements of GJ 436 (M2.5V) obtained at Keck Observatory. The velocities reveal a planetary companion with orbital period of 2.644 d, eccentricity of 0.12 (consistent with zero) and velocity semi-amplitude of $K =18.1$ \ms. The minimum mass (\msini) for the planet is 0.067 \mjup = 1.2 M$_{\rm NEP}$ = 21 M$_{\rm EARTH}$, making it the lowest mass exoplanet yet found around a main sequence star and the first candidate in the Neptune mass domain. GJ 436 (Mass = 0.41 \msune) is only the second M dwarf found to harbor a planet, joining the two--planet system around GJ 876. The low mass of the planet raises questions about its constitution, with possible compositions of primarily H and He gas, ice/rock, or rock--dominated. The implied semi--major axis is $a$ = 0.028 AU = 14 stellar radii, raising issues of planet formation, migration, and tidal coupling with the star. GJ 436 is $>3$ Gyr old, based on both kinematic and chromospheric diagnostics. The star exhibits no photometric variability on the 2.644-day Doppler period to a limiting amplitude of 0.0004 mag, supporting the planetary interpretation of the Doppler periodicity. Photometric transits of the planet across the star are ruled out for gas giant compositions and are also unlikely for solid compositions. As the third closest known planetary system, GJ 436 warrants follow--up observations by high resolution optical and IR

The N2K consortium is carrying out a distributed observing campaign with the Keck, Magellan and Subaru telescopes, as well as the automatic photometric telescopes of Fairborn Observatory, in order to search for short-period gas giant planets around metal-rich stars. We have established a reservoir of more than 14,000 main sequence and subgiant stars, closer than 110 pc, brighter than V=10.5 and with $0.4 < B-V < 1.2$. Because the fraction of stars with planets is a sensitive function of stellar metallicity, a broadband photometric calibration has been developed to identify a subset of 2000 stars with [Fe/H] $> 0.1$ dex for this survey. We outline the strategy and report the detection of a planet orbiting the metal-rich G5IV star HD 88133 with a period of 3.41 days, semi-velocity amplitude, K=35.7 \mse and \msini = 0.29 \mjupe. Photometric observations reveal that HD 88133 is constant on the 3.415-day radial velocity period to a limit of 0.0005 mag. Despite a transit probability of 15.6%, our photometry rules out the shallow transits predicted by the large stellar radius.

We present an updated analysis of the GJ 876 planetary system based on an augmented data set that incorporates 65 new high-precision radial velocities obtained with the Keck telescope from 2001 to 2004. These new radial velocities permit a more accurate characterization of the planet-planet interactions exhibited by the system. Self-consistent three-body orbital fits (which incorporate both the estimated instrumental uncertainties and 6 m/s Gaussian stellar jitter) continue to show that GJ 876 b and GJ 876 c are participating in a stable and symmetric 2:1 resonance condition in which the lowest order, eccentricity type mean-motion resonance variables are all librating. The planets are also locked in a secular resonance which causes them to librate about apsidal alignment. The small libration widths of all three resonances likely point to a dissipative history of differential migration for the two planets in the system. (for full abstract see paper).

The favored theoretical explanation for giant planet formation -- in both our solar system and others -- is the core accretion model (although it still has some serious difficulties). In this scenario, planetesimals accumulate to build up planetary cores, which then accrete nebular gas. With current opacity estimates for protoplanetary envelopes, this model predicts the formation of Jupiter-mass planets in 2--3 Myr at 5 AU around solar-mass stars, provided that the surface density of solids is enhanced over that of the minimum-mass solar nebula (by a factor of a few). Working within the core-accretion paradigm, this paper presents theoretical calculations which show that the formation of Jupiter-mass planets orbiting M dwarf stars is seriously inhibited at all radial locations (in sharp contrast to solar-type stars). Planet detection programs sensitive to companions of M dwarfs will test this prediction in the near future.

When stars form in small groups (N = 100 - 500 members), their circumstellar disks are exposed to little EUV radiation but a great deal of FUV radiation from massive stars in the group. This paper calculates mass loss rates for circumstellar disks exposed to external FUV radiation. Previous work treated large disks and/or intense radiation fields in which the disk radius exceeds the critical radius (supercritical disks) where the sound speed in the FUV heated layer exceeds the escape speed. This paper shows that significant mass loss still takes place for subcritical systems. Some of the gas extends beyond the disk edge (above the disk surface) to larger distances where the temperature is higher, the escape speed is lower, and an outflow develops. The evaporation rate is a sensitive function of the stellar mass and disk radius, which determine the escape speed, and the external FUV flux, which determines the temperature structure of the flow. Disks around red dwarfs are readily evaporated and shrink to disk radii of 15 AU on short time scales (10 Myr) when exposed to moderate FUV fields with $G_0$ = 3000. Although disks around solar type stars are more durable, these disks shrink to 15 AU in 10 Myr for intense FUV radiation fields with $G_0$ = 30,000; such fields exist in the central 0.7 pc of a cluster with N = 4000 stars. If our solar system formed in the presence of such strong FUV radiation fields, this mechanism could explain why Neptune and Uranus in our solar system are gas poor, whereas Jupiter and Saturn are gas rich. This mechanism for photoevaporation can also limit the production of Kuiper belt objects and can suppress giant planet formation in sufficiently large clusters, such as the Hyades, especially for disks associated with low mass stars.

We describe a project (transitsearch.org) currently attempting to discover transiting intermediate-period planets orbiting bright parent stars, and we simulate that project's performance. The discovery of such a transit would be an important astronomical advance, bridging the critical gap in understanding between HD209458b and Jupiter. However, the task is made difficult by intrinsically low transit probabilities and small transit duty cycles. This project's efficient and economical strategy is to photometrically monitor stars that are known (from radial velocity surveys) to bear planets, using a network of widely-spaced observers with small telescopes. These observers, each individually capable of precision (1%) differential photometry, monitor candidates during the time windows in which the radial velocity solution predicts a transit if the orbital inclination is close to 90 degrees. We use Monte Carlo techniques to simulate the performance of this network, performing simulations with different configurations of observers in order to optimize coordination of an actual campaign. Our results indicate that transitsearch.org can reliably rule out or detect planetary transits within the current catalog of known planet-bearing stars. A distributed network of skilled amateur astronomers and small college observatories is a cost-effective method for discovering the small number of transiting planets with periods in the range 10d < P < 200d that orbit bright (V < 11) stars.

The observational properties of brown dwarfs pose challenges to the theory of star formation. Because their mass is much smaller than the typical Jeans mass of interstellar clouds, brown dwarfs are most likely formed through secondary fragmentation processes, rather than through the direct collapse of a molecular cloud core. In order to prevent substantial post-formation mass accretion, young brown dwarfs must leave the high density formation regions in which they form. We propose here that brown dwarfs are formed in the circumbinary disks. Through post-formation dynamical interaction with their host binary stars, young brown dwarfs are either scattered to large distance or removed, with modest speed, from their cradles.

This paper examines how type I planet migration is affected by the presence of turbulent density fluctuations in the circumstellar disk. For type I migration, the planet does not clear a gap in the disk and its secular motion is driven by torques generated by the wakes it creates in the surrounding disk fluid. MHD turbulence creates additional density perturbations that gravitationally interact with the planet and can dominate the torques produced by the migration mechanism itself. This paper shows that conventional type I migration can be readily overwhelmed by turbulent perturbations and hence the usual description of type I migration should be modified in locations where the magnetorotational instability is active. In general, the migrating planet does not follow a smooth inward trned, but rather exhibits a random walk through phase space. Our main conclusion is that MHD turbulence will alter the time scales for type I planet migration and -- because of chaos -- requires the time scales to be described by a distribution of values.

We use hydrodynamical simulations to investigate the response of geometrically thin, self-gravitating, singular isothermal disks of gas to imposed rigidly rotating spiral potentials. By minimizing reflection-induced feedback from boundaries, and by restricting our attention to models where the swing parameter $X \sim 10$, we minimize the swing amplification of global normal modes even in models where Toomre's $Q_g \sim 1-2$ in the gas disk. We perform two classes of simulations: short-term ones over a few galactic revolutions where the background spiral forcing is large, and long-term ones over many galactic revolutions where the spiral forcing is considerably smaller. In both classes of simulations, the initial response of the gas disk is smooth and mimics the driving spiral field. At late times, many of the models evince substructure akin to the so-called branches, spurs, and feathers observed in real spiral galaxies. We comment on the parts played respectively by ultraharmonic resonances, reflection off internal features produced by nonlinear dredging, and local, transient, gravitational instabilites within spiral arms in the generation of such features. Our simulations reinforce the idea that spiral structure in the gaseous component becomes increasingly flocculent and disordered with the passage of time, even when the background population of old disk stars is a grand-design spiral. We speculate that truly chaotic behavior arises when many overlapping ultraharmonic resonances develop in reaction to an imposed spiral forcing that has itself a nonlinear, yet smooth, wave profile.

We have computed evolutionary models for extrasolar planets which range in mass from 0.1 to 3.0 Jovian Masses, and which range in equilibrium temperature from 113 K to 2000 K. We present four sequences of models, designed to show the structural effects of a solid core and of internal heating due to the conversion of kinetic to thermal energy at pressures of tens of bars. The model planetary radii are intended for comparisons with radii derived from observations of transiting extrasolar planets. To provide such comparisons, we expect that of order 10 transiting planets with orbital periods less than 200 days can be detected around bright (V<10) main-sequence stars for which accurate, well-sampled radial velocity measurements can be readily accumulated. Through these observations, structural properties of the planets will be derivable, particularly for low-mass, high-temperature planets. Implications regarding the transiting companion to OGLE-TR-56 recently announced by Konacki et al. are discussed. With regard to the confirmed transiting planet, HD 209458b, we find, in accordance with other recent calculations, that models without internal heating predict a radius that is 35 percent smaller than the observed radius. We explore the possibility that HD 209458b owes its large size to dissipation of energy arising from ongoing tidal circularization of the orbit. We show that residual scatter in the current radial velocity data set for HD 209458b is consistent with the presence of an as-of-yet undetected second companion, and that further radial velocity monitoring of the star is indicated.

This paper explores the intermediate-time dynamics of newly formed solar systems with a focus on possible mechanisms for planetary migration. We consider two limiting corners of the available parameter space -- crowded systems containing N=10 giant planets in the outer solar system, and solar systems with N=2 planets that are tidally interacting with a circumstellar disk. For a given set of initial conditions, dynamical relaxation leads to a well-defined distribution of possible solar systems. For each class of initial conditions, we perform large numbers of N-body simulations to obtain a statistical description of the possible outcomes. For N=10 planet systems, we consider several different planetary mass distributions; we also perform secondary sets of simulations to explore chaotic behavior and longer term dynamical evolution. For systems with 10 planets initially populating the radial range 5 - 30 AU, these scattering processes naturally produce planetary orbits with $a\sim1$ AU and the full range of possible eccentricity, but shorter period orbits are difficult to achieve. To account for the observed eccentric giant planets, we also explore a mechanism that combines dynamical scattering and tidal interactions with a circumstellar disk. This combined model naturally produces the observed range of semi-major axis and eccentricity. We discuss the relative merits of the different migration mechanisms for producing the observed eccentric giant planets.

(Abridged) We consider the expansion of an initially self-gravitating, static, singular, isothermal cloud core. For t>0, the gas is ionized and heated to a higher uniform temperature by the formation of a luminous, but massless, star in its center. The approximation that the mass and gravity of the central star is negligible for the subsequent motion of the HII region holds for distances much greater than about 100 AU and for the massive cloud cores that give rise to high-mass stars. If the initial ionization and heating is approximated to occur instantaneously at t=0, then the subsequent flow (for r >> 100 AU) caused by the resulting imbalance between self-gravity and thermal pressure is self-similar. Because of the steep density profile, pressure gradients produce a shock front that travels into the cloud, accelerating the gas to supersonic velocities in what has been called the ``champagne phase.'' We then study the self-similar solutions of the expansion of HII regions embedded in molecular clouds characterized by more general power-law density distributions with exponent 3/2 < n < 3. In these cases, the shock velocity is an increasing function of the exponent n, and diverges as n tends to 3. We show that this happens because the model includes an origin, where the pressure driving the shock diverges because the enclosed heated mass is infinite. Our results imply that the continued photoevaporation of massive reservoirs of neutral gas (e.g., surrounding disks and/or globules) nearby to the embedded ionizing source is required in order to maintain over a significant timescale the emission measure observed in champagne flows.

We report precise Doppler shift measurements of 55 Cancri (G8V) obtained from 1989 to 2002 at Lick Observatory. The velocities reveal evidence for an outer planetary companion to 55 Cancri orbiting at 5.5 AU. The velocities also confirm a second, inner planet at 0.11 AU. The outer planet is the first extrasolar planet found that orbits near or beyond the orbit of Jupiter. It was drawn from a sample of ~50 stars observed with sufficient duration and quality to detect a giant planet at 5 AU, implying that such planets are not rare. The properties of this jupiter analog may be compared directly to those of the Jovian planets in our Solar System. Its eccentricity is modest, e=0.16, compared with e=0.05 for both Jupiter and Saturn. Its mass is at least 4.0 jupiter masses (M sin i). The two planets do not perturb each other significantly. Moreover, a third planet of sub-Jupiter mass could easily survive in between these two known planets. Indeed a third periodicity remains in the velocity measurements with P = 44.3 d and a semi-amplitude of 13 m/s. This periodicity is caused either by a third planet at a=0.24 AU or by inhomogeneities on the stellar surface that rotates with period 42 d. The planet interpretation is more likely, as the stellar surface is quiet, exhibiting log(R'_{HK}) = -5.0 and brightness variations less than 1 millimag, and any hypothetical surface inhomogeneity would have to persist in longitude for 14 yr. Even with all three planets, an additional planet of terrestrial--mass could orbit stably at ~1 AU. The star 55 Cancri is apparently a normal, middle-aged main sequence star with a mass of 0.95 solar masses, rich in heavy elements ([Fe/H] = +0.27). This high metallicity raises the issue of the relationship between its age, rotation, and chromosphere.

We present an improved uvby-metallicity relation calibrated for F, G, and early K dwarfs, and an analogous uvby-T(eff) relation, both derived using a Levenberg-Marquardt minimization scheme. Our calibrations are based on 1533 stars which appear in both the Cayrel de Strobel (2001) metallicity compilation, and in the Hauck-Mermilliod (1998) catalog of uvby photometry. We also examine the speculative possibility of using uvby photometry to produce a uvby-planeticity calibration. We conclude that while there is likely no strong photometric indicator of the presence or absence of short-period planets, the possibility of a spectroscopic indicator of planeticity is well worth examining.

The mass and period ratios of the two planets orbiting 47 UMa suggest a possible kinship to the Jupiter-Saturn pair in our solar system. We explore the current dynamical state of the 47 UMa system with numerical integrations, and compare the results with analytic secular theory. We find that the planets in the system are likely participating in a secular resonance in which the difference in the longitudes of pericenter librates around zero. Alternately, it is possible that the system is participating in the 7:3 mean motion resonance. We show that stability considerations restrict the mutual inclination between the two planets to 40 degrees or less, and that this result is relatively insensitive to the total mass of the two planets. We present hydrodynamical simulations which measure the torques exerted on the planets by a hypothesized external protoplanetary disk. We show that planetary migration in response to torques from the disk may have led to capture of the system into a 7:3 mean-motion resonance, although it is unclear how the eccentricities of the planets would have been damped after capture occured. We show that Earth-mass planets can survive for long periods in some regions of the habitable zone of the nominal co-planar system. A set of planetary accretion calculations, however, shows that it is unlikely that large terrestrial planets can form in the 47 UMa habitable zone.

We explore the possibility that extrasolar planets might be found in the 1:1 mean-motion resonance. There are a variety of stable co-orbtial configurations, and we specifically examine three different versions of the 1:1 resonance. These include tadpole and horseshoe type orbits, as well as a more exotic configuration which occurs when one planet has a highly eccentric orbit while the other planet moves on a nearly circular orbit. We show that pairs of planets in 1:1 resonance yield characteristic radial velocity signatures which are not prone to the sin(i) degeneracy. Indeed, Keplerian fits to the radial velocities cannot reveal the presence of two planets in the 1:1 resonance. We discuss a dynamical fitting method for such systems, and illustrate its use with a simulated data set. Finally, we argue that hydrodynamical simulations and torqued three-body calculations indicate that 1:1 resonant pairs might readily form and migrate within protostellar disks.

Ten new low mass companions have emerged from the Keck precision Doppler velocity survey, with minimum (msini) masses ranging from 0.8 mjup to 0.34 msun. Five of these are planet candidates with msini < 12 mjup, two are brown dwarf candidates with msini ~30 mjup, and three are low mass stellar companions. Hipparcos astrometry reveals the orbital inclinations and masses for three of the (more massive) companions, and it provides upper limits to the masses for the rest. A new class of extrasolar planet is emerging, characterized by nearly circular orbits and orbital radii greater than 1 AU. The planet HD 4208b appears to be a member of this new class. The mass distribution of extrasolar planets continues to exhibit a rapid rise from 10 mjup toward the lowest detectable masses near 1 msat.

We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are executing independent Keplerian motions. We discuss two alternate fitting methods which can lead to an improved dynamical description of multiple planet systems. In the first method, the osculating orbital elements are determined via a Levenberg-Marquardt minimization scheme driving an N-body integrator. The second method is an improved analytic model in which orbital elements are allowed to vary according to a simple model for resonant interactions between the planets. Both of these methods can determine the true masses for the planets by eliminating the sin(i) degeneracy inherent in fits that assume independent Keplerian motions. We apply our fitting methods to the GJ876 radial velocity data (Marcy et al. 2001), and argue that the mass factors for the two planets are likely in the 1.25-2.0 range

The Sun's gradual brightening will seriously compromise the Earth's biosphere within ~ 1E9 years. If Earth's orbit migrates outward, however, the biosphere could remain intact over the entire main-sequence lifetime of the Sun. In this paper, we explore the feasibility of engineering such a migration over a long time period. The basic mechanism uses gravitational assists to (in effect) transfer orbital energy from Jupiter to the Earth, and thereby enlarges the orbital radius of Earth. This transfer is accomplished by a suitable intermediate body, either a Kuiper Belt object or a main belt asteroid. The object first encounters Earth during an inward pass on its initial highly elliptical orbit of large (~ 300 AU) semimajor axis. The encounter transfers energy from the object to the Earth in standard gravity-assist fashion by passing close to the leading limb of the planet. The resulting outbound trajectory of the object must cross the orbit of Jupiter; with proper timing, the outbound object encounters Jupiter and picks up the energy it lost to Earth. With small corrections to the trajectory, or additional planetary encounters (e.g., with Saturn), the object can repeat this process over many encounters. To maintain its present flux of solar energy, the Earth must experience roughly one encounter every 6000 years (for an object mass of 1E22 g). We develop the details of this scheme and discuss its ramifications.

We review the difficulties of the classical fission and fragmentation hypotheses for the formation of binary and multiple stars. A crucial missing ingredient in previous theoretical studies is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore formulate and solve the problem of the equilibria of isopedically magnetized, singular isothermal disks, without the assumption of axial symmetry. Considerable analytical progress can be made if we restrict our attention to models that are scale-free, i.e., that have surface densities that vary inversely with distance from the rotation axis of the system. In agreement with earlier analysis by Syer and Tremaine, we find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into M=2, 3, 4, ... separate bodies. On the basis of our experience in this paper, and the preceding Paper I, we advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

A crucial missing ingredient in previous theoretical studies of fragmentation is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore consider the equilibrium configuration of isopedically magnetized, scale-free, singular isothermal disks, without the assumption of axial symmetry. We find that lopsided (M=1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into separate bodies. We advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

Using the observed properties of our solar system, in particular the isotopic compositions of meteorites and the regularity of the planetary orbits, we constrain the star formation environment of the Sun within the scenario of (external) radioactive enrichment by a massive star. This calculation yields a probability distribution for the number of stars in the solar birth aggregate. The Sun is most likely to have formed within a stellar group containing $N$ = 2000 $\pm$ 1100 members. The {\it a priori} probability of a star forming in this type of environment is ${\cal P} \approx$ 0.0085, i.e., only about 1 out of 120 solar systems are expected to form under similar conditions. We discuss additional implications of this scenario, including possible effects from the radiation fields provided by the putative cluster environment and dynamical disruption of the Kuiper Belt. The constraints of this paper place tight restrictions on the properties of the solar birth aggregate for the scenario of external enrichment by a massive star; alternately, these tight constraints slightly favor a self-enrichment scenario for the short-lived radioactive species.

We examine the correlation between stellar metallicity and short period planets. It appears that approximately 1% of dwarf stars in the solar neighborhood harbor short-period planets characterized by near-circular orbits and orbital periods P<20 days. However, among the most metal-rich stars (defined as having [Fe/H]>0.2 dex), it appears that the fraction increases to 10%. Using the Hipparcos database and the Hauck & Mermilliod (1998) compilation of Stromgren uvby photometry, we identify a sample of 206 metal-rich stars of spectral type K, G, and F which have an enhanced probability of harboring short-period planets. Many of these stars would be excellent candidates for addition to radial velocity surveys. We have searched the Hipparcos epoch photometry for transiting planets within our 206 star catalog. We find that the quality of the Hipparcos data is not high enough to permit unambiguous transit detections. It is, however, possible to identify candidate transit periods. We then discuss various ramifications of the stellar metallicity - planet connection. First, we show that there is preliminary evidence for increasing metallicity with increasing stellar mass among known planet-bearing stars. This trend can be explained by a scenario in which planet-bearing stars accrete an average of 30 Earth Masses of rocky material after the gaseous protoplanetary disk phase has ended. We present dynamical calculation which suggest that a survey of metallicities of spectroscopic binary stars can be used to understand the root cause of the stellar metallicity - planet connection.

We explore possible effects of vacuum energy on the evolution of black holes. If the universe contains a cosmological constant, and if black holes can absorb energy from the vacuum, then black hole evaporation could be greatly suppressed. For the magnitude of the cosmological constant suggested by current observations, black holes larger than $\sim 4 \times 10^{24}$ g would accrete energy rather than evaporate. In this scenario, all stellar and supermassive black holes would grow with time until they reach a maximum mass scale of $\sim 6 \times 10^{55}$ g, comparable to the mass contained within the present day cosmological horizon.

We consider the long term fate and evolution of cold degenerate stars under the action of gravity alone. Although such stars cannot emit radiation through the Hawking mechanism, the wave function of the star will contain a small admixture of black hole states. These black hole states will emit radiation and hence the star can lose its mass energy in the long term. We discuss the allowed range of possible degenerate stellar evolution within this framework.

(Abridged Abstract) A favored interpretation of recent microlensing measurements towards the Large Magellanic Cloud implies that a large fraction (i.e. 10--50%) of the mass of the galactic halo is composed of white dwarfs. We compare model white dwarf luminosity functions to the data from the observational surveys in order to determine a lower bound on the age of any substantial white dwarf halo population (and hence possibly on the age of the Universe). We compare various theoretical white dwarf luminosity functions, in which we vary hese three parameters, with the abovementioned survey results. From this comparison, we conclude that if white dwarfs do indeed constitute more than 10% of the local halo mass density, then the Universe must be at least 10 Gyr old for our most extreme allowed values of the parameters. When we use cooling curves that account for chemical fractionation and more likely values of the IMF and the bolometric correction, we find tighter limits: a white dwarf MACHO fraction of 10% (30%) requires a minimum age of 14 Gyr (15.5 Gyr). Our analysis also indicates that the halo white dwarfs almost certainly have helium-dominated atmospheres.

A number of recently discovered extrasolar planet candidates have surprisingly small orbits, which may indicate that considerable orbital migration takes place in protoplanetary systems. A natural consequence of orbital migration is for a series of planets to be accreted, destroyed, and then thoroughly mixed into the convective envelope of the central star. We study the ramifications of planet accretion for the final main sequence metallicity of the star. If maximum disk lifetimes are on the order of 10 Myr, stars with masses near 1 solar mass are predicted to have virtually no metallicity enhancement. On the other hand, early F and late A type stars with masses of 1.5--2.0 solar masses can experience significant metallicity enhancements due to their considerably smaller convection zones during the first 10 Myr of pre-main-sequence evolution. We show that the metallicities of an aggregate of unevolved F stars are consistent with an average star accreting about 2 Jupiter-mass planets from a protoplanetary disk having a 10 Myr dispersal time.

Motivated by recent experimental proposals to search for extragalactic cosmic rays (including anti-matter from distant galaxies), we study particle propagation through the intergalactic medium (IGM). We first use estimates of the magnetic field strength between galaxies to constrain the mean free path for diffusion of particles through the IGM. We then develop a simple analytic model to describe the diffusion of cosmic rays. Given the current age of galaxies, our results indicate that, in reasonable models, a completely negligible number of particles can enter our Galaxy from distances greater than $\sim 100$ Mpc for relatively low energies ($E$ $< 10^6$ GeV/n). We also find that particle destruction in galaxies along the diffusion path produces an exponential suppression of the possible flux of extragalactic cosmic rays. Finally, we use gamma ray constraints to argue that the distance to any hypothetical domains of anti-matter must be roughly comparable to the horizon scale.

This paper outlines astrophysical issues related to the long term fate of the universe. We consider the evolution of planets, stars, stellar populations, galaxies, and the universe itself over time scales which greatly exceed the current age of the universe. (shortened version of abstract).

Motivated by recent measurements which suggest that roughly half the mass of the galactic halo may be in the form of white dwarfs, we study the implications of such a halo. We first use current limits on the infrared background light and the galactic metallicity to constrain the allowed initial mass function (IMF) of the stellar population that produced the white dwarfs. The IMF must be sharply peaked about a characteristic mass scale $M_C \approx 2.3 M_\odot$. Since only a fraction of the initial mass of a star is incorporated into the remnant white dwarf, we argue that the mass fraction of white dwarfs in the halo is likely to be 25\% or less, and that 50\% is an extreme upper limit. We use the IMF results to place corresponding constraints on the primordial initial conditions for star formation. The initial conditions must be much more homogeneous and skewed toward higher temperatures ($T_{\rm gas} \sim$ 200 K) than the conditions which lead to the present day IMF. Next we determine the luminosity function of white dwarfs. By comparing this result with the observed luminosity function, we find that the age of the halo population must be greater than $\sim 16$ Gyr. Finally, we calculate the radiative signature of a white dwarf halo. This infrared background is very faint, but is potentially detectable with future observations.