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.

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.

The coagulation of cosmic dust grains is a fundamental process which takes place in astrophysical environments, such as presolar nebulae and circumstellar and protoplanetary disks. Cosmic dust grains can become charged through interaction with their plasma environment or other processes, and the resultant electrostatic force between dust grains can strongly affect their coagulation rate. Since ions and electrons are collected on the surface of the dust grain at random time intervals, the electrical charge of a dust grain experiences stochastic fluctuations. In this study, a set of stochastic differential equations is developed to model these fluctuations over the surface of an irregularly-shaped aggregate. Then, employing the data produced, the influence of the charge fluctuations on the coagulation process and the physical characteristics of the aggregates formed is examined. It is shown that dust with small charges (due to the small size of the dust grains or a tenuous plasma environment) are affected most strongly.

The charging of dust grains in astrophysical environments has been investigated with the assumption these grains are homogeneous spheres. However, there is evidence which suggests many grains in astrophysical environments are irregularly-shaped aggregates. Recent studies have shown that aggregates acquire higher charge-to-mass ratios due to their complex structures, which in turn may alter their subsequent dynamics and evolution. In this paper, the charging of aggregates is examined including secondary electron emission and photoemission in addition to primary plasma currents. The results show that the equilibrium charge on aggregates can differ markedly from spherical grains with the same mass, but that the charge can be estimated for a given environment based on structural characteristics of the grain. The "small particle effect" due to secondary electron emission is also important for determining the charge of micron-sized aggregates consisting of nano-sized particles.

Titan is one of the more distinctive bodies in our solar system. In addition to being the largest of Saturn's moons, its thick atmosphere gene-rates interest because of its similarities and differences with Earth [1, 2]. Like Earth, Titan's lower atmosphere contains clouds which precipitate as rain [2]. This rain forms lakes and rivers of liquid methane and ethane which erode and shape the surface, much like water does on Earth, before evaporating into the atmosphere [2]. In Titan's atmospheric system, a single dominating factor controls the weather, the concentration of the organic molecule tholin.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a concept for an 8-meter to 16-meter UVOIR space observatory for launch in the 2025-2030 era. ATLAST will allow astronomers to answer fundamental questions at the forefront of modern astronphysics, including "Is there life elsewhere in the Galaxy?" We present a range of science drivers that define the main performance requirements for ATLAST (8 to 16 milliarcsec angular resolution, diffraction limited imaging at 0.5 {\mu}m wavelength, minimum collecting area of 45 square meters, high sensitivity to light wavelengths from 0.1 {\mu}m to 2.4 {\mu}m, high stability in wavefront sensing and control). We will also discuss the synergy between ATLAST and other anticipated future facilities (e.g., TMT, EELT, ALMA) and the priorities for technology development that will enable the construction for a cost that is comparable to current generation observatory-class space missions.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation of UVOIR space observatory with a primary aperture diameter in the 8-m to 16-m range that will allow us to perform some of the most challenging observations to answer some of our most compelling questions, including "Is there life elsewhere in the Galaxy?" We have identified two different telescope architectures, but with similar optical designs, that span the range in viable technologies. The architectures are a telescope with a monolithic primary mirror and two variations of a telescope with a large segmented primary mirror. This approach provides us with several pathways to realizing the mission, which will be narrowed to one as our technology development progresses. The concepts invoke heritage from HST and JWST design, but also take significant departures from these designs to minimize complexity, mass, or both. Our report provides details on the mission concepts, shows the extraordinary scientific progress they would enable, and describes the most important technology development items. These are the mirrors, the detectors, and the high-contrast imaging technologies, whether internal to the observatory, or using an external occulter. Experience with JWST has shown that determined competitors, motivated by the development contracts and flight opportunities of the new observatory, are capable of achieving huge advances in technical and operational performance while keeping construction costs on the same scale as prior great observatories.

The structure and evolution of aggregate grains formed within a plasma environment are dependent upon the charge acquired by the micron-sized dust grains during the coagulation process. The manner in which the charge is arranged on developing irregular structures can affect the fractal dimension of aggregates formed during collisions, which in turn influences the coagulation rate and size evolution of the dust within the plasma cloud. This paper presents preliminary models for the charge and size evolution of fractal aggregates immersed in a plasma environment calculated using a modification to the orbital-motion-limited (OML) theory. Primary electron and ion currents incident on points on the aggregate surface are determined using a line-of-sight (LOS) approximation: only those electron or ion trajectories which are not blocked by another grain within the aggregate contribute to the charging current. Using a self-consistent iterative approach, the equilibrium charge and dipole moment are calculated for the dust aggregate. The charges are then used to develop a heuristic charging scheme which can be implemented in coagulation models. While most coagulation theories assume that it is difficult for like-charged grains to coagulate, the OML_LOS approximation indicates that the electric potentials of aggregate structures are often reduced enough to allow significant coagulation to occur.

The charge on micron-sized dust grains plays a crucial role in the structure and evolution of forming aggregates within the dust population during the coagulation process. The manner in which the charge is arranged on developing irregular structures can affect the fractal dimension of aggregates formed during collisions, which in turn influences the coagulation rate and size evolution of the dust cloud. Preliminary models for the charge evolution on fractal aggregates immersed in a plasma environment calculated using a modification to the orbital-motion-limited (OML) theory are presented in this paper. The model calculates currents to each point on the aggregate surface using a line-of-sight (LOS) approximation: only those electron or ion trajectories which are not blocked by another grain within the aggregate contribute to the charging current. Both the total charge and the dipole moment are calculated for the dust aggregate. While most coagulation theories assume that it is difficult for like-charged grains to coagulate, the OML_LOS approximation indicates that the electric potentials of aggregate structures are often reduced enough to allow significant coagulation.

Dusty plasmas have generated a large amount of interest since the discovery of ordered structure (crystal) formation in experimentally generated complex plasmas in 1994. Dust particles within the plasma sheath generated in these complex plasmas can form vertical chains due to the streaming ion wakefield. For the simplest of these configurations (a two particle chain), the particle closest to the lower electrode will generally remain in the shadow of the particle farthest from the lower electrode. These results in the two particles feeling differing ion drag forces: the top particle is acted on by the ion drag force directed from the plasma to the lower electrode, while the bottom particle is acted upon by the resulting wakefield produced by the interaction of the upper particle with the ion drag force. This dynamic situation currently provides the best known experimental environment for examining the physics behind the ion drag force and its interaction with the plasma sheath. An experimental method for investigating the interaction between pair-particle chains based on modulating the bias on the lower electrode employing a DC bias modulation technique will be presented.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

Planetary formation is an efficient process now thought to take place on a relatively short astronomical time scale. Recent observations have shown that the dust surrounding a protostar emits more efficiently at longer wavelengths as the protoplanetary disk evolves, suggesting that the dust particles are coagulating into fluffy aggregates, "much as dust bunnies form under a bed." One poorly understood problem in this coagulation process is the manner in which micron-sized, charged grains form the fractal aggregate structures now thought to be the precursors of protoplanetary disk evolution. This study examines the characteristics of such fractal aggregates formed by the collision of spherical monomers and aggregates where the charge is distributed over the aggregate structure. The aggregates are free to rotate due to collisions and dipole-dipole electrostatic interactions. Comparisons are made for different precursor size distributions and like-charged, oppositelycharged, and neutral grains.

The levitation height of a dust particle layer within a RF discharge plasma sheath is known to be related to the DC bias, the background pressure, and the Debye length. In this paper, a new experimental technique for measurement of the Debye length is introduced. This technique is based on the relationship between an externally applied DC bias and the particle levitation height and shows that under appropriate conditions, the addition of an externally applied DC bias provides a mechanism for evaluation of the Debye length. When compared with existing techniques, this new method appears to be simpler to implement in some cases.

The formation of three-dimensional (3D) dust clusters within a complex plasma modeled as a spatially confined Yukawa system is simulated using the box_tree code. Similar to unscreened Coulomb clusters, the occurrence of concentric shells with characteristic occupation numbers was observed. Both the occupation numbers and radii were found to depend on the Debye length. Ground and low energy meta-stable states of the shielded 3D Coulomb clusters were determined for 4<N<20. The structure and energy of the clusters in different states was analyzed for various Debye lengths. Structural phase transitions, including inter-shell structural phase transitions and intra-shell structural phase transitions, were observed for varying Debye length and the critical value for transitions calculated.

In semiconductor manufacturing, contamination due to particulates significantly decreases the yield and quality of device fabrication, therefore increasing the cost of production. Dust particle clouds can be found in almost all plasma processing environments including both plasma etching devices and in plasma deposition processes. Dust particles suspended within such plasmas will acquire an electric charge from collisions with free electrons in the plasma. If the ratio of inter-particle potential energy to the average kinetic energy is sufficient, the particles will form either a liquid structure with short range ordering or a crystalline structure with long range ordering. Otherwise, the dust particle system will remain in a gaseous state. Many experiments have been conducted over the past decade on such colloidal plasmas to discover the character of the systems formed, but more work is needed to fully understand these structures. The preponderance of previous experiments used monodisperse spheres to form complex plasma systems.

The vertical as well as horizontal oscillation modes in a thermally excited two-dimensional (2D) dust coulomb cluster were investigated for particle numbers between and using both a box_tree simulation and an analytical method. The horizontal mode spectra is shown to agree with published results while the vertical mode spectra obtained from the box_tree simulation and the analytical method are shown to agree with one another. The maximum frequency of the vertical modes is shown to be the vertical oscillation frequency of the whole system acting as a solid plane and the frequency is shown to decrease as the mode number l decreases. The oscillation patterns of the vertical modes are also investigated. For clusters with large numbers of particles, the high frequency modes are shown to have oscillation patterns similar in shape to Bessel-Fourier functions with various index m and n. For specified values of n, modes with higher m (except for) have lower frequencies and for specified m, modes with higher n have lower frequencies. For low frequency modes, the largest amplitude particle motion is concentrated in a few inner rings with the outer rings remaining almost motionless. This is in contrast to the horizontal modes where the strongest motion of the particles is concentrated in the inner rings at high frequency. For clusters with small numbers of particles (), these modes where the strongest vertical motion is concentrated in the inner rings does not exist resulting in oscillation magnitudes for all the modes appearing in the shape of Bessel-Fourier functions.

It has been observed that binary Kuiper Belt Objects (KBOs) exist contrary to theoretical expectations. Their creation presents problems to most current models. However, the inclusion of a third body (for example, one of the outer planets) may provide the conditions necessary for the formation of these objects. The presence of a third massive body not only helps to clear the primordial Kuiper Belt but can also result in long lived binary Kuiper belt objects. The gravitational interaction between the KBOs and the third body causes one of four effects; scattering into the Oort cloud, collisions with the growing protoplanets, formation of binary pairs, or creation of a single Kuiper belt object. Additionally, the initial location of the progenitors of the Kuiper belt objects also has a significant effect on binary formation.

A self-consistent three-dimensional model for a complex (dusty) plasma is used to study the effects of multiple-sized dust grains in a dust crystal. In addition to the interparticle forces, which interact through a Yukawa potential, the model includes the effects of gravity, the variation of the sheath potential above the powered electrode, and a radial confining potential. Simulations studied various ratios of a mix of 6.5- and 8.9-micron monodisperse particles and compared their correlation functions, electric potential energy of the crystal formations, and the dispersion relations for in-plane and out-of-plane dust lattice wave (DLW) modes for two different sheath thicknesses. In the 7 mm sheath, the particles formed two layers in the vertical direction by size, and acted as a two-layer crystal with weak correlation between the layers. In the 3 mm sheath, the particles formed an essentially monolayer crystal; however the crystal dynamics showed some characteristics of a bilayer crystal.