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.

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 \pi^5 G^2 \Sigma_p^3 / \Omega^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..

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.

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.

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.

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.

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.

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 \alpha-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 \alpha-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.

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.

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.

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.

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.

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.