Dust ring formation due to sublimation of dust grains drifting radially inward by the Poynting-Robertson drag: An analytical model

Dust particles exposed to the stellar radiation and wind drift radially inward by the Poynting-Robertson (P-R) drag and pile up at the zone where they begin to sublime substantially. The reason they pile up or form a ring is that their inward drifts due to the P-R drag are suppressed by stellar radiation pressure when the ratio of radiation pressure to stellar gravity on them increases during their sublimation phases. We present analytic solutions to the orbital and mass evolution of such subliming dust particles, and find their drift velocities at the pileup zone are almost independent of their initial semimajor axes and masses. We derive analytically an enhancement factor of the number density of the particles at the outer edge of the sublimation zone from the solutions. We show that the formula of the enhancement factor reproduces well numerical simulations in the previous studies. The enhancement factor for spherical dust particles of silicate and carbon extends from 3 to more than 20 at stellar luminosities L_?=0.8-500L_⊙, where L_⊙ is solar luminosity. Although the enhancement factor for fluffy dust particles is smaller than that for spherical particles, sublimating particles inevitably form a dust ring as long as their masses decrease faster than their surface areas during sublimation. The formulation is applicable to dust ring formation for arbitrary shape and material of dust in dust-debris disks as well as in the Solar System.     

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.     

On the uncertainty of the transmission function of the optically thick AGB dust shells

The uncertainty of the spectral transmission function due to the nonsphericity of cosmic dust particles is analysed for optically thick C-rich and O-rich dust shells. The transmission function directly prescribes the intensity of radiation that passes through a dust shell. It is shown that nonspherical particles affect the stellar spectra in a different way than a system of equally sized spherical particles. Discrepancies in the stellar spectra for both morphological models (spheres and irregular targets) grow with optical thickness of the dust shell and the spectral behaviour of the optical thickness depends on the particle size distribution function. In particular, two most frequently used distrubutions, a power law and Dirac’s delta function, are considered in the presented numerical runs. Light transmission through C-rich dust shells is strongly influenced by absorption which dominates in case of carbonaceous particles. Irregularly shaped carbonaceous particles reduce the intensity of stellar spectra more efficiently than spherical particles of the same composition (the difference is about 10–30%). MgO particles which may be present in O-rich dust shells are almost pure scatterers, thus the ratio of transmission functions for irregularly shaped particles and spherical ones show specific resonant features (especially at wavelengths below 0.4 μm).     

Radiative acceleration in Schwarzschild space–times

We examine the radial motion of a material particle in the intense radiation field of a static spherically symmetric compact object with spherical emitting surface outside the Schwarzschild radius. This paper generalizes previous work which dealt with radial motion in the Thomson limit, where the radiation force is simply proportional to the radiative flux. In the general case the average time component of the 4-momentum transferred to the particle is not negligible compared with its rest mass. Consequently, we find that the frequency dependence of the radiation force owing to Compton scattering for highly energetic photons gives rise to an increase in the effective mass of the test particle. In this work we outline the effects of this frequency dependence and compare these with the results in the Thomson limit. We present the frequency dependent saturation velocity curves for a range of stellar luminosities and radiation frequencies and present the resulting phase-space diagrams corresponding to the radial test particle trajectories. In particular, the stable equilibrium points which exist in the Thomson limit are found to be absent in the general case.     

Orbital Evolution of Dust Particles in the Sublimation Zone near the Sun

We have performed the calculations of the orbital evolution of dust particles from volcanic glass (p-obsidian), basalt, astrosilicate, olivine, and pyroxene in the sublimation zone near the Sun. The sublimation (evaporation) rate is determined by the temperature of dust particles depending on their radius, material, and distance to the Sun. All practically important parameters that characterize the interaction of spherical dust particles with the radiation are calculated using the Mie theory. The influence of radiation and solar wind pressure, as well as the Poynting–Robertson drag force effects on the dust dynamics, are also taken into account. According to the observations (Shestakova and Demchenko, 2016), the boundary of the dust-free zone is 7.0–7.6 solar radii for standard particles of the zodiacal cloud and 9.1–9.2 solar radii for cometary particles. The closest agreement is obtained for basalt particles and certain kinds of olivine, pyroxene, and volcanic glass.     

Electromagnetic Radiation and Motion of a Particle

We consider the motion of uncharged dust grains of arbitrary shape including the effects of electromagnetic radiation and thermal emission. The resulting relativistically covariant equation of motion is expressed in terms of standard optical parameters. Explicit expressions for secular changes of osculating orbital elements are derived in detail for the special case of the Poynting-Robertson effect. Two subcases are considered: (i) central acceleration due to gravity and the radial component of radiation pressure independent of the particle velocity, (ii) central acceleration given by gravity and the radiation force as the disturbing force. The latter case yields results which may be compared with secular orbital evolution in terms of orbital elements for an arbitrarily shaped dust particle. The effects of solar wind are also presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamics of Dust Particles of Different Structure: Application to the Modeling of Dust Motion in the Vicinity of the Nucleus of Comet 67P/Churyumov–Gerasimenko

We consider the estimates of the main forces acting on dust particles near a cometary nucleus. On the basis of these estimates, the motion of dust particles of different structure and mass is analyzed. We consider the following forces: (1) the cometary nucleus gravity, (2) the solar radiation pressure, and (3) the drag on dust particles by a flow of gas produced in the sublimation of cometary ice. These forces are important for modeling the motion of dust particles relative to the cometary nucleus and may substantially influence the dust transfer over its surface. In the simulations, solid silicate spheres and homogeneous ballistic aggregates are used as model particles. Moreover, we propose a technique to build hierarchic aggregates—a new model of quasi-spherical porous particles. A hierarchic type of aggregates makes it possible to model rather large dust particles, up to a millimeter in size and larger, while no important requirements for computer resources are imposed. We have shown that the properties of such particles differ from those of classical porous ballistic aggregates, which are usually considered in the cometary physics problems, and considering the microscopic structure of particles is of crucial significance for the analysis of the observational data. With the described models, we study the dust dynamics near the nucleus of comet 67P/Churyumov–Gerasimenko at an early stage of the Rosetta probe observations when the comet was approximately at 3.2 AU from the Sun. The interrelations between the main forces acting on dust aggregates at difference distances from the nucleus have been obtained. The dependence of the velocity of dust aggregates on their mass has been found. The numerical modeling results and the data of spaceborne observations with the Grain Impact Analyzer and Dust Accumulator (GIADA) and the Cometary Secondary Ion Mass Analyzer (COSIMA) onboard the Rosetta probe are compared at a quantitative level.     

Scattered light models of the dust shell around HD 179821

The dust shell around the evolved star HD 179821 has been detected in scattered light in near-IR imaging polarimetry observations. Here, we subtract the contribution of the unpolarized stellar light to obtain an intrinsic linear polarization of between 30 and 40 per cent in the shell that seems to increase with radial offset from the star. The J - and K -band data are modelled using a scattering code to determine the shell parameters and dust properties. We find that the observations are well described by a spherically-symmetric distribution of dust with an r ⁻² density law, indicating that when mass-loss was occurring, the mass-loss rate was constant. The models predict that the detached nature of a spherically-symmetric, optically-thin dust shell, with a distinct inner boundary, will only be apparent in polarized flux. This is in accordance with the observations of this and other optically-thin circumstellar shells, such as IRAS 17436+5003. By fitting the shell brightness we derive an optical depth to the star that is consistent with V -band observations and that, assuming a distance of 6 kpc, gives an inner-shell radius of     , a dust number density of     at r _in and a dust mass of     . We have explored axisymmetric shell models but conclude that any deviations from spherical symmetry in the shell must be slight, with an equator-to-pole density contrast of less than 2:1. We have not been able to fit simultaneously the high linear polarizations and the small     colour excess of the shell and we attribute this to the unusual scattering properties of the dust. We suggest that the dust grains around HD 179821 either are highly elongated or consist of aggregates of smaller particles.     

Substances of cometary grains estimated from evaporation and radiation pressure mechanisms

The shape and intensity distribution of tails for several large comets are estimated on the basis of grain properties in the solar radiation field. The following results are obtained: (1) The ratio of the maximum radiation pressure force to the gravitaional force acting on dust grains in cometary tails is found to be less than 2.5. This means that grains such as graphite particles in the size range 0.02–0.2 μm do not exist in them, because such particles would allow forces greater than 2.5 (2) Tail substances supplied near the time of perihelion passage for the Sun-grazing comet Ikeya-Seki (1965 VIII) and Comet Seki-Lines (1962 III) were composed of particular grains which had values of radiation pressure ratio less than 1.0. Therefore, it is concluded that the material was composed of silicate grains only, since iron grains had sublimated and there were no graphite particles.     

Temperature-influenced dynamics of small dust particles

The motion of spherical dust particles under the action of gravity, electromagnetic radiation force and Lorentz force (LF) is studied theoretically for materials with temperature-dependent dielectric functions in the visible (VIS) spectral range. Even a weak variation of the optical constants with heliocentric distance may influence predominately a long-term dynamical behaviour of submicron-sized and small micron-sized dust grains. It is shown that the lifetime of carbonaceous or Si particles may change by several tens of per cent because of the temperature dependence of particle refractive indices. The orbital inclination is the most evident difference between the evolution of a dust particle with temperature-dependent optical properties and one without. While carbonaceous 2-μm-sized particles with optical constants independent of temperature may evolve in orbits with inclinations greater than an initial value, grains of the same size with variable refractive indices will be spread along orbits characterized with inclinations lower than the initial one. Here the temperature works as a separation factor for particles having slightly different temperature dependences of the optical constants.     

Light Scattering by Aggregate Particles Comparable in Size to Wavelength: Application to Cometary Dust

Many naturally occurring particles (including, most likely, cometary dust) have an aggregate structure. We study the scattering properties of polydisperse independent aggregate particles (clusters) comparable in size to visible wavelengths. The sizes of the monomers constituting a cluster play a significant role in forming the angular dependences of intensity and linear polarization of the scattered light. Irregularly structured aggregates composed of a moderate number of spheres (<50) with size parameters 1.3–1.65 exhibit properties typical of cometary dust particles: a slight increase in backscattering intensity, a negative polarization at small phase angles, an inversion phase angle close to the observed one, an increase in brightness, and a linear polarization with increasing wavelength. In this case, the imaginary part of the refractive index for particles can increase with decreasing wavelength in the visible spectral range, which is typical of silicates with an admixture of iron or organic material. The spectral dependence of extinction efficiency for aggregates is less steep than that for equivalent spherical particles, and its maximum is shifted to larger size parameters. Therefore, when analyzing extinction measurements, the scatterer shape must be taken into account to avoid underestimation of the scattering-particle sizes.     

On the Dynamics of Dust Grains in a Hierarchical Environment

Most main sequence stars are binaries or higher multiplicity Systems and it appears that at birth most stars have circumstellar disks. It is commonly accepted that planetary systems arise from the material of these disks; consequently, binary and multiple systems may have a main role in planet formation. In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimetre-to-metre sized primordial aggregates. We investigate the response of the particles, in a protoplanetary disk with radius R_D = 100 AU around a solar-like star, to the gravitational field of bound perturbing companions in a moderately wide (300–1600 AU) orbit. For this purpose, we have carried out a series of simulations of coplanar hierarchical configurations using a direct integration code that models gravitational and viscous forces. The massive protoplanetary disk is around one of the components of the binary. The evolution in time of the dust sub-disk depends mainly on the nature (prograde or retrograde) of the relative revolution of the stellar companion, and on the temperature and mass of the circumstellar disk. Our results show that for binary companions near the limit of tidal truncation of the disk, the perturbation leads to an enhanced accretion rate onto the primary, decreasing the lifetime of the particles in the protoplanetary disk with respect to the case of a single star. As a consequence of an enhanced accretion rate the mass of the disk decreases faster, which leads to a longer resultant lifetime for particles in the disk. On the other hand, binary companions may induce tidal arms in the dust phase of protoplanetary disks. Spiral perturbations with m = 1 may increase in a factor 10 or more the dust surface density in the neighbourhood of the arm, facilitating the growth of the particles. Moreover, in a massive disk (0.01M⊙) the survival time of particles is significantly shorter than in a less massive nebula (0.001M⊙) and the temperature of the disk severely influences the spiral-in time of particles. The rapid evolution of the dust component found in post T Tauri stars can be explained as a result of their binary nature. Binarity may also influence the evolution of circumpulsar disks. This revised version was published online in July 2006 with corrections to the Cover Date.     

Monte-Carlo and multifluid modelling of the circumnuclear dust coma II. Aspherical-homogeneous, and spherical-inhomogeneous nuclei

We use our newly developed Dust Monte-Carlo (DMC) simulation technique [Crifo, J.F., Lukianov, G.A., Rodionov, A.V., Zakharov, V.V., 2005. Icarus 176, 192-219] to study the dynamics of dust grains in the vicinity of some of the benchmark aspherical, homogeneous cometary nuclei and of the benchmark spherical, inhomogeneous nuclei studied by us precedingly. We use the interim unrealistic simplifying assumptions of grain sphericity, negligible nucleus rotation rate, and negligible tidal force, but take accurately into account the nucleus gravitational force, gas coma aerodynamic force, and solar radiation pressure force, and consider the full mass range of ejectable spherical grains. The resulting complicated grain motions are described in detail, as well as the resulting complicated and often counter-intuitive dust coma structure. The results are used to answer several important questions: (1) When computing coma dust distributions, (a) is it acceptable to take into consideration only one or two of the above mentioned forces (as currently done)? (b) to which accuracy must these forces be known, in particular is it acceptable to represent the gravity of an aspherical nucleus by a spherically symmetric gravity? (c) how do the more efficient but less general Dust Multi-Fluid (DMF) computations compare with the DMC results? (2) Are there simple structural relationships between the dust coma of a nucleus at small heliocentric distance rh, and that of the same nucleus at large rh? (3) Are there similarities between the gas coma structures and the associated dust coma structures? (4) Are there dust coma signatures revealing non-ambiguously a spherical nucleus inhomogeneity or an homogeneous nucleus asphericity? (5) What are the implications of the apparently quite general process of grain fall-backs for the evolution of the nucleus surface, and for the survival of a landed probe?     

Dynamics of charged dust particles in protoplanetary discs

We study the effect of an imposed magnetic field on the motion of charged dust particles in magnetically active regions of a protoplanetary disc. Assuming a power law structure for the vertical and the toroidal components of the magnetic field for the regions beyond magnetically dead region of the disc, the radial and the vertical velocities of the charged particles, in the asymptotic case of small particles, are calculated analytically. While grains with radii smaller than a critical radius significantly are affected by the magnetic force, motion of the particles with larger radii is independent of the magnetic field. The critical radius depends on the magnetic geometry and the charge of the grains. Assuming that a grain particle has one elementary charge and the physical properties of the disc correspond to a minimum-mass solar nebula, we show that only micron-sized grains are affected by the magnetic force. Also, charge polarity determines direction of the radial velocity. For such small particles, both the radial and the vertical velocities increase due to the magnetic force.     

Hamiltonian Formulation and Perturbations for Dust Motion Around Cometary Nuclei

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equilibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equilibria or the stable periodic orbits won’t escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric force, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we calculated the strengths of different forces. The motion of the dust grain depends on the area-mass ratio, the charge, and the distance relative to the comet’s mass center. For a large dust grain (> 1 mm) close to the cometary nucleus which has a small value of area-mass ratio, the comet gravity is the largest force acting on the dust grain. For a small dust grain (< 1 mm) close to the cometary nucleus with large value of area-mass ratio, both the solar radiation pressure and the comet gravity are two major forces. If the a small dust grain which is close to the cometary nucleus have the large value of charge, the magnetic force, the solar radiation pressure, and the electric force are all major forces. When the large dust grain is far away from the cometary nucleus, the solar gravity and solar radiation pressure are both major forces.     

Physical properties of cometary and interplanetary dust

In the absence of numerous in situ studies, physical properties of cosmic dust may be derived from observations of their light scattering and thermal properties, through numerical simulations making use of realistic assumptions. Estimations about cometary and interplanetary dust composition, structure, size, as well as about their light scattering and thermal properties, are first summarized. We then present and discuss the numerical simulations we have performed with different types of particles: core-mantle submicron-sized elongated grains (having contributed to the formation of cometary dust), fractal aggregates of such grains (found in cometary comae and in the interplanetary dust cloud), and fractal aggregates of large dust grains (found in cometary dust trails).A very satisfactory fit to the numerous polarimetric observations of comet Hale-Bopp is obtained for a mixture with about 33–60% of organics in mass, with a power law size distribution with an index of (−3) and a radius of 20 μm for the upper cut-off. For the less-constrained polarimetric observations of interplanetary dust near 1 AU, a fit is obtained for a mixture with about 40% of organics in mass, with a similar size distribution and a radius of about 50 μm for the upper cut-off. The ensemble of results obtained for the interplanetary dust strongly suggest that its light scattering and thermal properties stem from the presence of compact and fluffy particles, with compositions ranging from silicates to more absorbing materials, whose contribution decreases with decreasing distance to the Sun.     

Experiments on the photophoretic motion of chondrules and dust aggregates—Indications for the transport of matter in protoplanetary disks

In a set of 16 drop tower experiments the motion of sub-millimeter to millimeter-sized particles under microgravity was observed. Illumination by a halogen lamp induced acceleration of the particles due to photophoresis. Photophoresis on dust-free chondrules, on chondrules, glass spheres and metal spheres covered with SiC dust and on pure SiC dust aggregates was studied. This is the first time that photophoretic motion of millimeter-sized particles has been studied experimentally. The absolute values for the photophoretic force are consistent with theoretical expectations for spherical particles. The strength of the photophoretic force varies for chondrules, dust covered particles and pure dust from low to strong, respectively. The measurements support the idea that photophoresis in the early Solar System can be efficient to transport solid particles outward.     

The temperature of nonspherical circumstellar dust grains

The temperatures of prolate and oblate spheroidal dust grains in the envelopes of stars of various spectral types are calculated. Homogeneous particles with aspect ratios a/b≤10 composed of amorphous carbon, iron, dirty ice, various silicates, and other materials are considered. The temperatures of spherical and spheroidal particles were found to vary similarly with particle size, distance to the star, and stellar temperature. The temperature ratio T _d(spheroid)/T _d(sphere) depends most strongly on the grain chemical composition and shape. Spheroidal grains are generally colder than spherical particles of the same volume; only iron spheroids can be slightly hotter than iron spheres. At a/b≈2, the temperature differences do not exceed 10%. If a/b≥4, the temperatures can differ by 30–40%. For a fixed dust mass in the medium, the fluxes at wavelengths λ≥100 are higher if the grains are nonspherical, which gives overestimated dust masses from millimeter observations. The effect of grain shape should also be taken into account when modeling Galactic-dust emission properties, which are calculated when searching for fluctuations of the cosmic microwave background radiation in its Wien wing.     

On the Formation of Planetesimals: Radial Contraction of the Dust Layer Interacting with the Protoplanetary Disk Gas

Radial contraction of the dust layer in the midplane of a gas–dust protoplanetary disk that consists of large dust aggregates is modeled. Sizes of aggregates vary from centimeters to meters assuming the monodispersion of the layer. The highly nonlinear continuity equation for the solid phase of the dust layer is solved numerically. The purpose of the study is to identify the conditions under which the solid matter is accumulated in the layer, which contributes to the formation of planetesimals as a result of gravitational instability of the dust phase of the layer. We consider the collective interaction of the layer with the surrounding gas of the protoplanetary disk: shear stresses act on the gas in the dust layer that has a higher orbital velocity than the gas outside the layer, this leads to a loss of angular momentum and a radial drift of the layer. The stress magnitude is determined by the turbulent viscosity, which is represented as the sum of the α-viscosity associated with global turbulence in the disk and the viscosity associated with turbulence that is localized in a thin equatorial region comprising the dust layer and is caused by the Kelvin–Helmholtz instability. The evaporation of water ice and the continuity of the mass flux of the nonvolatile component on the ice line is also taken into account. It is shown that the accumulation of solid matter on either side of the ice line and in other regions of the disk is determined primarily by the ratio of the radii of dust aggregates on either side of the ice line. If after the ice evaporation the sizes (or density) of dust aggregates decrease by an order of magnitude or more, the density of the solid phase of the layer’s matter in the annular zone adjacent to the ice line from the inside increases sharply. If, however, the sizes of the aggregates on the inner side of the ice line are only a few times smaller than behind the ice line, then in the same zone there is a deficit of mass at the place of the modern asteroid belt. We have obtained constraints on the parameters at which the layer compaction is possible: the global turbulence viscosity parameter (α < 10^?5), the initial radial distribution of the surface density of the dust layer, and the distribution of the gas surface density in the disk. Restrictions on the surface density depend on the size of dust aggregates. It is shown that the timescale of radial contraction of a dust layer consisting of meter-sized bodies is two orders of magnitude and that of decimeter ones, an order of magnitude greater than the timescale of the radial drift of individual particles if there is no dust layer.