The streaming instability for solid particles in protoplanetary disks is re-examined assuming the familiar alpha ($\alpha$) model for isotropic turbulence. Turbulence always reduces the growth rates of the streaming instability relative to values calculated for globally laminar disks. While for small values of the turbulence parameter, $\alpha < 10^{-5}$, the wavelengths of the fastest-growing disturbances are small fractions of the local gas vertical scale height $H$, we find that for moderate values of the turbulence parameter, i.e., $\alpha \sim 10^{-5}-10^{-3}$, the lengthscales of maximally growing disturbances shift toward larger scales, approaching $H$. At these moderate turbulent intensities and for local particle to gas mass density ratios $\epsilon < 0.5$, the vertical scales of the most unstable modes begin to exceed the corresponding radial scales so that the instability appears in the form of vertically oriented sheets extending well beyond the particle scale height. We find that for hydrodynamical turbulent disk models reported in the literature, with $\alpha = 4\times 10^{-5} - 5\times 10^{-4}$, together with state of the art global evolution models of particle growth, the streaming instability is predicted to be viable within a narrow triangular patch of $\alpha$--$\tau_s$ parameter space centered on Stokes numbers, $\tau_s \sim 0.01$ and $\alpha \sim 4\times 10^{-5}$ and, further, exhibits growth rates on the order of several hundred to thousands of orbit times for disks with 1 percent ($Z= 0.01$) cosmic solids abundance or metallicity. Our results are consistent with, and place in context, published numerical studies of streaming instabilities.

We use a newly developed cascade model of turbulent concentration of particles in protoplanetary nebulae to calculate several properties of interest to the formation of primitive planetesimals and to the meteorite record. The model follows, and corrects, calculations of the primary planetesimal Initial Mass Function (IMF) by Cuzzi et al. (2010), in which an incorrect cascade model was used. Here we use the model of Hartlep et al. (2017), which has been validated against several published numerical simulations of particle concentration in turbulence. We find that, for a range of nebula and particle properties, planetesimals may be "born big", formed as sandpiles with diameters in the 10-100 km range, directly from freely floating particles. The IMFs have a modal nature, with a well-defined peak rather than a powerlaw size dependence. Predictions for the inner and outer nebula behave similarly in these regards, and observations of inner and outer nebula primitive bodies support such modal IMFs. Also, we present predictions of local particle concentrations on several lengthscales in which particles "commonly" find themselves, which have significance for meteoritical observations of the redox state and isotopic fractionation in regions of chondrule formation. An important difference between these results, and those of Cuzzi et al. (2010), is that particle growth-by-sticking must proceed to at least the 1-few cm radius range for the IMF and meteoritical properties to be most plausibly satisfied. That is, as far as the inner nebula goes, the predominant "particles" must be aggregates of chondrules (or chondrule-size precursors) rather than individual chondrules themselves.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

Quantification of the radial transport of gaseous species and solid particles is important to many applications in protoplanetary disk evolution. An especially important example is determining the location of the water snow lines in a disk, which requires computing the rates of outward radial diffusion of water vapor and the inward radial drift of icy particles; however, the application is generalized to evaporation fronts of all volatiles. We review the relevant formulas using a uniform formalism. This uniform treatment is necessary because the literature currently contains at least six mutually exclusive treatments of radial diffusion of gas, only one of which is correct. We derive the radial diffusion equations from first principles, using Fick's law. For completeness, we also present the equations for radial transport of particles. These equations may be applied to studies of diffusion of gases and particles in protoplanetary and other accretion disks.

A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability.

We model particle growth in a turbulent, viscously evolving protoplanetary nebula, incorporating sticking, bouncing, fragmentation, and mass transfer at high speeds. We treat small particles using a moments method and large particles using a traditional histogram binning, including a probability distribution function of collisional velocities. The fragmentation strength of the particles depends on their composition (icy aggregates are stronger than silicate aggregates). The particle opacity, which controls the nebula thermal structure, evolves as particles grow and mass redistributes. While growing, particles drift radially due to nebula headwind drag. Particles of different compositions evaporate at "evaporation fronts" (EFs) where the midplane temperature exceeds their respective evaporation temperatures. We track the vapor and solid phases of each component, accounting for advection and radial and vertical diffusion. We present characteristic results in evolutions lasting $2 \times 10^5$ years. In general, (a) mass is transferred from the outer to inner nebula in significant amounts, creating radial concentrations of solids at EFs, (b) particle sizes are limited by a combination of fragmentation, bouncing, and drift, (c) "lucky" large particles never represent a significant amount of mass, and (d) restricted radial zones just outside each EF become compositionally enriched in the associated volatiles. We point out implications for mm-submm SEDs and inference of nebula mass, radial banding, the role of opacity on new mechanisms for generating turbulence, enrichment of meteorites in heavy oxygen isotopes, variable and nonsolar redox conditions, primary accretion of silicate and icy planetesimals, and the makeup of Jupiter's core.

The James Webb Space Telescope (JWST) will provide unprecedented opportunities to observe the rings and small satellites in our solar system, accomplishing three primary objectives: 1) discovering new rings and moons, 2) unprecedented spectroscopy, and 3) time-domain observations. We give details on these science objectives and describe requirements that JWST must fulfill in order to accomplish the science objectives

Asteroids and meteorites provide key evidence on the formation of planetesimals in the Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation within a population of initially km-scale planetesimals. However, new models challenge this idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the gravitational collapse of small particles which have organised themselves in dense filaments and clusters in the turbulent gas. Particles concentrate passively between eddies down to the smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The streaming instability causes particles to take an active role in the concentration, by piling up in dense filaments whose friction on the gas reduces the radial drift compared to that of isolated particles. In this chapter we review new paradigms for asteroid formation and compare critically against the observed properties of asteroids as well as constraints from meteorites. Chondrules of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. We review recent progress on understanding the incorporation of chondrules into the asteroids, including layered accretion models where chondrules are accreted onto asteroids over millions of years. We highlight in the end ten unsolved questions in asteroid formation where we expect that progress will be made over the next decade.

Several mechanisms have been identified that create dense particle clumps in the solar nebula. The present work is concerned with the gravitational collapse of such clumps, idealized as being spherically symmetric. Calculations using the two-fluid model are performed (almost) up to the time when a central density singularity forms. The end result of the study is a parametrization for this time, in order that it may be compared with timescales for various disruptive effects to which clumps may be subject. An important effect is that as the clump compresses, it also compresses the gas due to drag. This increases gas pressure which retards particle collapse and leads to oscillation in the size and density of the clump. The ratio of gravitational force to gas pressure gives a two-phase Jeans parameter, $J_t$, which is the classical Jeans parameter with the sound speed replaced by an the wave speed in a coupled two-fluid medium. Its use makes the results insensitive to the initial density ratio of particles to gas as a separate parameter. An ordinary differential equation model is developed which takes the form of two coupled non-linear oscillators and reproduces key features of the simulations. Finally, a parametric study of the time to collapse is performed and a formula (fit to the simulations) is developed. In the incompressible limit $J_t \to 0$, collapse time equals sedimentation time. As $J_t$ increases, the collapse time decreases roughly linearly with $J_t$ until $J_t \gtrsim 0.4$ when it becomes approximately equal to the dynamical time.

The Goldreich and Ward (1973) (axisymmetric) gravitational instability of a razor thin particle layer occurs when the Toomre parameter $Q_T \equiv c_p \Omega_0 / \pi G \Sigma_p < 1$ ($c_p$ being the particle dispersion velocity). Ward(1976,2000) extended this analysis by adding the effect of gas drag upon particles and found that even when $Q_T > 1$, sufficiently long waves were always unstable. Youdin (2005a,b) carried out a detailed analysis and showed that the instability allows chondrule-sized ($\sim 1 $ mm) particles to undergo radial clumping with reasonable growth times even in the presence of a moderate amount of turbulent stirring. The analysis of Youdin includes the role of turbulence in setting the thickness of the dust layer and in creating a turbulent particle pressure in the momentum equation. However, he ignores the effect of turbulent mass diffusivity on the disturbance wave. Here we show that including this effect reduces the growth-rate significantly, by an amount that depends on the level of turbulence, and reduces the maximum intensity of turbulence the instability can withstand by 1 to 3 orders of magnitude. The instability is viable only when turbulence is extremely weak and the solid to gas surface density of the particle layer is considerably enhanced over minimum-mass-nebula values. A simple mechanistic explanation of the instability shows how the azimuthal component of drag promotes instability while the radial component hinders it. A gravito-diffusive overstability is also possible but never realized in the nebula models.

We present an analytical model for the transmission spectrum of a transiting exoplanet, showing that a cloud base can produce an observable inflection point in the spectrum. The wavelength and magnitude of the inflection can be used to break the degeneracy between the atmospheric pressure and the abundance of the main cloud material, however, the abundance still depends on cloud particle size. An observed inflection also provides a specific point on the atmospheric P-T profile, giving us a "thermometer" to directly validate or rule out postulated cloud species. We apply the model to the transit spectrum of HD 189733b.

As small solid grains grow into larger ones in protoplanetary nebulae, or in the cloudy atmospheres of exoplanets, they generally form porous aggregates rather than solid spheres. A number of previous studies have used highly sophisticated schemes to calculate opacity models for irregular, porous particles with size much smaller than a wavelength. However, mere growth itself can affect the opacity of the medium in far more significant ways than the detailed compositional and/or structural differences between grain constituents once aggregate particle sizes exceed the relevant wavelengths. This physics is not new; our goal here is to provide a model that provides physical insight and is simple to use in the increasing number of protoplanetary nebula evolution, and exoplanet atmosphere, models appearing in recent years, yet quantitatively captures the main radiative properties of mixtures of particles of arbitrary size, porosity, and composition. The model is a simple combination of effective medium theory with small-particle closed-form expressions, combined with suitably chosen transitions to geometric optics behavior. Calculations of wavelength-dependent emission and Rosseland mean opacity are shown and compared with Mie theory. The model's fidelity is very good in all comparisons we have made, except in cases involving pure metal particles or monochromatic opacities for solid particles with size comparable to the wavelength.

Over the last eight years, the Visual and Infrared Mapping Spectrometer (VIMS) aboard the Cassini orbiter has returned hyperspectral images in the 0.35-5.1 micron range of the icy satellites and rings of Saturn. These very different objects show significant variations in surface composition, roughness and regolith grain size as a result of their evolutionary histories, endogenic processes and interactions with exogenic particles. The distributions of surface water ice and chromophores, i.e. organic and non-icy materials, across the saturnian system, are traced using specific spectral indicators (spectral slopes and absorption band depths) obtained from rings mosaics and disk-integrated satellites observations by VIMS.

Clouds and hazes are commonplace in the atmospheres of solar system planets and are likely ubiquitous in the atmospheres of extrasolar planets as well. Clouds affect every aspect of a planetary atmosphere, from the transport of radiation, to atmospheric chemistry, to dynamics and they influence - if not control - aspects such as surface temperature and habitability. In this review we aim to provide an introduction to the role and properties of clouds in exoplanetary atmospheres. We consider the role clouds play in influencing the spectra of planets as well as their habitability and detectability. We briefly summarize how clouds are treated in terrestrial climate models and consider the far simpler approaches that have been taken so far to model exoplanet clouds, the evidence for which we also review. Since clouds play a major role in the atmospheres of certain classes of brown dwarfs we briefly discuss brown dwarf cloud modeling as well. We also review how the scattering and extinction efficiencies of cloud particles may be approximated in certain limiting cases of small and large particles in order to facilitate physical understanding. Since clouds play such important roles in planetary atmospheres, cloud modeling may well prove to be the limiting factor in our ability to interpret future observations of extrasolar planets.

Saturn's main rings exhibit variations in both their opacity and spectral properties on a broad range of spatial scales, and the correlations between these parameters can provide insights into the processes that shape the composition and dynamics of the rings. The Visual and Infrared Mapping Spectrometer (VIMS) instrument onboard the Cassini Spacecraft has obtained spectra of the rings between 0.35 and 5.2 microns with sufficient spatial resolution to discern variations on scales below 200 km. These relatively high-resolution spectral data reveal that both the depths of the near-infrared water-ice absorption bands and the visible spectral slopes are often correlated with structural parameters such as the rings' optical depth. Using a simplified model for the ring-particles' regolith properties, we have begun to disentangle the trends due to changes in the gross composition of the ring particles from those that may be due to shifts in the texture of the ring particles' regolith. Consistent with previous studies, this analysis finds that the C ring and the Cassini Division possess enhanced concentrations of a contaminant that absorbs light over a broad range of wavelengths. On the other hand, a second contaminant that preferentially absorbs at short visible and near-ultraviolet wavelengths is found to be more evenly distributed throughout the rings. The optical activity of this short-wavelength absorber increases in the inner B ring inwards of 100,000 km from Saturn center, which may provide clues to the origin of this contaminant. The spectral variations identified as shifts in the regolith texture are in some places clearly correlated with the ring's optical depth, and in other locations they appear to be associated with the disturbances generated by strong mean-motion resonances with Saturn's various moons.

In the last few years Cassini-VIMS, the Visible and Infared Mapping Spectrometer, returned to us a comprehensive view of the Saturn's icy satellites and rings. After having analyzed the satellites' spectral properties (Filacchione et al. (2007a)) and their distribution across the satellites' hemispheres (Filacchione et al. (2010)), we proceed in this paper to investigate the radial variability of icy satellites (principal and minor) and main rings average spectral properties. This analysis is done by using 2,264 disk-integrated observations of the satellites and a 12x700 pixels-wide rings radial mosaic acquired with a spatial resolution of about 125 km/pixel. The comparative analysis of these data allows us to retrieve the amount of both water ice and red contaminant materials distributed across Saturn's system and the typical surface regolith grain sizes. These measurements highlight very striking differences in the population here analyzed, which vary from the almost uncontaminated and water ice-rich surfaces of Enceladus and Calypso to the metal/organic-rich and red surfaces of Iapetus' leading hemisphere and Phoebe. Rings spectra appear more red than the icy satellites in the visible range but show more intense 1.5-2.0 micron band depths. The correlations among spectral slopes, band depths, visual albedo and phase permit us to cluster the saturnian population in different spectral classes which are detected not only among the principal satellites and rings but among co-orbital minor moons as well. Finally, we have applied Hapke's theory to retrieve the best spectral fits to Saturn's inner regular satellites using the same methodology applied previously for Rhea data discussed in Ciarniello et al. (2011).

The variation in sizes of chondrules from one chondrite to the next is thought to be due to some sorting process in the early solar nebula. Hypotheses for the sorting process include chondrule sorting by mass and sorting by some aerodynamic mechanism; one such aerodynamic mechanism is the process of turbulent concentration (TC). We present the results of a series of statistical tests of chondrule data from several different chondrites. The data do not clearly distinguish between various options for the sorting parameter, but we find that the data are inconsistent with being drawn from lognormal or (three-parameter) Weibull distributions in chondrule radius. We also find that all but one of the chondrule data sets tested are consistent with being drawn from the TC distribution.

We have conducted an intensive search for any material that may orbit Rhea, using images obtained by the Cassini ISS narrow-angle camera. We find no evidence for any such material, contradicting an earlier and surprising inference that Rhea, the second-largest moon of Saturn, possesses a system of narrow rings embedded in a broad circum-satellite disk or cloud (Jones et al. 2008, Science). If one accepts the calculations of the previous authors, our results would require any narrow rings around Rhea to be composed of objects no smaller than 10 meters in radius, in order to be massive enough to account for the observed charged-particle absorptions but diffuse enough to evade detection in our images. Such a particle-size distribution is unrealistic given that erosion processes will break down 10-meter objects to smaller sizes, which would then have been seen in our images. Furthermore, Jones et al. (2008) assumed that the absorption of electrons by hypothetical circum-Rhea particles is proportional to the volume of the particles, even when they are much larger than the electron penetration depth (van Allen 1983, 1987, JGR). If we assume instead that large particles can only absorb electrons efficiently near their surfaces, then our results rule out any solid material orbiting Rhea as the explanation for the absorptions detected by Jones et al. (2008), regardless of particle size, by a margin of four (for narrow rings) or two (for a broad diffuse cloud) orders of magnitude. We conclude that the anomalous electron absorptions reported by Jones et al. (2008) are the signature of a new and hitherto unknown magnetospheric process.

Our goal is to understand primary accretion of the first planetesimals. The primitive meteorite record suggests that sizeable planetesimals formed in the asteroid belt over a period longer than a million years, each composed entirely of an unusual, but homogeneous, mixture of mm-size particles. We sketch a scenario in which primary accretion of 10-100km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted mm-size particles (generically "chondrules"). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100km diameter, which has been identified as a "fossil" property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer solar system planetesimals, which may also have a preferred size in the 100km diameter range. We estimate formation rates of planetesimals and assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For nebula parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the "fossil" sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion (the model naturally helps explain the peculiar size distribution of chondrules within such objects). The scenario also produces 10-100km diameter primary KBOs. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between keplerian and pressure-supported orbital velocities, than "canonical" models of the solar nebula. We discuss several potential explanations for these differences.

An updated Science Vision for the SOFIA project is presented, including an overview of the characteristics and capabilities of the observatory and first generation instruments. A primary focus is placed on four science themes: 'The Formation of Stars and Planets', 'The Interstellar Medium of the Milky Way', 'Galaxies and the Galactic Center' and 'Planetary Science'.

We demonstrate an approach to solving the coagulation equation that involves using a finite number of moments of the particle size distribution. This approach is particularly useful when only general properties of the distribution, and their time evolution, are needed. The numerical solution to the integro-differential Smoluchowski coagulation equation at every time step, for every particle size, and at every spatial location is computationally expensive, and serves as the primary bottleneck in running evolutionary models over long periods of time. The advantage of using the moments method comes in the computational time savings gained from only tracking the time rate of change of the moments, as opposed to tracking the entire mass histogram which can contain hundreds or thousands of bins depending on the desired accuracy. The collision kernels of the coagulation equation contain all the necessary information about particle relative velocities, cross-sections, and sticking coefficients. We show how arbitrary collision kernels may be treated. We discuss particle relative velocities in both turbulent and non-turbulent regimes. We present examples of this approach that utilize different collision kernels and find good agreement between the moment solutions and the moments as calculated from direct integration of the coagulation equation. As practical applications, we demonstrate how the moments method can be used to track the evolving opacity, and also indicate how one may incorporate porous particles.

We outline a scenario which traces a direct path from freely-floating nebula particles to the first 10-100km-sized bodies in the terrestrial planet region, producing planetesimals which have properties matching those of primitive meteorite parent bodies. We call this "primary accretion". The scenario draws on elements of previous work, and introduces a new critical threshold for planetesimal formation. We presume the nebula to be weakly turbulent, which leads to dense concentrations of aerodynamically size-sorted particles having properties like those observed in chondrites. The fractional volume of the nebula occupied by these dense zones or clumps obeys a probability distribution as a function of their density, and the densest concentrations have particle mass density 100 times that of the gas. However, even these densest clumps are prevented by gas pressure from undergoing gravitational instability in the traditional sense (on a dynamical timescale). While in this state of arrested development, they are susceptible to disruption by the ram pressure of the differentially orbiting nebula gas. However, self-gravity can preserve sufficiently large and dense clumps from ram pressure disruption, allowing their entrained particles to sediment gently but inexorably towards their centers, producing 10-100 km "sandpile" planetesimals. Localized radial pressure fluctuations in the nebula, and interactions between differentially moving dense clumps, will also play a role that must be allowed for in future studies. The scenario is readily extended from meteorite parent bodies to primary accretion throughout the solar system.

We present a mechanism for chondrules to stick together by means of compaction of a porous dust rim they sweep up as they move through the dusty nebula gas. It is shown that dust aggregates formed out of micron-sized grains stick to chondrules, forming a porous dust rim. When chondrules collide, this dust can be compacted by means of rolling motions within the porous dust layer. This mechanism dissipates the collisional energy, compacting the rim and allowing chondrules to stick. The structure of the obtained chondrule-dust agglomerates (referred to as compounds) then consists of three phases: chondrules, porous dust, and dust that has been compacted by collisions. Subsequently, these compounds accrete their own dust and collide with other compounds. The evolution of the compound size distribution and the relative importance of the phases is calculated by a Monte Carlo code. Growth ends, and a simulation is terminated when all the dust in the compounds has been compacted. Numerous runs are performed, reflecting the uncertainty in the physical conditions at the chondrule formation time. It is found that compounds can grow by 1-2 orders of magnitudes in radius, upto dm-sizes when turbulence levels are low. However, relative velocities associated with radial drift form a barrier for further growth. Earlier findings that the dust sweep-up by chondrules is proportional to their sizes are confirmed. We contrast two scenarios regarding how this dust evolved further towards the densely packed rims seen in chondrites.

A cascade model is described based on multiplier distributions determined from 3D direct numerical simulations (DNS) of turbulent particle laden flows, which include two-way coupling between the phases at global mass loadings equal to unity. The governing Eulerian equations are solved using pseudo-spectral methods on up to 512**3 computional grid points. DNS results for particle concentration and enstrophy at Taylor microscale Reynolds numbers in the range 34 - 170 were used to directly determine multiplier distributions (PDFs) on spatial scales 3 times the Kolmogorov length scale. The width of the PDFs, which is a measure of intermittency, decreases with increasing mass loading within the local region where the multipliers are measured. The functional form of this dependence is not sensitive to Reynolds numbers in the range considered. A partition correlation probability is included in the cascade model to account for the observed spatial anticorrelation between particle concentration and enstrophy. Joint probability distribution functions of concentration and enstrophy generated using the cascade model are shown to be in excellent agreement with those derived directly from our 3D simulations. Probabilities predicted by the cascade model are presented at Reynolds numbers well beyond what is achievable by direct simulation. These results clearly indicate that particle mass loading significantly reduces the probabilities of high particle concentration and enstrophy relative to those resulting from unloaded runs. Particle mass density appears to reach a limit at around 100 times the gas density. This approach has promise for significant computational savings in certain applications.

In this note we present complete, closed-form expressions for random relative velocities between colliding particles of arbitrary size in nebula turbulence. These results are exact for very small particles (those with stopping times much shorter than the large eddy overturn time) and are also surprisingly accurate in complete generality (that is, also apply for particles with stopping times comparable to, or much longer than, the large eddy overturn time). We note that some previous studies may have adopted previous simple expressions, which we find to be in error regarding the size dependence in the large particle regime.

We discuss the results of laboratory measurements and theoretical models concerning the aggregation of dust in protoplanetary disks, as the initial step toward planet formation. Small particles easily stick when they collide and form aggregates with an open, often fractal structure, depending on the growth process. Larger particles are still expected to grow at collision velocities of about 1m/s. Experiments also show that, after an intermezzo of destructive velocities, high collision velocities above 10m/s on porous materials again lead to net growth of the target. Considerations of dust-gas interactions show that collision velocities for particles not too different in surface-to-mass ratio remain limited up to sizes about 1m, and growth seems to be guaranteed to reach these sizes quickly and easily. For meter sizes, coupling to nebula turbulence makes destructive processes more likely. Global aggregation models show that in a turbulent nebula, small particles are swept up too fast to be consistent with observations of disks. An extended phase may therefore exist in the nebula during which the small particle component is kept alive through collisions driven by turbulence which frustrates growth to planetesimals until conditions are more favorable for one or more reasons.

(Abridged) Astronomical observations have shown that protoplanetary disks are dynamic objects through which mass is transported and accreted by the central star. Age dating of meteorite constituents shows that their creation, evolution, and accumulation occupied several Myr, and over this time disk properties would evolve significantly. Moreover, on this timescale, solid particles decouple from the gas in the disk and their evolution follows a different path. Here we present a model which tracks how the distribution of water changes in an evolving disk as the water-bearing species experience condensation, accretion, transport, collisional destruction, and vaporization. Because solids are transported in a disk at different rates depending on their sizes, the motions will lead to water being concentrated in some regions of a disk and depleted in others. These enhancements and depletions are consistent with the conditions needed to explain some aspects of the chemistry of chondritic meteorites and formation of giant planets. The levels of concentration and depletion, as well as their locations, depend strongly on the combined effects of the gaseous disk evolution, the formation of rapidly migrating rubble, and the growth of immobile planetesimals. We present examples of evolution under a range of plausible assumptions and demonstrate how the chemical evolution of the inner region of a protoplanetary disk is intimately connected to the physical processes which occur in the outer regions.

Changes in the chemical and isotopic composition of the solar nebula with time are reflected in the properties of different constituents that are preserved in chondritic meteorites. CR carbonaceous chondrites are among the most primitive of all chondrite types and must have preserved solar nebula records largely unchanged. We have analyzed the oxygen and magnesium isotopes in a range of the CR constituents of different formation temperatures and ages, including refractory inclusions and chondrules of various types. The results provide new constraints on the time variation of the oxygen isotopic composition of the inner (<5 AU) solar nebula - the region where refractory inclusions and chondrules most likely formed. A chronology based on the decay of short-lived 26Al (t1/2 ~ 0.73 Ma) indicates that the inner solar nebula gas was 16O-rich when refractory inclusions formed, but less than 0.8 Ma later, gas in the inner solar nebula became 16O-poor and this state persisted at least until CR chondrules formed ~1-2 Myr later. We suggest that the inner solar nebula became 16O-poor because meter-size icy bodies, which were enriched in 17,18O due to isotopic self-shielding during the ultraviolet photo dissociation of CO in the protosolar molecular cloud or protoplanetary disk, agglomerated outside the snowline, drifted rapidly towards the Sun, and evaporated at the snowline. This led to significant enrichment in 16O-depleted water, which then spread through the inner solar system. Astronomical studies of the spatial and/or temporal variations of water abundance in protoplanetary disks may clarify these processes.

Solid material in a protoplanetary nebula is subject to vigorous redistribution processes relative to the nebula gas. Meter-sized particles drift rapidly inwards near the nebula midplane, and material evaporates when the particles cross a condensation/evaporation boundary. The material cannot be removed as fast in its vapor form as it is being supplied in solid form, so its concentration increases locally by a large factor (more than an order of magnitude under nominal conditions). As time goes on, the vapor phase enhancement propagates for long distances inside the evaporation boundary (potentially all the way in to the star). Meanwhile, material is enhanced in its solid form over a characteristic lengthscale outside the evaporation boundary. This effect is applicable to any condensible (water, silicates, {\it etc.}). Three distinct radial enhancement/depletion regimes can be discerned by use of a simple model. Meteoritics applications include oxygen fugacity and isotopic variations, as well as isotopic homogenization in silicates. Planetary system applications include more robust enhancement of solids in Jupiter's core formation region than previously suggested. Astrophysical applications include differential, time-dependent enhancement of vapor phase CO and H$_2$O in the terrestrial planet regions of actively accreting protoplanetary disks.

Size-selective concentration of particles in a weakly turbulent protoplanetary nebula may be responsible for the initial collection of chondrules and other constituents into primitive body precursors. This paper presents the main elements of this process of turbulent concentration. In the terrestrial planet region, both the characteristic size and size distribution of chondrules are explained. "Fluffier" particles would be concentrated in nebula regions which were at a lower gas density and/or more intensely turbulent. The spatial distribution of concentrated particle density obeys multifractal scaling}, suggesting a close tie to the turbulent cascade process. This scaling behavior allows predictions of the probability distributions for concentration in the protoplanetary nebula to be made. Large concentration factors (>10^5) are readily obtained, implying that numerous zones of particle density significantly exceeding the gas density could exist. If most of the available solids were actually in chondrule sized particles, the ensuing particle mass density would become so large that the feedback effects on gas turbulence due to mass loading could no longer be neglected. This paper describes the process, presenting its basic elements and some implications, without including the effects of mass loading.