Dust Density Distribution Research Articles

A Dust-Scattering Model for M1-92: A Revised Estimate of the Mass Distribution and Inclination

Preplanetary nebulae (PPNe) are formed from mass-ejecting late-stage AGB stars. Much of the light from the star gets scattered or absorbed by dust particles, giving rise to the observed reflection nebula seen at visible and near-IR wavelengths. Precursors to planetary nebulae (PNe), PPNe generally have not yet undergone any ionization by UV radiation from the still-buried stellar core. Bipolar PPNe are a common form of observed PPNe. This study lays the groundwork for future dynamical studies by reconstructing the dust density distribution of a particularly symmetric bipolar PPN, M1-92 (Minkowski’s Footprint, IRAS 19343+2926). For this purpose, we develop an efficient single-scattering radiative transfer model with corrections for double-scattering. Using a V-band image from the Hubble Space Telescope (HST), we infer the dust density profile and orientation of M1-92. These results indicate that M1-92’s slowly expanding equatorial torus exhibits an outer radial cutoff in its density, which implicates the influence of a binary companion during the formation of the nebula.

A deep-learning approach to the 3D reconstruction of dust density and temperature in star-forming regions

Aims. We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2 pc). Methods. We constructed a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 µm. We simplified the task by reconstructing the cloud structure along individual lines of sight (LoSs) and trained a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and it is able to predict full posterior distributions for the target dust properties. We tested different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluated the predictive performance of these models on synthetic test data. Results. We report an excellent reconstruction performance for the 23-wavelength cINN model, achieving median absolute relative errors of about 1.8% in log(n/m−3) and 1% in log(Tdust/K), respectively. We identify trends towards an overestimation at the low end of the density range and towards an underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. After limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 2.8% and 1.7% in log(n/m−3) and log(Tdust/K). Conclusions. This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and it is even compatible with a more realistically constrained wavelength coverage.

A uniform analysis of debris discs with the Gemini Planet Imager II: constraints on dust density distribution using empirically informed scattering phase functions

ABSTRACT Spatially resolved images of debris discs are necessary to determine disc morphological properties and the scattering phase function (SPF) thatantifies the brightness of scattered light as a function of phase angle. Current high-contrast imaging instruments have successfully resolved several dozens of debris discs around other stars, but few studies have investigated trends in the scattered-light, resolved population of debris discs in a uniform and consistent manner. We have combined Karhunen-Loeve Image Projection (KLIP) with radiative-transfer disc forward modelling in order to obtain the highest-quality image reductions and constrain disc morphological properties of eight debris discs imaged by the Gemini Planet Imager at H-band with a consistent and uniformly applied approach. In describing the scattering properties of our models, we assume a common SPF informed from solar system dust scattering measurements and apply it to all systems. We identify a diverse range of dust density properties among the sample, including critical radius, radial width, and vertical width. We also identify radially narrow and vertically extended discs that may have resulted from substellar companion perturbations, along with a tentative positive trend in disc eccentricity with relative disc width. We also find that using a common SPF can achieve reasonable model fits for discs that are axisymmetric and asymmetric when fitting models to each side of the disc independently, suggesting that scattering behaviour from debris discs may be similar to Solar system dust.

Single fluid versus multifluid: comparison between single-fluid and multifluid dust models for disc–planet interactions

ABSTRACT Recent observations of substructures such as dust gaps and dust rings in protoplanetary discs have highlighted the importance of including dust into purely gaseous disc models. At the same time, computational difficulties arise with the standard models of simulating the dust and gas separately. These include the cost of accurately simulating the interactions between well-coupled dust and gas and issues of dust concentration in areas below resolution of the gas phase. We test a single-fluid approach that incorporates the terminal velocity approximation valid for small particles, which can overcome these difficulties, through modification of FARGO3D. We compare this single-fluid model with a multifluid model for a variety of planet masses. We find differences in the dust density distribution in all cases. For high-mass, gap-opening planets, we find differences in the amplitude of the resulting dust rings, which we attribute to the failure of the terminal velocity approximation around shocks. For low-mass planets, both models agree everywhere except in the corotation region, where the terminal velocity approximation shows overdense dust lobes. We tentatively interpret these as dusty equivalents of thermal lobes seen in non-isothermal simulations with thermal diffusion, but more work is necessary to confirm this. At the same resolution, the computational time for the terminal velocity approximation model is significantly less than a two-fluid model. We conclude that the terminal velocity approximation is a valuable tool for modelling a protoplanetary disc, but care should be taken when shocks are involved.

Mid-infrared and Maser Flux Variability Correlation in Massive Young Stellar Object G036.70+00.09

We present the discovery of the simultaneous flux variation of a massive young stellar object (MYSO) G036.70+00.09 (G036.70) both in the maser emission and mid-infrared (MIR; λ = 3–5 μm) bands. Using the ALLWISE and NEOWISE archival databases that cover a long time span of approximately 10 yr with a cadence of 6 months, we confirm that G036.70 indicates a stochastic year-long MIR variability with no signs of a WISE band color change of W1 (3.4 μm) −W2 (4.6 μm). Cross-matching the MIR data set with the high-cadence 6.7 GHz class II methanol maser flux using the Hitachi 32 m radio telescope that discovered its periodicity in the methanol maser of 53.0–53.2 days, we also determine the flux correlations between the two bands at two different timescales, year-long and day-long, both of which have never been reported in MYSOs, except when they are in the accretion burst phase. The results of our study support the scenario that a class II methanol maser is pumped up by infrared emission from accreting disks of MYSOs. We also discuss the possible origins of MIR and maser variability. To explain the two observed phenomena, a stochastic year-long MIR variability with no signs of significant color change and maser-MIR variability correlation or a change in mass accretion rate and line-of-sight extinction because of the nonaxisymmetric dust density distribution in a rotating accretion disk are possible origins. Observations through spectroscopic monitoring of accretion-related emission lines are essential for determining the origin of the observed variability in G036.70.

Student Dust Counter Status Report: The First 50 au

Abstract The Venetia Burney Student Dust Counter (SDC) is an in situ dust detector on board the New Horizons spacecraft measuring the interplanetary dust particle (IDP) distribution for grains with mass m > 10−12 g. SDC provides a near-continuous measure of the interplanetary dust environment, with recent results spanning beyond 50 au. This coverage includes the Edgeworth–Kuiper Belt (EKB), suggested by numerical models to be the dominant source of IDP in the outer solar system. Here we present the updated dust density distribution to 50 au and compare estimated flux values to existing theoretical models. SDC observes peak dust flux and densities near 42 au, and we expect a decay with increasing heliocentric distance. Based on SDC measurements, we also discuss the effects of IDP generation, transport, and loss on the evolution of the surfaces of EKB objects, the continual intermixing of their surface material, and the general tendency to homogenize their spectral properties. Continued SDC measurements remain critical for revealing the large-scale structure of the EKB and to guide the interpretation of dust disks around other stars. Additionally, we consider the potential of an interstellar dust (ISD) and “outer” Kuiper Belt contribution to SDC measurements and its effect on anticipated SDC flux values beyond 50 au, and we show that the inclusion of either source to the predicted model results in a noticeable deviation in anticipated SDC measurements beyond 50 au. Current and future SDC measurements also serve to constrain the relative contribution of ISD to SDC’s flux and density estimates.

Eccentric Dust-ring Formation in Kozai–Lidov Gas Disks

Abstract A highly misaligned gas disk around one component of a binary star system can undergo global Kozai–Lidov (KL) oscillations for which the disk inclination and eccentricity are exchanged. With hydrodynamical simulations of a gas and dust disk we explore the effects of these oscillations on the dust density distribution. For dust that is marginally coupled to the gas (St ≈ 1), we find that the dust undergoes dynamical behavior similar to that of the gas disk but the radial distribution of dust may be very different from that of the gas. The inward radial drift of the dust is faster in an eccentric disk leading to a smaller outer dust disk radius. The dust breaks into multiple narrow eccentric rings during the highly eccentric disk phase. Eccentric dust-ring formation may have significant implications for the formation of planets in misaligned disks. We suggest that multiple dust rings may generally occur within gas disks that have sufficiently strong eccentricity peaks at intermediate radii.

Dust entrainment in photoevaporative winds: Densities and imaging

Context. X-ray- and extreme-ultraviolet- (together: XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. Aims. We constrain the dust densities in a typical XEUV-driven outflow, and determine whether these winds can be observed at μm-wavelengths. Methods. We used dust trajectories modelled atop a 2D hydrodynamical gas model of a protoplanetary disk irradiated by a central T-Tauri star. With these and two different prescriptions for the dust distribution in the underlying disk, we constructed wind density maps for individual grain sizes. We used the dust density distributions obtained to synthesise observations in scattered and polarised light. Results. For an XEUV-driven outflow around a M* = 0.7 M⊙ T-Tauri star with LX = 2 × 1030 erg s−1, we find a dust mass-loss rate Ṁdust ≲ 4.1 × 10−11 M⊙ yr−1 for an optimistic estimate of dust densities in the wind (compared to Ṁgas ≈ 3.7 × 10−8 M⊙ yr−1). The synthesised scattered-light images suggest a distinct chimney structure emerging at intensities I∕Imax < 10−4.5 (10−3.5) at λobs = 1.6 (0.4) μm, while the features in the polarised-light images are even fainter. Observations synthesised from our model do not exhibit clear features for SPHERE IRDIS, but show a faint wind signature for JWST NIRCam under optimal conditions. Conclusions. Unambiguous detections of photoevaporative XEUV winds launched from primordial disks are at least challenging with current instrumentation; this provides a possible explanation as to why disk winds are not routinely detected in scattered or polarised light. Our calculations show that disk scale heights retrieved from scattered-light observations should be only marginally affected by the presence of an XEUV wind.

Streaming instability with multiple dust species – II. Turbulence and dust–gas dynamics at non-linear saturation

ABSTRACT The streaming instability is a fundamental process that can drive dust–gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio ϵ. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into non-linear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time τs,max = 0.1 and ϵ = 2 drives turbulent vertical dust–gas vortices, while the other with τs,max = 2 and ϵ = 0.2 leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cut-off. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multiwavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs.

Revealing asymmetrical dust distribution in the inner regions of HD 141569

Context. The combination of high-contrast imaging with spectroscopy and polarimetry offers a pathway to studying the grain distribution and properties of debris disks in exquisite detail. Here, we focus on the case of a gas-rich debris disk around HD 141569A, which features a multiple-ring morphology first identified with SPHERE in the near-infrared. Aims. We obtained polarimetric differential imaging with SPHERE in the H-band to compare the scattering properties of the innermost ring at 44 au with former observations in total intensity with the same instrument. In polarimetric imaging, we observed that the intensity of the ring peaks in the south-east, mostly in the forward direction, whereas in total intensity imaging, the ring is detected only at the south. This noticeable characteristic suggests a non-uniform dust density in the ring. With these two sets of images, we aim to study the distribution of the dust to solve for the actual dust distribution. Methods. We implemented a density function varying azimuthally along the ring and generated synthetic images both in polarimetry and in total intensity, which are then compared to the actual data. The search for the best-fit model was performed both with a grid-based and an MCMC approach. Using the outcome of this modelization, we further measured the polarized scattering phase function for the observed scattering angle between 33° and 147° as well as the spectral reflectance of the southern part of the ring between 0.98 and 2.1 μm. We tentatively derived the grain properties by comparing these quantities with MCFOST models and assuming Mie scattering. Results. We find that the dust density peaks in the south-west at an azimuthal angle of 220°~238° with a rather broad width of 61°~127°. The difference in the intensity distributions observed in polarimetry and total intensity is the result of this particular morphology. Although there are still uncertainties that remain in the determination of the anisotropic scattering factor, the implementation of an azimuthal density variation to fit the data proved to be robust. Upon elaborating on the origin of this dust density distribution, we conclude that it could be the result of a massive collision when we account for the effect of the high gas mass that is present in the system on the dynamics of grains. In terms of grain composition, our preliminary interpretation indicates a mixture of porous sub-micron sized astro-silicate and carbonaceous grains. Conclusions. The SPHERE observations have allowed, for the first time, for meaningful constraints to be placed on the dust distribution beyond the standard picture of a uniform ring-like debris disk. However, future studies with a multiwavelength approach and additional detailed modeling would be required to better characterize the grain properties in the HD 141569 system.

Galactic magnetic field reconstruction using the polarized diffuse Galactic emission: formalism and application to Planck data

The polarized Galactic synchrotron and thermal dust emission constitutes a major tool in the study of the Galactic magnetic field (GMF) and in constraining its strength and geometry for the regular and turbulent components. In this paper, we review the modeling of these two components of the polarized Galactic emission and present our strategy for optimally exploiting the currently existing data sets. We investigate a Markov chain Monte Carlo (MCMC) method to constrain the model parameter space through maximum-likelihood analysis, focusing mainly on dust polarized emission. Relying on simulations, we demonstrate that our methodology can be used to constrain the regular GMF geometry. Fitting for the reduced Stokes parameters, this reconstruction is only marginally dependent of the accuracy of the reconstruction of the Galactic dust grain density distribution. However, the reconstruction degrades, apart from the pitch angle, when including a turbulent component on the order of the regular one as suggested by current observational constraints. Finally, we applied this methodology to a set of Planck polarization maps at 353 GHz to obtain the first MCMC based constrains on the large-scale regular-component of the GMF from the polarized diffuse Galactic thermal dust emission. By testing various models of the dust density distribution and of the GMF geometry, we prove that it is possible to infer the large-scale geometrical properties of the GMF. We obtain coherent three-dimensional views of the GMF, from which we infer a mean pitch angle of 27 degrees with 14% scatter, which is in agreement with results obtained in the literature from synchrotron emission.

First 3D grid-based gas-dust simulations of circumstellar discs with an embedded planet

ABSTRACT Substructures are ubiquitous in high resolution (sub-)millimeter continuum observations of circumstellar discs. They are possibly caused by forming planets embedded in their disc. To investigate the relation between observed substructures and young planets, we perform novel 3D two-fluid (gas+1-mm-dust) hydrodynamic simulations of circumstellar discs with embedded planets (Neptune-, Saturn-, Jupiter-, 5 Jupiter-mass) at different orbital distances from the star (5.2 AU, 30 AU, 50 AU). We turn these simulations into synthetic (sub-)millimeter ALMA images. We find that all but the Neptune-mass planet open annular gaps in both the gas and the dust component of the disc. We find that the temporal evolution of the dust density distribution is distinctly different from the gas’. For example, the planets cause significant vertical stirring of the dust in the circumstellar disc which opposes the vertical settling. This creates a thicker dust disc than discs without a planet. We find that this effect greatly influences the dust masses derived from the synthetic ALMA images. Comparing the dust disc masses in the 3D simulations to the disc masses derived from the 2D ALMA synthetic images using the optically thin approximation, we find the former to be a factor of a few (up to 10) larger, pointing to the conclusion that real discs are significantly more massive than previously thought based on ALMA continuum images. Finally, we analyse the synthetic ALMA images and provide an empirical relationship between the planet mass and the width of the gap in the ALMA images, including the effects of the beam size.

Gravitoviscous protoplanetary discs with a dust component – IV. Disc outer edges, spectral indices, and opacity gaps

ABSTRACT The crucial initial step in planet formation is the agglomeration of micron-sized dust into macroscopic aggregates. This phase is likely to happen very early during the protostellar disc formation, which is characterized by active gas dynamics. We present numerical simulations of protostellar/protoplanetary disc long-term evolution, which includes gas dynamics with self-gravity in the thin-disc limit, and bidisperse dust grain evolution due to coagulation, fragmentation, and drift through the gas. We show that the decrease of the grain size to the disc periphery leads to sharp outer edges in dust millimetre emission, which are explained by a drop in dust opacity coefficient rather than by dust surface density variations. These visible outer edges are at the location where average grain size ≈λ/2$\pi$, where λ is the observational wavelength, so discs typically look more compact at longer wavelengths if dust size decreases outwards. This allows a simple recipe for reconstructing grain sizes in disc outer regions. Discs may look larger at longer wavelengths if grain size does not reach λ/2$\pi$ for some wavelength. Disc visible sizes evolve non-monotonically over the first million years and differ from dust and gas physical sizes by factor of a few. We compare our model with recent observation data on gas and dust disc sizes, far-infrared fluxes, and spectral indices of protoplanetary discs in Lupus. We also show that non-monotonic variations of the grain size in radial direction can cause wavelength-dependent opacity gaps, which are not associated with any physical gaps in the dust density distribution.

Multiwavelength modelling of the circumstellar environment of the massive protostar AFGL 2591 VLA 3

ABSTRACT We have studied the dust density, temperature, and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2\times 10^5\, \mbox{L$_{\odot}$}$) and distance (d = 3.3 kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multiwavelength data and radiative transfer to obtain a well-constrained 2D axisymmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70- μm observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH3CN observations with unrealistically low stellar masses (<15 M⊙), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our three-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M⊙) and radius (2200 au). A combination of SED and near-infrared observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles.

Inner dusty envelope of the AGB stars W Hydrae, SW Virginis, and R Crateris using SPHERE/ZIMPOL

Context. The asymptotic giant branch (AGB) marks the final evolutionary stage of stars with initial masses between ~0.8 and 8 M⊙. During this phase, stars undergo copious mass loss, which contributes significantly to the enrichment of the interstellar medium. The well-accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the density-enhanced and extended AGB stellar atmospheres. The details of the mass-loss process are not yet well understood, however. For oxygen-rich AGB stars, which are the focus of this study, the dust grains that drive the wind are expected to scatter visible light very efficiently because their sizes are relative large. Aims. We study the distribution of dust in the inner wind of oxygen-rich AGB stars to advance our understanding of the wind-driving process. Methods. We observed light scattered off dust grains that form around three oxygen-rich AGB stars (W Hya, SW Vir, and R Crt) with mass-loss rates between 10−7 and 10−6 M⊙ yr−1 using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL with three filters centred at 0.65, 0.75 and 0.82 μm. We compared the observed morphologies and the spectral dependence of the scattered light between the three sources and determined the radial profile, per image octant, of the dust density distribution around the closest target, W Hya. Results. We find the distribution of dust to be asymmetric for the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be inferred from the ZIMPOL data, but we speculate that it might be the consequence of a circumstellar disc or of the action of strong magnetic fields. The dust grains polarise light more efficiently at 0.65 μm for R Crt and SW Vir and at 0.82 μm for W Hya. This indicates that at the time of the observations, the grains around SW Vir and R Crt had sizes <0.1 μm, while those around W Hya were larger, with sizes ≳0.1 μm. The asymmetric distribution of dust around R Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. We fit radiative transfer models to the radial profile of the polarised light observed around W Hya and find a steep dust density profile, with steepness varying considerably with direction. We find the wind-acceleration region of W Hya to extend to at least ~7 R⋆. This is in agreement with theoretical predictions of wind acceleration up to ~12 R⋆, and highlights that ZIMPOL observations probe the crucial region around AGB stars where dust forms and is accelerated.

Predicting the Observational Signature of Migrating Neptune-sized Planets in Low-viscosity Disks

Abstract The migration of planetary cores embedded in a protoplanetary disk is an important mechanism within planet-formation theory, relevant for the architecture of planetary systems. Consequently, planet migration is actively discussed, yet often results of independent theoretical or numerical studies are unconstrained due to the lack of observational diagnostics designed in light of planet migration. In this work we follow the idea of inferring the migration behavior of embedded planets by means of the characteristic radial structures that they imprint in the disk’s dust density distribution. We run hydrodynamical multifluid simulations of gas and several dust species in a locally isothermal α-disk in the low-viscosity regime (α = 10−5) and investigate the obtained dust structures. In this framework, a planet of roughly Neptune mass can create three (or more) rings in which dust accumulates. We find that the relative spacing of these rings depends on the planet’s migration speed and direction. By performing subsequent radiative transfer calculations and image synthesis we show that—always under the condition of a near-inviscid disk—different migration scenarios are, in principle, distinguishable by long-baseline, state-of-the-art Atacama Large Millimeter/submillimeter Array observations.

Quantitative Detection of Dust Storms with the Millimeter Wave Radar in the Taklimakan Desert

In order to conduct real-time quantitative monitoring of dust storms, Ka-band millimeter wave radar (MMWR) was utilized for the consecutive detection of dust storms over the Taklimakan Desert from April to June 2018. The retrievals of the reflectivity factor, dust spectrum distribution and dust mass concentration were carried out with the power spectrum data detected by MMWR for three dust storm processes. The analysis shows that: The probability density distribution of dust conforms to the lognormal distribution. During the dust storm processes, the effective detection height of the reflectivity factor was within 2000 m and the range of the reflectivity factors was between-25 dBZ and 25 dBZ. During the floating dust period, the effective height of the dust spectrum distribution was lower than 300 m and the values of dust mass concentration were less than 31.62 μg·m-3,at a height of 200 m. Furthermore, during the blowing sand stage, the effective height of the dust spectrum distribution was normally lower than 600 m and the values of dust mass concentration were mainly less than 316.23 μg·m-3,at a height of 200 m. During the dust storm period, the effective height of the dust spectrum distributionexceeded1000 m; when the height was 100 m, the values of dust mass concentration were between 1220 μg∙m-3 and 42,146 μg∙m-3 and the average mass concentration was 9287 μg·m-3; whereas, the values of dust mass concentration were between 2 μg∙m-3 and 820 μg∙m-3 when the height was 1200m and the average mass concentration was 24 μg∙m-3. The relationship between the reflectivity factor Z and the dust mass concentration M isdefined as Z=651.6M0.796. Compared with the observational data from Grimm180 particle detector, the data of the retrieved dust mass concentration are basically accurate and this retrieved method proves to be feasible. Thus, the MMWR cans be used as a new device for quantitative monitoring of dust storms.

Characterization of Ring Substructures in the Protoplanetary Disk of HD 169142 from Multiwavelength Atacama Large Millimeter/submillimeter Array Observations

Abstract We present a detailed multiwavelength characterization of the multi-ring disk of HD 169142. We report new Atacama Large Millimeter/submillimeter Array (ALMA) observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is ∼1 cm, with slightly flatter power-law size distributions than the interstellar medium-like size distribution (p ∼ 3.5) found in the gaps. In particular, the inner ring (∼26 au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet–disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multiwavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes.

Interplanetary Dust, Meteoroids, Meteors and Meteorites

Interplanetary dust particles and meteoroids mostly originate from comets and asteroids. Understanding their distribution in the Solar system, their dynamical behavior and their properties, sheds light on the current state and the dynamical behavior of the Solar system. Dust particles can endanger Earth-orbiting satellites and deep-space probes, and a good understanding of the spatial density and velocity distribution of dust and meteoroids in the Solar system is important for designing proper spacecraft shielding. The study of interplanetary dust and meteoroids provides clues to the formation of the Solar system. Particles having formed 4.5 billion years ago can survive planetary accretion and those that survived until now did not evolve significantly since then. Meteoroids and interplanetary dust can be observed by measuring the intensity and polarization of the zodiacal light, by observing meteors entering the Earth’s atmosphere, by collecting them in the upper atmosphere, polar ices and snow, and by detecting them with in-situ detectors on space probes.

Infrared dust arcs around stars – I. Effect of the radiation pressure

In this paper we consider the distribution of the interstellar dust in the vicinity of the star under an influence of the stellar gravitation and radiation pressure. This study is applicable to the stars with relatively weak stellar wind and strong radiation, when the stellar radiation swept out the interstellar dust much further from the star than the position of the bow shock created by an interaction of the stellar and interstellar plasma flows. In this case the number density of dust for a certain dust grains radius can be calculated analytically based on the classical `cold model'. The dust density distribution for the mixture of dust grains with different radii is calculated. We also calculated intensity maps of the thermal infrared emission at 24~$\mu$m from dust due to heating by the stellar radiation. It is shown that the obtained maps of the infrared emission extremely depends on the model parameters: material of dust grains, dust size distribution assumed in the interstellar medium, the approach used to calculate the dust temperature. The bright distinct arc at the intensity maps is seen for graphite dust grains and almost disappears for silicates. Absolute values of intensity in the case of graphite are several order of magnitudes larger than for silicates due to more intensive heating of graphite. Possible application of the presented theory is proposed as an algorithm for analysis of the observational images of the infrared arc around the star.