Showing 1–18 of 18 results for author: Kretke, K

Planet Formation Theory in the Era of ALMA and Kepler: from Pebbles to Exoplanets

Authors: Joanna Drazkowska, Bertram Bitsch, Michiel Lambrechts, Gijs D. Mulders, Daniel Harsono, Allona Vazan, Beibei Liu, Chris W. Ormel, Katherine Kretke, Alessandro Morbidelli

Abstract: Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constr… ▽ More Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constraints on the planet formation process. In this chapter, we summarize the new information derived from the exoplanets population and the circumstellar disks observations. We present the new developments in planet formation theory, from dust evolution to the growth of planetary cores by accretion of planetesimals, pebbles, and gas. We review the state-of-the-art models for the formation of diverse planetary systems, including the population synthesis approach which is necessary to compare theoretical model outcomes to the exoplanet population. We emphasize that the planet formation process may not be spatially uniform in the disk and there are preferential locations for the formation of planetesimals and planets. Outside of these locations, a significant fraction of solids is not growing past the pebble-sizes. The reservoir of pebbles plays an important role in the growth of planetary cores in the pebble accretion process. The timescale of the emergence of massive planetary cores is an important aspect of the present models and it is likely that the cores within one disk form at different times. In addition, there is growing evidence that the first planetary cores start forming early, during the circumstellar disk buildup process. △ Less

Submitted 24 May, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 38 pages, 13 figures, review chapter accepted to Protostars and Planets VII, editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura, corrected label in Fig. 1

Resilience of Planetesimal Formation in Weakly-Reinforced Pressure Bumps

Authors: Daniel Carrera, Andrew Thomas, Jacob B. Simon, Matthew A. Small, Katherine A. Kretke, Hubert Klahr

Abstract: The discovery that axisymmetric dust rings are ubiquitous in protoplanetary disks has provoked a flurry of research on the role of pressure bumps in planet formation. High-resolution simulations by our group have shown that even a modest bump can collect enough dust to trigger planetesimal formation by the streaming instability. In this work, we probe the limits of planetesimal formation when the… ▽ More The discovery that axisymmetric dust rings are ubiquitous in protoplanetary disks has provoked a flurry of research on the role of pressure bumps in planet formation. High-resolution simulations by our group have shown that even a modest bump can collect enough dust to trigger planetesimal formation by the streaming instability. In this work, we probe the limits of planetesimal formation when the external source of pressure bump reinforcement is extremely weak. We conduct simulations of radially elongated shearing boxes to capture the entire bump, which is generated and maintained over some timescale $t_{\rm reinf}$ by a Newtonian relaxation scheme. We find that planetesimal formation is extremely resilient for cm-sized grains. We reduced the strength of reinforcement by up to a factor of 100 and the location and initial masses of planetesimals were essentially unaffected. However, we do find that strong reinforcement causes much faster pebble drift compared to the the standard pebble drift rates. The resulting larger pebble flux enhances the planetesimal growth rate by pebble accretion. We hypothesize that to sustain the bump, our code has to extract angular momentum (the strength of this negative torque depends on $t_{\rm reinf}$), and some of this torque is transferred to the particles, causing them to drift faster for a stronger torque (i.e., smaller $t_{\rm reinf}$). Since any physical process that sustains a pressure bump must do so by torquing the gas, we conjecture that the effect on pebble drift is a real phenomenon, motivating further work with physically realistic sources to generate the bump. △ Less

Submitted 1 February, 2022; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: 12 pages, 8 figures, accepted to ApJ

Linking planetary embryo formation to planetesimal formation II: The impact of pebble accretion in the terrestrial planet zone

Authors: Oliver Voelkel, Rogerio Deienno, Katherine Kretke, Hubert Klahr

Abstract: The accretion of pebbles on planetary cores has been widely studied in recent years and is found to be a highly effective mechanism for planetary growth. While most studies assume planetary cores as an initial condition in their simulation, the question how, where and when these cores form is often neglected. We study the impact of pebble accretion during the formation phase and subsequent evoluti… ▽ More The accretion of pebbles on planetary cores has been widely studied in recent years and is found to be a highly effective mechanism for planetary growth. While most studies assume planetary cores as an initial condition in their simulation, the question how, where and when these cores form is often neglected. We study the impact of pebble accretion during the formation phase and subsequent evolution of planetary embryos in the early stages of circumstellar disk evolution. In doing so we aim to quantify the timescales and local dependency of planetary embryo formation, based on the solid evolution of the disk. We connect a one dimensional two population model for solid evolution and pebble flux regulated planetesimal formation to the N-body code LIPAD. In our study we focus on the growth of planetesimals with an initial size of 100 km in diameter by planetesimal collisions and pebble accretion for the first 1 million years of a viscously evolving disk. We compare 18 different N-body simulations in which we vary the total planetesimal mass after 1 million years, the surface density profile of the planetesimal disk, the radial pebble flux and the possibility of pebble accretion. Pebble accretion leads to the formation of fewer, but substantially more massive embryos. The area of possible embryo formation is weakly influenced by the accretion of pebbles and the innermost embryos tend to form slightly earlier compared to the simulations in which pebble accretion is neglected. Pebble accretion strongly enhances the formation of super earths in the terrestrial planet region, but it does not enhance the formation of embryos at larger distances. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Linking planetary embryo formation to planetesimal formation I: The impact of the planetesimal surface density in the terrestrial planet zone

Authors: Oliver Voelkel, Rogerio Deienno, Katherine Kretke, Hubert Klahr

Abstract: The growth time scales of planetary embryos and their formation process are imperative for our understanding on how planetary systems form and develop. They determine the subsequent growth mechanisms during the life stages of a circumstellar disk. We quantify the timescales and spatial distribution of planetary embryos via collisional growth and fragmentation of dynamically forming 100km sized pla… ▽ More The growth time scales of planetary embryos and their formation process are imperative for our understanding on how planetary systems form and develop. They determine the subsequent growth mechanisms during the life stages of a circumstellar disk. We quantify the timescales and spatial distribution of planetary embryos via collisional growth and fragmentation of dynamically forming 100km sized planetesimals. In our study, the formation timescales of viscous disk evolution and planetesimal formation are linked to the formation of planetary embryos in the terrestrial planet zone. We connect a one dimensional model for viscous gas evolution, dust and pebble dynamics and pebble flux regulated planetesimal formation to the N-body code LIPAD. Our framework enables us to study the formation, growth, fragmentation and evolution of planetesimals with an initial size of 100km in diameter for the first million years of a viscous disk. Our study shows the effect of the planetesimal surface density evolution on the preferential location and timescales of planetary embryo formation. A one dimensional analytically derived model for embryo formation based on the local planetesimal surface density evolution is presented. This model manages to reproduce the spatial distribution, formation rate and total number of planetary embryos at a fraction of the computational cost of the N-body simulations. The formation of planetary embryos in the terrestrial planet zone occurs simultaneously to the formation of planetesimals. The local planetesimal surface density evolution and the orbital spacing of planetary embryos in the oligarchic regime serve well as constraints to model planetary embryo formation analytically. Our embryo formation model will be a valuable asset in future studies regarding planet formation. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Protoplanetary Disk Rings as Sites for Planetesimal Formation

Authors: Daniel Carrera, Jacob B. Simon, Rixin Li, Katherine A. Kretke, Hubert Klahr

Abstract: Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI… ▽ More Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those observed by ALMA. In particular, we place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of $Z = 0.01$, and we include both particle back-reaction and particle self-gravity. We consider both mm-sized and cm-sized particles separately. For simulations with cm-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the streaming instability; a pressure bump does {\it not} need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the streaming instability) is not responsible for planetesimal formation. For mm-sized particles, we find tentative evidence that planetesimal formation does not occur. This result, if it holds up at higher resolution and for a broader range of parameters, could put strong constraints on where in protoplanetary disks planetesimals can form. Ultimately, however, our results suggest that for cm-sized particles, planetesimal formation in pressure bumps is an extremely robust process. △ Less

Submitted 12 January, 2021; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: accepted to ApJ; 22 pages, 16 figures

Constraining the parameter space for the Solar Nebula

Authors: Christian T. Lenz, Hubert Klahr, Tilman Birnstiel, Katherine Kretke, Sebastian Stammler

Abstract: If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)pl… ▽ More If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)planets, once the gas disk is gone. Protoplanetary disks provide a plethora of different parameters to explore. The question is whether this parameter space can be constrained, allowing simulations to reproduce the Solar System. Models and observations of planet formation provide constraints on the initial planetesimal mass in certain regions of the Solar Nebula. By making use of pebble flux-regulated planetesimal formation, we perform a parameter study with nine different disk parameters like the initial disk mass, initial disk size, initial dust-to-gas ratio, turbulence level, and more. We find that the distribution of mass in planetesimals in the disk depends on the planetesimal formation timescale and the pebbles' drift timescale. Multiple disk parameters can influence pebble properties and thus planetesimal formation. However, it is still possible to draw some conclusions on potential parameter ranges. Pebble flux-regulated planetesimal formation seems to be very robust, allowing simulations with a wide range of parameters to meet the initial planetesimal constraints for the Solar Nebula. I.e., it does not require a lot of fine tuning. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: A&A accepted

The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. II. The Effect of the Radial Gas Pressure Gradient

Authors: Charles P. Abod, Jacob B. Simon, Rixin Li, Philip J. Armitage, Andrew N. Youdin, Katherine A. Kretke

Abstract: The streaming instability concentrates solid particles in protoplanetary disks, leading to gravitational collapse into planetesimals. Despite its key role in producing particle clumping and determining critical length scales in the instability's linear regime, the influence of the disk's radial pressure gradient on planetesimal properties has not been examined in detail. Here, we use streaming ins… ▽ More The streaming instability concentrates solid particles in protoplanetary disks, leading to gravitational collapse into planetesimals. Despite its key role in producing particle clumping and determining critical length scales in the instability's linear regime, the influence of the disk's radial pressure gradient on planetesimal properties has not been examined in detail. Here, we use streaming instability simulations that include particle self-gravity to study how the planetesimal initial mass function depends on the radial pressure gradient. Fitting our results to a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, we find a single value $p \approx 1.6$ describes simulations in which the pressure gradient varies by $\gtrsim 2$. An exponentially truncated power-law provides a significantly better fit, with a low mass slope of $p^\prime \approx 1.3$ that weakly depends on the pressure gradient. The characteristic truncation mass is found to be $\sim M_G = 4 π^5 G^2 Σ_p^3 / Ω^4$. We exclude the cubic dependence of the characteristic mass with pressure gradient suggested by linear considerations, finding instead a linear scaling. These results strengthen the case for a streaming-derived initial mass function that depends at most weakly on the aerodynamic properties of the disk and participating solids. A simulation initialized with zero pressure gradient---which is {\em not} subject to the streaming instability---also yields a top-heavy mass function but with modest evidence for a different shape. We discuss the consistency of the theoretically predicted mass function with observations of Kuiper Belt planetesimals, and describe implications for models of early stage planet formation.. △ Less

Submitted 29 August, 2019; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: 18 pages, 10 figures, 3 tables, accepted to ApJ

Growing the terrestrial planets from the gradual accumulation of sub-meter sized objects

Authors: Harold F. Levison, Katherine A. Kretke, Kevin Walsh, William Bottke

Abstract: Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants.… ▽ More Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants. This model envisions that the cores of the giant planets formed from 100 to 1000 km bodies that directly accreted a population of pebbles --- sub-meter sized objects that slowly grew in the protoplanetary disk. Here we apply this model to the terrestrial planet region and find that it can reproduce the basic structure of the inner Solar System, including a small Mars and a low-mass asteroid belt. Our models show that for an initial population of planetesimals with sizes similar to those of the main belt asteroids, VSPA becomes inefficient beyond $\sim\!$1.5 AU. As a result, Mars's growth is stunted and nothing large in the asteroid belt can accumulate. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Comments: Accepted for PNAS

Growing the gas-giant planets by the gradual accumulation of pebbles

Authors: Harold F. Levison, Katherine A. Kretke, Martin J. Duncan

Abstract: It is widely held that the first step in forming the gas giant planets, such as Jupiter and Saturn, is to form solid `cores' of roughly 10 M$_\oplus$. Getting the cores to form before the solar nebula dissipates ($\sim\!1-10\,$Myr) has been a major challenge for planet formation models. Recently models have emerged in which `pebbles' (centimeter- to meter-size objects) are first concentrated by ae… ▽ More It is widely held that the first step in forming the gas giant planets, such as Jupiter and Saturn, is to form solid `cores' of roughly 10 M$_\oplus$. Getting the cores to form before the solar nebula dissipates ($\sim\!1-10\,$Myr) has been a major challenge for planet formation models. Recently models have emerged in which `pebbles' (centimeter- to meter-size objects) are first concentrated by aerodynamic drag and then gravitationally collapse to form 100 --- 1000 km objects. These `planetesimals' can then efficiently accrete leftover pebbles and directly form the cores of giant planets. This model known as `pebble accretion', theoretically, can produce 10 M$_\oplus$ cores in only a few thousand years. Unfortunately, full simulations of this process show that, rather than creating a few 10 M$_\oplus$ cores, it produces a population of hundreds of Earth-mass objects that are inconsistent with the structure of the Solar System. Here we report that this difficulty can be overcome if pebbles form slowly enough to allow the planetesimals to gravitationally interact with one another. In this situation the largest planetesimals have time to scatter their smaller siblings out of the disk of pebbles, thereby stifling their growth. Our models show that, for a large, and physically reasonable region of parameter space, this typically leads to the formation of one to four gas giants between 5 and 15 AU in agreement with the observed structure of the Solar System. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Journal ref: Nature, Volume 524, Issue 7565, pp. 322-324 (2015)

Evidence for Pebbles in Comets

Authors: K. A. Kretke, H. F. Levison

Abstract: When the EPOXI spacecraft flew by Comet 103P/Hartley 2, it observed large particles floating around the comet nucleus. These particles are likely low-density, centimeter- to decimeter-sized clumps of ice and dust. While the origin of these objects remains somewhat mysterious, it is possible that they are giving us important information about the earliest stages of our Solar System's formation. Rec… ▽ More When the EPOXI spacecraft flew by Comet 103P/Hartley 2, it observed large particles floating around the comet nucleus. These particles are likely low-density, centimeter- to decimeter-sized clumps of ice and dust. While the origin of these objects remains somewhat mysterious, it is possible that they are giving us important information about the earliest stages of our Solar System's formation. Recent advancements in planet formation theory suggest that planetesimals (or cometestimals) may grow directly from the gravitational collapse of aerodynamically concentrated small particles, often referred to as "pebbles." Here we show that the particles observed in the coma of 103P are consistent with the sizes of pebbles expected to efficiently form planetesimals in the region that this comet likely formed, while smaller pebbles are may be expected in the majority of comets, whose chemistry is often indicative of formation in the colder, outer regions of the protoplanetary disk. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Journal ref: Icarus, Volume 262, December 2015, pp 9-13

Challenges in Forming the Solar System's Giant Planet Cores via Pebble Accretion

Authors: K. A. Kretke, H. F. Levison

Abstract: Though ~10 Earth mass rocky/icy cores are commonly held as a prerequisite for the formation of gas giants, theoretical models still struggle to explain how these embryos can form within the lifetimes of gaseous circumstellar disks. In recent years, aerodynamic-aided accretion of "pebbles," objects ranging from centimeters to meters in size, has been suggested as a potential solution to this long-s… ▽ More Though ~10 Earth mass rocky/icy cores are commonly held as a prerequisite for the formation of gas giants, theoretical models still struggle to explain how these embryos can form within the lifetimes of gaseous circumstellar disks. In recent years, aerodynamic-aided accretion of "pebbles," objects ranging from centimeters to meters in size, has been suggested as a potential solution to this long-standing problem. While pebble accretion has been demonstrated to be extremely effective in local simulations that look at the detailed behavior of these pebbles in the vicinity of a single planetary embryo, to date there have been no global simulations demonstrating the effectiveness of pebble accretion in a more complicated, multi-planet environment. Therefore, we have incorporated the aerodynamic-aided accretion physics into LIPAD, a Lagrangian code that can follow the collisional / accretional / dynamical evolution of a protoplanetary system, to investigate how pebble accretion manifests itself in the larger planet formation picture. We find that under generic circumstances, pebble accretion naturally leads to an "oligarchic" type of growth in which a large number of planetesimals grow to similar-sized planets. In particular, our simulations tend to form hundreds of Mars- and Earth-mass objects between 4 and 10 AU. While the merging of some oligarchs may allow them to grow massive enough to form giant planet cores, leftover oligarchs lead to planetary systems that cannot be consistent with our own solar system. We investigate various ideas presented in the literature (including evaporation fronts and planet traps) and find that none easily overcome this tendency toward oligarchic growth. △ Less

Submitted 13 October, 2014; v1 submitted 15 September, 2014; originally announced September 2014.

Comments: 20 pages, 33 figures, accepted in AJ

Stellar irradiated discs and implications on migration of embedded planets III: viscosity transitions

Authors: Bertram Bitsch, Alessandro Morbidelli, Elena Lega, Kathrine Kretke, Aurélien Crida

Abstract: The migration strength and direction of embedded low-mass planets depends on the disc structure. In discs with an efficient radiative transport, the migration can be directed outwards for planets with more than 3-5 Earth masses. This is due to the entropy driven corotation torque, a process that extends the lifetimes of growing planetary embryos. We investigate the influence on the disc structure… ▽ More The migration strength and direction of embedded low-mass planets depends on the disc structure. In discs with an efficient radiative transport, the migration can be directed outwards for planets with more than 3-5 Earth masses. This is due to the entropy driven corotation torque, a process that extends the lifetimes of growing planetary embryos. We investigate the influence on the disc structure caused by a jump in the alpha parameter of the viscosity to model a dead-zone structure in the disc. We focus on M-dot discs, which have a constant net mass flux. Using the resulting disc structure, we investigate the consequences for the formation of planetesimals and determine the regions of outward migration for proto-planets. We performed numerical hydrosimulations of M-dot discs in the r-z-plane. We used the explicit/implicit hydrodynamical code FARGOCA that includes a full tensor viscosity and stellar irradiation as well as a two-temperature solver that includes radiation transport in the flux-limited diffusion approximation. Viscosity transitions inside the disc create transitions in density that stop inward migration for small planets through the so-called "planet trap" mechanism. This mechanism also works for planets down to M_P > 0.5M_E. Additionally, the viscosity transitions change the pressure gradient in the disc, which facilitates planetesimal formation via the streaming instability. However, a very steep transition in viscosity is needed to achieve in a pressure bump in the disc. The transition in viscosity facilitates planetesimal formation and can stop the migration of small-mass planets (M_P>0.5M_E), but still does not halt inward migration of smaller planets and planetesimals that are affected by gas drag. A very steep, probably unrealistic viscosity gradient is needed to trap planets of smaller masses and halt gas-drag-driven planetesimal migration at a pressure bump. △ Less

Submitted 5 August, 2014; originally announced August 2014.

Comments: 12 pages, accepted by A&A

Journal ref: A&A 570, A75 (2014)

The Importance of Disk Structure in Stalling Type I Migration

Authors: Katherine A. Kretke, D. N. C. Lin

Abstract: As planets form they tidally interact with their natal disks. Though the tidal perturbation induced by Earth and super-Earth mass planets is generally too weak to significantly modify the structure of the disk, the interaction is potentially strong enough to cause the planets to undergo rapid type I migration. This physical process may provide a source of short-period super-Earths, though it may a… ▽ More As planets form they tidally interact with their natal disks. Though the tidal perturbation induced by Earth and super-Earth mass planets is generally too weak to significantly modify the structure of the disk, the interaction is potentially strong enough to cause the planets to undergo rapid type I migration. This physical process may provide a source of short-period super-Earths, though it may also pose a challenge to the emergence and retention of cores on long-period orbits with sufficient mass to evolve into gas giants. Previous numerical simulations have shown that the type I migration rate sensitively depends upon the circumstellar disk's properties, particularly the temperature and surface density gradients. Here, we derive these structure parameters for 1) a self-consistent viscous-disk model based on a constant α-prescription, 2) an irradiated disk model that takes into account heating due to the absorption of stellar photons, and 3) a layered-accretion disk model with variable α-parameter. We show that in the inner viscously-heated regions of typical protostellar disks, the horseshoe and corotation torques of super-Earths can exceed their differential Lindblad torque and cause them to undergo outward migration. However, the temperature profile due to passive stellar irradiation causes type I migration to be inwards throughout much of the disk. For disks in which there is outwards migration, we show that location and the mass range of the "planet traps" depends on some uncertain assumptions adopted for these disk models. Competing physical effects may lead to dispersion in super-Earths' mass-period distribution. △ Less

Submitted 17 May, 2012; originally announced May 2012.

Comments: 12 pages, Submitted to ApJ

A Method to Constrain the Size of the Protosolar Nebula

Authors: K. A. Kretke, H. F. Levison, M. W. Buie, A. Morbidelli

Abstract: Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own pro… ▽ More Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk. △ Less

Submitted 10 February, 2012; originally announced February 2012.

Comments: 10 pages, Accepted to AJ

Teaching Astronomy with an Inquiry Activity on Stellar Populations

Authors: Marc Rafelski, Michael Foley, Genevieve J. Graves, Katherine A. Kretke, Elisabeth Mills, Michael Nassir, Shannon Patel

Abstract: We describe a new inquiry design aimed at teaching advanced high-school to senior college students the basics of stellar populations. The inquiry is designed to have students come up with their own version of the Hertzsprung-Russell diagram as a tool to understand how stars evolve based on their color, mass, and luminosity. The inquiry makes use of pictures and spectra of stars, which the students… ▽ More We describe a new inquiry design aimed at teaching advanced high-school to senior college students the basics of stellar populations. The inquiry is designed to have students come up with their own version of the Hertzsprung-Russell diagram as a tool to understand how stars evolve based on their color, mass, and luminosity. The inquiry makes use of pictures and spectra of stars, which the students analyze and interpret to answer the questions they come up with at the beginning. The students undergo a similar experience to real astronomers, using the same tools and methods to figure out the phenomena they are trying to understand. Specifically, they use images and spectra of stars, and organize the data via tables and plots to find trends that will then enable them to answer their questions. The inquiry also includes a "thinking tool" to help connect the trends students observe to the larger picture of stellar evolution. We include a description of the goals of the inquiry, the activity description, the motivations and thoughts that went into the design of the inquiry, and reflections on how the inquiry activity worked in practice. △ Less

Submitted 27 September, 2010; originally announced September 2010.

Comments: 11 pages, 2 figure, 2 tables. ASP Conference Series 436, Learning from Inquiry in Practice, in press

Journal ref: ASP Conference Series 436 (2010) 108

Assembling the Building Blocks of Giant Planets around Intermediate Mass Stars

Authors: K. A. Kretke, D. N. C. Lin, P. Garaud, N. J. Turner

Abstract: We examine a physical process that leads to the efficient formation of gas giant planets around intermediate mass stars. In the gaseous protoplanetary disks surrounding rapidly-accreting intermediate-mass stars we show that the midplane temperature (heated primarily by turbulent dissipation) can reach > 1000 K out to 1 AU. Thermal ionization of this hot gas couples the disk to the magnetic field… ▽ More We examine a physical process that leads to the efficient formation of gas giant planets around intermediate mass stars. In the gaseous protoplanetary disks surrounding rapidly-accreting intermediate-mass stars we show that the midplane temperature (heated primarily by turbulent dissipation) can reach > 1000 K out to 1 AU. Thermal ionization of this hot gas couples the disk to the magnetic field, allowing the magneto-rotational instability (MRI) to generate turbulence and transport angular momentum. Further from the central star the ionization fraction decreases, decoupling the disk from the magnetic field and reducing the efficiency of angular momentum transport. As the disk evolves towards a quasi-steady state, a local maximum in the surface density and in the midplane pressure both develop at the inner edge of the MRI-dead zone, trapping inwardly migrating solid bodies. Small particles accumulate and coagulate into planetesimals which grow rapidly until they reach isolation mass. In contrast to the situation around solar type stars, we show that the isolation mass for cores at this critical radius around the more massive stars is large enough to promote the accretion of significant amounts of gas prior to disk depletion. Through this process, we anticipate a prolific production of gas giants at ~1 AU around intermediate-mass stars. △ Less

Submitted 9 June, 2008; originally announced June 2008.

Comments: 24 pages with 5 figures, preprint format. Submitted to ApJ

Journal ref: Astrophys.J.690:407-415,2009

Grain Retention and Formation of Planetesimals near the Snow Line in MRI-driven Turbulent Protoplanetary Disks

Authors: Katherine A. Kretke, D. N. C. Lin

Abstract: The first challenge in the formation of both terrestrial planets and the cores of gas giants is the retention of grains in protoplanetary disks. In most regions of these disks, gas attains sub-Keplerian speeds as a consequence of a negative pressure gradient. Hydrodynamic drag leads to orbital decay and depletion of the solid material in the disk, with characteristic timescales as short as only… ▽ More The first challenge in the formation of both terrestrial planets and the cores of gas giants is the retention of grains in protoplanetary disks. In most regions of these disks, gas attains sub-Keplerian speeds as a consequence of a negative pressure gradient. Hydrodynamic drag leads to orbital decay and depletion of the solid material in the disk, with characteristic timescales as short as only a few hundred years for meter-sized objects at 1 AU. In this paper, we suggest a particle retention mechanism which promotes the accumulation of grains and the formation of planetesimals near the water sublimation front or ``snow line.'' This model is based on the assumption that, in the regions most interesting for planet formation, the viscous evolution of the disk is due to turbulence driven by the magneto-rotational instability (MRI) in the surface layers of the disk. The depth to which MRI effectively generates turbulence is a strong function of the grain size and abundance. A sharp increase in the grain-to-gas density ratio across the snow line reduces the column depth of the active layer. As the disk evolves towards a quasi-steady-state, this change in the active layer creates a local maximum in radial distribution of the gas surface density and pressure, causing the gas to rotate at super-Keplerian speed and halting the inward migration of grains. This senario presents a robust process for grain retention which may aid in the formation of proto-gas-giant cores preferentially near the snow line. △ Less

Submitted 8 June, 2007; originally announced June 2007.

Comments: 12 pages in preprint format, 2 figures, accepted ApJ Letters

The N2K Consortium. VII. Atmospheric Parameters of 1907 Metal-Rich Stars: Finding Planet-Search Targets

Authors: Sarah E. Robinson, S. Mark Ammons, Katherine A. Kretke, Jay Strader, Jeremy G. Wertheimer, Debra A. Fischer, Gregory Laughlin

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

Submitted 22 January, 2007; originally announced January 2007.

Comments: Accepted for publication in Astrophysical Journal Supplements. 24 pages, including 8 figures and one table. Two more tables will appear in ApJS online edition