Coagulation of dust aggregates plays an important role in the formation of planets and is of key importance to the evolution of protoplanetary disks (PPDs). Characteristics of dust, such as the diversity of particle size, porosity, charge, and the manner in which dust couples to turbulent gas, affect the collision outcome and the rate of dust growth. Here we present a numerical model of the evolution of the dust population within a PPD which incorporates all of these effects. The probability that any two particles collide depends on the particle charge, cross-sectional area and their relative velocity. The actual collision outcome is determined by a detailed collision model which takes into account the aggregate morphology, trajectory, orientation, and electrostatic forces acting between charged grains. The data obtained in this research reveal the characteristics of dust populations in different environments at the end of the hit-and-stick growth, which establishes the foundation for the onset of the next growth stage where bouncing, mass transfer and fragmentation become important. For a given level of turbulence, neutral and weakly charged particles collide more frequently and grow faster than highly charged particles. However, highly charged particles grow to a larger size before reaching the bouncing barrier, and exhibit a "Runaway" growth, in which a few large particles grow quickly by accreting smaller particles while the rest of the population grows very slowly. In general, highly charged aggregates have a more compact structure and are comprised of larger monomers than neutral/weakly charged aggregates. The differences in the particle structure/composition not only affect the threshold velocities for bouncing and fragmentation, but also change the scattering and absorption opacity of dust, influencing the appearance of PPDs.

In order to characterize the early growth of fine-grained dust rims (FGRs) that commonly surround chondrules, we simulate the growth of FGRs through direct accretion of monomers of various sizes onto the chondrule surfaces. Dust becomes charged to varying degrees in the radiative plasma environment of the solar nebula (SN), and the resulting electrostatic force alters the trajectories of colliding dust grains, influencing the structure of the dust rim as well as the time scale of rim formation. We compare the growth of FGRs in protoplanetary disks (PPD) with different turbulence strengths and plasma conditions to previous models which assumed neutral dust grains. We use a combination of a Monte Carlo method and an N-body code to simulate the collision of dust monomers with chondrules: a Monte Carlo algorithm is used to randomly select dust particles that will collide with the chondrule as well as determine the elapsed time interval between collisions; at close approach, the detailed collision process is modeled using an N-body algorithm, Aggregate Builder (AB), to determine the collision outcome, as well as any restructuring of the chondrule rim. The collisions are driven by Brownian motion and coupling to turbulent gas motion in the protoplanetary disk. The charge distribution of the dust rim is modeled, used to calculate the trajectories of dust grains, and then analyze the resulting morphology of the dust rim. In a weakly turbulent region, the decreased relative velocity between charged particles causes small grains to be repelled from the chondrule, causing dust rims to grow more slowly and be composed of larger monomers, which results in a more porous structure. In a highly turbulent region, the presence of charge mainly affects the porosity of the rim by causing dust particles to deviate from the extremities of the rim and reducing the amount of restructuring.

In order to characterize the early growth of fine-grained dust rims (FGRs) that commonly surround chondrules, we perform numerical simulations of dust accretion onto chondrule surfaces. We employ a Monte Carlo algorithm to simulate the collision of dust monomers having radii between 0.5 and 10 $\mu$m with chondrules whose radii are between 500 and 1000 $\mu$m, in 100-$\mu$m increments. The collisions are driven by Brownian motion and solar nebula turbulence. After each collision, the colliding particles either stick at the point of contact, roll or bounce. We limit accretion of dust monomers (and in some cases, dust aggregates) to a small patch of the chondrule surface, for computational expediency. We model the morphology of the dust rim and the trajectory of the dust particle, which are not considered in most of the previous works. Radial profiles of FGR porosity show that rims formed in weak turbulence are more porous (with a porosity of 60-74\%) than rims formed in stronger turbulence (with a porosity of 52-60\%). The lower end of each range corresponds to large chondrules and the upper end to small chondrules, meaning that the chondrule size also has an impact on FGR porosity. The thickness of FGRs depends linearly on chondrule radius, and the slope of this linear dependency increases with time, and decreases with the turbulence strength. The porosity of FGRs formed by dust aggregates is $\sim 20\%$ on average greater than that of FGRs formed by single monomers. In general, the relatively high porosities that we obtain are consistent with those calculated by previous authors from numerical simulations, as well as with initial FGR porosities inferred from laboratory measurements of rimmed chondrule samples and rimmed chondrule analogs.

(Abridged) We perform the first study of magnetohydrodynamic processes in the protolunar disk (PLD). With the use of published data on the chemical composition of the PLD, along with existing analytical models of the disk structure, we show that the high temperatures that were prevalent in the disk would have led to ionization of Na, K, SiO, Zn and, to a lesser extent, O$_2$. We assume that the disk has a vapor structure. The resulting ionization fractions, together with a relatively weak magnetic field, would have been sufficient to trigger the magneto-rotational instability, or MRI. We calculate the intensity of the resulting magnetohydrodynamic turbulence, as parameterized by the dimensionless ratio $\alpha$ of turbulent stresses to gas pressure, and obtain maximum values $\alpha\sim 10^{-2}$ along most of the vertical extent of the disk, and at different orbital radii. Under these conditions, turbulent mixing within the PLD due to the MRI was likely capable of transporting isotopic and chemical species efficiently. To test these results in a conservative manner, we carry out a numerical magnetohydrodynamic simulation of a small, rectangular patch of the PLD, located at 4 Earth radii ($r_{\rm E}$) from the center of the Earth, and assuming once again that the disk is completely gaseous. We obtain turbulence with an average intensity $\alpha\sim 7\times 10^{-6}$ over the course of 280 orbital periods (133 days at 4$r_{\rm{E}}$). Despite this relatively low value of $\alpha$, the effective turbulent diffusivity of a passive tracer introduced in the flow is large enough to allow the tracer to spread across a radial distance of $10r_{\rm{E}}$ in $\sim13-129$ yr, less than the estimated cooling time of the PLD of $\sim250$ yr. Further improvements to our model will need to incorporate the energy balance in the disk, a complete two-phase structure, and a more realistic equation of state.

We analyze the coagulation of dust in and around a gap opened by a Jupiter-mass planet. To this end, we carry out a high-resolution magnetohydrodynamic (MHD) simulation of the gap environment, which is turbulent due to the magnetorotational instability. From the MHD simulation, we obtain values of the gas velocities, densities and turbulent stresses a) close to the gap edge, b) in one of the two gas streams that accrete onto the planet, c) inside the low-density gap, and d) outside the gap. The MHD values are then supplied to a Monte Carlo dust coagulation algorithm, which models grain sticking and compaction. We consider two dust populations for each region: one whose initial size distribution is monodisperse, with monomer radius equal to 1 $\mu$m, and another one whose initial size distribution follows the Mathis-Rumpl-Nordsieck distribution for interstellar dust grains, with an initial range of monomer radii between 0.5 and 10 $\mu$m. Our Monte Carlo calculations show initial growth of dust aggregates followed by compaction in all cases but one, that of aggregates belonging to the initially monodisperse population subject to gas conditions outside the gap. In this latter case, the mass-weighted (MW) average porosity of the population reaches extremely high final values of 98\%. The final MW porosities in all other cases range between 30\% and 82\%. The efficiency of compaction is due to high turbulent relative speeds between dust particles. Future studies will need to explore the effect of different planet masses and electric charge on grains.

(Abridged) Numerical magnetohydrodynamic (MHD) simulations of a turbulent solar nebula are used to study the growth of dust mantles swept up by chondrules. A small neighborhood of the solar nebula is represented by an orbiting patch of gas at a radius of 3 AU, and includes vertical stratification of the gas density. The differential rotation of the nebular gas is replaced by a shear flow. Turbulence is driven by destabilization of the flow as a result of the magnetorotational instability (MRI), whereby magnetic field lines anchored to the gas are continuously stretched by the shearing motion. A passive contaminant mimics small dust grains that are aerodynamically well coupled to the gas, and chondrules are modeled by Lagrangian particles that interact with the gas through drag. Whenever a chondrule enters a region permeated by dust, its radius grows at a rate that depends on the local dust density and the relative velocity between itself and the dust. The local dust abundance decreases accordingly. Different chondrule volume densities lead to varying depletion and rimmed-chondrule size growth times. Most of the dust sweep-up occurs within 1 gas scale height of the nebula midplane. Chondrules can reach their asymptotic radius in 10 to 800 years. The vertical variation of nebula turbulent intensity results in a moderate dependence of mean rimmed-chondrule size with nebula height. The technique used here could be combined with Monte Carlo (MC) methods that include the physics of dust compaction, in a self-consistent MHD-MC model of dust rim growth around chondrules in the solar nebula.

The central engine of Gamma Ray Bursts is hidden from direct probing with photons mainly due to the high densities involved. Inferences on their properties are thus made from their cosmological setting, energetics, low-energy counterparts and variability. If GRBs are powered by hypercritical accretion onto compact objects, on small spatial scales the flow will exhibit fluctuations, which could in principle be reflected in the power output of the central engine and ultimately in the high energy prompt emission. Here we address this issue by characterizing the variability in neutrino cooled accretion flows through local shearing box simulations with magnetic fields, and then convolving them on a global scale with large scale dynamical simulations of accretion disks. The resulting signature is characteristic, and sensitive to the details of the cooling mechanism, providing in principle a discriminant for GRB central engine properties.

We apply ionization balance and MHD calculations to investigate whether magnetic activity moderated by recombination on dust can account for the mass accretion rates and the mid-infrared spectra and variability of protostellar disks. The MHD calculations use the stratified shearing-box approach and include grain settling and the feedback from the changing dust abundance on the resistivity of the gas. The two-decade spread in accretion rates among T Tauri stars is too large to result solely from variety in the grain size and stellar X-ray luminosity, but can be produced by varying these together with the disk magnetic flux. The diversity in the silicate bands can come from the coupling of grain settling to the distribution of the magneto-rotational turbulence, through three effects: (1) Recombination on grains yields a magnetically inactive dead zone extending above two scale heights, while turbulence in the magnetically active disk atmosphere overshoots the dead zone boundary by only about one scale height. (2) Grains deep in the dead zone oscillate vertically in waves driven by the turbulent layer above, but on average settle at the laminar rates, so the interior of the dead zone is a particle sink and the disk atmosphere becomes dust-depleted. (3) With sufficient depletion, the dead zone is thinner and mixing dredges grains off the midplane. The MHD results also show that the magnetic activity intermittently lifts clouds of dust into the atmosphere. The photosphere height changes by up to one-third over a few orbits, while the extinction along lines of sight grazing the disk surface varies by factors of two over times down to 0.1 orbit. We suggest that the changing shadows cast by the dust clouds on the outer disk are a cause of the daily to monthly mid-infrared variability in some young stars. (Abridged.)

We perform numerical simulations of solid particle motion in a shearing box model of a protoplanetary disc. The accretion flow is turbulent due to the action of the magnetorotational instability. Aerodynamic drag on the particles is modelled using the Epstein law with the gas velocity interpolated to the particle position. The effect of the magnetohydrodynamic turbulence on particle velocity dispersions is quantified for solids of different stopping times t_s, or equivalently, different sizes. The anisotropy of the turbulence is reflected upon the dispersions of the particle velocity components, with the radial component larger than both the azimuthal and vertical components for particles larger than ~ 10 cm (assuming minimum-mass solar nebula conditions at 5 AU). The dispersion of the particle velocity magnitude, as well as that of the radial and azimuthal components, as functions of stopping time, agree with previous analytical results for isotropic turbulence. The relative speed between pairs of particles with the same value of t_s decays faster with decreasing separation than in the case of solids with different stopping time. Correlations in the particle number density introduce a non-uniform spatial distribution of solids in the 10 to 100 cm size range. Any clump of particles is disrupted by the turbulence in less than one tenth on an orbital period, and the maximally concentrated clumps are stable against self-gravitational collapse.

We study the vertical settling of solid bodies in a turbulent protoplanetary disc. We consider the situation when the coupling to the gas is weak or equivalently when the particle stopping time tau_{st} due to friction with the gas is long compared to the orbital timescale Omega^{-1}. An analytical model, which takes into account the stochastic nature of the sedimentation process using a Fokker-Planck equation for the particle distribution function in phase space, is used to obtain the vertical scale height of the solid layer as a function of the vertical component of the turbulent gas velocity correlation function and the particle stopping time. This is found to be of the same form as the relation obtained for strongly coupled particles in previous work. We compare the predictions of this model with results obtained from local shearing box MHD simulations of solid particles embedded in a vertically stratified disc in which there is turbulence driven by the MRI. We find that the ratio of the dust disc thickness to the gas disc thickness satifies H_d/H=0.08 (Omega tau_{st})^{-1/2}, which is in very good agreement with the analytical model. By discussing the conditions for gravitational instability in the outer regions of protoplanetary discs in which there is a similar level of turbulence, we find that bodies in the size range 50 to 600 metres can aggregate to form Kuiper belt-like objects with characteristic radii ranging from tens to hundreds of kilometres.

We calculate the radial diffusion coefficient for a passive contaminant in an accretion disc which is turbulent due to the action of the magnetorotational instability. Numerical MHD simulations are used to follow the evolution of a local patch of the disc using the shearing box formalism. A separate continuity equation for the mass fraction of contaminant is integrated along with the MHD system, and radial profiles of this fraction are obtained as a function of time. Solutions of a linear diffusion equation are fitted to the numerical measured profiles of the contaminant, treating the diffusion coefficient D as the fitting parameter. At early times, the value of D is found to vary, however once the contaminant is spread over scales comparable to the box size, it saturates at a steady value. The ratio of D to the transport coefficient of angular momentum due to shear stress is small. If D can be used as a proxy for the turbulent magnetic diffusivity, the effective magnetic Prandtl number P_eff=\nu/D (where \nu is the coefficient of "efective viscosity" associated with shear stress) would be large.