Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

Near-infrared polarimetry of the GG Tauri A binary system

A high angular resolution near-infrared image that shows the intensity of polarization for the GG Tau A binary system was obtained with the Subaru Telescope.The image shows a circumbinary disk scattering the light from the central binary. The azimuthal profile of the intensity of polarization for the circumbinary disk is roughly reproduced by a simple disk model with the Henyey-Greenstein phase function and the Rayleigh function, indicating there are small dust grains at the surface of the disk.Combined with a previous observation of the circumbinary disk, our image indicates that the gap structure in the circumbinary disk orbits counterclockwise, but material in the disk orbits clockwise. We propose that there is a shadow caused by material located between the central binary and the circumbinary disk. The separations and position angles of the stellar components of the binary in the past 20 yr are consistent with the binary orbit with a = 33.4 AU and e = 0.34.     

Dust-induced chemical differentiation in dense regions

Dust grains respond to the physical and chemical conditions of the interstellar region in which they are embedded. The interaction produces an extinction curve which depends on the local environment and on the past history of the dust grains. In this work we present a theoretical study of the effects of radial variations of dust extinction properties on gas-phase chemistry in spherical core–halo clouds. We use observational constraints on the variation range of the extinction curve, and we analyse if the degree of dust environmental processing could be reflected by chemical signatures in the gas-phase molecular concentrations. The results of this work show that significant variations in the photodestruction rates and in the thermal profile of the cloud might induce chemical patterns otherwise excluded in the standard dense-cloud chemistry. Some discrepancies between observations and theoretical provisions are discussed in the light of the present results.     

The local oxygen abundance evolution in models with galactic dust expulsion

The abundance evolution of oxygen in the local galactic disk is discussed. The age-metallicity relation of nearby stars is confronted with predictions from simple evolution models where infall of primordial gas and outflow of disk gas and dust are allowed. Gas infall and expulsion of dust grains have considerable effects on the age-metallicity relation. Dust outflow may be important in order to solve the discrepancy between the observed rate of infall of primordial gas and rates predicted by simple evolution models.     

Self-magnetization of the early Solar System matter

The remnant magnetization of chondrite type meteorite matter indicates the existence of 10⁻⁵-10⁻³ T magnetic fields in the early Solar System accretion disk. Taking into account parameter regimes being typical for this evolutionary stage of Sun and planets we consider the protosolar disk matter as partially ionized dusty plasma consisting of massive charged dust grains, neutral gas, electrons and ions. Results of systematic multifluid neutral gas-plasma-dust simulations show that shear flow driven collisional interactions yield a self-magnetization of the early Solar System matter which is able to explain the measured remnant magnetization of meteorite material.     

Sweeping-out of dust from spiral galaxy as a result of “magnetic sling” effect

A new mechanism of sweeping out of dust grains beyond galactic disks both in the radial direction along the galactic plane and in the vertical, cross-disk direction is proposed. The mechanism is driven by the interaction of dust grains with the bisymmetric nonstationary magnetic field of the galaxy, whose lines are curved and corotate with the stellar spiral density wave responsible for the arms. We attribute the radial transfer of interstellar dust grains in the plane of galactic disks to the fact that charged dust grains are “glued” to magnetic field lines and are therefore pushed outward because of the rotation of magnetic field lines and their tilt with respect to the radial direction parallel to the disk plane. In addition, dust is swept out vertically in the cross-disk direction because of the drift motion in crossed magnetic and gravitational fields (both are parallel to the galactic plane). Numerical computations of the motion of dust grains in real magneto-gravitational fields with the allowance for the drag force from interstellar gas show that the time scale of dust grain transport beyond galactic disks is on the order of 1 Gyr or shorter.     

Dust in the disk winds from young stars as a source of the circumstellar extinction

We consider the problem of dust grain survival in the disk winds from T Tauri and Herbig Ae stars. For our analysis, we have chosen a disk wind model in which the gas component of the wind is heated through ambipolar diffusion to a temperature of ～10⁴ K. We show that the heating of dust grains through their collisions with gas atoms is inefficient compared to their heating by stellar radiation and, hence, the grains survive even in the hot wind component. As a result, the disk wind can be opaque to the ultraviolet and optical stellar radiation and is capable of absorbing an appreciable fraction of it. Calculations show that the fraction of the wind-absorbed radiation for T Tauri stars can be from 20 to 40% of the total stellar luminosity at an accretion rate ? _a = 10^?8-10^?6 M _⊙yr^?1. This means that the disk winds from T Tauri stars can play the same role as the puffed-up inner rim in current accretion disk models. In Herbig Ae stars, the inner layers of the disk wind (r ≤ 0.5 AU) are dust-free, since the dust in this region sublimates under the effect of stellar radiation. Therefore, the fraction of the radiation absorbed by the disk wind in this case is considerably smaller and can be comparable to the effect from the puffed-up inner rim only at an accretion rate of the order of or higher than 10^?6 M _⊙yr^?1. Since the disk wind is structurally inhomogeneous, its optical depth toward the observer can be variable, which should be reflected in the photometric activity of young stars. For the same reason, moving shadows from gas and dust streams with a spiral-like shape can be observed in high-angular-resolution circumstellar disk images.     

Astronomical evidence for the rapid growth of millimeter‐sized particles in protoplanetary disks

Abstract– I summarize recent surveys of protoplanetary disks at millimeter wavelengths and show that the distribution of luminosity, equivalent to the mass in small dust grains, declines rapidly. This contrasts with statistics on the lifetime of disks from infrared observations and the high occurrence of planets from radial velocity and transit surveys. I suggest that these disparate results can be reconciled if most of the dust in a disk is locked up in millimeter and larger‐sized particles within about 2 Myr. This statistical result on disk evolution agrees with detailed modeling of a small number of individual disks and with cosmochemical measurements of rapid planetesimal formation.     

Chemical evolution of the gas in C-type shocks in dark clouds

A magnetohydrodynamic model of a steady, transverse C-type shock in a dense molecular cloud is presented. A complete gas–grain chemical network is taken into account: the gas-phase chemistry, the adsorption of gas species on dust grains, various desorption mechanisms, the grain surface chemistry, the ion neutralization on dust grains, the sputtering of grain mantles. The population densities of energy levels of ions CI, CII and OI and molecules H₂, CO, H₂O are computed in parallel with the dynamical and chemical rate equations. The large velocity gradient approximation is used in the line radiative transfer calculations. The simulations consist of two steps: (i) modelling of the chemical and thermal evolution of a static molecular cloud and (ii) shock simulations. A comparison is made with the results of publicly available models of similar physical systems.The focus of the paper is on the chemical processing of gas material and ice mantles of dust grains by the shock. Sputtering of ice mantles takes place in the shock region close to the temperature peak of the neutral gas. At high shock speeds, molecules ejected from ice mantles are effectively destroyed in hot gas, and their survival time is low—of the order of dozens of years. After a passage of high-speed C-type shock, a zone of high abundance of atomic hydrogen appears in the cooling postshock gas that triggers formation of complex organic species such as methanol. It is shown that abundances of some complex organic molecules (COMs) in the postshock region can be much higher than in the preshock gas. These results are important for interpretation of observations of COMs in protostellar outflows.     

On the diversity of giant planets—Simulating the evolution of solids in protoplanetary disks

We describe a model designed to track simultaneously the evolution of gas and solids in protoplanetary disks from an early stage, when all solids are in the dust form, to the stage when most solids are in the form of a planetesimal swarm. The model is computationally efficient and allows for a global, comprehensive approach to the evolution of solid particles due to gas–solid coupling, coagulation, sedimentation, and evaporation/condensation. We have used it to calculate the co-evolution of gas and solids starting from a comprehensive domain of initial conditions. Then based on the core accretion-gas capture scenario, we have estimated the planet-bearing capability of the environment defined by the final planetesimal swarm and the still evolving gaseous component of the disk. We describe how the disk's capability of formation of giant planets depends on the initial mass and size of a protoplanetary disk, its thermal structure, mass of the central star and properties of the material forming solid grains.     

Infrared photometry of the X-ray binary XTE J1118+480 in April 2000

We present the IR photometry of the X-ray binary XTE J1118+480 performed during seven nights in April and two nights in May–June 2000. A significant IR excess has been detected in the object, which may be due to the thermal radiation from a dust envelope/cloud. The observed energy distribution in the range 1.25–3.5 μm can be interpreted in terms of the sum of the fluxes from an accretion disk with a temperature of ～20 000 K and a dust envelope with grains heated to ～900 K. The distance to the X-ray binary estimated from the total flux from the dust envelope is no less than 0.6–3 kpc. The mean optical depth of the dust envelope for the accretion-disk radiation is about 0.06.     

The thermal structure and the location of the snow line in the protosolar nebula: Axisymmetric models with full 3-D radiative transfer

The precise location of the water ice condensation front (‘snow line’) in the protosolar nebula has been a debate for a long time. Its importance stems from the expected substantial jump in the abundance of solids beyond the snow line, which is conducive to planet formation, and from the higher ‘stickiness’ in collisions of ice-coated dust grains, which may help the process of coagulation of dust and the formation of planetesimals. In an optically thin nebula, the location of the snow line is easily calculated to be around 3 AU, subject to brightness variations of the young Sun. However, in its first 5-10 myr, the solar nebula was optically thick, implying a smaller snowline radius due to shielding from direct sunlight, but also a larger radius because of viscous heating. Several models have attempted to treat these opposing effects. However, until recently treatments beyond an approximate 1 + 1D radiative transfer were unfeasible. We revisit the problem with a fully self-consistent 3D treatment in an axisymmetric disk model, including a density-dependent treatment of the dust and ice sublimation. We find that the location of the snow line is very sensitive to the opacities of the dust grains and the mass accretion rate of the disk. We show that previous approximate treatments are quite efficient at determining the location of the snow line if the energy budget is locally dominated by viscous accretion. Using this result we derive an analytic estimate of the location of the snow line that compares very well with results from this and previous studies. Using solar abundances of the elements we compute the abundance of dust and ice and find that the expected jump in solid surface density at the snow line is smaller than previously assumed. We further show that in the inner few AU the refractory species are also partly evaporated, leading to a significantly smaller solid state surface density in the regions where the rocky planets were formed.     

Hydrodynamic Aspects of Modeling of the Mass Transfer and Coagulation Processes in Turbulent Accretion Disks

This paper considers, in the context of modeling the evolution of a protoplanetary cloud, the hydrodynamic aspects of the theory of concurrent processes of mass transfer and coagulation in a two-phase medium in the presence of shear turbulence in a differentially rotating gas–dust disk and of polydisperse solid particles suspended in a carrying flow of solid particles. The defining relations are derived for diffuse fluxes of particles of different sizes in the equations of turbulent diffusion in the gravitational field, which describe the convective transfer, turbulent mixing, and sedimentation of disperse dust grains onto the central plane of the disk, as well as their coagulation growth. A semiempirical method is developed for calculating the coefficients of turbulent viscosity and turbulent diffusion for particles of different kinds. This method takes into account the inverse effects of dust transfer on the turbulence evolution in the disk and the inertial differences between disperse solid particles. To solve rigorously the problem of the mutual influence of the turbulent mixing and coagulation kinetics in forming the gas–dust subdisk, the possible mechanisms of gravitational, turbulent, and electric coagulation in a protoplanetary disk are explored and the parametric method of moments for solving the Smoluchowski integro-differential coagulation equation for the particles' size distribution function is considered. This method takes into account the fact that this distribution belongs to a definite parametric class of distributions.     

Radiation Pressure and the Poynting-Robertson Effect for Fluffy Dust Particles

A correct understanding of the dynamical effect of solar radiation exerted on fluffy dust particles can be achieved with assistance of a light scattering theory as well as the equation of motion. We reformulate the equation of motion so that the radiation pressure and the Poynting-Robertson effect on fluffy grains are given in both radial and nonradial directions from the center of the Sun. This allows numerical estimates of these radiation forces on fluffy dust aggregates in the framework of the discrete dipole approximation, in which the first term of the scattering coefficients in Mie theory determines the polarizability of homogeneous spheres forming the aggregates.The nonsphericity in shape turns out to play a key role in the dynamical evolution of dust particles, while its consequence depends on the rotation rate and axis of the grains. Unless a fluffy dust particle rapidly revolves on its randomly oriented axis, the nonradial radiation forces may prevent, apart from the orbital eccentricity and semimajor axis, the orbital inclination of the particle from being preserved in orbit around the Sun. However, a change in the inclination is most probably controlled by the Lorentz force as a consequence of the interaction between electric charges on the grains and the solar magnetic field. Although rapidly and randomly rotating grains spiral into the Sun under the Poynting-Robertson effect in spite of their shapes and structures, fluffy grains drift inward on time scales longer at submicrometer sizes and shorter at much larger sizes than spherical grains of the same sizes. Numerical calculations reveal that the dynamical lifetimes of fluffy particles are determined by the material composition of the grains rather than by their morphological structures and sizes. The Poynting-Robertson effect alone is nevertheless insufficient for giving a satisfactory estimate of lifetimes for fluffy dust grains since their large ratios of cross section to mass would reduce the lifetimes by enhancing the collisional probabilities. We also show that the radiation pressure on a dust particle varies with the orbital velocity of the particle but that this effect is negligibly small for dust grains in the Solar System.     

A Study on the Characteristics of the Structure of Vega's Debris Disk

The clumpy structure in the Vega's debris disk was reported at millimeter wavelengths previously, and attributed to the concentration of dust grains trapped in resonances with a potential high-eccentricity planet. However, current imaging at multi-wavelengths with higher sensitivity indicates that the Vega's debris disk has a smooth structure. But a planet orbiting Vega could not be neglected, and the present-day observations may place a severe constraint on the orbital parameters for the potential planet. Herein, we utilize the modi- fied MERCURY codes to numerically simulate the Vega system, which consists of a debris disk and a planet. In our simulations, the initial inner and outer boundaries of the debris disk are assumed to be 80 AU and 120 AU, respectively. The dust grains in the disk have the sizes from 10 μm to 100 μm, and the nearly coplanar orbits. From the outcomes, we show that the evolution of debris disk is consistent with recent observations, if there is no planet orbiting Vega. However, if Vega owns a planet with a high eccentricity (e.g., e = 0.6), the planet's semi- major axis cannot be larger than 60 AU, otherwise, an aggregation phenomenon will occur in the debris disk due to the existence of the postulated planet. In addition, the 2:1 mean motion resonances may play a significant role in forming the structure of debris disk.     

Salten Skov I: A Martian magnetic dust analogue

A fine grained magnetic iron oxide precipitate found in Denmark has been studied with regard to grain size, magnetic properties, aerosol transport, grain electrification, aggregation and optical reflectance. It has shown itself to be a good Martian dust analogue. The fraction of the Salten Skov I soil sample <63 μm was separated from the natural sample by dry sieving. This fraction could be dispersed by ultrasonic treatment into grains of diameter ～1 μm, in reasonable agreement with suspended dust grains in the Martian atmosphere estimated from the Viking, Pathfinder and Mars Exploration Rover missions. Though mineralogical and chemical differences exist between this analogue and Martian dust material, in wind tunnel experiments many of the physical properties of the atmospheric dust aerosol are reproduced.     

The effect of grain size distribution on H2 formation rate in the interstellar medium

The formation of molecular hydrogen  (H₂)  in the interstellar medium takes place on the surfaces of dust grains. Hydrogen molecules play a role in gas-phase reactions that produce other molecules, some of which serve as coolants during gravitational collapse and star formation. Thus, the evaluation of the production rate of hydrogen molecules and its dependence on the physical conditions in the cloud are of great importance. Interstellar dust grains exhibit a broad size distribution in which the small grains capture most of the surface area. Recent studies have shown that the production efficiency strongly depends on the grain composition and temperature as well as on its size. In this paper, we present a formula that provides the total production rate of  H₂  per unit volume in the cloud, taking into account the grain composition and temperature as well as the grain size distribution. The formula agrees very well with the master equation results. It shows that for a physically relevant range of grain temperatures, the production rate of  H₂  is significantly enhanced due to their broad size distribution.     

The role of sticky interstellar organic material in the formation of asteroids

Abstract— Collision experiments and measurements of viscoelastic properties were performed involving an interstellar organic material analogue to investigate the growth of organic grains in the protosolar nebula. The organic material was found to be stickiest at a radius of between 2.3 and 3.0 AU, with a maximum sticking velocity of 5 m s^?1 for millimeter‐size organic grains. This stickiness is considered to have resulted in the very rapid coagulation of organic grain aggregates and subsequent formation of planetesimals in the early stage of the turbulent accretion disk. The planetesimals formed in this region appear to be represent achondrite parent bodies. In contrast, the formation of planetesimals at <2.1 and >3.0 AU begins with the establishment of a passive disk because silicate and ice grains are not as sticky as organic grains.     

Scattering Properties and Composition of Cometary Dust

Composition of the Comet dust obtained by the dust impact analyzer on the Halley probes indicated that the comet dust is a mixture of silicate and carbonaceous material. The collected interplanetary dust particles (IDP's) are fluffy and composite, having grains of several different types stuck together. Using discrete dipole approximation (DDA) we study the scattering properties of composite grains. In particular, we study the angular distribution of the scattered intensity and linear polarization of composite grains. We assume that the composite grains are made up of a host silicate sphere/spheroid with the inclusions of graphite. Results of our calculations on the composite grains show that the angle of maximum polarization shifts, and the degree of polarization varies with the volume fraction of the inclusions. We use these results on the composite grains to interpret the observed scattering in cometary dust.     

Surface thermoelastic erosion of atmosphereless solar system bodies under bombardment by multicharge cosmic ray ions

The mechanism of ion-stimulated erosion of atmosphereless solar system bodies is suggested and investigated. A theoretical model for the brittle surface erosion resulting under the effect of multicharge ion cosmic rays is analyzed. It is shown that the thermoelastic waves originated in the energetic track of a very heavy ion can result in the near-surface stresses exceeding the dynamic tensile strength of the surface material for any atmosphereless solar system body. The thermoelastic wave surface arrival yields brittle erosion of the material and ejection of this latter fragments (the track-breaking process). Thus ejected dust grains have plano-oblong shape, average mass on the order of 10^–17 g and velocity up to 400 m/sec providing the surface erosion rate of 10^–1 ÷ 3 · 10² »/year (near the Earth orbit) which depends upon the surface material (rock or ice). Possible track-breaking consequences, in particular, presence of the dust fraction of ultramicron grains and their aggregates on the lunar surface are discussed. Near the bodies with the radii from 10 to 300 km predicted is the existence of extended dust cocoons consisting of ultramicron and submicron grains. Smaller objects (asteroids, comets, smallest satellites of planets, meteoroids, etc.) can serve sources of permanent dust wind of ultramicron and submicron sized grains escaping from their surfaces. The interplanetary dust yield owing to the ion-stimulated erosion of these bodies is not less than 10¹² g/year. Possible interpreting in the frames of track-breaking process some observational data and effects, including existence of dust grains with the mass of 10^–18 ÷ 10^–17 g near the Halley's comet and the nature of 2060 Chiron dust coma is discussed. To prove the theory, observational identification and investigation of dust phenomena complex related to the ion-stimulated erosion of atmosphereless bodies, suggested is employing extreme ultraviolet and far infrared/submillimeter wavelengths, as well as polarimetric methods.