Search

Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing GEMS both through confirmed planets, as well as protoplanetary disk observations, and a combination of tests to reconcile these with theoretical predictions. We then introduce the \textit{Searching for GEMS} survey, where we utilize multi-dimensional nonparameteric statistics to simulate hypothetical survey scenarios to predict the required sample size of transiting GEMS with mass measurements to robustly compare their bulk-density with canonical hot-Jupiters orbiting FGK stars. Our Monte-Carlo simulations predict that a robust comparison requires about 40 transiting GEMS (compared to the existing sample of $\sim$ 15) with 5-$\sigma$ mass measurements. Furthermore, we discuss the limitations of existing occurrence estimates for GEMS, and provide a brief description of our planned systematic search to improve the occurrence rate estimates for GEMS., Comment: 16 pages + references, including 7 figures. Accepted in AAS Journals

Recent discoveries of gas giant exoplanets around M-dwarfs (GEMS) from transiting and radial velocity (RV) surveys are difficult to explain with core-accretion models. We present here a homogeneous suite of 162 models of gravitationally unstable gaseous disks. These models represent an existence proof for gas giants more massive than 0.1 Jupiter masses to form by the gas disk gravitational instability (GDGI) mechanism around M-dwarfs for comparison with observed exoplanet demographics and protoplanetary disk mass estimates for M-dwarf stars. We use the Enzo 2.6 adaptive mesh refinement (AMR) 3D hydrodynamics code to follow the formation and initial orbital evolution of gas giant protoplanets in gravitationally unstable gaseous disks in orbit around M-dwarfs with stellar masses ranging from 0.1 $M_\odot$ to 0.5 $M_\odot$. The gas disk masses are varied over a range from disks that are too low in mass to form gas giants rapidly to those where numerous gas giants are formed, therefore revealing the critical disk mass necessary for gas giants to form by the GDGI mechanism around M-dwarfs. The disk masses vary from 0.01 $M_\odot$ to 0.05 $M_\odot$ while the disk to star mass ratios explored range from 0.04 to 0.3. The models have varied initial outer disk temperatures (10 K to 60 K) and varied levels of AMR grid spatial resolution, producing a sample of expected gas giant protoplanets for each star mass. Broadly speaking, disk masses of at least 0.02 $M_\odot$ are needed for the GDGI mechanism to form gas giant protoplanets around M-dwarfs., Comment: 34 pages, 8 figures, in press, ApJ

Goal 1 of the National Academies of Science, Engineering and Mathematics Exoplanet Science Strategy is "to understand the formation and evolution of planetary systems as products of the process of star formation, and characterize and explain the diversity of planetary system architectures, planetary compositions, and planetary environments produced by these processes", with the finding that "Current knowledge of the demographics and characteristics of planets and their systems is substantially incomplete." One significant roadblock to our ongoing efforts to improve our demographics analyses is the lack of comprehensive meta-data accompanying published exoplanet surveys. The Exoplanet Program Analysis Group (ExoPAG) Science Interest Group 2: Exoplanet Demographics has prepared this document to provide guidance to survey architects, authors, referees and funding agencies as to the most valuable such data products for five different exoplanet detection techniques - transit, radial velocity, direct imaging, microlensing and astrometry. We find that making these additional data easily available would greatly enhance the community's ability to perform robust, reproducible demographics analyses, and make progress on achieving the most important goals identified by the exoplanet and wider astronomical community., Comment: 21 pages, final report after community feedback addressed

Protoplanets formed in a marginally gravitationally unstable (MGU) disk by either core accretion or disk instability will be subject to dynamical interactions with massive spiral arms, possibly resulting in inward or outward orbital migration, mergers with each other, or even outright ejection from the protoplanetary system. The latter process has been hypothesized as a possible formation scenario for the unexpectedly high frequency of unbound gas giant exoplanets (free floating planets, FFP). Previous calculations with the EDTONS fixed grid three dimensional (3D) hydrodynamics code found that protoplanets with masses from 0.01 $M_\oplus$ to 3 $M_{Jup}$ could undergo chaotic orbital evolutions in MGU disks for $\sim$ 1000 yrs without undergoing monotonic inward or outward migration. Here the Enzo 2.5 adaptive mesh refinement (AMR) 3D hydrodynamics code is used to follow the formation and orbital evolution of protoplanets in MGU disks for up to 2000 yrs. The Enzo results confirm the basic disk fragmentation results of the EDTONS code, as well as the absence of monotonic inward or outward orbital migration. In addition, Enzo allows protoplanet mergers to occur, unlike EDTONS, resulting in a significant decrease in the number of protoplanets that survive for 1000 to 2000 yrs in the Enzo models. These models also imply that gas giants should be ejected frequently in MGU disks that fragment into large numbers of protoplanets, supporting ejection as a possible source mechanism for the observed FFPs., Comment: 27 pages, 10 figures, 2 tables, in press, Astrophysical Journal

We present the discovery of TOI-5205~b, a transiting Jovian planet orbiting a solar metallicity M4V star, which was discovered using Transiting Exoplanet Survey Satellite photometry and then confirmed using a combination of precise radial velocities, ground-based photometry, spectra, and speckle imaging. TOI-5205~b has one of the highest mass ratios for M dwarf planets with a mass ratio of almost 0.3$\%$, as it orbits a host star that is just $0.392 \pm 0.015$ \solmass{}. Its planetary radius is $1.03 \pm 0.03~R_J$, while the mass is $1.08 \pm 0.06~M_J$. Additionally, the large size of the planet orbiting a small star results in a transit depth of $\sim 7\%$, making it one of the deepest transits of a confirmed exoplanet orbiting a main-sequence star. The large transit depth makes TOI-5205~b a compelling target to probe its atmospheric properties, as a means of tracing the potential formation pathways. While there have been radial-velocity-only discoveries of giant planets around mid-M dwarfs, this is the first transiting Jupiter with a mass measurement discovered around such a low-mass host star. The high mass of TOI-5205~b stretches conventional theories of planet formation and disk scaling relations that cannot easily recreate the conditions required to form such planets., Comment: Accepted in AJ. arXiv admin note: text overlap with arXiv:2203.07178

Cook et al. (2021) found that iron meteorites have an initial abundance ratio of the short-lived isotope $^{60}$Fe to the stable isotope $^{56}$Fe of $^{60}$Fe/$^{56}$Fe $\sim$ $(6.4 \pm 2.0) \times 10^{-7}$. This appears to require the injection of live $^{60}$Fe from a Type II supernova (SN II) into the presolar molecular cloud core, as the observed ratio is over a factor of ten times higher than would be expected to be found in the ambient interstellar medium (ISM) as a result of galactic chemical evolution. The supernova triggering and injection scenario offers a ready explanation for an elevated initial $^{60}$Fe level, and in addition provides a physical mechanism for explaining the non-carbonaceous -- carbonaceous (NC-CC) dichotomy of meteorites. The NC-CC scenario hypothesizes the solar nebula first accreted material that was enriched in supernova-derived nuclides, and then later accreted material depleted in supernova-derived nuclides. While the NC-CC dichotomy refers to stable nuclides, not short-lived isotopes like $^{60}$Fe, the SN II triggering hypothesis provides an explanation for the otherwise unexplained change in nuclides being accreted by the solar nebula. Three dimensional hydrodynamical models of SN II shock-triggered collapse show that after triggering collapse of the presolar cloud core, the shock front sweeps away the local ISM while accelerating the resulting protostar/disk to a speed of several km/s, sufficient for the protostar/disk system to encounter within $\sim$ 1 Myr the more distant regions of a giant molecular cloud complex that might be expected to have a depleted inventory of supernova-derived nuclides., Comment: 26 pages, 14 figures, accepted for ApJ

While collisional accumulation is nearly universally accepted as the formation mechanism of rock and ice worlds, the situation regarding gas giant planet formation is more nuanced. Gas accretion by solid cores formed by collisional accumulation is the generally favored mechanism, but observations increasingly suggest that gas disk gravitational instability might explain the formation of at least the massive or wide-orbit gas giant exoplanets. This paper continues a series aimed at refining three-dimensional (3D) hydrodynamical models of disk instabilities, where the handling of the gas thermodynamics is a crucial factor. Boss (2017, 2019, 2021) used the $\beta$ cooling approximation (Gammie 2001) to calculate 3D models of disks with initial masses of 0.091 $M_\odot$ extending from 4 to 20 au around 1 $M_\odot$ protostars. Here we employ 3D flux-limited diffusion (FLD) approximation models of the same disks, in order to provide a superior treatment of disk gas thermodynamics. The new models have quadrupled spatial resolution compared to previous 3D FLD models (Boss 2008, 2012), in both the radial and azimuthal spherical coordinates, resulting in the highest spatial resolution 3D FLD models to date. The new models continue to support the hypothesis that such disks can form self-gravitating, dense clumps capable of contracting to form gas giant protoplanets, and suggest that the FLD models yield similar numbers of clumps as $\beta$ cooling models with $\beta \sim$ 1 to $\sim$ 10, including the critical value of $\beta$ = 3 for fragmentation proposed by Gammie (2001)., Comment: 26 pages, 7 figures, 3 tables, accepted by ApJ. arXiv admin note: text overlap with arXiv:2103.02566

To fully constrain the orbits of low mass circumstellar companions, we conduct combined analyses of the radial velocity data as well as the Gaia and Hipparcos astrometric data for eight nearby systems. Our study shows that companion-induced position and proper motion differences between Gaia and Hipparcos are significant enough to constrain orbits of low mass companions to a precision comparable with previous combined analyses of direct imaging and radial velocity data. We find that our method is robust to whether we use Gaia DR2 or Gaia EDR3, as well as whether we use all of the data, or just proper motion differences. In particular, we fully characterize the orbits of HD 190360 b and HD 16160 C for the first time. With a mass of 1.8$\pm$0.2$m_{\rm Jup}$ and an effective temperature of 123-176 K and orbiting around a Sun-like star, HD 190360 b is the smallest Jupiter-like planet with well-constrained mass and orbit, belonging to a small sample of fully characterized Jupiter analogs. It is separated from its primary star by 0.25$''$ and thus may be suitable for direct imaging by the CGI instrument of the Roman Space Telescope., Comment: Accepted for publication in MNRAS; 15 pages, 4 tables, 5 figures

Observations support the hypothesis that gas disk gravitational instability might explain the formation of massive or wide-orbit gas giant exoplanets. The situation with regard to Jupiter-mass exoplanets orbiting within $\sim$ 20 au is more uncertain. Theoretical models yield divergent assessments often attributed to the numerical handling of the gas thermodynamics. Boss (2019) used the $\beta$ cooling approximation to calculate three dimensional hydrodynamical models of the evolution of disks with initial masses of 0.091 $M_\odot$ extending from 4 to 20 au around 1 $M_\odot$ protostars. The models considered a wide range (1 to 100) of $\beta$ cooling parameters and started from an initial minimum Toomre stability parameter of $Q_i = 2.7$ (gravitationally stable). The disks cooled down from initial outer disk temperatures of 180 K to as low as 40 K as a result of the $\beta$ cooling, leading to fragmentation into dense clumps, which were then replaced by virtual protoplanets (VPs) and evolved for up to $\sim$ 500 yr. The present models test the viability of replacing dense clumps with VPs by quadrupling the spatial resolution of the grid once dense clumps form, sidestepping in most cases VP insertion. After at least $\sim$ 200 yr of evolution, the new results compare favorably with those of Boss (2019): similar numbers of VPs and dense clumps form by the same time for the two approaches. The results imply that VP insertion can greatly speed disk instability calculations without sacrificing accuracy., Comment: 25 pages, 3 tables, 7 figures, in press, ApJ

We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $\eta_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orbiting stars with effective temperatures between 4800 K and 6300 K. We find that $\eta_\oplus$ for the conservative HZ is between $0.37^{+0.48}_{-0.21}$ (errors reflect 68\% credible intervals) and $0.60^{+0.90}_{-0.36}$ planets per star, while the optimistic HZ occurrence is between $0.58^{+0.73}_{-0.33}$ and $0.88^{+1.28}_{-0.51}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates using both a Poisson likelihood Bayesian analysis and Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with $95\%$ confidence that, on average, the nearest HZ planet around G and K dwarfs is about 6 pc away, and there are about 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun., Comment: To appear in The Astronomical Journal

Several observations suggest that the Solar system has been located in a region affected by massive stellar feedback for at least a few Myr; these include detection of live $^{60}\text{Fe}$ in deep-sea archives and Antarctic snow, the broad angular distribution of $^{26}\text{Al}$ around the Galactic plane seen in all-sky $\gamma$-ray maps, and the all-sky soft X-ray background. However, our position inside the Galactic disc makes it difficult to fully characterise this environment, and our limited time baseline provides no information about its formation history or relation to large-scale Galactic dynamics. We explore these questions by using an $N$-body+hydrodynamics simulation of a Milky-Way-like galaxy to identify stars on Sun-like orbits whose environments would produce conditions consistent with those we observe. We find that such stars are uncommon but not exceptionally rare. These stars are found predominantly near the edges of spiral arms, and lie inside kpc-scale bubbles that are created by multiple generations of star formation in the arm. We investigate the stars' trajectories and find that the duration of the stay in the bubble ranges from 20 Myr to 90 Myr. The duration is governed by the crossing time of stars across the spiral arm. This is generally shorter than the bubble lifetime, which is $\sim 100$ Myr as a result of the continuous gas supply provided by the arm environment., Comment: 7 pages, 5 figures. Accepted for publication in MNRAS

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)., Comment: 54 pages, 2 tables, 19 figures; PEXO is available at 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., Comment: 24 pages, 2 tables, 8 figures, in press, ApJ

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.", Comment: 8 pages, 2 figures, Astro2020 decadal submission

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., Comment: Astro2020 Science White paper

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., Comment: Astro2020 Science White Paper

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., Comment: Astro2020 Science White Paper

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$., Comment: White paper for the Astro2020 decadal survey

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey

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., Comment: 14 pages including 3 figures

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., Comment: 31 pages, 2 tables, 13 figures, accepted by ApJ

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., Comment: Submitted to ApJ. Cleared by statistics referee. Incorporates response to science referee. Includes link to MCMC code

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., Comment: Accepted to The Astrophysical Journal

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., Comment: For the report with higher-resolution images and for updates and supporting documentations please visit http://eos-nexus.org/sag15/

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., Comment: 16 pages, 5 figures, 1 table, in press, Astronomical Journal

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., Comment: 24 pages, 8 figures, 1 table, accepted by ApJ

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., Comment: 40 pages, 10 figures, 1 table, accepted by ApJ

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., Comment: Accepted for publication in ApJ. 24 pages, 8 figures, 5 tables

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)., Comment: 22 pages, 2 tables, no figures

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., Comment: Corrected planet mass error (7.8 Mj --> 1.09 Mj, in agreement with previous studies)

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., Comment: 24 pages, 12 figures, 1 table, accepted by ApJ

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., Comment: 13 pages, 7 figures, to appear in Proc. SPIE 9148

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., Comment: 21 pages, 8 figures, 1 table, accepted, ApJ

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., Comment: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning

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)., Comment: 94 pages, 55 figures, 25 tables. Accepted by ApJS

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., Comment: Accepted for publication as a chapter in Protostars and Planets VI, to be published in 2014 by University of Arizona Press

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., Comment: 96 pages, 23 figures, accepted for publication in ApJ

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., Comment: Accepted for publication in ApJ

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., Comment: Accepted to ApJ on June 24, 2013. 67 pages, 17 figures, 12 tables

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., Comment: 28 pages, 4 tables, 12 figures, in press, ApJ

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., Comment: 84 pages, 26 figures, 10 tables. ApJ in press

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., Comment: 28 pages, 5 figures, 8 tables, in press, ApJ

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., Comment: 22 pages, 6 figures, Accepted for publication in ApJ

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., Comment: 31 pages, 10 figures, 4 tables, in press, ApJ

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., Comment: 27 pages, 7 figures, 3 tables - in press, ApJ

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., Comment: 23 pages, 5 B/W and 3 color figures, accepted to ApJ

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., Comment: A&A accepted. 10 pages, 7 figures, 8 tables - A few more astrometric epochs added. Trigonometric parallax refined

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., Comment: To appear on Science Express August 28, 11 pages, 3 figures, one table (main text), 56 pages, 28 figures, 10 tables