The ability to make independent detections of the signatures of exoplanets with complementary telescopes and instruments brings a new potential for robust identification of exoplanets and precision characterization. We introduce PEXO, a package for Precise EXOplanetology to facilitate the efficient modeling of timing, astrometry, and radial velocity data, which will benefit not only exoplanet science but also various astrophysical studies in general. PEXO is general enough to account for binary motion and stellar reflex motions induced by planetary companions and is precise enough to treat various relativistic effects both in the solar system and in the target system. We also model the post-Newtonian barycentric motion for future tests of general relativity in extrasolar systems. We benchmark PEXO with the pulsar timing package TEMPO2 and find that PEXO produces numerically similar results with timing precision of about 1 ns, space-based astrometry to a precision of 1 $\mu$as, and radial velocity of 1 $\mu$m/s and improves on TEMPO2 for decade-long timing data of nearby targets, due to its consideration of third-order terms of Roemer delay. PEXO is able to avoid the bias introduced by decoupling the target system and the solar system and to account for the atmospheric effects which set a practical limit for ground-based radial velocities close to 1 cm/s. Considering the various caveats in barycentric correction and ancillary data required to realize cm/s modeling, we recommend the preservation of original observational data. The PEXO modeling package is available at GitHub (https://github.com/phillippro/pexo).

Observational evidence suggests that gas disk instability may be responsible for the formation of at least some gas giant exoplanets, particularly massive or distant gas giants. With regard to close-in gas giants, Boss (2017) used the $\beta$ cooling approximation to calculate hydrodynamical models of inner gas disk instability, finding that provided disks with low values of the initial minimum Toomre stability parameter (i.e., $Q_i < 2$ inside 20 au) form, fragmentation into self-gravitating clumps could occur even for $\beta$ as high as 100 (i.e., extremely slow cooling). Those results implied that the evolution of disks toward low $Q_i$ must be taken into account. This paper presents such models: initial disk masses of 0.091 $M_\odot$ extending from 4 to 20 au around a 1 $M_\odot$ protostar, with a range (1 to 100) of $\beta$ cooling parameters, the same as in Boss (2017), but with all the disks starting with $Q_i = 2.7$, i.e., gravitationally stable, and allowed to cool from their initial outer disk temperature of 180 K to as low as 40 K. All the disks eventually fragment into at least one dense clump. The clumps were again replaced by virtual protoplanets (VPs) and the masses and orbits of the resulting ensemble of VPs compare favorably with those of Boss (2017), supporting the claim that disk instability can form gas giants rapidly inside 20 au, provided that sufficiently massive protoplanetary disks exist.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

Microlensing is the only method that can detect and measure mass of wide orbit, low mass, solar system analog exoplanets. Mass measurements of such planets would yield massive science on planet formation, exoplanet demographics, free floating planets, planet frequencies towards the galaxy. High res follow-up observations of past microlens targets provide a mass measurement of microlens planets and hosts at an uncertainty of <20%. This will be primary method for mass measurement with WFIRST. We advocate for the fact that high resolution observations with AO, HST and JWST(in future) remain necessary in coming decade to develop the methods, to determine the field and filter selection, understand the systematics and to develop a robust pipeline to release high quality data products from WFIRST microlensing survey such that the astronomy community can promptly engage in the science. We also support future high res obs with US ELTs with advanced Laser AO systems in context of enhancing the science return of WFIRST microlensing survey. We endorse the 2018 Exoplanet Science Strategy report published by the National Academy. This white paper extends and complements the material presented therein. In particular, this white paper supports the recommendation of the National Academy Exoplanet Science Strategy report that: NASA should launch WFIRST to conduct its microlensing survey of distant planets and to demonstrate the technique of coronagraphic spectroscopy on exoplanet targets. This white paper also supports to the finding from that report which states "A number of activities, including precursor and concurrent observations using ground- and space-based facilities, would optimize the scientific yield of the WFIRST microlensing survey."

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

Gas giants are the most readily detectable exoplanets but fundamental questions about their system architectures, formation, migration, and atmospheres have been unanswerable with the current generation of ground- and space-based facilities. The dominant techniques to detect and characterize giant planets $-$ radial velocities, transits, direct imaging, microlensing, and astrometry $-$ are each isolated to a limited range of planet masses, separations, ages, and temperatures. These windows into the arrangement and physical properties of giant planets have spawned new questions about the timescale and location of their assembly; the distributions of planet mass and orbital separation at young and old ages; the composition and structure of their atmospheres; and their orbital and rotational angular momentum architectures. The ELTs will address these questions by building bridges between these islands of mass, orbital distance, and age. The angular resolution, collecting area, all-sky coverage, and novel instrumentation suite of these facilities are needed to provide a complete map of the orbits and atmospheric evolution of gas giant planets (0.3$-$10 $M_\mathrm{Jup}$) across space (0.1$-$100 AU) and time (1 Myr to 10 Gyr). This white paper highlights the scientific potential of the GMT and TMT to address these outstanding questions, with a particular focus on the role of direct imaging and spectroscopy of large samples of giant planets that will soon be made available with $Gaia$.

The transfer of circumstellar disk mass and momentum onto the protostar and out into the environment occurs via a variety of mechanisms including magnetospheric accretion, jets, outflows, and disk winds. The interplay of these processes determine both the conditions under which planet formation occurs and the lifetime of the disk. Metallic emission lines, along with the Balmer series of hydrogen, probe the kinematics of gas within the planet-forming and central regions of circumstellar disks. High-spectral resolution study of these emission lines provides critical information on mass and momentum loss, turbulence, and disk wind origins.

We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Cosmochemical evaluations of the initial meteoritical abundance of the short-lived radioisotope (SLRI) $^{26}$Al have remained fairly constant since 1976, while estimates for the initial abundance of the SLRI $^{60}$Fe have varied widely recently. At the high end of this range, $^{60}$Fe initial abundances have seemed to require $^{60}$Fe nucleosynthesis in a core collapse supernova, followed by incorporation into primitive meteoritical components within $\sim$ 1 Myr. This paper continues the detailed exploration of this classical scenario, using models of the self-gravitational collapse of molecular cloud cores that have been struck by suitable shock fronts, leading to the injection of shock front gas into the collapsing cloud through Rayleigh-Taylor fingers formed at the shock-cloud interface. As before, these models are calculated using the FLASH three dimensional, adaptive mesh refinement (AMR), gravitational hydrodynamical code. While the previous models used FLASH 2.5, the new models employ FLASH 4.3, which allows sink particles to be introduced to represent the newly formed protostellar object. Sink particles permit the models to be pushed forward farther in time to the phase where a $\sim 1 M_\odot$ protostar has formed, orbited by a rotating protoplanetary disk. These models are thus able to define what type of target cloud core is necessary for the supernova triggering scenario to produce a plausible scheme for the injection of SLRIs into the presolar cloud core: a $\sim 3 M_\odot$ cloud core rotating at a rate of $\sim 3 \times 10^{-14}$ rad s$^{-1}$ or higher.

We report individual dynamical masses for the brown dwarfs Epsilon Indi B and C, which have spectral types of T1.5 and T6, respectively, measured from astrometric orbit mapping. Our measurements are based on a joint analysis of astrometric data from the Carnegie Astrometric Planet Search and the Cerro Tololo Inter-American Observatory Parallax Investigation as well as archival high resolution imaging, and use a Markov Chain Monte Carlo method. We find dynamical masses of 75.0+-0.82 Mjup for the T1.5 B component and 70.1+-0.68 Mjup for the T6 C component. These masses are surprisingly high for substellar objects and challenge our understanding of substellar structure and evolution. We discuss several evolutionary scenarios proposed in the literature and find that while none of them can provide conclusive explanations for the high substellar masses, evolutionary models incorporating lower atmospheric opacities come closer to approximating our results. We discuss the details of our astrometric model, its algorithm implementation, and how we determine parameter values via Markov Chain Monte Carlo Bayesian inference.

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

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

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

Recent meteoritical analyses support an initial abundance of the short-lived radioisotope $^{60}$Fe that may be high enough to require nucleosynthesis in a core collapse supernova, followed by rapid incorporation into primitive meteoritical components, rather than a scenario where such isotopes were inherited from a well-mixed region of a giant molecular cloud polluted by a variety of supernovae remnants and massive star winds. This paper continues to explore the former scenario, by calculating three dimensional, adaptive mesh refinement, hydrodynamical code (FLASH 2.5) models of the self-gravitational, dynamical collapse of a molecular cloud core that has been struck by a thin shock front with a speed of 40 km/sec, leading to the injection of shock front matter into the collapsing cloud through the formation of Rayleigh-Taylor fingers at the shock-cloud intersection. These models extend the previous work into the nonisothermal collapse regime using a polytropic approximation to represent compressional heating in the optically thick protostar. The models show that the injection efficiencies of shock front material are enhanced compared to previous models, which were not carried into the nonisothermal regime and so did not reach such high densities. The new models, combined with the recent estimates of initial $^{60}$Fe abundances, imply that the supernova triggering and injection scenario remains as a plausible explanation for the origin of the short-lived radioisotopes involved in the formation of our solar system.

Observational evidence exists for the formation of gas giant planets on wide orbits around young stars by disk gravitational instability, but the roles of disk instability and core accretion for forming gas giants on shorter period orbits are less clear. The controversy extends to population synthesis models of exoplanet demographics and to hydrodynamical models of the fragmentation process. The latter refers largely to the handling of radiative transfer in three dimensional (3D) hydrodynamical models, which controls heating and cooling processes in gravitationally unstable disks, and hence dense clump formation. A suite of models using the $\beta$ cooling approximation is presented here. The initial disks have masses of 0.091 $M_\odot$ and extend from 4 to 20 AU around a 1 $M_\odot$ protostar. The initial minimum Toomre $Q_i$ values range from 1.3 to 2.7, while $\beta$ ranges from 1 to 100. We show that the choice of $Q_i$ is equal in importance to the $\beta$ value assumed: high $Q_i$ disks can be stable for small $\beta$, when the initial disk temperature is taken as a lower bound, while low $Q_i$ disks can fragment for high $\beta$. These results imply that the evolution of disks toward low $Q_i$ must be taken into account in assessing disk fragmentation possibilities, at least in the inner disk, i.e., inside about 20 AU. The models suggest that if low $Q_i$ disks can form, there should be an as yet largely undetected population of gas giants orbiting G dwarfs between about 6 AU and 16 AU.

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

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

We present high-contrast Magellan adaptive optics (MagAO) images of HD 7449, a Sun-like star with one planet and a long-term radial velocity (RV) trend. We unambiguously detect the source of the long-term trend from 0.6-2.15 \microns ~at a separation of \about 0\fasec 54. We use the object's colors and spectral energy distribution to show that it is most likely an M4-M5 dwarf (mass \about 0.1-0.2 \msun) at the same distance as the primary and is therefore likely bound. We also present new RVs measured with the Magellan/MIKE and PFS spectrometers and compile these with archival data from CORALIE and HARPS. We use a new Markov chain Monte Carlo procedure to constrain both the mass ($> 0.17$ \msun ~at 99$\%$ confidence) and semimajor axis (\about 18 AU) of the M dwarf companion (HD 7449B). We also refine the parameters of the known massive planet (HD 7449Ab), finding that its minimum mass is $1.09^{+0.52}_{-0.19}$ \mj, its semimajor axis is $2.33^{+0.01}_{-0.02}$ AU, and its eccentricity is $0.8^{+0.08}_{-0.06}$. We use N-body simulations to constrain the eccentricity of HD 7449B to $\lesssim$ 0.5. The M dwarf may be inducing Kozai oscillations on the planet, explaining its high eccentricity. If this is the case and its orbit was initially circular, the mass of the planet would need to be $\lesssim$ 1.5 \mj. This demonstrates that strong constraints on known planets can be made using direct observations of otherwise undetectable long-period companions.

This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed).

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

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

One of the primary goals of exoplanet science is to find and characterize habitable planets, and direct imaging will play a key role in this effort. Though imaging a true Earth analog is likely out of reach from the ground, the coming generation of giant telescopes will find and characterize many planets in and near the habitable zones (HZs) of nearby stars. Radial velocity and transit searches indicate that such planets are common, but imaging them will require achieving extreme contrasts at very small angular separations, posing many challenges for adaptive optics (AO) system design. Giant planets in the HZ may even be within reach with the latest generation of high-contrast imagers for a handful of very nearby stars. Here we will review the definition of the HZ, and the characteristics of detectable planets there. We then review some of the ways that direct imaging in the HZ will be different from the typical exoplanet imaging survey today. Finally, we present preliminary results from our observations of the HZ of {\alpha} Centauri A with the Magellan AO system's VisAO and Clio2 cameras.

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

Observations from optical to centimeter wavelengths have demonstrated that multiple systems of two or more bodies is the norm at all stellar evolutionary stages. Multiple systems are widely agreed to result from the collapse and fragmentation of cloud cores, despite the inhibiting influence of magnetic fields. Surveys of Class 0 protostars with mm interferometers have revealed a very high multiplicity frequency of about 2/3, even though there are observational difficulties in resolving close protobinaries, thus supporting the possibility that all stars could be born in multiple systems. Near-infrared adaptive optics observations of Class I protostars show a lower binary frequency relative to the Class 0 phase, a declining trend that continues through the Class II/III stages to the field population. This loss of companions is a natural consequence of dynamical interplay in small multiple systems, leading to ejection of members. We discuss observational consequences of this dynamical evolution, and its influence on circumstellar disks, and we review the evolution of circumbinary disks and their role in defining binary mass ratios. Special attention is paid to eclipsing PMS binaries, which allow for observational tests of evolutionary models of early stellar evolution. Many stars are born in clusters and small groups, and we discuss how interactions in dense stellar environments can significantly alter the distribution of binary separations through dissolution of wider binaries. The binaries and multiples we find in the field are the survivors of these internal and external destructive processes, and we provide a detailed overview of the multiplicity statistics of the field, which form a boundary condition for all models of binary evolution. Finally we discuss various formation mechanisms for massive binaries, and the properties of massive trapezia.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

We report results of a direct imaging survey for giant planets around 80 members of the Beta Pic, TW Hya, Tucana-Horologium, AB Dor, and Hercules-Lyra moving groups, observed as part of the Gemini NICI Planet-Finding Campaign. For this sample, we obtained median contrasts of \Delta H=13.9 mag at 1" in combined CH4 narrowband ADI+SDI mode and median contrasts of \Delta H=15.1 mag at 2" in H-band ADI mode. We found numerous (>70) candidate companions in our survey images. Some of these candidates were rejected as common-proper motion companions using archival data; we reobserved with NICI all other candidates that lay within 400 AU of the star and were not in dense stellar fields. The vast majority of candidate companions were confirmed as background objects from archival observations and/or dedicated NICI campaign followup. Four co-moving companions of brown dwarf or stellar mass were discovered in this moving group sample: PZ Tel B (36+-6 MJup, 16.4+-1.0 AU, Biller et al. 2010), CD -35 2722B (31+-8 MJup, 67+-4 AU, Wahhaj et al. 2011), HD 12894B (0.46+-0.08 MSun, 15.7+-1.0 AU), and BD+07 1919C (0.20+-0.03 MSun, 12.5+-1.4 AU). From a Bayesian analysis of the achieved H band ADI and ASDI contrasts, using power-law models of planet distributions and hot-start evolutionary models, we restrict the frequency of 1--20 MJup companions at semi-major axes from 10--150 AU to <18% at a 95.4% confidence level using DUSTY models and to <6% at a 95.4% using COND models.

We announce the discovery of a ~2 Jupiter-mass planet in an eccentric 11-year orbit around the K7/M0 dwarf GJ 328. Our result is based on 10 years' worth of radial velocity (RV) data from the Hobby-Eberly and Harlan J. Smith telescopes at McDonald Observatory, and from the Keck Telescope at Mauna Kea. Our analysis of GJ 328's magnetic activity via the Na I D features reveals a long-period stellar activity cycle, which creates an additional signal in the star's RV curve with amplitude 6-10 m/s. After correcting for this stellar RV contribution, we see that the orbit of the planet is more eccentric than suggested by the raw RV data. GJ 328b is currently the most massive, longest-period planet discovered around a low-mass dwarf.

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.5" and 14.1 mag at 1" separation. Follow-up observations of the 66 candidates with projected separation < 500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known Beta Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a >5MJup planet beyond 80 AU, and <21% of debris disk stars have a >3MJup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly-imaged planets as d^2N/dMda ~ m^alpha a^beta, where m is planet mass and a is orbital semi-major axis (with a maximum value of amax). We find that beta < -0.8 and/or alpha > 1.7. Likewise, we find that beta < -0.8 and/or amax < 200 AU. If we ignore the Beta Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that < 20% of debris disk stars have a > 3MJup planet beyond 10 AU, and beta < -0.8 and/or alpha < -1.5. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation > 40 AU and planet masses > 3 MJup do not carve the central holes in these disks.

Analyses of primitive meteorites and cometary samples have shown that the solar nebula must have experienced a phase of large-scale outward transport of small refractory grains as well as homogenization of initially spatially heterogeneous short-lived isotopes. The stable oxygen isotopes, however, were able to remain spatially heterogenous at the $\sim$ 6% level. One promising mechanism for achieving these disparate goals is the mixing and transport associated with a marginally gravitationally unstable (MGU) disk, a likely cause of FU Orionis events in young low-mass stars. Several new sets of MGU models are presented that explore mixing and transport in disks with varied masses (0.016 to 0.13 $M_\odot$) around stars with varied masses (0.1 to 1 $M_\odot$) and varied initial $Q$ stability minima (1.8 to 3.1). The results show that MGU disks are able to rapidly (within $\sim 10^4$ yr) achieve large-scale transport and homogenization of initially spatially heterogeneous distributions of disk grains or gas. In addition, the models show that while single-shot injection heterogeneity is reduced to a relatively low level ($\sim$ 1%), as required for early solar system chronometry, continuous injection of the sort associated with the generation of stable oxygen isotope fractionations by UV photolysis leads to a sustained, relatively high level ($\sim$ 10%) of heterogeneity, in agreement with the oxygen isotope data. These models support the suggestion that the protosun may have experienced at least one FU Orionis-like outburst, which produced several of the signatures left behind in primitive chondrites and comets.

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (~1.5-2.5 M_sun) conducted to date and includes the planet hosts beta Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58 (+21, -20) M_Jup and 55 (+20, -19) M_Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M_sun stars can have giant planets greater than 4 M_Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M_Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M_Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

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

We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.

Planets embedded within dust disks may drive the formation of large scale clumpy dust structures by trapping dust into resonant orbits. Detection and subsequent modeling of the dust structures would help constrain the mass and orbit of the planet and the disk architecture, give clues to the history of the planetary system, and provide a statistical estimate of disk asymmetry for future exoEarth-imaging missions. Here we present the first search for these resonant structures in the inner regions of planetary systems by analyzing the light curves of hot Jupiter planetary candidates identified by the Kepler mission. We detect only one candidate disk structure associated with KOI 838.01 at the 3-sigma confidence level, but subsequent radial velocity measurements reveal that KOI 838.01 is a grazing eclipsing binary and the candidate disk structure is a false positive. Using our null result, we place an upper limit on the frequency of dense exozodi structures created by hot Jupiters. We find that at the 90% confidence level, less than 21% of Kepler hot Jupiters create resonant dust clumps that lead and trail the planet by ~90 degrees with optical depths >~5*10^-6, which corresponds to the resonant structure expected for a lone hot Jupiter perturbing a dynamically cold dust disk 50 times as dense as the zodiacal cloud.

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

Core accretion and disk instability require giant protoplanets to form in the presence of disk gas. Protoplanet migration models generally assume disk masses low enough that the disk's self-gravity can be neglected. However, disk instability requires a disk massive enough to be marginally gravitationally unstable (MGU). Even for core accretion, a FU Orionis outburst may require a brief MGU disk phase. We present a new set of three dimensional, gravitational radiation hydrodynamics models of MGU disks with multiple protoplanets, which interact gravitationally with the disk and with each other, including disk gas mass accretion. Initial protoplanet masses are 0.01 to 10 $M_\oplus$ for core accretion models, and 0.1 to 3 $M_{Jup}$ for Nice scenario models, starting on circular orbits with radii of 6, 8, 10, or 12 AU, inside a 0.091 $M_\odot$ disk extending from 4 to 20 AU around a $1 M_\odot$ protostar. Evolutions are followed for up to $\sim$ 4000 yr and involve phases of relative stability ($e \sim$ 0.1) interspersed with chaotic phases ($e \sim$ 0.4) of orbital interchanges. The 0.01 to 10 $M_\oplus$ cores can orbit stably for $\sim$ 1000 yr: monotonic inward or outward orbital migration of the type seen in low mass disks does not occur. A system with giant planet masses similar to our Solar System (1.0, 0.33, 0.1, 0.1 $M_{Jup}$) was stable for over 1000 yr, and a Jupiter-Saturn-like system was stable for over 3800 yr, implying that our giant planets might well survive a MGU disk phase.

GJ 1214 is orbited by a transiting super-Earth-mass planet. It is a primary target for ongoing efforts to understand the emerging population of super-Earth-mass planets around M dwarfs. We present new precision astrometric measurements, a re-analysis of HARPS radial velocity measurements, and medium-resolution infrared spectroscopy of GJ 1214. We combine these measurements with recent transit follow-up observations and new catalog photometry to provide a comprehensive update of the star-planet properties. The distance is obtained with 0.6% relative uncertainty using CAPScam astrometry. The new value increases the nominal distance to the star by ~10% and is significantly more precise than previous measurements. Updated Doppler measurements combined with published transit observations significantly refine the constraints on the orbital solution. The analysis of the infrared spectrum and photometry confirm that the star is enriched in metals compared to the Sun. Using all this information, combined with empirical mass-luminosity relations for low mass stars, we derive updated values for the bulk properties of the star-planet system. We also use infrared absolute fluxes to estimate the stellar radius and to re-derive the star-planet properties. Both approaches provide very consistent values for the system. Our analysis shows indicates that the favoured mean density of GJ 1214b is 1.6 +/-0.6 g cm^{-3}. We illustrate how fundamental properties of M dwarfs are of paramount importance in the proper characterization of the low mass planetary candidates orbiting them. Given that the distance is now known to better than 1%, interferometric measurements of the stellar radius, additional high precision Doppler observations, and/or or detection of the secondary transit (occultation), are necessary to further improve the constraints on the GJ 1214 star-planet properties.

From common proper motion and signatures of youth, researchers have identified about 30 members of a putative TW Hydrae Association. Only four of these had parallactic distances from Hipparcos. We have measured parallaxes and proper motions for 14 primary members. We combine these with literature values of radial velocities to show that the Galactic space motions of the stars, with the exception of TWA 9 and 22, are parallel and do not indicate convergence at a common formation point sometime in the last few million years. The space motions of TWA 9 and 22 do not agree with the others and indicate that they are not TWA members. The median parallax is 18 mas or 56 pc. We further analyze the stars' absolute magnitudes on pre-main sequence evolutionary tracks and find a range of ages with a median of 10.1 Myr and no correlation between age and Galactic location. The TWA stars may have formed from an extended and filamentary molecular cloud but are not necessarily precisely coeval.

We report the detection of Kepler-47, a system consisting of two planets orbiting around an eclipsing pair of stars. The inner and outer planets have radii 3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a Sun-like star and a companion roughly one-third its size, orbiting each other every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the inner planet have been observed, allowing a detailed characterization of its orbit and those of the stars. The outer planet's orbital period is 303.2 days, and although the planet is not Earth-like, it resides within the classical "habitable zone", where liquid water could exist on an Earth-like planet. With its two known planets, Kepler-47 establishes that close binary stars can host complete planetary systems.

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

We have conducted a kinematic study of 165 young M dwarfs with ages of <300 Myr. Our sample is composed of stars and brown dwarfs with spectral types ranging from K7 to L0, detected by ROSAT and with photometric distances of <25 pc assuming the stars are single and on the main-sequence. In order to find stars kinematically linked to known young moving groups (YMGs), we measured radial velocities for the complete sample with Keck and CFHT optical spectroscopy and trigonometric parallaxes for 75 of the M dwarfs with the CAPSCam instrument on the du Pont 2.5-m Telescope. Due to their youthful overluminosity and unresolved binarity, the original photometric distances for our sample underestimated the distances by 70% on average, excluding two extremely young (<3 Myr) objects found to have distances beyond a few hundred parsecs. We searched for kinematic matches to 14 reported YMGs and identified 9 new members of the AB Dor YMG and 2 of the Ursa Majoris group. Additional possible candidates include 6 Castor, 4 Ursa Majoris, 2 AB Dor members, and 1 member each of the Her-Lyr and beta Pic groups. Our sample also contains 27 young low-mass stars and 4 brown dwarfs with ages <150 Myr which are not associated with any known YMG. We identified an additional 15 stars which are kinematic matches to one of the YMGs, but the ages from spectroscopic diagnostics and/or the positions on the sky do not match. These warn against grouping stars together based only on kinematics and that a confluence of evidence is required to claim that a group of stars originated from the same star-forming event.

The solar nebula is thought to have undergone a number of episodes of FU Orionis outbursts during its early evolution. We present here the first calculations of the trajectories of particles in a marginally gravitationally unstable solar nebula during an FU Orionis outburst, which show that 0.1 to 10 cm-sized particles traverse radial distances of 10 AU or more, inward and outward, in less than 200 yrs, exposing the particles to temperatures from $\sim$ 60 K to $\sim$ 1500 K. Such trajectories can thus account for the discovery of refractory particles in comets. Refractory particles should acquire Wark-Lovering-like rims as they leave the highest temperature regions of the disk, and these rims should have significant variations in their stable oxygen isotope ratios. Particles are likely to be heavily modified or destroyed if they pass within 1 AU of the Sun, and so are only likely to survive if they formed in the final few FU Orionis outbursts, or were transported to the outer reaches of the solar system. Calcium, aluminum-rich inclusions (CAIs) from primitive meteorites are the oldest known solar system objects and have a very narrow age range. Most CAIs may have formed at the end of the FU Orionis outbursts phase, with an age range reflecting the period between the last few outbursts.

We present the results of a search for planetary companions orbiting near hot Jupiter planet candidates (Jupiter-size candidates with orbital periods near 3 days) identified in the Kepler data through its sixth quarter of science operations. Special emphasis is given to companions between the 2:1 interior and exterior mean-motion resonances. A photometric transit search excludes companions with sizes ranging from roughly 2/3 to 5 times the size of the Earth, depending upon the noise properties of the target star. A search for dynamically induced deviations from a constant period (transit timing variations or TTVs) also shows no significant signals. In contrast, comparison studies of warm Jupiters (with slightly larger orbits) and hot Neptune-size candidates do exhibit signatures of additional companions with these same tests. These differences between hot Jupiters and other planetary systems denote a distinctly different formation or dynamical history.

We report the discovery by the Gemini NICI Planet-Finding Campaign of two low-mass companions to the young A0V star HD 1160 at projected separations of 81 +/- 5 AU (HD 1160 B) and 533 +/- 25 AU (HD 1160 C). VLT images of the system taken over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm that both companions are physically associated. By comparing the system to members of young moving groups and open clusters with well-established ages, we estimate an age of 50 (+50,-40) Myr for HD 1160 ABC. While the UVW motion of the system does not match any known moving group, the small magnitude of the space velocity is consistent with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 +/- 0.5 star with an estimated mass of 0.22 (+0.03,-0.04) M_Sun, while NIR photometry of HD 1160 B suggests a brown dwarf with a mass of 33 (+12,-9) M_Jup. The very small mass ratio (0.014) between the A and B components of the system is rare for A star binaries, and would represent a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun.

New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.

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

Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star's rotation period is 35.1 +/- 1.0 days, and its projected obliquity with respect to the stellar binary orbit is 1.6 +/- 2.4 degrees. Therefore the three largest sources of angular momentum---the stellar orbit, the planetary orbit, and the primary's rotation---are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the "pseudosynchronous" period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.

Short-lived radioisotopes (SLRI) such as $^{60}$Fe and $^{26}$Al were likely injected into the solar nebula in a spatially and temporally heterogeneous manner. Marginally gravitationally unstable (MGU) disks, of the type required to form gas giant planets, are capable of rapid homogenization of isotopic heterogeneity as well as of rapid radial transport of dust grains and gases throughout a protoplanetary disk. Two different types of new models of a MGU disk in orbit around a solar-mass protostar are presented. The first set has variations in the number of terms in the spherical harmonic solution for the gravitational potential, effectively studying the effect of varying the spatial resolution of the gravitational torques responsible for MGU disk evolution. The second set explores the effects of varying the initial minimum value of the Toomre $Q$ stability parameter, from values of 1.4 to 2.5, i.e., toward increasingly less unstable disks. The new models show that the basic results are largely independent of both sets of variations. MGU disk models robustly result in rapid mixing of initially highly heterogeneous distributions of SLRIs to levels of $\sim$ 10% in both the inner ($<$ 5 AU) and outer ($>$ 10 AU) disk regions, and to even lower levels ($\sim$ 2%) in intermediate regions, where gravitational torques are most effective at mixing. These gradients should have cosmochemical implications for the distribution of SLRIs and stable oxygen isotopes contained in planetesimals (e.g., comets) formed in the giant planet region ($\sim$ 5 to $\sim$ 10 AU) compared to those formed elsewhere.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.

Doppler surveys have shown that more massive stars have significantly higher frequencies of giant planets inside $\sim$ 3 AU than lower mass stars, consistent with giant planet formation by core accretion. Direct imaging searches have begun to discover significant numbers of giant planet candidates around stars with masses of $\sim$ 1 $M_\odot$ to $\sim$ 2 $M_\odot$ at orbital distances of $\sim$ 20 AU to $\sim$ 120 AU. Given the inability of core accretion to form giant planets at such large distances, gravitational instabilities of the gas disk leading to clump formation have been suggested as the more likely formation mechanism. Here we present five new models of the evolution of disks with inner radii of 20 AU and outer radii of 60 AU, for central protostars with masses of 0.1, 0.5, 1.0, 1.5, and 2.0 $M_\odot$, in order to assess the likelihood of planet formation on wide orbits around stars with varied masses. The disk masses range from 0.028 $M_\odot$ to 0.21 $M_\odot$, with initial Toomre $Q$ stability values ranging from 1.1 in the inner disks to $\sim 1.6$ in the outer disks. These five models show that disk instability is capable of forming clumps on time scales of $\sim 10^3$ yr that, if they survive for longer times, could form giant planets initially on orbits with semimajor axes of $\sim$ 30 AU to $\sim$ 70 AU and eccenticities of $\sim$ 0 to $\sim$ 0.35, with initial masses of $\sim 1 M_{Jup}$ to $\sim 5 M_{Jup}$, around solar-type stars, with more protoplanets forming as the mass of the protostar (and protoplanetary disk) are increased. In particular, disk instability appears to be a likely formation mechanism for the HR 8799 gas giant planetary system.

We present the discovery of a wide (67 AU) substellar companion to the nearby (21 pc) young solar-metallicity M1 dwarf CD-35 2722, a member of the ~100 Myr AB Doradus association. Two epochs of astrometry from the NICI Planet-Finding Campaign confirm that CD-35 2722 B is physically associated with the primary star. Near-IR spectra indicate a spectral type of L4\pm1 with a moderately low surface gravity, making it one of the coolest young companions found to date. The absorption lines and near-IR continuum shape of CD-35 2722 B agree especially well the dusty field L4.5 dwarf 2MASS J22244381-0158521, while the near-IR colors and absolute magnitudes match those of the 5 Myr old L4 planetary-mass companion, 1RXS J160929.1-210524 b. Overall, CD-35 2722 B appears to be an intermediate-age benchmark for L-dwarfs, with a less peaked H-band continuum than the youngest objects and near-IR absorption lines comparable to field objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find T=1700-1900 K and log(g) =4.5\pm0.5. The spectra also show that the radial velocities of components A and B agree to within \pm10 km/s, further confirming their physical association. Using the age and bolometric luminosity of CD-35 2722 B, we derive a mass of 31\pm8 Mjup from the Lyon/Dusty evolutionary models. Altogether, young late-M to mid-L type companions appear to be over-luminous for their near-IR spectral type compared to field objects, in contrast to the under-luminosity of young late-L and early-T dwarfs.

Meteoritic data, especially regarding chondrules and calcium-rich, aluminum-rich inclusions (CAIs), and isotopic evidence for short-lived radionuclides (SLRs) in the solar nebula, potentially can constrain how planetary systems form. Intepretation of these data demands an astrophysical model, and the "X-wind" model of Shu et al. (1996) and collaborators has been advanced to explain the origin of chondrules, CAIs and SLRs. It posits that chondrules and CAIs were thermally processed < 0.1 AU from the protostar, then flung by a magnetocentrifugal outflow to the 2-3 AU region to be incorporated into chondrites. Here we critically examine key assumptions and predictions of the X-wind model. We find a number of internal inconsistencies: theory and observation show no solid material exists at 0.1 AU; particles at 0.1 AU cannot escape being accreted into the star; particles at 0.1 AU will collide at speeds high enough to destroy them; thermal sputtering will prevent growth of particles; and launching of particles in magnetocentrifugal outflows is not modeled, and may not be possible. We also identify a number of incorrect predictions of the X-wind model: the oxygen fugacity where CAIs form is orders of magnitude too oxidizing; chondrule cooling rates are orders of magnitude lower than those experienced by barred olivine chondrules; chondrule-matrix complementarity is not predicted; and the SLRs are not produced in their observed proportions. We conclude that the X-wind model is not relevant to chondrule and CAI formation and SLR production. We discuss more plausible models for chondrule and CAI formation and SLR production.

Forming giant planets by disk instability requires a gaseous disk that is massive enough to become gravitationally unstable and able to cool fast enough for self-gravitating clumps to form and survive. Models with simplified disk cooling have shown the critical importance of the ratio of the cooling to the orbital timescales. Uncertainties about the proper value of this ratio can be sidestepped by including radiative transfer. Three-dimensional radiative hydrodynamics models of a disk with a mass of $0.043 M_\odot$ from 4 to 20 AU in orbit around a $1 M_\odot$ protostar show that disk instabilities are considerably less successful in producing self-gravitating clumps than in a disk with twice this mass. The results are sensitive to the assumed initial outer disk ($T_o$) temperatures. Models with $T_o$ = 20 K are able to form a single self-gravitating clump, whereas models with $T_o$ = 25 K form clumps that are not quite self-gravitating. These models imply that disk instability requires a disk with a mass of at least $\sim 0.043 M_\odot$ inside 20 AU in order to form giant planets around solar-mass protostars with realistic disk cooling rates and outer disk temperatures. Lower mass disks around solar-mass protostars must rely upon core accretion to form inner giant planets.

The Palomar High-precision Astrometric Search for Exoplanet Systems monitored 51 subarcsecond binary systems to evaluate whether tertiary companions as small as Jovian planets orbited either the primary or secondary stars, perturbing their otherwise smooth Keplerian motions. Six binaries are presented that show evidence of substellar companions orbiting either the primary or secondary star. Of these six systems, the likelihoods of two of the detected perturbations to represent real objects are considered to be "high confidence", while the remaining four systems are less certain and will require continued observations for confirmation.

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

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

The Kepler Mission, launched on Mar 6, 2009 was designed with the explicit capability to detect Earth-size planets in the habitable zone of solar-like stars using the transit photometry method. Results from just forty-three days of data along with ground-based follow-up observations have identified five new transiting planets with measurements of their masses, radii, and orbital periods. Many aspects of stellar astrophysics also benefit from the unique, precise, extended and nearly continuous data set for a large number and variety of stars. Early results for classical variables and eclipsing stars show great promise. To fully understand the methodology, processes and eventually the results from the mission, we present the underlying rationale that ultimately led to the flight and ground system designs used to achieve the exquisite photometric performance. As an example of the initial photometric results, we present variability measurements that can be used to distinguish dwarf stars from red giants.

Five new planets orbiting G and K dwarfs have emerged from the Magellan velocity survey. These companions are jovian-mass planets in eccentric (e \geq 0.24) intermediate and long-period orbits. HD 86226b orbits a solar metallicity G2 dwarf. The MP sin i mass of the planet is 1.5 MJUP, the semi-major axis is 2.6 AU, and the eccentricity 0.73. HD 129445b orbits a metal rich G6 dwarf. The minimum mass of the planet is MP sin i =1.6 MJUP, the semi-major axis is 2.9 AU, and the eccentricity 0.70. HD 164604b orbits a K2 dwarf. The MP sin i mass is 2.7 MJUP, semi-major axis is 1.3 AU, and the eccentricity is 0.24. HD 175167b orbits a metal rich G5 star. The MP sin i mass is 7.8 MJUP, the semi-major axis is 2.4 AU, and the eccentricity 0.54. HD 152079b orbits a G6 dwarf. The MP sin i mass of the planet is 3 MJUP, the semi-major axis is 3.2 AU, and the eccentricity is 0.60.

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

We are undertaking an astrometric search for gas giant planets and brown dwarfs orbiting nearby low mass dwarf stars with the 2.5-m du Pont telescope at the Las Campanas Observatory in Chile. We have built two specialized astrometric cameras, the Carnegie Astrometric Planet Search Cameras (CAPSCam-S and CAPSCam-N), using two Teledyne Hawaii-2RG HyViSI arrays, with the cameras' design having been optimized for high accuracy astrometry of M dwarf stars. We describe two independent CAPSCam data reduction approaches and present a detailed analysis of the observations to date of one of our target stars, NLTT 48256. Observations of NLTT 48256 taken since July 2007 with CAPSCam-S imply that astrometric accuracies of around 0.3 milliarcsec per hour are achievable, sufficient to detect a Jupiter-mass companion orbiting 1 AU from a late M dwarf 10 pc away with a signal-to-noise ratio of about 4. We plan to follow about 100 nearby (primarily within about 10 pc) low mass stars, principally late M, L, and T dwarfs, for 10 years or more, in order to detect very low mass companions with orbital periods long enough to permit the existence of habitable, Earth-like planets on shorter-period orbits. These stars are generally too faint and red to be included in ground-based Doppler planet surveys, which are often optimized for FGK dwarfs. The smaller masses of late M dwarfs also yield correspondingly larger astrometric signals for a given mass planet. Our search will help to determine whether gas giant planets form primarily by core accretion or by disk instability around late M dwarf stars.

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

The collapse and fragmentation of initially prolate and oblate, magnetic molecular clouds is calculated in three dimensions with a gravitational, radiative hydrodynamics code. The code includes magnetic field effects in an approximate manner: magnetic pressure, tension, braking, and ambipolar diffusion are all modelled. The parameters varied for both the initially prolate and oblate clouds are the initial degree of central concentration of the radial density profile, the initial angular velocity, and the efficiency of magnetic braking (represented by a factor $f_{mb} = 10^{-4}$ or $10^{-3}$). The oblate cores all collapse to form rings that might be susceptible to fragmentation into multiple systems. The outcome of the collapse of the prolate cores depends strongly on the initial density profile. Prolate cores with central densities 20 times higher than their boundary densities collapse and fragment into binary or quadruple systems, whereas cores with central densities 100 times higher collapse to form single protostars embedded in bars. The inclusion of magnetic braking is able to stifle protostellar fragmentation in the latter set of models, as when identical models were calculated without magnetic braking (Boss 2002), those cores fragmented into binary protostars. These models demonstrate the importance of including magnetic fields in studies of protostellar collapse and fragmentation, and suggest that even when magnetic fields are included, fragmentation into binary and multiple systems remains as a possible outcome of protostellar collapse.

We present an analysis of three years of precision radial velocity measurements of 160 metal-poor stars observed with HIRES on the Keck 1 telescope. We report on variability and long-term velocity trends for each star in our sample. We identify several long-term, low-amplitude radial-velocity variables worthy of follow-up with direct imaging techniques. We place lower limits on the detectable companion mass as a function of orbital period. Our survey would have detected, with a 99.5% confidence level, over 95% of all companions on low-eccentricity orbits with velocity semi-amplitude K > 100 m/s, or M_p*sin(i) > 3.0 M_JUP*(P/yr)^(1/3), for orbital periods P< 3 yr. None of the stars in our sample exhibits radial-velocity variations compatible with the presence of Jovian planets with periods shorter than the survey duration. The resulting average frequency of gas giants orbiting metal-poor dwarfs with -2.0 < [Fe/H] < -0.6 is f_p<0.67% (at the 1-sigma confidence level). We examine the implications of this null result in the context of the observed correlation between the rate of occurrence of giant planets and the metallicity of their main-sequence solar-type stellar hosts. By combining our dataset with the Fischer & Valenti (2005) uniform sample, we confirm that the likelihood of a star to harbor a planet more massive than Jupiter within 2 AU is a steeply rising function of the host's metallicity. However, the data for stars with -1.0 < [Fe/H] < 0.0 are compatible, in a statistical sense, with a constant occurrence rate f_p~1%. Our results can usefully inform theoretical studies of the process of giant planet formation across two orders of magnitude in metallicity.

The detection of lower mass planets now being reported via radial velocity and microlensing surveys suggests that they may be ubiquitous. If missions such as Kepler are able to confirm this, the detection and study of rocky planets via direct imaging with ground-based telescopes of apertures > 20 m is feasible in the thermal infrared. We discuss two cases for detecting rocky planets, the first via detection of molten Earths formed though an Earth-Moon like impact event, and the second via detection of planets around very nearby stars. These observations have the potential to give us a first look at a rocky planet similar to the Earth.

The disk instability mechanism for giant planet formation is based on the formation of clumps in a marginally-gravitationally unstable protoplanetary disk, which must lose thermal energy through a combination of convection and radiative cooling if they are to survive and contract to become giant protoplanets. While there is good observational support for forming at least some giant planets by disk instability, the mechanism has become theoretically contentious, with different three dimensional radiative hydrodynamics codes often yielding different results. Rigorous code testing is required to make further progress. Here we present two new analytical solutions for radiative transfer in spherical coordinates, suitable for testing the code employed in all of the Boss disk instability calculations. The testing shows that the Boss code radiative transfer routines do an excellent job of relaxing to and maintaining the analytical results for the radial temperature and radiative flux profiles for a spherical cloud with high or moderate optical depths, including the transition from optically thick to optically thin regions. These radial test results are independent of whether the Eddington approximation, diffusion approximation, or flux-limited diffusion approximation routines are employed. The Boss code does an equally excellent job of relaxing to and maintaining the analytical results for the vertical (theta) temperature and radiative flux profiles for a disk with a height proportional to the radial distance. These tests strongly support the disk instability mechanism for forming giant planets.

We report low mass companions orbiting five Solar-type stars that have emerged from the Magellan precision Doppler velocity survey, with minimum (Msini) masses ranging from 1.2 to 25 Mjup. These nearby target stars range from mildly metal-poor to metal-rich, and appear to have low chromospheric activity. The companions to the brightest two of these stars have previously been reported from the CORALIE survey. Four of these companions (HD 48265-b, HD 143361-b, HD 28185-b, HD 111232-b) are low-mass Jupiter-like planets in eccentric intermediate and long-period orbits. On the other hand, the companion to HD 43848 appears to be a long period brown dwarf in a very eccentric orbit.

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

Both core accretion and disk instability appear to be required as formation mechanisms in order to explain the entire range of giant planets found in extrasolar planetary systems. Disk instability is based on the formation of clumps in a marginally-gravitationally unstable protoplanetary disk. These clumps can only be expected to contract and survive to become protoplanets if they are able to lose thermal energy through a combination of convection and radiative cooling. Here we present several new three dimensional, radiative hydrodynamics models of self-gravitating protoplanetary disks, where radiative transfer is handled in the flux-limited diffusion approximation. We show that while the flux-limited models lead to higher midplane temperatures than in a diffusion approximation model without the flux-limiter, the difference in temperatures does not appear to be sufficiently high to have any significant effect on the formation of self-gravitating clumps. Self-gravitating clumps form rapidly in the models both with and without the flux-limiter. These models suggest that the reason for the different outcomes of numerical models of disk instability by different groups cannot be attributed solely to the handling of radiative transfer, but rather appears to be caused by a range of numerical effects and assumptions. Given the observational imperative to have disk instability form at least some extrasolar planets, these models imply that disk instability remains as a viable giant planet formation mechanism.

The discovery of refractory grains amongst the particles collected from Comet 81P/Wild 2 by the Stardust spacecraft (Brownlee et al. 2006) provides the ground truth for large-scale transport of materials formed in high temperature regions close to the protosun outward to the comet-forming regions of the solar nebula. While accretion disk models driven by a generic turbulent viscosity have been invoked as a means to explain such large-scale transport, the detailed physics behind such an ``alpha'' viscosity remains unclear. We present here an alternative physical mechanism for large-scale transport in the solar nebula: gravitational torques associated with the transient spiral arms in a marginally gravitationally unstable disk, of the type that appears to be necessary to form gas giant planets. Three dimensional models are presented of the time evolution of self-gravitating disks, including radiative transfer and detailed equations of state, showing that small dust grains will be transported upstream and downstream (with respect to the mean inward flow of gas and dust being accreted by the central protostar) inside the disk on time scales of less than 1000 yr inside 10 AU. These models furthermore show that any initial spatial heterogeneities present (e.g., in short-lived isotopes such as 26Al) will be homogenized by disk mixing down to a level of ~10%, preserving the use of short-lived isotopes as accurate nebular chronometers, while simultaneously allowing for the spread of stable oxygen isotope ratios. This finite level of nebular spatial heterogeneity appears to be related to the coarse mixing achieved by spiral arms, with radial widths of order 1 AU, over time scales of ~1000 yrs.

We summarize the results of two experiments to address important issues related to the correlation between planet frequencies and properties and the metallicity of the hosts. Our results can usefully inform formation, structural, and evolutionary models of gas giant planets.

We report the detection of a planetary companion with a minimum mass of m sin i = 0.0771 M_Jup = 24.5 M_Earth to the nearby (d = 9.4 pc) M2.5V star GJ 176. The star was observed as part of our M dwarf planet search at the Hobby-Eberly Telescope (HET). The detection is based on 5 years of high-precision differential radial velocity (RV) measurements using the High-Resolution-Spectrograph (HRS). The orbital period of the planet is 10.24 d. GJ 176 thus joins the small (but increasing) sample of M dwarfs hosting short-periodic planets with minimum masses in the Neptune-mass range. Low mass planets could be relatively common around M dwarfs and the current detections might represent the tip of a rocky planet population.

We review the models and results of simulations of self-gravitating, gaseous protoplanetary disks in binary star systems. These models have been calculated by three different groups with three different computational methods, two particle-based and one grid-based. We show that interactions with the companion star can affect the temperature distribution and structural evolution of disks, and discuss the implications for giant planet formation by gravitational instability as well as by core accretion. Complete consensus has not been reached yet on whether the formation of giant planets is promoted or suppressed by tidal interactions with a companion star. While systems with binary separations of order 100 AU or larger appear to behave more or less as in isolation, systems with smaller separations exhibit an increased or decreased susceptibility to fragmentation, depending on the details of thermodynamics, in particular on the inclusion or absence of artificial viscosity, and on the initial conditions. While code comparisons on identical problems need to be carried out (some of which are already in progress) to decide which computer models are more realistic, it is already clear that relatively close binary systems, with separations of order 60 AU or less, should provide strong constraints on how giant planets form in these systems.

Disk instability is an attractive yet controversial means for the rapid formation of giant planets in our solar system and elsewhere. Recent concerns regarding the first adiabatic exponent of molecular hydrogen gas are addressed and shown not to lead to spurious clump formation in the author's disk instability models. A number of disk instability models have been calculated in order to further test the robustness of the mechanism, exploring the effects of changing the pressure equation of state, the vertical temperature profile, and other parameters affecting the temperature distribution. Possible reasons for differences in results obtained by other workers are discussed. Disk instability remains as a plausible formation mechanism for giant planets.

In this paper, we calculate simulated scattered light images of a circumstellar disk in which a planet is forming by gravitational instability. The simulated images bear no correlation to the vertically integrated surface density of the disk, but rather trace the density structure in the tenuous upper layers of the disk. Although the density at high altitudes does not bear a direct relation to activity at the midplane, the very existence of structure at high altitudes along with high time variability is an indicator of gravitational instability within the disk. The timescale for variations is much shorter than the orbital period of the planet, which facilitates observation of the phenomenon. Scattered light images may not necessarily be able to tell us where exactly a planet might be forming in a disk, but can still be a useful probe of active planet formation within a circumstellar disk. Although these phenomena are unlikely to be observable by current telescopes, future large telescopes, such as the Giant Magellan Telescope, may be able to detect them.

Isotopic abundances of short-lived radionuclides such as 26Al provide the most precise chronometers of events in the early solar system, provided that they were initially homogeneously distributed. On the other hand, the abundances of the three stable isotopes of oxygen in primitive meteorites show a mass-independent fractionation that survived homogenization in the solar nebula. As as result of this and other cosmochemical evidence, the degree of spatial heterogeneity of isotopes in the solar nebula has long been a puzzle. We show here that based on hydrodynamical models of the mixing and transport of isotopic anomalies formed at, or injected onto, the surface of the solar nebula, initially high levels of isotopic spatial heterogeneity are expected to fall to steady state levels (~10%) low enough to validate the use of 26Al for chronometry, but high enough to preserve the evidence for mass-independent fractionation of oxygen isotopes. The solution to this puzzle relies on the mixing being accomplished by the chaotic fluid motions in a marginally gravitationally unstable disk, as seems to be required for the formation of gas giant planets and by the inability of alternative physical processes to drive large-scale mixing and transport in the planet-forming midplane of the solar nebula. Such a disk is also capable of large-scale outward transport of the thermally annealed dust grains found in comets, and of driving the shock fronts that appear to be responsible for much of the thermal processing of the components of primitive meteorites, creating a self-consistent picture of the basic physical processes shaping the early solar nebula.

The collapse and fragmentation of initially filamentary, magnetic molecular clouds is calculated in three dimensions with a gravitational, radiative hydrodynamics code. The code includes magnetic field effects in an approximate manner: magnetic pressure, tension, braking, and ambipolar diffusion are all modelled. Three types of outcomes are observed: direct collapse and fragmentation into a multiple protostar system, periodic contraction and expansion without collapse, or periodic contraction and expansion leading eventually to collapse. While the models begin their evolution at rest except for the assumed solid-body rotation, they develop weakly supersonic velocity fields as a result of the rebounding prior to collapse. The models show that magnetically-supported clouds subject to magnetic braking can undergo dynamic collapse leading to protostellar fragmentation on scales of 10 AU to 100 AU, consistent with typical binary star separations.

Stable, hydrogen-burning, M dwarf stars comprise about 75% of all stars in the Galaxy. They are extremely long-lived and because they are much smaller in mass than the Sun (between 0.5 and 0.08 MSun), their temperature and stellar luminosity are low and peaked in the red. We have re-examined what is known at present about the potential for a terrestrial planet forming within, or migrating into, the classic liquid-surface-water habitable zone close to an M dwarf star. Observations of protoplanetary disks suggest that planet-building materials are common around M dwarfs, but N-body simulations differ in their estimations of the likelihood of potentially-habitable, wet planets residing within their habitable zones, which are only ~ 1/5 to 1/50 of the width of that for a G star. Particularly in light of the claimed detection of the planets with masses as small as 5.5 and 7.5 MEarth orbiting M stars, there seems no reason to exclude the possibility of terrestrial planets. Tidally locked synchronous rotation within the narrow habitable zone doesn't necessarily lead to atmospheric collapse, and active stellar flaring may not be as much of an evolutionarily disadvantageous factor as has previously been supposed. We conclude that M dwarf stars may indeed be viable hosts for planets on which the origin and evolution of life can occur. A number of planetary processes such as cessation of geothermal activity, or thermal and non-thermal atmospheric loss processes may limit the duration of planetary habitability to periods far shorter than the extreme lifetime of the M dwarf star. Nevertheless, it makes sense to include M dwarf stars in programs that seek to find habitable worlds and evidence of life.

We describe a high-precision Doppler search for giant planets orbiting a well-defined sample of metal-poor dwarfs in the field. This experiment constitutes a fundamental test of theoretical predictions which will help discriminate between proposed giant planet formation and migration models. We present here details on the survey as well as an overall assessment of the quality of our measurements, making use of the results for the stars that show no significant velocity variation.

In the era of high resolution submillimeter interferometers, it will soon be possible to observe the neutral circumstellar medium directly involved in gas giant planet (GGP) formation at physical scales previously unattainable. In order to explore possible signatures of gas giant planet formation via disk instabilities, we have combined a 3D, non-local thermodynamic equilibrium (LTE) radiative transfer code with a 3D, finite differences hydrodynamical code to model molecular emission lines from the vicinity of a 1.4 M_J self-gravitating proto-GGP. Here, we explore the properties of rotational transitions of the commonly observed dense gas tracer, HCO+. Our main results are the following: 1. Very high lying HCO+ transitions (e.g. HCO+ J=7-6) can trace dense planet forming clumps around circumstellar disks. Depending on the molecular abundance, the proto-GGP may be directly imageable by the Atacama Large Millimeter Array (ALMA). 2. HCO+ emission lines are heavily self-absorbed through the proto-GGP's dense molecular core. This signature is nearly ubiquitous, and only weakly dependent on assumed HCO+ abundances. The self-absorption features are most pronounced at higher angular resolutions. Dense clumps that are not self-gravitating only show minor self-absorption features. 3. Line temperatures are highest through the proto-GGP at all assumed abundances and inclination angles. Conversely, due to self-absorption in the line, the velocity-integrated intensity may not be. High angular resolution interferometers such as the Submillimeter Array (SMA) and ALMA may be able to differentiate between competing theories of gas giant planet formation.

While the recent microlensing discoveries of super-Earths orbiting two M dwarf stars have been taken as support for the core accretion mechanism of giant planet formation, we show here that these planets could also have been formed by the competing mechanism of disk instability, coupled with photoevaporative loss of their gaseous envelopes by a strong external source of UV radiation, i.e., an O star. M dwarfs that form in regions of future high-mass star formation would then be expected to have super-Earths orbiting at distances of several AU and beyond, while those that form in regions of low-mass star formation would be expected to have gas giants at those distances. Given that most stars are born in the former rather than in the latter regions, M dwarfs should have significantly more super-Earths than gas giants, as seems to be indicated by the microlensing surveys.

Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.

Extrasolar planet surveys have begun to detect gas giant planets in orbit around M dwarf stars. While the frequency of gas giant planets around M dwarfs so far appears to be lower than that around G dwarfs, it is clearly not zero. Previous work has shown that the core accretion mechanism does not seem to be able to form gas giant planets around M dwarfs, because the time required for core formation scales with the orbital period, which lengthens for lower mass stars, resulting in failed (gas-poor) cores unless the gaseous protoplanetary disk survives for > 10 Myr. Disk instability, on the other hand, is rapid enough (~ 1000 yrs) that it should be able to form gas giant protoplanets around even low mass stars well before the gaseous disk disappears. A new suite of three dimensional radiative, gravitational hydrodynamical models is presented that calculates the evolution of initially marginally gravitationally unstable disks with masses of 0.021 to 0.065 solar masses orbiting around stars with masses of 0.1 and 0.5 solar masses, respectively. The models show that gas giant planets are indeed likely to form by the disk instability mechanism in orbit around M dwarf stars, the opposite of the prediction for formation by the core accretion mechanism. This difference offers another observational test for discriminating between these two theoretical end members for giant planet formation. Ongoing and future extrasolar planet searches around M dwarfs by spectroscopy, microlensing, photometry, and astrometry offer the opportunity to help decide between the dominance of the two mechanisms.

A new suite of three dimensional radiative, gravitational hydrodynamical models is used to show that gas giant planets are unlikely to form by the disk instability mechanism at distances of ~100 AU to ~200 AU from young stars. A similar result seems to hold for the core accretion mechanism. These results appear to be consistent with the paucity of detections of gas giant planets on wide orbits by infrared imaging surveys, and also imply that if the object orbiting GQ Lupus is a gas giant planet, it most likely did not form at a separation of ~100 AU. Instead, a wide planet around GQ Lup must have undergone a close encounter with a third body that tossed the planet outward to its present distance from its protostar. If it exists, the third body may be detectable by NASA's Space Interferometry Mission.

We present a suite of three dimensional radiative gravitational hydrodynamics models suggesting that binary stars may be quite capable of forming planetary systems similar to our own. The new models with binary companions do not employ any explicit artificial viscosity, and also include the third (vertical) dimension in the hydrodynamic calculations, allowing for transient phases of convective cooling. The calculations of the evolution of initially marginally gravitationally stable disks show that the presence of a binary star companion may actually help to trigger the formation of dense clumps that could become giant planets. We also show that in models without binary companions, which begin their evolution as gravitationally stable disks, the disks evolve to form dense rings, which then break-up into self-gravitating clumps. These latter models suggest that the evolution of any self-gravitating disk with sufficient mass to form gas giant planets is likely to lead to a period of disk instability, even in the absence of a trigger such as a binary star companion.

Chondrules are mm-sized spherules found throughout primitive, chondritic meteorites. Flash heating by a shock front is the leading explanation of their formation. However, identifying a mechanism for creating shock fronts inside the solar nebula has been difficult. In a gaseous disk capable of forming Jupiter, the disk must have been marginally gravitationally unstable at and beyond Jupiter's orbit. We show that this instability can drive inward spiral shock fronts with shock speeds of up to about 10 km/s at asteroidal orbits, sufficient to account for chondrule formation. Mixing and transport of solids in such a disk, combined with the planet-forming tendencies of gravitational instabilities, results in a unified scenario linking chondrite production with gas giant planet formation.

We present results from our ongoing spectroscopic search for giant planets within 1 AU around a well-defined sample of metal-poor stars with HIRES on the Keck 1 telescope. We have achieved an rms radial velocity precision of $\sim 8$ m/s over a time-span of 1.5 years. The data collected so far build toward evidence of the absence of very short-period ($< 1$ month) giant planets. However, about 7% of the stars in our sample exhibits velocity trends indicative of the existence of companions. We place preliminary upper limits on the detectable companion mass as a function of orbital period, and compare them with the performance of ESA's future space-borne high-precision astrometric observatory Gaia.

We present preliminary results from our spectroscopic search for planets within 1 AU of metal-poor field dwarfs using NASA time with HIRES on Keck I. The core accretion model of gas giant planet formation is sensitive to the metallicity of the raw material, while the disk instability model is not. By observing metal-poor stars in the field we eliminate the role of dynamical interactions in dense stellar environments, such as a globular cluster. The results of our survey should allow us to distinguish the relative roles of the two competing giant planet formation scenarios.

We study the motions of small solids, ranging from micron-sized dust grains to meter-sized objects, in the vicinity of local pressure enhancements of a gaseous nebula. Integrating numerically, we show that as a result of the combined effect of gas drag and pressure gradients, solids tend to accumulate at the locations where the pressure of the gas maximizes. The rate of migration of solids varies with their sizes and densities and also with the physical properties of the gas. The results of our numerical simulations indicate that such migrations are most rapid for meter-sized objects. The applicability of the results to the enhancement of the collision and coagulation of solids and also to the growth-rate of planetesimals is discussed.

We study the motions of small solids, ranging from micron-sized dust grains to 100-m objects, in the vicinity of a local density enhancement of an isothermal gaseous solar nebula. Being interested in possible application of the results to the formation of clumps and spiral arms in a circumstellar disk, we numerically integrate the equations of motion of such solids and study their migration for different values of their sizes and masses and also for different physical properties of the gas, such as its density and temperature. We show that, considering the drag force of the gas and also the gravitational attraction of the nebula, it is possible for solids, within a certain range of size and mass, to migrate rapidly (i.e. within ~1000 years) toward the location of a local maximum density where collisions and coagulation may result in an accelerated rate of planetesimal formation.

The existence of Uranus and Neptune presents severe difficulties for the core accretion model for the formation of ice giant planets. We suggest an alternative mechanism, namely disk instability leading to the formation of gas giant protoplanets, coagulation and settling of dust grains to form ice/rock cores at their centers, and photoevaporation of their gaseous envelopes by a nearby OB star, as a possible means of forming ice giant planets.

Meteorite studies have revealed the presence of short-lived radioactivities in the early solar system. The current data suggests that the origin of at least some of the radioactivities requires contribution from recent nucleosynthesis at a stellar site. This sets a strict time limit on the time available for the formation of the solar system and argues for the theory of the triggered origin of the solar system. According to this scenario, the formation of our planetary system was initiated by the impact of an interstellar shock wave on a molecular cloud core. The shock wave originated from a nearby explosive stellar event and carried with it radioactivities produced in the stellar source. In addition to triggering the collapse of the molecular cloud core, the shock wave also deposited some of the freshly synthesized radioactivities into the collapsing system. The radioactivities were then incorporated into the first solar system solids, in this manner leaving a record of the event in the meteoritic material. The viability of the scenario can be investigated through numerical simulations studying the processes involved in mixing shock wave material into the collapsing system. The high-resolution calculations presented here show that injection occurs through Rayleigh-Taylor instabilities, the injection efficiency is approximately 10%, and temporal and spatial heterogeneities in the abundances of the radioactivities existed at the time of their arrival in the forming solar system.