The origin of close-in Jovian planets is still elusive. We examine the in-situ gas accretion scenario as a formation mechanism of these planets. We reconstruct natal disk properties from the occurrence rate distribution of close-in giant planets, under the assumption that the occurrence rate may reflect the gas accretion efficiency onto cores of these planets. We find that the resulting gas surface density profile becomes an increasing function of the distance from the central star with some structure at $r \simeq 0.1$ au. This profile is quite different from the standard minimum-mass solar nebula model, while our profile leads to better reproduction of the population of observed close-in super-Earths based on previous studies. We compute the resulting magnetic field profiles and find that our profiles can be fitted by stellar dipole fields ($\propto r^{-3}$) in the vicinity of the central star and large-scale fields ($\propto r^{-2}$) at the inner disk regions, either if the isothermal assumption breaks down or if nonideal MHD effects become important. For both cases, the transition between these two profiles occurs at $r \simeq 0.1$ au, which corresponds to the period valley of giant exoplanets. Our work provides an opportunity to test the in-situ gas accretion scenario against disk quantities, which may constrain the gas distribution of the minimum-mass {\it extra}solar nebula.

The TRAPPIST-1 system has 7 known terrestrial planets arranged compactly in a mean motion resonant chain around an ultra-cool central star, some within the estimated habitable zone. Given their short orbital periods of just a few days, it is often presumed that the planets are tidally locked such that the spin rate is equal to that of the orbital mean motion. However, the compact, and resonant, nature of the system implies that there can be significant variations in the mean motion of these planets due to their mutual interactions. We show that such fluctuations can then have significant effects on the spin states of these planets. In this paper, we analyze, using detailed numerical simulations, the mean motion histories of the three planets that are thought to lie within or close to the habitable zone of the system: planets d, e, and f. We demonstrate that, depending on the strength of the mutual interactions within the system, these planets can be pushed into spin states which are effectively non-synchronous. We find that it can produce significant libration of the spin state, if not complete circulation in the frame co-rotating with the orbit. We also show that these spin states are likely to be unable to sustain long-term stability, with many of our simulations suggesting that the spin evolves, under the influence of tidal synchronization forces, into quasi-stable attractor states, which last on timescales of thousands of years.

Using the updated exoplanet population parameters of our previous study, which includes the planetary radius updates from Gaia DR2 and an inferred multiplicity distribution, we provide a revised $\eta_{\oplus}$ calculation. This is achieved by sampling planets from our derived population model and determining which planets meet our criterion for habitability. To ensure robust results, we provide probabilities calculated over a range of upper radius limits. Our most optimistic criterion for habitability provides an $\eta_{\oplus}$ value of $0.34\pm 0.02 \frac{\rm planets}{\rm star}$. We also consider the effects of multiplicity and the number of habitable planets each system may contain. Our calculation indicates that $6.4\pm0.5\%$ of GK dwarfs have more than one planet within their habitable zone. This optimistic habitability criterion also suggests that $0.036\pm0.009\%$ of solar-like stars will harbor 5 or more habitable planets. These tightly packed highly habitable system should be extremely rare, but still possible. Even with our most pessimistic criterion we still expect that $1.8\pm0.2\%$ of solar-like stars harbor more than one habitable planet.

Identification of the main planet formation site is fundamental to understanding how planets form and migrate to the current locations. We consider the heavy-element content trend of observed exoplanets derived from improved measurements of mass and radius, and explore how this trend can be used as a tracer of their formation sites. Using gas accretion recipes obtained from detailed hydrodynamical simulations, we confirm that the disk-limited gas accretion regime is most important for reproducing the heavy-element content trend. Given that such a regime is specified by two characteristic masses of planets, we compute these masses as a function of the distance ($r$) from the central star, and then examine how the regime appears in the mass-semimajor axis diagram. Our results show that a plausible solid accretion region emerges at $r \simeq 0.6$ au and expands with increasing $r$, using the conventional disk model. Given that exoplanets that possess the heavy-element content trend distribute currently near their central stars, our results imply the importance of planetary migration that would occur after solid accretion onto planets might be nearly completed at $r \geq 0.6$ au. Self-consistent simulations would be needed to verify the predictions herein.

We present precision 4.5 $\mu$m Spitzer transit photometry of eight planet candidates discovered by the K2 mission: K2-52 b, K2-53 b, EPIC 205084841.01, K2-289 b, K2-174 b, K2-87 b, K2-90 b, and K2-124 b. The sample includes four sub-Neptunes and two sub-Saturns, with radii between 2.6 and 18 $R_\oplus$, and equilibrium temperatures between 440 and 2000 K. In this paper we identify several targets of potential interest for future characterization studies, demonstrate the utility of transit follow-up observations for planet validation and ephemeris refinement, and present new imaging and spectroscopy data. Our simultaneous analysis of the K2 and Spitzer light curves yields improved estimates of the planet radii, and multi-wavelength information which help validate their planetary nature, including the previously un-validated candidate EPIC 205686202.01 (K2-289 b). Our Spitzer observations yield an order of magnitude increase in ephemeris precision, thus paving the way for efficient future study of these interesting systems by reducing the typical transit timing uncertainty in mid-2021 from several hours to a dozen or so minutes. K2-53 b, K2-289 b, K2-174 b, K2-87 b, and K2-90 b are promising radial velocity (RV) targets given the performance of spectrographs available today or in development, and the M3V star K2-124 hosts a temperate sub-Neptune that is potentially a good target for both RV and atmospheric characterization studies.

We investigate the role that planet detection order plays in the Kepler planet detection pipeline. The Kepler pipeline typically detects planets in order of descending signal strength (MES). We find that the detectability of transits experiences an additional $5.5\%$ and $15.9\%$ efficiency loss, for periods $<200$ days and $>200$ days respectively, when detected after the strongest signal transit in a multiple-planet system. We provide a method for determining the transit probability for multiple-planet systems by marginalizing over the empirical Kepler dataset. Furthermore, because detection efficiency appears to be a function of detection order, we discuss the sorting statistics that affect the radius and period distributions of each detection order. Our occurrence rate dataset includes radius measurement updates from the California Kepler Survey (CKS), Gaia DR2, and asteroseismology. Our population model is consistent with the results of Burke et al. (2015), but now includes an improved estimate of the multiplicity distribution. From our obtained model parameters, we find that only $4.0\pm4.6\%$ of solar-like GK dwarfs harbor one planet. This excess is smaller than prior studies and can be well modeled with a modified Poisson distribution, suggesting that the Kepler Dichotomy can be accounted for by including the effects of multiplicity on detection efficiency. Using our modified Poisson model we expect the average number of planets is $5.86\pm0.18$ planets per GK dwarf within the radius and period parameter space of Kepler.

We show that a model in which Mars grows near Earth and Venus but is then scattered out of the terrestrial region yields a natural pathway to explain the low masses of the Martian moons Phobos & Deimos. In this scenario, the last giant impact experienced by Mars is followed by an extended period (tens to hundreds of Myr) of close passages by other planetary embryos. These close passages perturb and dynamically heat any system of forming satellites left over by the giant impact and can substantially reduce the mass in the satellite system (sometimes to zero). The close passage of massive perturbing bodies also offers the opportunity to capture small objects by three-body scattering. Both mechanisms lead to low mass moon systems with a substantially collisional history.

We show that the rate of pollution of white dwarfs by asteroidal material implies a concomitant rate of material ejection that can contribute significantly to the population of interstellar minor bodies. We note also that the irradiation during post main sequence evolution implies that much of this ejected material may lose volatiles, providing a rationale for the curious properties of the recently discovered interstellar object Oumuamua.

We investigate the constraints on the formation of, and tidal dissipation processes in, hot Jupiters (HJs) that can be inferred based on reliable knowledge of the age of a system or population. Particular attention is paid to the role of young systems (such as those in open clusters or star-forming regions) in such studies. For an ensemble of coeval HJ (or proto-HJ) systems, we quantify the effect of age on the distribution of orbital eccentricities with respect to orbital periods as well as the location of the observed "pile-up" feature. We expect the effects of pre-main-sequence stellar evolution to be important only if a substantial fraction of HJs approach their current orbits early in protostellar contraction (ages less than 10 Myr). Application to the HJs presently known in the cluster M67 yields constraints on the dissipation roughly consistent with those gleaned from planets in the field; for those in the Hyades and Praesepe, our results suggest a higher degree of dissipation at early times than that inferred from other populations.

One longstanding problem for the potential habitability of planets within M dwarf systems is their likelihood to be tidally locked in a synchronously rotating spin state. This problem thus far has largely been addressed only by considering two objects: the star and the planet itself. However, many systems have been found to harbor multiple planets, with some in or very near to mean-motion resonances. The presence of a planetary companion near a mean-motion resonance can induce oscillatory variations in the mean-motion of the planet, which we demonstrate can have significant effects on the spin-state of an otherwise synchronously rotating planet. In particular, we find that a planetary companion near a mean-motion resonance can excite the spin states of planets in the habitable zone of small, cool stars, pushing otherwise synchronously rotating planets into higher amplitude librations of the spin state, or even complete circulation resulting in effective stellar days with full surface coverage on the order of years or decades. This increase in illuminated area can have potentially dramatic influences on climate, and thus on habitability. We also find that the resultant spin state can be very sensitive to initial conditions due to the chaotic nature of the spin state at early times within certain regimes. We apply our model to two hypothetical planetary systems inspired by the K00255 and TRAPPIST-1 systems, which both have Earth-sized planets in mean-motion resonances orbiting cool stars.

HD 3167 is a bright (V = 8.9), nearby K0 star observed by the NASA K2 mission (EPIC 220383386), hosting two small, short-period transiting planets. Here we present the results of a multi-site, multi-instrument radial velocity campaign to characterize the HD 3167 system. The masses of the transiting planets are 5.02+/-0.38 MEarth for HD 3167 b, a hot super-Earth with a likely rocky composition (rho_b = 5.60+2.15-1.43 g/cm^3), and 9.80+1.30-1.24 MEarth for HD 3167 c, a warm sub-Neptune with a likely substantial volatile complement (rho_c = 1.97+0.94-0.59 g/cm^3). We explore the possibility of atmospheric composition analysis and determine that planet c is amenable to transmission spectroscopy measurements, and planet b is a potential thermal emission target. We detect a third, non-transiting planet, HD 3167 d, with a period of 8.509+/-0.045 d (between planets b and c) and a minimum mass of 6.90+/-0.71 MEarth. We are able to constrain the mutual inclination of planet d with planets b and c: we rule out mutual inclinations below 1.3 degrees as we do not observe transits of planet d. From 1.3-40 degrees, there are viewing geometries invoking special nodal configurations which result in planet d not transiting some fraction of the time. From 40-60 degrees, Kozai-Lidov oscillations increase the system's instability, but it can remain stable for up to 100Myr. Above 60 degrees, the system is unstable. HD 3167 promises to be a fruitful system for further study and a preview of the many exciting systems expected from the upcoming NASA TESS mission.

We report precise mass and density measurements of two extremely hot sub-Neptune-size planets from the K2 mission using radial velocities, K2 photometry, and adaptive optics imaging. K2-66 harbors a close-in sub-Neptune-sized (2.49$^{+0.34}_{-0.24} R_\oplus$) planet (K2-66b) with a mass of 21.3 $\pm$ 3.6 $M_\oplus$. Because the star is evolving up the sub-giant branch, K2-66b receives a high level of irradiation, roughly twice the main sequence value. K2-66b may reside within the so-called "photoevaporation desert", a domain of planet size and incident flux that is almost completely devoid of planets. Its mass and radius imply that K2-66b has, at most, a meager envelope fraction (< 5%) and perhaps no envelope at all, making it one of the largest planets without a significant envelope. K2-106 hosts an ultra-short-period planet ($P$ = 13.7 hrs) that is one of the hottest sub-Neptune-size planets discovered to date. Its radius (1.82$^{+0.20}_{-0.14} R_\oplus$) and mass (9.0 $\pm$ 1.6 $M_\oplus$) are consistent with a rocky composition, as are all other small ultra-short-period planets with well-measured masses. K2-106 also hosts a larger, longer-period planet (Rp = 2.77$^{+0.37}_{-0.23} R_\oplus$, $P$ = 13.3 days) with a mass less than 24.4 $M_\oplus$ at 99.7% confidence. K2-66b and K2-106b probe planetary physics in extreme radiation environments. Their high densities reflect the challenge of retaining a substantial gas envelope in such extreme environments.

We report precise radial velocity (RV) measurements of WASP-47, a G star that hosts three transiting planets in close proximity (a hot Jupiter, a super-Earth and a Neptune-sized planet) and a non-transiting planet at 1.4 AU. Through a joint analysis of previously published RVs and our own Keck-HIRES RVs, we significantly improve the planet mass and bulk density measurements. For the super-Earth WASP-47e ($P$ = 0.79 days), we measure a mass of 9.11 $\pm$ 1.17 $M_\oplus$, and a bulk density of 7.63 $\pm$ 1.90 g cm$^{-3}$, consistent with a rocky composition. For the hot Jupiter WASP-47b ($P$ = 4.2 days), we measure a mass of 356 $\pm$ 12 $M_\oplus$ (1.12 $\pm$ 0.04 $M_\rm{Jup}$) and constrain its eccentricity to $<0.021$ at 3-$\sigma$ confidence. For the Neptune-size planet WASP-47d ($P$ = 9.0 days), we measure a mass of 12.75 $\pm$ 2.70 $M_\oplus$, and a bulk density of 1.36 $\pm$ 0.42 g cm$^{-3}$, suggesting it has a thick H/He envelope. For the outer non-transiting planet, we measure a minimum mass of 411 $\pm$ 18 $M_\oplus$ (1.29 $\pm$ 0.06 $M_\rm{Jup}$), an orbital period of 595.7 $\pm$ 5.0 days, and an orbital eccentricity of 0.27 $\pm$ 0.04. Our new measurements are consistent with but 2$-$4$\times$ more precise than previous mass measurements.

We examine the effect of secular perturbations by giant planets on systems of multiple, lower mass planets orbiting Sun-like stars. We simulate the effects of forcing both eccentricity and inclination, separately and together. We compare our results to the statistics of the observed Kepler data and examine whether these results can be used to explain the observed excess of single transiting planets. We cannot explain the observed excess by pumping only inclination without driving most systems over the edge of dynamical instability. Thus, we expect the underlying planetary population for systems with a single transitting planet to contain an intrinsically low multiplicity population. We can explain the Kepler statistics and occurrence rates for R< 2 Rearth planets with a perturber population consistent with that inferred from radial velocity surveys, but require too many giant planets if we wish to explain all planets with R < 4 Rearth. These numbers can be brought into agreement if we posit the existence of an equivalent size population of planets below the RV detection limit (of characteristic mass 0.1 Mjupiter). This population would need to be dynamically hot to produce sufficiently strong perturbations and would leave the imprint of high obliquities amongst the surviving planets. Thus, an extensive sample of obliquities for low mass planets can help to indicate the presence of such a population. (...) Some of our simulations also produce planetary systems with planets that survive in the habitable zone but have no planets interior to them -- much as in the case of our Solar System. This suggests that such a configuration may not be altogether rare, but may occur around a few percent of FGK stars.

We present 197 planet candidates discovered using data from the first year of the NASA K2 mission (Campaigns 0-4), along with the results of an intensive program of photometric analyses, stellar spectroscopy, high-resolution imaging, and statistical validation. We distill these candidates into sets of 104 validated planets (57 in multi-planet systems), 30 false positives, and 63 remaining candidates. Our validated systems span a range of properties, with median values of R_P = 2.3 R_E, P=8.6 d, Tef = 5300 K, and Kp=12.7 mag. Stellar spectroscopy provides precise stellar and planetary parameters for most of these systems. We show that K2 has increased by 30% the number of small planets known to orbit moderately bright stars (1-4 R_E, Kp=9-13 mag). Of particular interest are 37 planets smaller than 2 R_E, 15 orbiting stars brighter than Kp=11.5, five receiving Earth-like irradiation levels, and several multi-planet systems -- including four planets orbiting the M dwarf K2-72 near mean-motion resonances. By quantifying the likelihood that each candidate is a planet we demonstrate that our candidate sample has an overall false positive rate of 15-30%, with rates substantially lower for small candidates (< 2 R_E) and larger for candidates with radii > 8 R_E and/or with P < 3 d. Extrapolation of the current planetary yield suggests that K2 will discover between 500-1000 planets in its planned four-year mission -- assuming sufficient follow-up resources are available. Efficient observing and analysis, together with an organized and coherent follow-up strategy, is essential to maximize the efficacy of planet-validation efforts for K2, TESS, and future large-scale surveys.

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4\sigma$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7\sigma$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

We present a catalog of 11 multi-planet systems from Campaigns 1 and 2 of the K2 mission. We report the sizes and orbits of 26 planets split between seven 2-planet systems and four 3-planet systems. These planets stem from a systematic search of the K2 photometry for all dwarf stars observed by K2 in these fields. We precisely characterized the host stars with adaptive optics imaging and analysis of high-resolution optical spectra from Keck/HIRES and medium-resolution spectra from IRTF/SpeX. We confirm two planet candidates by mass detection and validate the remaining 24 candidates to $>99\%$ confidence. Thirteen planets were previously validated or confirmed by other studies and 24 were previously identified as planet candidates. The planets are mostly smaller than Neptune (21/26 planets) as in the Kepler mission and all have short periods ($P < 50$ d) due to the duration of the K2 photometry. The host stars are relatively bright (most have $Kp < 12.5$ mag) and are amenable to follow-up characterization. For K2-38, we measured precise radial velocities using Keck/HIRES and provide initial estimates of the planet masses. K2-38b is a short-period super-Earth with a radius of $1.55 \pm 0.16~R_\oplus$, a mass of $12.0 \pm 2.9~M_\oplus$, and a high density consistent with an iron-rich composition. The outer planet K2-38c is a lower density sub-Neptune-size planet with a radius of $2.42 \pm 0.29~R_\oplus$ and a mass of $9.9 \pm 4.6~M_\oplus$ that likely has a substantial envelope. This new planet sample demonstrates the capability of K2 to discover numerous planetary systems around bright stars.

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

Warm jupiters are an unexpected population of extrasolar planets that are too near to their host to have formed in situ, but distant enough to retain a significant eccentricity in the face of tidal damping. These planets are curiously absent around stars larger than two solar radii. We hypothesize that the warm jupiters are migrating due to Kozai-Lidov oscillations, which leads to transient episodes of high eccentricity and a consequent tidal decay. As their host evolves, such planets would be rapidly dragged in or engulfed at minimum periapse, leading to a rapid depletion of the population with increasing stellar radius, as is observed. Using numerical simulations, we determine the relationship between periapse distance and orbital migration rate for planets 0.1 to 10 Jupiter masses and with orbital periods between 10 and 100 days. We find that Kozai-Lidov oscillations effectively result in planetary removal early in the evolution of the host star, possibly accounting for the observed deficit. While the observed eccentricity distribution is inconsistent with the simulated distribution for an oscillating and migrating warm jupiter population, observational biases may explain the discrepancy.

Discoveries from the prime Kepler mission demonstrated that small planets (< 3 Earth-radii) are common outcomes of planet formation. While Kepler detected many such planets, all but a handful orbit faint, distant stars and are not amenable to precise follow up measurements. Here, we report the discovery of two small planets transiting K2-21, a bright (K = 9.4) M0 dwarf located 65$\pm$6 pc from Earth. We detected the transiting planets in photometry collected during Campaign 3 of NASA's K2 mission. Analysis of transit light curves reveals that the planets have small radii compared to their host star, 2.60 $\pm$ 0.14% and 3.15 $\pm$ 0.20%, respectively. We obtained follow up NIR spectroscopy of K2-21 to constrain host star properties, which imply planet sizes of 1.59 $\pm$ 0.43 Earth-radii and 1.92 $\pm$ 0.53 Earth-radii, respectively, straddling the boundary between high-density, rocky planets and low-density planets with thick gaseous envelopes. The planets have orbital periods of 9.32414 days and 15.50120 days, respectively, and have a period ratio of 1.6624, very near to the 5:3 mean motion resonance, which may be a record of the system's formation history. Transit timing variations (TTVs) due to gravitational interactions between the planets may be detectable using ground-based telescopes. Finally, this system offers a convenient laboratory for studying the bulk composition and atmospheric properties of small planets with low equilibrium temperatures.

We present Hubble Space Telescope observations of the upper part (T_eff> 10 000 K) of the white dwarf cooling sequence in the globular cluster 47 Tucanae and measure a luminosity function of hot white dwarfs. Comparison with previous determinations from large scale field surveys indicates that the previously determined plateau at high effective temperatures is likely a selection effect, as no such feature is seen in this sample. Comparison with theoretical models suggests that the current estimates of white dwarf neutrino emission (primarily by the plasmon channel) are accurate, and variations are restricted to no more than a factor of two globally, at 95% confidence. We use these constraints to place limits on various proposed exotic emission mechanisms, including a non-zero neutrino magnetic moment, formation of axions, and emission of Kaluza-Klein modes into extra dimensions.

We demonstrate that tidal evolution of the inner planet (`e') of the system orbiting the star rho 55 Cancri could have led to passage through two secular resonances with other planets in the system. The consequence of this evolution is excitation of both the planetary eccentricity and inclination relative to the original orbital plane. The large mass ratio between the innermost planet and the others means that these excitations can be of substantial amplitude and can have dramatic consequences for the system organisation. Such evolution can potentially explain the large observed mutual inclination between the innermost and outermost planets in the system, and implies that tidal heating could have substantially modified the structure of planet e, and possibly reduced its mass by Roche lobe overflow. Similar inner secular resonances may be found in many multiple planet systems and suggest that many of the innermost planets in these systems could have suffered similar evolutions.

Small, cool planets represent the typical end-products of planetary formation. Studying the archi- tectures of these systems, measuring planet masses and radii, and observing these planets' atmospheres during transit directly informs theories of planet assembly, migration, and evolution. Here we report the discovery of three small planets orbiting a bright (Ks = 8.6 mag) M0 dwarf using data collected as part of K2, the new transit survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and K2 photometry indicate that the system hosts three transiting planets with radii 1.5-2.1 R_Earth, straddling the transition region between rocky and increasingly volatile-dominated compositions. With orbital periods of 10-45 days the planets receive just 1.5-10x the flux incident on Earth, making these some of the coolest small planets known orbiting a nearby star; planet d is located near the inner edge of the system's habitable zone. The bright, low-mass star makes this system an excellent laboratory to determine the planets' masses via Doppler spectroscopy and to constrain their atmospheric compositions via transit spectroscopy. This discovery demonstrates the ability of K2 and future space-based transit searches to find many fascinating objects of interest.

We describe the long-term evolution of compact systems of terrestrial planets, using a set of simulations that match the statistical properties of the observed exoplanet distribution. The evolution is driven by tidal dissipation in the planetary interiors, but the systems evolve as a whole due to secular gravitational interactions. We find that, for Earth-like dissipation levels, planetary orbits can be circularised out to periods of order 100 days, an order of magnitude larger than is possible for single planets. The resulting distribution of eccentricities is a qualitative match to that inferred from transit timing variations, with a minority of non-zero eccentricities maintained by particular secular configurations. The coupling of the tidal and secular processes enhance the inward migration of the innermost planets in these systems, and can drive them to short orbital periods. Resonant interactions of both the mean motion and secular variety are observed, although the interactions are not strong enough to drive systemic instability in most cases. However, we demonstrate that these systems can easily be driven unstable if coupled to giant planets on longer period orbits.

We present a model for the in situ assembly of planetary systems around a 0.5~$M_{\odot}$ star, and compare the resulting statistics to the observed sample of cool Kepler planet candidates from Dressing & Charbonneau (2013). We are able to reproduce the distribution of planetary periods and period ratios, although we once again find an underabundance of single transit systems relative to the observations. We also demonstrate that almost every planetary system assembled in this fashion contains at least one planet in the habitable zone, and that water delivery to these planets can potentially produce a water content comparable to that of Earth. Our results broadly support the notion that habitable planets are plentiful around M~dwarfs in the solar neighbourhood.

We examine the consequences of a model for the circulation of solids in a protoplanetary nebula in which aerodynamic drag is counterbalanced by the recycling of material to the outer disc by a protostellar outflow or a disc wind. This population of circulating dust eventually becomes unstable to the formation of planetesimals by gravitational instability, and results in the ultimate deposition of 30--50 earth masses in planetesimals on scales R< 1 AU. Such a model may provide an appropriate justification for the approximately power law initial conditions needed to reproduce observed planetary systems by in situ assembly.

A significant fraction of white dwarfs (WDs) are observed to be polluted with metals despite high surface gravities and short settling times. The current theoretical model for this pollution is accretion of rocky bodies delivered to the WD through perturbations by orbiting planets. Using N-body simulations, we examine the possibility of a single planet as the source of pollution. We determine the stability of test particles on circular orbits in systems with a single planet located at 4 au for a range of masses and eccentricities, comparing the fractions that are ejected and accreted. In particular, we compare the instabilities that develop before and after the star loses mass to form a WD, a process which causes orbiting bodies to migrate outward. We determine that a planet must be eccentric (e > 0.02) to deliver significant (> 0.5 per cent) amounts of material to the host and that the amount increases with the planetary eccentricity. This result is robust with respect to the initial eccentricities of the particles for planetary eccentricity above ~0.4 and for randomly-distributed particle long. of pericentre. We also find that the efficiency of pollution is enhanced as planetary mass is reduced. We demonstrate that a 0.03 M_Jup planet with substantial eccentricity (e > 0.4) can account for the observed levels of pollution for initial disc masses of order 1 M_Earth. Such discs are within the range estimated for initial planetesimals discs and below that estimated for the solar system. However, their survival to the WD stage is uncertain as estimates for the collisional evolution of planetesimal discs suggest they should be ground down below the required levels on Gyr timescales. Thus, planetary scattering by eccentric, sub-Jovian planets can explain the observed levels of WD pollution, but only if current estimates of the collisional erosion of planetesimal discs are in error.

We report our spectroscopic investigation of the transiting ice giant GJ 3470b's atmospheric transmission, and the first results of extrasolar planet observations from the new Keck/MOSFIRE spectrograph. We measure a planet/star radius ratio of Rp/Rs = 0.0789 +/- 0.0020 in a bandpass from 2.09-2.36 micron and in six narrower bands across this wavelength range. When combined with existing broadband photometry, these measurements rule out cloud-free atmospheres in chemical equilibrium assuming either solar abundances (5.4 sigma confidence) or a moderate level of metal enrichment (50x solar abundances, 3.8 sigma), confirming previous results that such models are not representative for cool, low-mass, externally irradiated extrasolar planets. Current measurements are consistent with a flat transmission spectrum, which suggests that the atmosphere is explained by high-altitude clouds and haze, disequilibrium chemistry, unexpected abundance patterns, or the atmosphere is extremely metal-rich (>200x solar). Because GJ 3470b's low bulk density sets an upper limit on the planet's atmospheric enrichment of <300x solar, the atmospheric mean molecular weight must be <9. Thus, if the atmosphere is cloud-free its spectral features should be detectable with future observations. Transit observations at shorter wavelengths will provide the best opportunity to discriminate between plausible scenarios. We obtained optical spectroscopy with the GMOS spectrograph, but these observations exhibit large systematic uncertainties owing to thin, persistent cirrus conditions. Finally, we also provide the first detailed look at the steps necessary for well-calibrated MOSFIRE observations, and provide advice for future observations with this instrument.

Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397.

We present a Monte Carlo model for the structure of low mass (total mass < 25 earth mass) planetary systems that form by the in situ gravitational assembly of planetary embryos into final planets. Our model includes distributions of mass, eccentricity, inclination and period spacing that are based on the simulation of a disk of 20 earth masses, forming planets around a solar mass star, and assuming a power law surface density distribution $\propto a^{-1.5}$. The output of the Monte Carlo model is then subjected to the selection effects that mimic the observations of a transiting planet search such as that performed by the Kepler satellite. The resulting comparison of the output to the properties of the observed sample yields an encouraging agreement in terms of the relative frequencies of multiple planet systems and the distribution of the mutual inclinations, when moderate tidal circularisation is taken into account. The broad features of the period distribution and radius distribution can also be matched within this framework, although the model underpredicts the distribution of small period ratios. This likely indicates that some dissipation is still required in the formation process. The most striking deviation between model and observations is in the ratio of single to multiple systems, in that there are roughly 50% more single planet candidates observed than are produced in any model population. This suggests that some systems must suffer additional attrition to reduce the number of planets or increase the range of inclinations.

The hot Jupiter WASP-12b is one of the largest, hottest, and best-studied extrasolar planets. We revisit our recent analysis of WASP-12b's emission spectrum in light of near-infrared spectroscopic measurements which have been claimed to support either a hydride-dominated or carbon-rich atmospheric composition. We show that this new spectrum is still consistent with a featureless blackbody, indicating a nearly isothermal photosphere on the planet's day side. Thus the ensemble of occultation measurements for WASP-12b is still insufficient to constrain the planet's atmospheric composition.

We present a new measurement of the stellar initial mass function (IMF) based on ultra-deep, high-resolution photometry of >5,000 stars in the outskirts of the Small Magellanic Cloud (SMC) galaxy. The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations reveal this rich, co-spatial population behind the foreground globular cluster 47 Tuc, which we targeted for 121 HST orbits. The stellar main sequence of the SMC is measured in the F606W, F814W color-magnitude diagram (CMD) down to ~30th magnitude, and is cleanly separated from the foreground star cluster population using proper motions. We simulate the SMC population by extracting stellar masses (single and unresolved binaries) from specific IMFs, and converting those masses to luminosities in our bandpasses. The corresponding photometry for these simulated stars is drawn directly from a rich cloud of 4 million artificial stars, thereby accounting for the real photometric scatter and completeness of the data. Over a continuous and well populated mass range of M = 0.37 - 0.93 Msun (i.e., down to a ~75% completeness limit at F606W = 28.7), we demonstrate that the IMF is well represented by a single power-law form with slope \alpha = -1.90 (^{+0.15}_{-0.10}) (3 sigma error) (i.e., dN/dM \propto M^{\alpha}). This is shallower than the Salpeter slope of \alpha = -2.35, which agrees with the observed stellar luminosity function at higher masses. Our results indicate that the IMF does {\it not} turn over to a more shallow power-law form within this mass range. We discuss implications of this result for the theory of star formation, the inferred masses of galaxies, and the (lack of a) variation of the IMF with metallicity.

We revisit the atmospheric properties of the extremely hot Jupiter WASP-12b in light of several new developments. First, new narrowband (2.315 micron) secondary eclipse photometry that we present here, which exhibits a planet/star flux ratio of 0.45% +/- 0.06 %, corresponding to a brightness temperature of 3640 K +/- 230 K; second, recent Spitzer/IRAC and Hubble/WFC3 observations; and third, a recently observed star only 1" from WASP-12, which has diluted previous observations and which we further characterize here. We correct past WASP-12b eclipse measurements for the presence of this object, and we revisit the interpretation of WASP-12b's dilution-corrected emission spectrum. The resulting planetary emission spectrum is well-approximated by a blackbody, and consequently our primary conclusion is that the planet's infrared photosphere is nearly isothermal. Thus secondary eclipse spectroscopy is relatively ill-suited to constrain WASP-12b's atmospheric abundances, and transmission spectroscopy may be necessary to achieve this goal.

We investigated the probability that an inelastic collision of planetesimals within the Hill sphere of the Jovian planets could explain the presence and orbits of observed irregular satellites. Capture of these satellites via this mechanism is highly dependent on not only the mass of the protoplanetary disk, but also the shape of the planetesimal size distribution. We performed 2000 simulations for integrated time intervals $\sim 2$ Myr and found that given the currently accepted value for the minimum mass solar nebula and planetesimal number density based upon the \citet{Bottke2010b} size distribution $dN \sim D^{-1.8} dD$, the collision rates for the different Jovian planets range between $\sim 60$ and $\gtrsim 10^3 \, \Myr^{-1}$ Additionally we found that the probability that these collisions remove enough orbital energy to yield a bound orbit was $\lesssim 10^{-5}$ and had very little dependence on the relative size of the planetesimals. Of these collisions, the collision energy between two masses were $\gtrsim 10^3$ times the gravitational binding energy for masses with radii $\sim 100$ km. Additionally, the combination that enough energy is dissipated by the collision is $\lesssim 10^{-6}$ and the ratio of collision energy to gravitational binding energy is sufficiently large enough to break apart, or completely obliterate, the colliding planetesimals, it is likely that for irregular satellites to be captured as a result of the collision of planetesimals, only a small fraction of the colliding material remained bound to the planet and would most likely be produced by a large mass with a radius $\sim 100$ km colliding with a smaller mass with radius $\lesssim 20$ km.

The spectral distribution of field white dwarfs shows a feature called the "non-DA gap". As defined by Bergeron et al., this is a temperature range (5100K--6100K) where relatively few non-DA stars are found, even though such stars are abundant on either side of the gap. It is usually viewed as an indication that a significant fraction of white dwarfs switch their atmospheric compositions back and forth between hydrogen-rich and helium-rich as they cool. In this paper, we present a Monte Carlo model of the Galactic disk white dwarf population, based on the spectral evolution model of Chen and Hansen. We find that the non-DA gap emerges naturally, even though our model only allows white dwarf atmospheres to evolve monotonically from hydrogen-rich to helium-rich through convective mixing. We conclude by discussing the effects of convective mixing on the white dwarf luminosity function and the use thereof for Cosmochronology.

We demonstrate that the observed distribution of `Hot Neptune'/`Super-Earth' systems is well reproduced by a model in which planet assembly occurs in situ, with no significant migration post-assembly. This is achieved only if the amount of mass in rocky material is $\sim 50$--$100 M_{\oplus}$ interior to 1 AU. Such a reservoir of material implies that significant radial migration of solid material takes place, and that it occur before the stage of final planet assembly. The model not only reproduces the general distribution of mass versus period, but also the detailed statistics of multiple planet systems in the sample. We furthermore demonstrate that cores of this size are also likely to meet the criterion to gravitationally capture gas from the nebula, although accretion is rapidly limited by the opening of gaps in the gas disk. If the mass growth is limited by this tidal truncation, then the scenario sketched here naturally produces Neptune-mass objects with substantial components of both rock and gas, as is observed. The quantitative expectations of this scenario are that most planets in the `Hot Neptune/Super-Earth' class inhabit multiple-planet systems, with characteristic orbital spacings. The model also provides a natural division into gas-rich (Hot Neptune) and gas-poor (Super-Earth) classes at fixed period. The dividing mass ranges from $\sim 3 M_{\oplus}$ at 10 day orbital periods to $\sim 10 M_{\oplus}$ at 100 day orbital periods. For orbital periods $< 10$ days, the division is less clear because a gas atmosphere may be significantly eroded by stellar radiation.

We present a new empirical calibration of equilibrium tidal theory for extrasolar planet systems, extending a prior study by incorporating detailed physical models for the internal structure of planets and host stars. The resulting strength of the stellar tide produces a coupling that is strong enough to reorient the spins of some host stars without causing catastrophic orbital evolution, thereby potentially explaining the observed trend in alignment between stellar spin and planetary orbital angular momentum. By isolating the sample whose spins should not have been altered in this model, we also show evidence for two different processes that contribute to the population of planets with short orbital periods. We apply our results to estimate the remaining lifetimes for short period planets, examine the survival of planets around evolving stars, and determine the limits for circularisation of planets with highly eccentric orbits. Our analysis suggests that the survival of circularised planets is strongly affected by the amount of heat dissipated, which is often large enough to lead to runaway orbital inflation and Roche lobe overflow.

We report the tentative detection of the near-infrared emission of the Hot Jupiter WASP-12b with the low-resolution prism on IRTF/SpeX. We find a K-H contrast color of 0.137% +/- 0.054%, corresponding to a blackbody of temperature 2400 (+1500/-500) K and consistent with previous, photometric observations. We also revisit WASP-12b's energy budget on the basis of secondary eclipse observations: the dayside luminosity is a relatively poorly constrained (2.0-4.3) x 10^30 erg/s, but this still allows us to predict a day/night effective temperature contrast of 200-1,000 K (assuming A_B=0). Thus we conclude that WASP-12b probably does not have both a low albedo and low recirculation efficiency. Our results show the promise and pitfalls of using single-slit spectrographs for characterization of extrasolar planet atmospheres, and we suggest future observing techniques and instruments which could lead to further progress. Limiting systematic effects include the use of a too-narrow slit on one night -- which observers could avoid in the future -- and chromatic slit losses (resulting from the variable size of the seeing disk) and variations in telluric transparency -- which observers cannot control. Single-slit observations of the type we present remain the best option for obtaining lambda > 1.7 micron spectra of transiting exoplanets in the brightest systems. Further and more precise spectroscopy is needed to better understand the atmospheric chemistry, structure, and energetics of this, and other, intensely irradiated planet.

In HST Cycle 17, we imaged the well known globular star cluster 47 Tucanae for 121 orbits using ACS and both the UVIS and IR channels of WFC3 (GO-11677, PI - H. Richer). This unique data set was obtained to address many scientific questions that demand a very deep, panchromatic, and panoramic view of the cluster's stellar populations. In total, the program obtained over 0.75 Ms of imaging exposure time with the three HST cameras, over a time span of 9 months in 2010. The primary ACS field was imaged in the two broadband filters F606W and F814W filters, at 13 orientations, for all 121 orbits. The parallel WFC3 imaging provides a panchromatic (0.4 - 1.7 micron) and contiguous imaging swath over a 250 degree azimuthal range at impact radii of 6.5 -- 17.9 pc in 47 Tuc. This imaging totals over 60 arcmin^2 in area and utilizes the F390W and F606W broadband filters on WFC3/UVIS and the F110W and F160W broadband filters on WFC3/IR. In this paper, we describe the observational design of the new survey and one of the methods used to analyze all of the imaging data. This analysis combines over 700 full-frame images taken with the three HST cameras into a handful of ultra-deep, well-sampled combined images in each of the six filters. The results reveal unprecedented color-magnitude diagrams (CMDs) of the cluster extending to >30th magnitude in the optical, 29th magnitude in the UV, and 27th magnitude in the IR. The data set provides a characterization of the complete stellar populations of 47 Tuc, extending from the faintest hydrogen burning dwarfs through the main-sequence and giant branches, down to very cool white dwarf remnants in the cluster. The imaging also provides the deepest probe of the stellar populations of the background Small Magellanic Cloud (SMC) galaxy, resolving low mass main-sequence dwarfs with M < 0.2 Msun.

We investigated the probability that an inelastic collision of planetesimals within the Hill sphere of the Jovian planets could explain the presence and orbits of observed irregular satellites. Capture of satellites via this mechanism is highly dependent on not only the mass of the protoplanetary disk, but also the shape of the planetesimal size distribution. We performed 2000 simulations for integrated time intervals $\sim 2$ Myr and found that, given the currently accepted value for the minimum mass solar nebula and planetesimal number density based upon the \citet{Nesvorny2003} and \citet{Charnoz2003} size distribution $dN \sim D^{-3.5} dD$, the collision rates for the different Jovian planets range between $\sim 0.6$ and $\gtrsim 170 \, \Myr^{-1}$ for objects with radii, $1 \, \km \le r \le 10 \, \km$. Additionally, we found that the probability that these collisions remove enough orbital energy to yield a bound orbit was $\lesssim 10^{-5}$ and had very little dependence on the relative size of the planetesimals. Of these collisions, the collision energy between two objects was $\gtrsim 10^3$ times the gravitational binding energy for objects with radii $\sim 100$ km. We find that, capturing irregular satellites via collisions between unbound objects can only account for $\sim 0.1%$ of the observed population, hence can this not be the sole method of producing irregular satellites.

The nearby star GJ 1214 hosts a planet intermediate in radius and mass between Earth and Neptune, resulting in some uncertainty as to its nature. We have observed this planet, GJ 1214b, during transit with the high-resolution, near-infrared NIRSPEC spectrograph on the Keck II telescope, in order to characterize the planet's atmosphere. By cross-correlating the spectral changes through transit with a suite of theoretical atmosphere models, we search for variations associated with absorption in the planet atmosphere. Our observations are sufficient to rule out tested model atmospheres with wavelength-dependent transit depth variations >5e-4 over the wavelength range 2.1 - 2.4 micron. Our sensitivity is limited by variable slit loss and telluric transmission effects. We find no positive signatures but successfully rule out a number of plausible atmospheric models, including the default assumption of a gaseous, H-dominated atmosphere in chemical equilibrium. Such an atmosphere can be made consistent if the absorption due to methane is reduced. Clouds can also render such an atmosphere consistent with our observations, but only if they lie higher in the atmosphere than indicated by recent optical and infrared measurements. When taken in concert with constraints from other groups, our results support a consensus model in which the atmosphere of GJ 1214b contains significant H and He, but where methane is depleted. If this depletion is the result of photochemical processes, it may also produce a haze that suppresses spectral features in the optical.

We present theoretical analysis on old white dwarf stars with \emph{some} hydrogen, that possess a mass of surface hydrogen from $1\times10^{-11}M_{\sun}$ to $1\times10^{-7}M_{\sun}$. The evolution of such objects is complicated by convective mixing from surface convection zone to the underlying helium layer. In this paper, we provide first self-consistent, quantitative investigation on the subject of convective mixing. Numerical cooling curves and chemical evolution curves are obtained as a function of white dwarf mass and hydrogen content. Such results will be applied to the investigation of the non-DA gap of \citet{1997ApJS..108..339B} in a later paper.

We provide an 'effective theory' of tidal dissipation in extrasolar planet systems by empirically calibrating a model for the equilibrium tide. The model is valid to high order in eccentricity and parameterised by two constants of bulk dissipation - one for dissipation in the planet and one for dissipation in the host star. We are able to consistently describe the distribution of extrasolar planetary systems in terms of period, eccentricity and mass (with a lower limit of a Saturn mass) with this simple model. Our model is consistent with the survival of short-period exoplanet systems, but not with the circularisation period of equal mass stellar binaries, suggesting that the latter systems experience a higher level of dissipation than exoplanet host stars. Our model is also not consistent with the explanation of inflated planetary radii as resulting from tidal dissipation. The paucity of short period planets around evolved A stars is explained as the result of enhanced tidal inspiral resulting from the increase in stellar radius with evolution.

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

We use signal enhancement techniques and a matched filter analysis to search for the K band spectroscopic absorption signature of the close orbiting extrasolar giant planet, HD 189733b. With timeseries observations taken with NIRSPEC at the Keck II telescope, we investigate the relative abundances of H2O and carbon bearing molecules, which have now been identified in the dayside spectrum of HD 189733b. We detect a candidate planet signature with a low level of significance, close to the ~153 km/s velocity amplitude of HD 189733b. However, some systematic variations, mainly due to imperfect telluric line removal, remain in the residual spectral timeseries in which we search for the planetary signal. The robustness of our candidate signature is assessed, enabling us to conclude that it is not possible to confirm the presence of any planetary signal which appears at Fp/F* contrasts deeper than the 95.4 per cent confidence level. Our search does not enable us to detect the planet at a contrast ratio of Fp/F* = 1/1920 with 99.9 per cent confidence. We also investigate the effect of model uncertainties on our ability to reliably recover a planetary signal. The use of incorrect temperature, model opacity wavelengths and model temperature-pressure profiles have important consequences for the least squares deconvolution procedure that we use to boost the S/N ratio in our spectral timeseries observations. We find that mismatches between the empirical and model planetary spectrum may weaken the significance of a detection by ~30-60 per cent, thereby potentially impairing our ability to recover a planetary signal with high confidence.

Globular star clusters are among the first stellar populations to have formed in the Milky Way, and thus only a small sliver of their initial spectrum of stellar types are still burning hydrogen on the main-sequence today. Almost all of the stars born with more mass than 0.8 M_sun have evolved to form the white dwarf cooling sequence of these systems, and the distribution and properties of these remnants uniquely holds clues related to the nature of the now evolved progenitor stars. With ultra-deep HST imaging observations, rich white dwarf populations of four nearby Milky Way globular clusters have recently been uncovered, and are found to extend an impressive 5 - 8 magnitudes in the faint-blue region of the H-R diagram. In this paper, we characterize the properties of these population II remnants by presenting the first direct mass measurements of individual white dwarfs near the tip of the cooling sequence in the nearest of the Milky Way globulars, M4. Based on Gemini/GMOS and Keck/LRIS multiobject spectroscopic observations, our results indicate that 0.8 M_sun population II main-sequence stars evolving today form 0.53 +/- 0.01 M_sun white dwarfs. We discuss the implications of this result as it relates to our understanding of stellar structure and evolution of population II stars and for the age of the Galactic halo, as measured with white dwarf cooling theory.

We show that the assembly of the Solar System terrestrial planets can be successfully modelled with all of the mass initially confined to a narrow annulus between 0.7 and 1.0 AU. With this configuration, analogues of Mercury and Mars often form from the collisional evolution of material diffusing out of the annulus under the scattering of the forming Earth and Venus analogues. The final systems also possess eccentricities and inclinations that match the observations, without recourse to dynamical friction from remnant small body populations. Finally, the characteristic assembly timescale for Earth analogues is rapid in this model, and consistent with cosmochemical models based on the $^{182}$Hf--$^{182}$W isotopes. The agreement between this model and the observations suggests that terrestrial planet systems may also be formed in `planet traps', as has been proposed recently for the cores of giant planets in our solar system and others.

We model the assembly of planets from planetary embryos under the conditions suggested by various scenarios for the formation of the planetary system around the millisecond pulsar B1257+12. We find that the most likely models fall at the low angular momentum end of the proposed range. Models that invoke supernova fallback produce such disks, although we find that a solar composition disk produces a more likely evolution than one composed primarily of heavy elements. Furthermore, we find that dust sedimentation must occur rapidly as the disk cools, in order that the solid material be confined to a sufficiently narrow range of radii. A quantitative comparison between the observations and the best-fit models shows that the simulations can reproduce the observed eccentricities and masses, but have difficulty reproducing the compactness of the pulsar planet system. Finally, we examine the results of similar studies of solar system terrestrial planet accumulation and discuss what can be learned from the comparison.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

The evidence is mounting that star formation necessarily involves planet formation. We clearly have a vested interest in finding other Earths but a true understanding of planet formation requires completing the census and mapping planetary architecture in all its grandeur and diversity. Here, we show that a 2000-star survey undertaken with SIM Lite will uniquely probe planets around B-A-F stars, bright and binary stars and white dwarfs. In addition, we show that the high precision of SIM Lite allows us to gain unique insights into planet formation via accurate measurements of mutual inclinations.

We present the results of three body simulations focused on understanding the fates of intermediate mass black holes (IBH) that drift within the central 0.5 pc of the Galaxy. In particular, we modeled the interactions between pairs of $4000 {\rm M}_{\odot}$ black holes as they orbit a central blac k hole of mass $4 \times 10^6 {\rm M}_{\odot}$. The simulations performed assume a Schwarzschild geometry and account for Chandrasekhar dynamical friction as well as acceleration resulting from energy lost due to gravitational radiation. We found the branching ratio for one of the orbiting IBHs to merge with the CBH was 0.95 and is independent of the inner IBH's initial eccentricity as well as the rate of sinking. This, coupled with an infall rate of $\sim 10^7 $ yrs for an IBH to drift into the Galactic center, results in an IBH-CBH merger every $\lesssim 11$ Myrs. Lastly we found that the IBH-IBH-CBH triple body system ``resets'' itself, in the sense that a system with an inner I BH with an initially circular orbit generally left behind an IBH with a large eccentricity, whereas a system in which the inner IBH had a high eccentricity ($e_0 \sim 0.9$) usually left a remnant with low eccentricity. Branching ratios for different outcomes are also similar in the two cases.

Turbulence is ubiquitous in Solar System planetary atmospheres. In hot Jupiter atmospheres, the combination of moderately slow rotation and thick pressure scale height may result in dynamical weather structures with unusually large, planetary-size scales. Using equivalent-barotropic, turbulent circulation models, we illustrate how such structures can generate a variety of features in the thermal phase curves of hot Jupiters, including phase shifts and deviations from periodicity. Such features may have been spotted in the recent infrared phase curve of HD 189733b. Despite inherent difficulties with the interpretation of disk-integrated quantities, phase curves promise to offer unique constraints on the nature of the circulation regime present on hot Jupiters.

We describe here the reduction methods that we developed to study the faintest red dwarfs and white dwarfs in an outer field of NGC6397, which was observed by \hst for 126 orbits in 2005. The particular challenge of this data set is that the faintest stars are not readily visible in individual exposures, so special care must be taken to combine the information in all the exposures in order to identify and measure them. Unfortunately, it is hard to find the faintest stars without also finding a large number of faint galaxies, so we developed specialized tools to distinguish between the point-like stars and the barely resolved galaxies. We found that artificial-star tests, while obviously necessary for completeness determination, can also play an important role in helping us optimize our finding and measuring algorithms. Although this paper focuses on this data set specifically, many of the techniques are new and might find application in other work, particularly when a large number of images is available for a single field.

We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.

We present a sequence of toy models for irradiated planet atmospheres, in which the effects of geometry and energy redistribution are modelled self-consistently. We use separate but coupled grey atmosphere models to treat the ingoing stellar irradiation and outgoing planetary reradiation. We investigate how observed quantities such as full phase secondary eclipses and orbital phase curves depend on various important parameters, such as the depth at which irradiation is absorbed and the depth at which energy is redistributed. We also compare our results to the more detailed radiative transfer models in the literature, in order to understand how those map onto the toy model parameter space. Such an approach can prove complementary to more detailed calculations, in that they demonstrate, in a simple way, how the solutions change depending on where, and how, energy redistribution occurs. As an example of the value of such models, we demonstrate how energy redistribution and temperature equilibration at moderate optical depths can lead to temperature inversions in the planetary atmosphere, which may be of some relevance to recent observational findings.

The most commonly accepted explanation for the origin of hypervelocity stars in the halo of the Milky Way is that they are the result of tidal disruption of binaries by the massive black hole at the center of the Galaxy. We show that, if this scenario is correct, and if the original binary properties are similar to those in the local stellar neighbourhood, then the hypervelocity stars should rotate with velocities measureably lower than those for field stars of similar spectral type. This may prove to be a more direct test of the model than trying to predict the position and velocity distributions.

We present results from a set of over 300 pseudospectral simulations of atmospheric circulation on extrasolar giant planets with circular orbits. The simulations are of high enough resolution (up to 341 total and sectoral modes) to resolve small-scale eddies and waves, required for reasonable physical accuracy. In this work, we focus on the global circulation pattern that emerges in a shallow, ``equivalent-barotropic'', turbulent atmosphere on both tidally synchronized and unsynchronized planets. A full exploration of the large physical and numerical parameter-space is performed to identify robust features of the circulation. For some validation, the model is first applied to Solar System giant planets. For extrasolar giant planets with physical parameters similar to HD209458b--a presumably synchronized extrasolar giant planet, representative in many dynamical respects--the circulation is characterized by the following features: 1) a coherent polar vortex that revolves around the pole in each hemisphere; 2) a low number--typically two or three--of slowly-varying, broad zonal (east-west) jets that form when the maximum jet speed is comparable to, or somewhat stronger than, those observed on the planets in the Solar System; and, 3) motion-associated temperature field, whose detectability and variability depend on the strength of the net heating rate and the global root mean square wind speed in the atmosphere. In many ways, the global circulation is Earth-like, rather than Jupiter-like. However, if extrasolar giant planets rotate faster and are not close-in (therefore not synchronized), their circulations become more Jupiter-like, for Jupiter-like rotation rates.

The initial-final mass relation represents a mapping between the mass of a white dwarf remnant and the mass that the hydrogen burning main-sequence star that created it once had. The relation thus far has been constrained using a sample of ~40 stars in young open clusters, ranging in initial mass from ~2.75 -- 7 Msun, and shows a general trend that connects higher mass main-sequence stars with higher mass white dwarfs. In this paper, we present CFHT/CFH12K photometric and Keck/LRIS multiobject spectroscopic observations of a sample of 22 white dwarfs in two older open clusters, NGC 7789 (t = 1.4 Gyr) and NGC 6819 (t = 2.5 Gyr). We measure masses for the highest S/N spectra by fitting the Balmer lines to atmosphere models and place the first direct constraints on the low mass end of the initial-final mass relation. Our results indicate that the observed general trend at higher masses continues down to low masses, with M_initial = 1.16 Msun main-sequence stars forming M_final = 0.53 Msun white dwarfs (including our data from the very old open cluster, NGC 6791). This extention of the relation represents a four fold increase in the total number of hydrogen burning stars for which the integrated mass loss can now be calculated, assuming a Salpeter initial mass function. The new leverage at the low mass end is used to derive a purely empirical initial-final mass relation without the need for any indirectly derived anchor points. The sample of white dwarfs in these clusters also shows several very interesting systems that we discuss further: a DB (helium atmosphere) white dwarf, a magnetic white dwarf, a DAB (mixed hydrogen/helium atmosphere or a double degenerate DA+DB) white dwarf(s), and two possible equal mass DA double degenerate binary systems.

We present the first detailed study of the properties (temperatures, gravities, and masses) of the NGC 6791 white dwarf population. This unique stellar system is both one of the oldest (8 Gyr) and most metal-rich ([Fe/H] ~ 0.4) open clusters in our Galaxy, and has a color-magnitude diagram (CMD) that exhibits both a red giant clump and a much hotter extreme horizontal branch. Fitting the Balmer lines of the white dwarfs in the cluster, using Keck/LRIS spectra, suggests that most of these stars are undermassive, <M> = 0.43 +/- 0.06 Msun, and therefore could not have formed from canonical stellar evolution involving the helium flash at the tip of the red giant branch. We show that at least 40% of NGC 6791's evolved stars must have lost enough mass on the red giant branch to avoid the flash, and therefore did not convert helium into carbon-oxygen in their core. Such increased mass loss in the evolution of the progenitors of these stars is consistent with the presence of the extreme horizontal branch in the CMD. This unique stellar evolutionary channel also naturally explains the recent finding of a very young age (2.4 Gyr) for NGC 6791 from white dwarf cooling theory; helium core white dwarfs in this cluster will cool ~3 times slower than carbon-oxygen core stars and therefore the corrected white dwarf cooling age is in fact ~7 Gyr, consistent with the well measured main-sequence turnoff age. These results provide direct empirical evidence that mass loss is much more efficient in high metallicity environments and therefore may be critical in interpreting the ultraviolet upturn in elliptical galaxies.

We present the CMD from deep HST imaging in the globular cluster NGC 6397. The ACS was used for 126 orbits to image a single field in two colors (F814W, F606W) 5 arcmin SE of the cluster center. The field observed overlaps that of archival WFPC2 data from 1994 and 1997 which were used to proper motion (PM) clean the data. Applying the PM corrections produces a remarkably clean CMD which reveals a number of features never seen before in a globular cluster CMD. In our field, the main sequence stars appeared to terminate close to the location in the CMD of the hydrogen-burning limit predicted by two independent sets of stellar evolution models. The faintest observed main sequence stars are about a magnitude fainter than the least luminous metal-poor field halo stars known, suggesting that the lowest luminosity halo stars still await discovery. At the bright end the data extend beyond the main sequence turnoff to well up the giant branch. A populous white dwarf cooling sequence is also seen in the cluster CMD. The most dramatic features of the cooling sequence are its turn to the blue at faint magnitudes as well as an apparent truncation near F814W = 28. The cluster luminosity and mass functions were derived, stretching from the turn off down to the hydrogen-burning limit. It was well modeled with either a very flat power-law or a lognormal function. In order to interpret these fits more fully we compared them with similar functions in the cluster core and with a full N-body model of NGC 6397 finding satisfactory agreement between the model predictions and the data. This exercise demonstrates the important role and the effect that dynamics has played in altering the cluster IMF.

We report the relative abundances of 17 elements in the atmosphere of the white dwarf star GD 362, material that, very probably, was contained previously in a large asteroid or asteroids with composition similar to the Earth/Moon system. The asteroid may have once been part of a larger parent body not unlike one of the terrestrial planets of our solar system.

Hydrogen atmosphere white dwarfs with metal lines, so-called DAZs, require external accretion of material to explain the presence of weak metal line absorption in their photospheres. The source of this material is currently unknown, but could come from the interstellar medium, unseen companions, or relic planetesimals from asteroid belt or Kuiper belt analogues. Accurate mid-infrared photometry of these white dwarfs provide additional information to solve the mystery of this accretion and to look for evidence of planetary systems that have survived post main sequence evolution. We present {\em Spitzer} IRAC photometry accurate to $\sim$3% for four DAZs and one DA with circumstellar absorption lines in the UV. We search for excesses due to unseen companions or circumstellar dust disks. We use {\em Hubble Space Telescope} NICMOS imaging of these white dwarfs to gauge the level of background contamination to our targets as well as rule out common proper motion companions to WD 1620-391. All of our targets show no excesses due to companions $>$20 M$_{J}$, ruling out all but very low mass companions to these white dwarfs at all separations. No excesses due to circumstellar disks are observed, and we place limits on what types of disks may still be present.

We identify two classes of transiting planet, based on their equilibrium temperatures and Safronov numbers. We examine various possible explanations for the dichotomy. It may reflect the influence of planet or planetesimal scattering in determining when planetary migration stops. Another possibility is that some planets lose more mass to evaporation than others. If this evaporation process preferentially removes Helium from the planet, the consequent reduction in the mean molecular weight may explain why some planets have anomalously large radii.

We model the evolution of protoplanetary disks surrounding millisecond pulsars, using PSR 1257+12 as a test case. Initial conditions were chosen to correspond to initial angular momenta expected for supernova-fallback disks and disks formed from the tidal disruption of a companion star. Models were run under two models for the viscous evolution of disks: fully viscous and layered accretion disk models. Supernova-fallback disks result in a distribution of solids confined to within 1-2 AU and produce the requisite material to form the three known planets surrounding PSR 1257+12. Tidal disruption disks tend to slightly underproduce solids interior to 1 AU, required for forming the pulsar planets, while overproducing the amount of solids where no body, lunar mass or greater, exists. Disks evolving under 'layered' accretion spread somewhat less and deposit a higher column density of solids into the disk. In all cases, circumpulsar gas dissipates on $\lesssim 10^{5}$ year timescales, making formation of gas giant planets highly unlikely.

Physical conditions in the atmospheres of tidally-locked, slowly-rotating hot Jupiters correspond to dynamical circulation regimes with Rhines scales and Rossby deformation radii comparable to the planetary radii. Consequently, the large spatial scales of moving atmospheric structures could generate significant photospheric variability. Here, we estimate the level of thermal infrared variability expected in successive secondary eclipse depths, according to hot Jupiter turbulent ``shallow-layer'' models. The variability, at the few percent level or more in models with strong enough winds, is within the reach of Spitzer measurements. Eclipse depth variability is thus a valuable tool to constrain the circulation regime and global wind speeds in hot Jupiter atmospheres.

Recent Spitzer infrared measurements of hot Jupiter eclipses suggest that eclipse mapping techniques could be used to spatially resolve the day-side photospheric emission of these planets using partial occultations. As a first step in this direction, we simulate ingress/egress lightcurves for the three brightest known eclipsing hot Jupiters and evaluate the degree to which parameterized photospheric emission models can be distinguished from each other with repeated, noisy eclipse measurements. We find that the photometric accuracy of Spitzer is insufficient to use this tool effectively. On the other hand, the level of photospheric details that could be probed with a few JWST eclipse measurements could greatly inform hot Jupiter atmospheric modeling efforts. A JWST program focused on non-parametric eclipse map inversions for hot Jupiters should be actively considered.

We present the results of a deep Hubble Space Telescope (HST) exposure of the nearby globular cluster NGC6397, focussing attention on the cluster's white dwarf cooling sequence. This sequence is shown to extend over 5 magnitudes in depth, with an apparent cutoff at magnitude F814W=27.6. We demonstrate, using both artificial star tests and the detectability of background galaxies at fainter magnitudes, that the cutoff is real and represents the truncation of the white dwarf luminosity function in this cluster. We perform a detailed comparison between cooling models and the observed distribution of white dwarfs in colour and magnitude, taking into account uncertainties in distance, extinction, white dwarf mass, progenitor lifetimes, binarity and cooling model uncertainties. After marginalising over these variables, we obtain values for the cluster distance modulus and age of \mu_0 = 12.02 \pm 0.06 and T_c = 11.47 \pm 0.47Gyr (95% confidence limits). Our inferred distance and white dwarf initial-final mass relations are in good agreement with other independent determinations, and the cluster age is consistent with, but more precise than, prior determinations made using the main sequence turnoff method. In particular, within the context of the currently accepted \Lambda CDM cosmological model, this age places the formation of NGC6397 at a redshift z=3, at a time when the cosmological star formation rate was approaching its peak.

NGC 6397 is the second closest globular star cluster to the Sun. Using 5 days of time on the Hubble Space Telescope, we have constructed the deepest ever color-magnitude diagram for this cluster. We see a clear truncation in each of its two major stellar sequences. Faint red main sequence stars run out well above our observational limit and near to the theoretical prediction for the lowest mass stars capable of stable hydrogen-burning in their cores. We also see a truncation in the number counts of faint blue stars, namely white dwarfs. This reflects the limit to which the bulk of the white dwarfs can cool over the lifetime of the cluster. There is also a turn towards bluer colors in the least luminous of these objects. This was predicted for the very coolest white dwarfs with hydrogen-rich atmospheres as the formation of H2 causes their atmospheres to become largely opaque to infrared radiation due to collision-induced absorption.

As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.

The star upsilon Andromeda is orbited by three known planets, the innermost of which has an orbital period of 4.617 days and a mass at least 0.69 that of Jupiter. This planet is close enough to its host star that the radiation it absorbs overwhelms its internal heat losses. Here we present the 24 micron light curve of this system, obtained with the Spitzer Space Telescope. It shows a clear variation in phase with the orbital motion of the innermost planet. This is the first demonstration that such planets possess distinct hot substellar (day) and cold antistellar (night) faces.

We present results of N-body simulations aimed at understanding the dynamics of young stars near the Galactic center. Specifically, we model the inspiral of a cluster core containing an intermediate mass black hole and $N \sim 50$ cluster stars in the gravitational potential of a supermassive black hole. We first study the elliptic three-body problem to isolate issues of tidal stripping and subsequent scattering, followed by full N-body simulations to treat the internal dynamics consistently. We find that our simulations reproduce several dynamical features of the observed population. These include the observed inner edge of the claimed clockwise disk, as well as the thickness of said disk. We find that high density clumps, such as that claimed for IRS13E, also result generically from our simulations. However, not all features of the observations are reproduced. In particular, the surface density profile of the simulated disk scales as $\Sigma \propto r^{-0.75}$, which is considerably shallower than that observed. Further, at no point is any significant counter-rotating population formed.

We examine 14 hot white dwarfs for signs of infrared excess using the Spitzer Space Telescope. Twelve of the objects are massive white dwarfs which have been suggested to be the result of binary mergers. The remaining two objects are undermassive white dwarfs which again may be the result of mergers or the inspiral of a substellar companion. In no case do we find any evidence for significant infrared excesses out to wavelengths of 8 microns. This places strong constraints on the presence of orbiting dust and weaker constraints on the presence of close substellar companions.

We propose that the anomalously bright white dwarf luminosity function observed in NGC 6791 (Bedin et al 2005) is the consequence of the formation of 0.5 Msun white dwarfs with Helium cores instead of Carbon cores. This may happen if mass loss during the ascent of the Red Giant Branch is strong enough to prevent a star from reaching the Helium flash. Such a model can explain the slower white dwarf cooling (relative to standard models) and fits naturally with scenarios advanced to explain Extreme Horizontal Branch stars, a population of which are also found in this cluster.

We discuss atmosphere models of HD209458b in light of the recent day-side flux measurement of HD209458b's secondary eclipse by Spitzer-MIPS at 24 microns. In addition, we present a revised secondary eclipse IRTF upper limit at 2.2 microns which places a stringent constraint on the adjacent H2O absorption band depths. These two measurements are complementary because they are both shaped by H2O absorption and because the former is on the Wien tail of the planet's thermal emission spectrum and the latter is near the thermal emission peak. A wide range of models fit the observational data, confirming our basic understanding of hot Jupiter atmospheric physics. Although a range of models are viable, some models at the hot and cold end of the plausible temperature range can be ruled out. One class of previously unconsidered hot Jupiter atmospheric models that fit the data are those with C/O >~ 1 (as Jupiter may have), which have a significant paucity of H2O compared to solar abundance models with C/O = 0.5. The models indicate that HD209458b is in a situation intermediate between pure in situ reradiation and very efficient redistribution of heat; one which will require a careful treatment of atmospheric circulation. We discuss how future wavelength-dependent and phase-dependent observations will further constrain the atmospheric circulation regime. In the shorter term, additional planned measurements for HD209458b, especially Spitzer IRAC photometry, should lift many of the model degeneracies. Multiwavelength IR observations constrain the atmospheric structure and circulation properties of hot Jupiters and thus open a new chapter in quantitative extrasolar planetology.

We present new observations of very faint white dwarfs (WDs) in the rich open star cluster NGC 2099 (M37). Following deep, wide field imaging of the cluster using CFHT, we have now obtained spectroscopic observations of candidate WDs using both GMOS on Gemini and LRIS on Keck. Of our 24 WD candidates (all fainter than V = 22.4), 21 are spectroscopically confirmed to be bona fide WDs, 4-5 of which are most likely field objects. Fitting 18 of the 21 WD spectra with model atmospheres, we find that most WDs in this cluster are quite massive (0.7-0.9 Msun), as expected given the cluster's young age (650 Myr) and, hence, high turnoff mass (~2.4 Msun). We determine a new initial-final mass relationship and almost double the number of existing data points from previous studies. The results indicate that stars with initial masses between 2.8 and 3.4 Msun lose 70-75% of their mass through stellar evolution. For the first time, we find some evidence of a metallicity dependence on the initial-final mass relationship.

Spectra have been obtained of 21 white dwarfs (WDs) in the direction of the young, rich open star cluster NGC 2099. This represents an appreciable fraction (>30%) of the cluster's total WD population. The mean derived mass of the sample is 0.8 Msun - about 0.2 Msun larger than the mean seen among field WDs. A surprising result is that all of the NGC 2099 WDs have hydrogen-rich atmospheres (DAs) and none exhibit helium-rich ones (DBs), or any other spectral class. The number ratio in the field at the temperatures of the NGC 2099 WDs is DA/DB ~ 3.5. While the probability of seeing no DB WDs in NGC 2099 solely by chance is ~2%, if we include WDs in other open clusters of similar age it then becomes highly unlikely that the dearth of DB WDs in young open clusters is just a statistical fluctuation. We explore possible reasons for the lack of DBs in these clusters and conclude that the most promising scenario for the DA/DB number ratio discrepancy in young clusters is that hot, high-mass WDs do not develop large enough helium convection zones to allow helium to be brought to the surface and turn a hydrogen-rich WD into a helium-rich one.

We respond to the paper by De Marchi, Paresce, Straniero and Moroni (2003) on the white dwarf cooling age of M4. The authors question the data analysis and interpretation that led to the conclusions in Hansen et al. (2002). In their paper, De Marchi et al. are unable to obtain photometry as deep as ours from the same data set and therefore assert that only a lower limit to the white dwarf cooling age for this cluster of approximately 9 Gyr can be obtained. In this short contribution we show that shortcomings in the data analysis and reduction techniques of De Marchi et al. are responsible for their inability to reach the photometry limits that our study reports. In a forthcoming paper in which the complete techniques for age determination with white dwarfs are laid out, we demonstrate that their method of fitting the luminosity function gives a spuriously low white dwarf cooling age.

We investigate in detail the white dwarf cooling sequence of the globular cluster Messier 4. In particular we study the influence of various systematic uncertainties, both observational and theoretical, on the determination of the cluster age from the white dwarf cooling sequence. These include uncertainties in the distance to the cluster and the extinction along the line of sight, as well as the white dwarf mass, envelope and core compositions and the white dwarf --main sequence mass relation. We find that fitting to the full two-dimensional colour-magnitude diagram offers a more robust method for age determination than the traditional method of fitting the one-dimensional white dwarf luminosity function. After taking into account the various uncertainties, we find a best fit age of 12.1 Gyr, with a 95% lower limit of 10.3 Gyr. We also perform fits using two other sets of cooling models from the literature. The models of Chabrier et al (2000) yield an encouragingly similar result, although the models of Salaris et al (2000) do not provide as good a fit. Our results support our previous determination of a delay between the formation of the Galactic halo and the onset of star formation in the Galactic disk.

This paper presents new results from a photometric study of the main-sequence stars in M4, the closest globular cluster to the sun. Multi-field, multi-epoch observations at approximately 1, 2, and 6 core-radii were obtained with the WFPC2 camera on the HST. The multi-epoch observations allowed us to clean the data on the basis of proper motion. In all the fields the cluster main sequence can be traced to at least V = 27.0 but there remains a trail of stars to the limit of the data near V = 30 in the deepest outer field. There is no evidence that we have reached the end of the hydrogen burning main sequence in any of our fields, however, there is some indication that very few stars remain to be detected in the deepest data. The scatter about the cluster main sequence yields a surprisingly small and variable binary fraction; f_b = 2% in the inner parts of the cluster falling to the 1% range outside. For the currently visible main sequence stars, the cluster mass function is very flat (alpha =0.1) in the outer field and flattens further in the inner fields suggesting well developed mass segregation. The observed variation in the mass function is broadly consistent with isotropic, multi-mass Michie-King models. Because we have a large sample of white dwarfs in the outer field we are able to show that the cluster IMF above 0.8Msun was considerably steeper than the present day mass function for low mass stars. Two appendicies are included in this contribution; a detailed discussion of the techniques used to reduce the data and a direct comparison between the cluster stars and those belonging to the inner spheroid of the Galaxy yielding a Galactocentric distance of 7.5+/-0.6 kpc.

In this paper we present new and improved data on 38 cool white dwarfs identified by Oppenheimer et al. 2001 (OHDHS) as candidate dark halo objects. Using the high-res spectra obtained with LRIS, we measure radial velocities for 13 WDs that show an H alpha line. We show that the knowledge of RVs decreases the UV-plane velocities by only 6%. The radial velocity sample has a W-velocity dispersion of sig_W = 59 km/s--in between the values associated with the thick disk and the stellar halo. We also see indications for the presence of two populations by analyzing the velocities in the UV plane. In addition, we present CCD photometry for half of the sample, and with it recalibrate the photographic photometry of the remaining WDs. Using the new photometry in standard bands, and by applying the appropriate color-magnitude relations for H and He atmospheres, we obtain new distance estimates. New distances of the WDs that were not originally selected as halo candidates yield 13 new candidates. On average, new distances produce velocities in the UV plane that are larger by 10%, with already fast objects gaining more. Using the new data, while applying the same UV-velocity cut (94 km/s) as in OHDHS, we find a density of cool WDs of 1.7e-4 pc^-3, confirming the value of OHDHS. In addition, we derive the density as a function of the UV-velocity cutoff. The density (corrected for losses due to higher UV cuts) starts to flatten out at 150 km/s (0.4e-4 pc^-3), and is minimized (thus minimizing a possible non-halo contamination) at 190 km/s (0.3e-4 pc^-3). These densities are in a rough agreement with the estimates for the stellar halo WDs, corresponding to a factor of 1.9 and 1.4 higher values.

We establish an extragalactic, zero-motion frame of reference within the deepest optical image of a globular star cluster, an HST 123-orbit exposure of M4 (GO 8679, cycle 9). The line of sight beyond M4 (l,b (deg) = 351,16) intersects the inner halo (spheroid) of our Galaxy at a tangent-point distance of 7.6 kpc (for Ro = 8 kpc). We isolate these spheroid stars from the cluster based on their proper motions over the 6-year baseline between these and previous epoch HST data (GO 5461, cycle 4). Distant background galaxies are also found on the same sight line using image-morphology techniques. This fixed reference frame allows us to independently determine the fundamental Galactic constant, Vo/Ro = 25.3 +/- 2.6 km/s/kpc, thus providing a velocity of the Local Standard of Rest, v = 202.7 +/- 24.7 km/s for Ro = 8.0 +/- 0.5 kpc. Secondly, the galaxies allow a direct measurement of M4's absolute proper motion, mu_total = 22.57 +/- 0.76 mas/yr, in excellent agreement with recent studies. The clear separation of galaxies from stars in these deep data also allow us to search for inner-halo white dwarfs. We model the conventional Galactic contributions of white dwarfs along our line of sight and predict 7.9 (thin disk), 6.3 (thick disk) and 2.2 (spheroid) objects to the limiting magnitude at which we can clearly delineate stars from galaxies (V = 29). An additional 2.5 objects are expected from a 20% white dwarf dark halo consisting of 0.5 Mo objects, 70% of which are of the DA type. After considering the kinematics and morphology of the objects in our data set, we find the number of white dwarfs to be consistent with the predictions for each of the conventional populations. However, we do not find any evidence for dark halo white dwarfs.

Time-series data taken with the Hubble Space Telescope of a field six core radii (~5') from the center of the globular cluster Messier 4, and covering a period of about 10 weeks in early 2001, have been analyzed in search of variable objects. Various criteria were employed to select candidate variable stars. Period searches were performed on the selected candidates using phase dispersion minimization (PDM). The reliability of the PDM search results was tested using synthetic light curves of eclipsing binary stars and sinusoidal light curves of different periods. Results from this analysis showed that there are probably no eclipsing binary stars or periodic variables in our field with periods on the order of a few hours to a few days, down to limiting magnitudes of V~25 and I~24, which is consistent with the absence of contact binaries such as W Ursae Majoris systems. However, one candidate variable star does show an increase in brightness of ~0.1 magnitudes in both bandpasses, which seems to last for a few days. Possible explanations concerning the nature of this object include a binary system with a white dwarf primary and a low-mass main sequence secondary, or a BY Draconis variable star. We are able to set an upper limit to the observed fraction of photometric variability in this dataset of 0.05%.

The pulsar B1620-26 has two companions, one of stellar mass and one of planetary mass. We detected the stellar companion using Hubble Space Telescope observations. The color and magnitude of the stellar companion indicate that it is an undermassive white dwarf ($0.34 \pm 0.04 M_{\odot}$) of age $480\pm 140$ Myr. This places a constraint on the recent history of this triple system and supports a scenario in which the current configuration arose through a dynamical exchange interaction in the cluster core. This implies that planets may be relatively common in low-metallicity globular clusters and that planet formation is more widespread and happened earlier than previously believed.

Deep infra-red observations and long-term monitoring programs have provided dynamical evidence for a supermassive black hole of mass 3.e6 solar masses associated with the radio source Sagitarrius A* at the center of our Galaxy. The brightest stars orbiting within 0.1 parsecs of the black hole appear to be young, massive main sequence stars, n spite of an environment near the black hole that is hostile to star formation. We discuss mechanisms by which stars born outside the central parsec can sink towards the black hole and conclude that the drag coming from plausible stellar populations does not operate on the short timescales required by the stellar ages. We propose that these stars were dragged in by a second black hole of mass of 1.e3-1.e4 solar masses, which would be classified as an intermediate-mass black hole. We discuss the implications for the stellar populations and the kinematics in the Galactic center. Finally we note that continued astrometric monitoring of the central radio source offers the prospect for a direct detection of such objects.

High-resolution atmospheric flow simulations of the tidally-locked extrasolar giant planet, HD 209458b, show large-scale spatio-temporal variability. This is in contrast to the simple, permanent day/night (i.e., hot/cold) picture. The planet's global circulation is characterized by a polar vortex in motion around each pole and a banded structure corresponding to ~3 broad zonal (east-west) jets. For very strong jets, the circulation-induced temperature difference between moving hot and cold regions can reach up to ~1000 K, suggesting that atmospheric variability could be observed in the planet's spectral and photometric signatures.

Shallow-water numerical simulations show that the atmospheric circulation of the close-in extrasolar giant planet (EGP) HD 209458b is characterized by moving circumpolar vortices and few bands/jets (in contrast with ~10 bands/jets and absence of polar vortices on cloud-top Jupiter and Saturn). The large spatial scales of moving circulation structures on HD 209458b may generate detectable variability of the planet's atmospheric signatures. In this Letter, we generalize these results to other close-in EGPs, by noting that shallow-water dynamics is essentially specified by the values of the Rossby (Ro) and Burger (Bu) dimensionless numbers. The range of likely values of Ro (~ 0.01 - 10) and Bu (~ 1 - 200) for the atmospheric flow of known close-in EGPs indicates that their circulation should be qualitatively similar to that of HD 209458b. This results mostly from the slow rotation of these tidally-synchronized planets.

The deepest optical image ever in a globular star cluster, a Hubble Space Telescope (HST) 123 orbit exposure in a single field of Messier 4, was obtained in 2 filters (F606W, F814W) over a 10 week period in early 2001. A somewhat shallower image obtained in 1995 allowed us to select out cluster and field objects via their proper motion displacement resulting in remarkably clean color-magnitude diagrams that reach to V = 30, I = 28. The cluster main sequence luminosity function contains very few stars fainter than MV = 15.0, MI = 11.8 which, in both filters, is more than 2 magnitudes brighter than our limit. This is about the faintest luminosity seen among field Population II subdwarfs of the same metallicity. However, there remains a sprinkling of potential cluster stars to lower luminosity all the way down to our limiting magnitudes. These latter objects are significantly redder than any known metal-poor field subdwarf. Comparison with the current generation of theoretical stellar models implies that the masses of the lowest luminosity cluster stars observed are near 0.09 Msun. We derive the mass function of the cluster in our field and find that it is very slowly rising towards the lowest masses with no convincing evidence of a turnover even below 0.1 Msun. The formal slope between 0.65 and 0.09 Msun is alpha = 0.75 (Salpeter of 2.35) with a 99% confidence interval 0.55 - 1.05. A consistency check between these slopes and the number of observed cluster white dwarfs yields a range of possible conclusions, one of which is that we have indeed seen the termination of the white dwarf cooling sequence in M4.

I examine the hypothesis that many of the high velocity white dwarfs observed by Oppenheimer et al. (2001) are the remnants of donor stars from binaries that produced type Ia supernovae via the `single degenerate' channel. If this channel is a significant contributor to the Galactic SN Ia supernova rate, then the local density of such remnants with V>100 km/s could be as high as 2.e-4 per cubic parsec, comparable to the densities found by Oppenheimer et al. This white dwarf population differs from others in that it is composed exclusively of single stars.

We examine the theoretical implications of a population of low-mass helium-core white dwarfs in globular clusters. In particular, we focus on the observed population in the core of NGC 6397, where several low-mass white dwarf canditates have been identified as ``non-flickerers'' by Cool and collaborators. Age and mass estimates from cooling models, combined with dynamical and evolutionary considerations, lead us to infer that the dark binary companions are C/O white dwarfs rather than neutron stars. Furthermore, we find that the progenitor binaries very likely underwent an exchange interaction within the last 1.e9 years. We examine the prospects for detecting a similar population in other globular clusters, with particular attention to the case of 47 Tuc.

We present the white dwarf sequence of the globular cluster M4, based on a 123 orbit Hubble Space Telescope exposure, with limiting magnitude V = 30, I = 28. The white dwarf luminosity function rises sharply for I >25.5, consistent with the behaviour expected for a burst population. The white dwarfs of M4 extend to approximately 2.5 magnitudes fainter than the peak of the local Galactic disk white dwarf luminosity function. This demonstrates a clear and significant age difference between the Galactic disk and the halo globular cluster M4. Using the same standard white dwarf models (Hansen 1999) to fit each luminosity function yields ages of 7.3 +/- 1.5 Gyr for the disk and 12.7 +/- 0.7 Gyr for M4 (2-sigma statistical errors).

The discovery of a population of high proper motion white dwarfs by Oppenheimer et al (2001) has caused a lot of speculation as to the origin of these stars. I show that the age distribution of the white dwarfs offers a kind of sanity check in these discussions. In particular, this population appears to have a similar age distribution to those in the standard, thin disk white dwarf population. This is not what is expected for either the halo or thick disk, which are thought to be old populations. It may indicate a different dynamical origin entirely, or it may be just be the tail of a larger distribution, implying a significantly larger total density in even older white dwarfs.

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.

We study spectroscopically determined iron abundances of 642 solar-type stars to search for the signature of accreted iron-rich material. We find that the metallicity [Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte Carlo models we find that, on average, these stars have accreted about 0.4 Earth masses of iron while on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, which are slightly evolved and whose convection zones are significantly more massive, have lower average [Fe/H], and their metallicity shows no clear variation with stellar mass. These findings suggest that terrestrial-type material is common around solar type stars.

I describe models for the formation of planetary systems surrounding the remnants of stellar mergers and collisions. I focus primarily on models for the viscous evolution of disks suitable for the formation of the planets surrounding the pulsar B1257+12. I show that the adaptation of models for traditional protoplanetary disks which invoke quiescent or `dead' zones are quite successful in producing disks appropriate for the formation of the pulsar planets. I also briefly describe some even more exotic possibilities that may arise from compact object mergers.

A popular interpretation of recent microlensing studies of the line of sight towards the Large Magellanic Cloud invokes a population of old white dwarf stars in the Galactic halo. Below I review the basic properties of old white dwarf stars and the ongoing efforts to detect this population directly.

Planet formation is accompanied by the formation of comet clouds. In systems where planets migrate on rapid timescales, the diffusive evolution of comet orbits may stall, resulting in a comet cloud intermediate between a flattened Kuiper Belt and a spherical Oort cloud. These `failed Oort clouds' may provide a `smoking gun', indicating that planetary migration has taken place. If some fraction of the scattered component consists of planetary embryos, it may be possible to observe transits of such bodies even when the planetary system is not edge-on to the line of sight.

We describe the possible electromagnetic signals expected from the magnetospheric interactions of a neutron star binary prior to merger. We find that both radio and X-ray signals of detectable strength are possible. We discuss possible links with the phenomenon of gamma-ray bursts and describe the prospects for direct detection of these signals in searches for radio and X-ray transients.

Recent discoveries of extrasolar planets at small orbital radii, or with significant eccentricities, indicate that interactions between massive planets and the disks of gas and dust from which they formed are vital for determining the final shape of planetary systems. We show that if this interaction occurs at an early epoch, when the protoplanetary disc was still massive, then rapid planet growth through accretion causes an otherwise stable disc to fragment into additional planetary mass bodies when the planetary mass reaches 4-5 Jupiter masses. We suggest that such catastrophic planet formation could account for apparent differences in the mass function of massive planets and brown dwarfs, and the existence of young stars that appear to have dissipated their discs at an early epoch. Subsequent gravitational interactions will lead to planetary systems comprising a small number of massive planets in eccentric orbits.

Using a new grid of models of cooling white dwarfs, we calculate isochrones and luminosity functions in the Johnson-Kron/Cousins and HST filter sets for systems containing old white dwarfs. These new models incorporate a non-grey atmosphere which is necessary to properly describe the effects of molecular opacity at the cool temperatures of old white dwarfs. The various functions calculated and extensively tabulated and plotted are meant to be as utilitarian as possible for observers so all results are listed in quantities that observers will obtain. The tables and plots developed should eventually prove critical in interpreting the results of HST's Advanced Camera observations of the oldest white dwarfs in nearby globular clusters, in understanding the results of searches for old white dwarfs in the Galactic halo, and in determining ages for star clusters of all ages using white dwarfs. As a practical application we demonstrate the use of these results by deriving the white dwarf cooling age of the old Galactic cluster M67.

I present a scenario for the production of low mass, degenerate dwarfs of mass $>0.1 M_{\odot}$ via the mechanism of Lenzuni, Chernoff & Salpeter (1992). Such objects meet the mass limit requirements for halo dark matter from microlensing surveys while circumventing the chemical evolution constraints on normal white dwarf stars. I describe methods to observationally constrain this scenario and suggest that such objects may originate in small clusters formed from the thermal instability of shocked, heated gas in dark matter haloes, such as suggested by Fall & Rees (1985) for globular clusters.

We present new white dwarf cooling models which incorporate an accurate outer boundary condition based on new opacity and detailed radiative transfer calculations. We find that helium atmosphere dwarfs cool considerably faster than has previously been claimed, while old hydrogen atmosphere dwarfs will deviate significantly from black body appearance. We use our new models to derive age limits for the Galactic disk. We find that the Liebert, Dahn & Monet (1988) luminosity function yields an age of only 6 Gyr if it is complete to stated limits. However, age estimates of individual dwarfs and the luminosity function of Oswalt et al (1995) are both consistent with disk ages as large as \sim 11 Gyr. We have also used our models to place constraints on white dwarf dark matter in Galactic halos. We find that previous attempts using inadequate cooling models were too severe and that direct detection limits allow a halo that is 11 Gyr old. If the halo is composed solely of helium atmosphere dwarfs, the lower age limit is only 7.5 Gyr. We also demonstrate the importance of studying the cooling sequences of white dwarfs in Globular clusters.

Under the hypothesis that some fraction of massive stellar core collapses give rise to unusually energetic events, termed hypernovae, I examine the required rates assuming some fraction of such events yield gamma ray bursts. I then discuss evidence from studies of pulsars and r-process nucleosynthesis that independently suggests the existence of a class of unusually energetic events. Finally I describe a scenario which links these different lines of evidence as supporting the hypernova hypothesis.

The recent discovery of microlensing of stars in the Large Magellanic Cloud has excited much interest in the nature of the lensing population. Detailed analyses indicate that the mass of these objects ranges from 0.3-0.8 solar masses, suggesting that they might be white dwarfs, the faint remnants of stellar evolution. The confirmation of such an hypothesis would yield profound insights into the early history of our galaxy and the early generations of stars in the universe. Previous attempts have been made to place theoretical constraints on this scenario, but were unduly pessimistic because they relied on inadequate evolutionary models. Here we present the first results from detailed evolutionary models appropriate for the study of white dwarfs of truly cosmological vintage. We find that the commonly held notion that old white dwarfs are red to hold only for helium atmosphere dwarfs and that hydrogen atmosphere dwarfs will be blue, with colours similar to those of the faint point sources found in the Hubble Deep Field. Thus, any direct observational search for the microlensing population should search for faint blue objects rather than faint red ones.

We propose that the cosmological gamma ray bursts arise from the collapse of neutron stars to black holes triggered by collisions or mergers with main sequence stars. This scenario represents a cosmological history qualitatively different from most previous theories because it contains a significant contribution from an old stellar population, namely the globular clusters. Furthermore, the gas poor central regions of globular clusters provide an ideal environment for the generation of the recently confirmed afterglows via the fireball scenario. Collisions resulting from neutron star birth kicks in close binaries may also contribute to the overall rate and should lead to associations between some gamma ray bursts and supernovae of type Ib/c.

A planet orbiting in a disk of planetesimals can experience an instability in which it migrates to smaller orbital radii. Resonant interactions between the planet and planetesimals remove angular momentum from the planetesimals, increasing their eccentricities. Subsequently, the planetesimals either collide with or are ejected by the planet, reducing the semimajor axis of the planet. If the surface density of planetesimals exceeds a critical value, corresponding to 0.03 solar masses of gas inside the orbit of Jupiter, the planet will migrate inward a large distance. This instability may explain the presence of Jupiter-mass objects in small orbits around nearby stars.

We use the grid of models described in paper~I to analyse those millisecond pulsar binaries whose secondaries have been studied optically. In particular, we find cooling ages for these binary systems that range from $< 1 \rm Gyr$ to $\sim \rm 15 Gyr$. Comparison of cooling ages and characteristic spin down ages allows us to constrain the initial spin periods and spin-up histories for individual systems, showing that at least some millisecond pulsars had sub-Eddington accretion rates and long magnetic field decay times.

The presence of low mass, degenerate secondaries in millisecond pulsar binaries offers the opportunity to determine an age for the binary system independent of the rotational properties of the pulsar. To this end, we present here a detailed calculation of the evolution of a grid of low mass ($<0.5 M_{\odot}$) helium core white dwarfs. We investigate the effects of different Hydrogen layer masses and provide results for well-known optical band-passes. We supplement the OPAL opacity calculations with our own calculations for low effective temperatures ($ T_{eff} < 6000 K$) and also provide fitting formulae for the gravity as a function of mass and effective temperature. In paper II we shall apply these results to individual cases.

We analyse the sample of pulsar proper motions, taking detailed account of the selection effects of the original surveys. We treat censored data using survival statistics. From a comparison of our results with Monte Carlo simulations, we find that the mean birth speed of a pulsar is 250-300 km/s, rather than the 450 km/s foundby Lyne & Lorimer (1994). The resultant distribution is consistent with a maxwellian with dispersion $ \sigma_v = 190 km/s$. Despite the large birth velocities, we find that the pulsars with long characteristic ages show the asymmetric drift, indicating that they are dynamically old. These pulsars may result from the low velocity tail of the younger population, although modified by their origin in binaries and by evolution in the galactic potential.