Dynamics Of Dust Grains Research Articles

Thermodynamics of Giant Molecular Clouds: The Effects of Dust Grain Size

The dust grain size distribution (GSD) likely varies significantly across star-forming environments in the Universe, but its impact on star formation remains unclear. This ambiguity arises because the GSD interacts nonlinearly with processes like heating, cooling, radiation, and chemistry, which have competing effects and varying environmental dependencies. Processes such as grain coagulation, expected to be efficient in dense star-forming regions, reduce the abundance of small grains and increase that of larger grains. Motivated by this, we investigate the effects of similar GSD variations on the thermochemistry and evolution of giant molecular clouds (GMCs) using magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that grain size variations significantly alter GMC thermochemistry: the leading-order effect is that larger grains, under fixed dust mass, GSD dynamic range, and dust-to-gas ratio, result in lower dust opacities. This reduced opacity permits interstellar radiation field and internal radiation photons to penetrate more deeply. This leads to rapid gas heating and inhibited star formation. Star formation efficiency is highly sensitive to grain size, with an order-of-magnitude reduction when grain size dynamic range increases from 10−3–0.1 μm to 0.1–10 μm. Additionally, warmer gas suppresses low-mass star formation, and decreased opacities result in a greater proportion of gas in diffuse ionized structures.

Dust-evacuated Zones near Massive Stars: Consequences of Dust Dynamics on Star-forming Regions

Stars form within dense cores composed of both gas and dust within molecular clouds. However, despite the crucial role that dust plays in the star formation process, its dynamics is frequently overlooked, with the common assumption being a constant, spatially uniform dust-to-gas ratio and grain size spectrum. In this study, we introduce a set of radiation-dust-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project. These simulations expand upon the earlier radiation MHD models, which included cooling, individual star formation, and feedback. Notably, they explicitly address the dynamics of dust grains, considering radiation, drag, and Lorentz forces acting on a diverse size spectrum of live dust grains. We find that once stars exceed a certain mass threshold (∼2 M ⊙), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size ∼100 au. This removal of dust, which interacts with gas through cooling, chemistry, drag, and radiative transfer, alters the gas properties in the region. Commencing during the early accretion stages and preceding the main-sequence phase, this process results in a mass-dependent depletion in the accreted dust-to-gas (ADG) mass ratio within both the circumstellar disk and the star. We predict that massive stars (≳10 M ⊙) would exhibit ADG ratios that are approximately 1 order of magnitude lower than that of their parent clouds. Consequently, stars, their disks, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen.

Quantifying Solar Radiation Pressure Effects on Dust at Small Heliocentric Distances

The dynamics of dust grains has traditionally been modeled under the assumption that the ratio of the opposing solar radiation pressure and gravitational forces is fixed for an individual, unchanging grain. However, at small heliocentric distances, this assumption is inaccurate, as the variation in radiation pressure force with heliocentric distance does not follow an inverse square law there, as the Sun cannot be approximated as a point source of light. To correct for this, we have calculated the required correction to assumed radiation pressure force values, taking into account the decrease in both the visible solar area and in the radial solar radiation force components. We present here an empirical fit to this correction, for heliocentric distances 0.028–1.00 au. This correction may be used to more accurately simulate dust motion near the Sun.

Dynamics of dust grains in turbulent molecular clouds

Context. Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which interstellar dust grains decouple from the dynamics of gas in molecular clouds remain poorly constrained. Aims. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of two fundamentally different numerical implementations for the evolution of dust and gas mixtures in turbulent molecular clouds. Methods. We carried out a suite of numerical experiments using two different schemes to integrate the dust grain equation of motion within the same framework. First, we used a monofluid formalism (or often referred to as single fluid) in the terminal velocity approximation. This scheme follows the evolution of the barycentre of mass between the gas and the dust on a Eulerian grid. Second, we used a two-fluid scheme, in which the dust dynamics is handled with Lagrangian super-particles, and the gas dynamics on a Eulerian grid. Results. The monofluid results are in good agreement with the theoretical critical size for decoupling. We report dust dynamics decoupling for Stokes number St > 0.1, that is, dust grains of s > 4 μm in size. We find that the terminal velocity approximation is well suited for grain sizes of 10 μm in molecular clouds, in particular in the densest regions. However, the maximum dust enrichment measured in the low-density material - where St > 1 - is questionable. In the Lagrangian dust experiments, we show that the results are affected by the numerics for all dust grain sizes. At St ≪ 1, the dust dynamics is largely affected by artificial trapping in the high-density regions, leading to spurious variations of the dust concentration. At St > 1 , the maximum dust enrichment is regulated by the grid resolution used for the gas dynamics. Conclusions. Dust enrichment of submicron dust grains is unlikely to occur in the densest parts of molecular clouds. Two fluid implementations using a mixture of Eulerian and Lagrangian descriptions for the dust and gas mixture dynamics lead to spurious dust concentration variations in the strongly and weakly coupled regimes. Conversely, the monofluid implementation using the terminalvelocity approximation does not accurately capture dust dynamics in the low-density regions, that is, where St > 1 . The results of previous similar numerical work should therefore be revisited with respect to the limitations we highlight in this study.

Self-Rotational Dynamics of Dust Grains in a Magnetic Field

The self-rotation of dust particles in a magnetic field up to 0.02 T is researched. The stratified glow discharge and a coaxial magneticfield are used. The dust particles are hollow glass microspheres, whose horizontal deflection is observed. The dependences of the rotation frequency and the direction of dust grains on the magnitude and the direction of a superimposed magnetic field are studied. The hysteresis effect of a particle self-rotation direction at the variation of the magnetic field intensity is examined.

Numerical Simulations of Dust Dynamics Around Small Asteroids

This paper presents a simulation study of the dynamics of dust grains around a small asteroid. Plasma-asteroid interactions and asteroid charging are obtained from fully kinetic particle-in-cell simulations. Dust transport simulation takes into account the electrostatic force, solar radiation pressure, and the gravitational force. Both the result derived from laboratory measurements of regolith surface charging and the single dust charging model are considered for dust charging. For the small asteroid considered in this paper, we find that dust transport is mostly determined by the competition between the electrostatic force and the solar pressure and dust distribution is sensitive to both dust charging and dust grain size. For dust grains in the tens of micrometer size range, the solar radiation pressure dominates dust transport at small dust charge-to-mass ratio, whereas the electrostatic force dictates dust transport at large dust charge-to-mass ratio.

Macroscopic Dynamics of Dust Grains and Formation of Multiple Rings in Protoplanetary Disks

Macroscopic Dynamics of Dust Grains and Formation of Multiple Rings in Protoplanetary Disks

Nonlinear Wave Structures and Plasma−Dust Effects in the Earth’s Atmosphere

The interaction of charged dust grains with nonlinear vortical structures in the Earth’s atmosphere is analyzed. Certain aspects of the atmosphere−ionosphere interaction, in particular, mechanisms for the appearance of dust grains at ionospheric altitudes, are discussed. It is shown that, at certain altitudes, there are regions in the wavenumber space in which conditions leading to the excitation of acoustic−gravity waves are satisfied. The interaction of nonlinear acoustic−gravity waves with dust grains of meteoric origin at ionospheric altitudes, which leads to the mixing and redistribution of dust grains over the region where vortices exist, is investigated. The possibility of formation of vertical and horizontal dust flows in dusty ionospheric plasma as a result of modulational instability is analyzed. The dynamics of dust grains in dust devils frequently arising in the atmosphere above well-heated surfaces is modeled. The vortical structure of such a dust devil is characterized by a reduced pressure in the center, which facilitates the lifting of small dust grains from the surface. The formulated model is used to calculate the trajectories of dust grains in dust devils with allowance for the influence of the electric field generated in the vortex by colliding dust grains. The calculations show that dust devils play an important role in the transport of dust grains.

Enforcing dust mass conservation in 3D simulations of tightly coupled grains with the Phantom SPH code

We describe a new implementation of the one-fluid method in the SPH code Phantom to simulate the dynamics of dust grains in gas protoplanetary discs. We revise and extend previously developed algorithms by computing the evolution of a new fluid quantity that produces a more accurate and numerically controlled evolution of the dust dynamics. Moreover, by limiting the stopping time of uncoupled grains that violate the assumptions of the terminal velocity approximation, we avoid fatal numerical errors in mass conservation. We test and validate our new algorithm by running 3D SPH simulations of a large range of disc models with tightly- and marginally-coupled grains.

Plasma Effects in Dust Devils near the Martian Surface

Plasma−dust effects in the martian atmosphere are discussed. A specific feature of the martian atmosphere is the presence of dust grains in a wide range of altitudes. Taking into account the presence of the martian ionosphere and the high conductivity of the medium at lower altitudes, the appearance of plasma systems in the martian atmosphere can be considered quite a common phenomenon. Special attention is paid to dust devils that frequently form in the martian atmosphere and can efficiently lift dust grains. The processes of dust grain charging as a result of triboelectric effect and generation of electric fields in a dust devil are discussed. The dynamics of dust grains in such a vortex is simulated with allowance for their charging and the generated electric field.

DUST DYNAMICS IN PROTOPLANETARY DISK WINDS DRIVEN BY MAGNETOROTATIONAL TURBULENCE: A MECHANISM FOR FLOATING DUST GRAINS WITH CHARACTERISTIC SIZES

ABSTRACT We investigate the dynamics of dust grains of various sizes in protoplanetary disk winds driven by magnetorotational turbulence, by simulating the time evolution of the dust grain distribution in the vertical direction. Small dust grains, which are well-coupled to the gas, are dragged upward with the upflowing gas, while large grains remain near the midplane of a disk. Intermediate-size grains float near the sonic point of the disk wind located at several scale heights from the midplane, where the grains are loosely coupled to the background gas. For the minimum mass solar nebula at 1 au, dust grains with size of 25–45 μm float around 4 scale heights from the midplane. Considering the dependence on the distance from the central star, smaller-size grains remain only in an outer region of the disk, while larger-size grains are distributed in a broader region. We also discuss the implications of our result for observations of dusty material around young stellar objects.

Localized structures in complex plasmas in the presence of a magnetic field

In this work, the general framework in which fits our investigation is that of modeling the dynamics of dust grains therein dusty plasma (complex plasma) in the presence of electromagnetic fields. The generalized discrete complex Ginzburg-Landau equation (DCGLE) is thus obtained to model discrete dynamical structure in dusty plasma with Epstein friction. In the collisionless limit, the equation reduces to the modified discrete nonlinear Schrodinger equation (MDNLSE). The modulational instability phenomenon is studied and we present the criterion of instability in both cases and it is shown that high values of damping extend the instability region. Equations thus obtained highlight the presence of soliton-like excitation in dusty plasma. We studied the generation of soliton in a dusty plasma taking in account the effects of interaction between dust grains and theirs neighbours. Numerical simulations are carried out to show the validity of analytical approach.

Dust grains from the heart of supernovae

Dust grains are classically thought to form in the winds of asymptotic giant branch (AGB) stars. However, there is increasing evidence today for dust formation in supernovae (SNe). To establish the relative importance of these two classes of stellar sources of dust, it is important to know the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows following the dynamics of dust grains in the shocked SN ejecta and computing the time evolution of the mass, composition, and size distribution of the grains. We considered four well-studied SNe in the Milky Way and Large Magellanic Cloud: SN 1987A, CasA, the Crab nebula, and N49. These sources have been observed with both Spitzer and Herschel, and the multiwavelength data allow a better assessment the mass of warm and cold dust associated with the ejecta. For each SN, we first identified the best explosion model, using the mass and metallicity of the progenitor star, the mass of 56Ni, the explosion energy, and the circumstellar medium density inferred from the data. We then ran a recently developed dust formation model to compute the properties of freshly formed dust. Starting from these input models, GRASH_Rev self-consistently follows the dynamics of the grains, considering the effects of the forward and reverse shock, and allows predicting the time evolution of the dust mass, composition, and size distribution in the shocked and unshocked regions of the ejecta. All the simulated models aagree well with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Hence the observed dust mass of 0.7−0.9 M⊙ in this source can be safely considered as indicative of the mass of freshly formed dust in SN ejecta. Conversely, in the other three SNe, the reverse shock has already destroyed between 10−40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 103 and 105 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1−8% of the currently observed mass will survive, resulting in an average SN effective dust yield of (1.55 ± 1.48) × 10-2M⊙. This agrees well with the values adopted in chemical evolution models that consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.

Dust trapping by spiral arms in gravitationally unstable protostellar discs

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

The magnetized sheath of a dusty plasma with grains size distribution

The structure of a plasma sheath in the presence of dust grains size distribution (DGSD) is investigated in the multi-fluid framework. It is shown that effect of the dust grains with different sizes on the sheath structure is a collective behavior. The spatial distributions of electric potential, the electron and ion densities and velocities, and the dust grains surface potential are strongly affected by DGSD. The dynamics of dust grains with different sizes in the sheath depend on not only DGSD but also their radius. By comparison of the sheath structure, it is found that under the same expected value of DGSD condition, the sheath length is longer in the case of lognormal distribution than that in the case of uniform distribution. In two cases of normal and lognormal distributions, the sheath length is almost equal for the small variance of DGSD, and then the difference of sheath length increases gradually with increase in the variance.

On the theory of dynamics of dust grain in plasma

The dynamics of rotationally symmetric dust grains in plasma embedded in a magnetic field are of concern. The general expressions for forces and torques acting on dust are found. It is shown that dust spinning is determined by torques related to both the Lorentz force (dominant for relatively small grains) and the gyro-motion of plasma particles impinging the grain (which prevails for large grains). The stability of grain spinning is analyzed and it is shown that, for some cases (e.g., oblate spheroid), there is no stable dynamic equilibrium of grain spinning.

Influence of fast interstellar gas flow on the dynamics of dust grains

The orbital evolution of a dust particle under the action of a fast interstellar gas flow is investigated. The secular time derivatives of Keplerian orbital elements and the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit are derived. The secular time derivatives of the semi-major axis, eccentricity, and of the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle’s orbit constitute a system of equations that determines the evolution of the particle’s orbit in space with respect to the gas flow velocity vector. This system of differential equations can be easily solved analytically. From the solution of the system we found the evolution of the Keplerian orbital elements in the special case when the orbital elements are determined with respect to a plane perpendicular to the gas flow velocity vector. Transformation of the Keplerian orbital elements determined for this special case into orbital elements determined with respect to an arbitrary oriented plane is presented. The orbital elements of the dust particle change periodically with a constant oscillation period or remain constant. Planar, perpendicular and stationary solutions are discussed. The applicability of this solution in the Solar System is also investigated. We consider icy particles with radii from 1 to 10 μm. The presented solution is valid for these particles in orbits with semi-major axes from 200 to 3000 AU and eccentricities smaller than 0.8, approximately. The oscillation periods for these orbits range from 105 to 2 × 106 years, approximately.

Characterisation of SPH noise in simulations of protoplanetary discs

AbstractThe effects of turbulence on the dynamics of dust grains in protoplanetary discs is of relevant importance in the study of pre-planetesimal formation. The complex interplay between gas and dust and the modelling of turbulence require numerical simulations.A statistical study of the noise in SPH simulations of gas-only protoplanetary accretion discs is performed in order to determine if it could mimic turbulence and to what extent.

Experiments on dust transport in plasma to investigate the origin of the lunar horizon glow

Dust grains on the lunar surface are exposed to UV radiation and solar wind plasma and can collect electrical charges, leading to their possible lift‐off and transport in the presence of near‐surface electric fields. Motivated by the long‐standing open questions about the physics of electrostatic lunar dust transport, we investigated the dynamics of dust grains on a conducting surface in a laboratory plasma. The dust used in these experiments was a nonconducting JSC‐Mars‐1 sample with particle size of less than 25 microns. We found that dust grains placed on a conducting surface, which is biased more negatively than its floating potential, charge positively, and an initial pile spreads to form a dust ring. Dust particles were observed to land on insulating blocks, indicating the height of their hopping motion. The measured electrostatic potential distribution above the dust pile shows that an outward pointing electric field near the edge of the pile is responsible for spreading the positively charged grains. A nonmonotonic potential dip was measured in the sheath above an insulating patch, indicating a localized upward electric field causing the dust lift‐off from the surface. Faraday cup measurements showed that the grains near the boundary of the dust pile collect more charge than those closer to the center of the dust pile and can be more readily lifted and moved in the radial direction, leading to the formation of a spreading ring.

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 r h , and that of the same nucleus at large r h ? (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?