Explore the vast range of titles available.

We present a new model of interstellar dust in which large grains are a single composite material, “astrodust,” and nanoparticle-sized grains come in distinct varieties including polycyclic aromatic hydrocarbons (PAHs). We argue that a single-composition model for grains larger than ∼0.02 μm most naturally explains the lack of frequency dependence in the far-infrared (FIR) polarization fraction and the characteristic ratio of optical to FIR polarization. We derive a size distribution and alignment function for 1.4:1 oblate astrodust grains that, with PAHs, reproduce the mean wavelength dependence and polarization of Galactic extinction and emission from the diffuse interstellar medium while respecting constraints on solid-phase abundances. All model data and Python-based interfaces are made publicly available.

Dust particle sizes constrained from dust continuum and polarization observations by radio interferometry are inconsistent by at least an order of magnitude. Motivated by porous dust observed in small solar system bodies (e.g., from the Rosetta mission), we explore how the dust particle’s porosity affects the estimated particle sizes from these two methods. Porous particles have lower refractive indices, which affect both opacity and polarization fraction. With weaker Mie interference patterns, the porous particles have lower opacity at millimeter wavelengths than the compact particles if the particle size exceeds several hundred microns. Consequently, the inferred dust mass using porous particles can be up to a factor of six higher. The most significant difference between compact and porous particles is their scattering properties. The porous particles have a wider range of particle sizes with high linear polarization from dust self-scattering, allowing millimeter- to centimeter-sized particles to explain polarization observations. With a Bayesian approach, we use porous particles to fit HL Tau disk’s multiwavelength continuum and millimeter-polarization observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Very Large Array (VLA). The moderately porous particles with sizes from 1 mm–1 m can explain both continuum and polarization observations, especially in the region between 20 and 60 au. If the particles in HL Tau are porous, the porosity should be from 70%–97% from current polarization observations. We also predict that future observations of the self-scattering linear polarization at longer wavelengths (e.g., ALMA B1 and ngVLA) have the potential to further constrain the particle’s porosity and size.

The dust extinction curve is typically parameterized by a single variable, R(V), in optical and near-infrared wavelengths. R(V) controls the slope of the extinction-versus-wavelength curve, and is thought to reflect the grain-size distribution and composition of dust. Low-resolution, flux-calibrated BP/RP spectra from Gaia have allowed the determination of the extinction curve along sightlines to 130 million stars in the Milky Way and Magellanic Clouds. We show that these extinction curves contain more than a single degree of freedom—that is, that they are not simply described by R(V). We identify a number of components that are orthogonal to R(V) variation, and we show that these components vary across the sky in coherent patterns that resemble interstellar medium (ISM) structure. These components encode variation in the 770 nm extinction feature, intermediate-scale and very broad structure, and a newly identified feature at 850 nm, and they likely trace both dust composition and local conditions in the ISM. Correlations of the 770 and 850 nm features with R(V) suggest that their carriers become more abundant as the carrier of the 2175 Å feature is destroyed. Our 24 million extinction-curve decompositions and feature equivalent-width measurements are publicly available at doi:10.5281/zenodo.14005028.

Despite rapidly growing disk observations, it remains a mystery what primordial dust aggregates look like and what the physical and chemical properties of their constituent grains (monomers) are in young planet-forming disks. Confrontation of models with observations to answer this mystery has been a notorious task because we have to abandon a commonly used assumption, perfectly spherical grains, and take into account particles with complex morphology. In this Letter, we present the first thorough comparison between near-infrared scattered light of the young planet-forming disk around IM Lup and the light-scattering properties of complex-shaped dust particles. The availability of scattering observations at multiple wavelengths and over a significant range of scattering angles allows for the first determination of the monomer size, fractal dimension, and size of dust aggregates in a planet-forming disk. We show that the observations are best explained by fractal aggregates with a fractal dimension of 1.5 and a characteristic radius larger than ∼2 μm. We also determined the radius of the monomer to be ∼200 nm, and monomers much smaller than this size can be ruled out on the premise that the fractal dimension is less than 2. Furthermore, dust composition comprising amorphous carbon is found to be favorable to simultaneously account for the faint scattered light and the flared disk morphology. Our results support that planet formation begins with fractal coagulation of submicron-sized grains. All the optical properties of complex dust particles computed in this study are publicly available.

Superthermal gas giant planets or their progenitor cores are known to open deep gaps in protoplanetary disks, which stop large, drifting dust particles on their way to the inner disk. The possible separation of the disk into distinct reservoirs and the resulting dust depletion interior to the gap have important implications for planetesimal formation and the chemical and isotopic composition of the inner regions of protoplanetary disks. Dust fragmentation, however, maintains a reservoir of small grains that can traverse the gap. Dust evolution models are thus instrumental for studies of a gap’s filtration efficiency. We present 2D multifluid hydrodynamic simulations of planet–disk systems with dust coagulation and fragmentation. For the first time, we evolve a series of 2D simulations with dust coagulation over 45,000 planetary orbits and track the dust’s size evolution and origin by using the TriPoD dust coagulation method. We investigate the effects of different planetary masses, fragmentation velocities, and viscosities on the inner disk’s dust mass budget and composition, and highlight the advantages of multidimensional simulations over 1D models. Filtering can only be efficient for high planetary masses, high fragmentation velocities, and low diffusivities. Clear compositional distinctions between the inner and outer disk could not have been maintained by Jupiter’s core if the fragmentation velocity was low, even if α ≲ 5 × 10−4. Significant “contamination” of the inner disk by outer-disk dust occurs in much less than 2 × 105 yr in this case and even for more-massive objects. This either places tight constraints on the physical conditions in the solar nebula or mandates consideration of alternative explanations for the dichotomy between non-carbonaceous chondrites (NC) and carbonaceous chondrites (CC). Astrophysical constraints on the parameters could discriminate between these possibilities.

The large uncertainty in the mineral dust direct radiative effect (DRE) hinders projections of future climate change due to anthropogenic activity. Resolving modeled dust mineral speciation allows for spatially and temporally varying refractive indices consistent with dust aerosol composition. Here, for the first time, we quantify the range in dust DRE at the top of the atmosphere (TOA) due to current uncertainties in the surface soil mineralogical content using a dust mineral-resolving climate model. We propagate observed uncertainties in soil mineral abundances from two soil mineralogy atlases along with the optical properties of each mineral into the DRE and compare the resultant range with other sources of uncertainty across six climate models. The shortwave DRE responds region-specifically to the dust burden depending on the mineral speciation and underlying shortwave surface albedo: positively when the regionally averaged annual surface albedo is larger than 0.28 and negatively otherwise. Among all minerals examined, the shortwave TOA DRE and single scattering albedo at the 0.44–0.63 µm band are most sensitive to the fractional contribution of iron oxides to the total dust composition. The global net (shortwave plus longwave) TOA DRE is estimated to be within −0.23 to +0.35 W/sq. m. Approximately 97 % of this range relates to uncertainty in the soil abundance of iron oxides. Representing iron oxide with solely hematite optical properties leads to an overestimation of shortwave DRE by +0.10 W/sq. m at the TOA, as goethite is not as absorbing as hematite in the shortwave spectrum range. Our study highlights the importance of iron oxides to the shortwave DRE: they have a disproportionally large impact on climate considering their small atmospheric mineral mass fractional burden (∼2 %). An improved description of iron oxides, such as those planned in the Earth Surface Mineral Dust Source Investigation (EMIT), is thus essential for more accurate estimates of the dust DRE.

Supernova remnants (SNRs) can strongly affect the chemical composition of the interstellar dust. In this paper, we investigate to what degree the dust and ices are modified by observing four stars expected to be absorbed by a giant molecular cloud interacting with SNR W44, using medium-resolution spectroscopy in 2–5 μm. Absorption from H2O ice around 3.0 μm and aliphatic hydrocarbon dust around 3.4 μm were detected toward two stars, while probable CO ice at 4.67 μm was detected toward one of them. Millimeter gas-phase CO J = 1–0 lines and three-dimensional (3D) dust extinction maps show that the dense molecular gas associated with W44 dominates (≳60%) the total interstellar extinction (AK ∼ 2.6) along these two sight lines. The H2O ice column densities are a factor of 1.5–3 lower than nearby molecular clouds at similar extinctions, possibly because of the destruction of ice by shocks and cosmic rays (CRs) from W44, consistent with the low CO ice abundance relative to H2O (≲12%). One of the sight lines shows an unusually strong 3.4 μm aliphatic hydrocarbon absorption. If the carriers are located in diffuse dust along the sight line, unrelated to W44, their strength is ∼4 times larger than that typically observed for diffuse dust clouds. Alternatively, the carriers may be enhanced in the W44 environment. We discuss several possible explanations, including shock formation of aliphatic hydrocarbons in diffuse clouds associated with W44, contribution from aliphatic hydrocarbons in shocked and CR-bombarded molecular clouds, and changes in the extinction law due to the SNR interaction.

Polluted white dwarfs provide unique constraints on the elemental compositions of planetary bodies. The tidal disruption of accreting bodies is thought to form circumstellar dusty disks, whose emission spectra could offer additional insights into the mineral phases of the accreted solid material. Silicates are detected in the mid-infrared spectra of several disks, but cannot fully account for the near-infrared excess in the disks’ spectra. Conductive materials, such as metallic iron, are potential sources of near-infrared emissivity. We investigate the role of metallic iron within silicate dust in the observed spectra of the white dwarfs G29-38 and GD56. Using calculations of thermal emission spectra, we analyze the abundance of metallic iron in the dust and the disk structure parameters that best fit the observed spectra. We find that metallic-iron-bearing dust enhances the near-infrared opacity, thereby providing a better fit to the G29-38 spectrum for various silicate compositions than metallic-iron-free dust. The best-fit metal-to-silicate mixing ratio is approximately unity, and for Mg-rich pyroxenes, this value is also consistent with G29-38’s stellar atmospheric composition within 1σ observational uncertainties. Based on the spectral fitting and compositional consistency, Fe-rich silicates without metallic iron cannot be ruled out. The observed GD56 spectrum also favors iron-bearing dust. However, the large observational uncertainties of GD56’s stellar elemental abundances hinder a precise comparison between the stellar and dust iron abundances. Upcoming high-precision JWST observations will provide a larger sample, enabling a statistical analysis of the correlation between the iron abundances in the atmospheres and circumstellar dust of polluted white dwarfs.

In the current era of high-z galaxy discovery with JWST and the Atacama Large Millimeter/submillimeter Array, our ability to study the stellar populations and interstellar medium conditions in a diverse range of galaxies at Cosmic Dawn has rapidly improved. At the same time, the need to understand the current limitations in modeling galaxy formation processes and physical properties in order to interpret these observations is critical. Here, we study the challenges in modeling galaxy dust temperatures, both in the context of forward modeling galaxy spectral properties from a hydrodynamical simulation and via backwards modeling galaxy physical properties from mock observations of far-infrared dust emission. Using the simba model for galaxy formation combined with powderday radiative transfer, we can accurately predict the evolution of dust at high redshift, though several aspects of the model are essentially free parameters (dust composition, subresolution dust in star-forming regions) that dull the predictive power of the model dust temperature distributions. We also highlight the uncertainties in the backwards modeling methods, where we find the commonly used models and assumptions to fit far-infrared spectral energy distributions and infer dust temperatures (e.g., single temperature, optically thin modified blackbody) largely fail to capture the complexity of high-z dusty galaxies. We caution that conclusions inferred from both simulations—limited by resolution and post-processing techniques—and observations—limited by sparse data and simplistic model parameterizations—are susceptible to unique and nuanced uncertainties that can limit the usefulness of current high-z dust measurements.

Observing in six frequency bands from 27 to 280 GHz over a large sky area, the Simons Observatory (SO) is poised to address many questions in Galactic astrophysics in addition to its principal cosmological goals. In this work, we provide quantitative forecasts on astrophysical parameters of interest for a range of Galactic science cases. We find that SO can: constrain the frequency spectrum of polarized dust emission at a level of Δβ d ≲ 0.01 and thus test models of dust composition that predict that β d in polarization differs from that measured in total intensity; measure the correlation coefficient between polarized dust and synchrotron emission with a factor of two greater precision than current constraints; exclude the nonexistence of exo-Oort clouds at roughly 2.9σ if the true fraction is similar to the detection rate of giant planets; map more than 850 molecular clouds with at least 50 independent polarization measurements at 1 pc resolution; detect or place upper limits on the polarization fractions of CO(2–1) emission and anomalous microwave emission at the 0.1% level in select regions; and measure the correlation coefficient between optical starlight polarization and microwave polarized dust emission in 1° patches for all lines of sight with N H ≳ 2 × 1020 cm−2. The goals and forecasts outlined here provide a roadmap for other microwave polarization experiments to expand their scientific scope via Milky Way astrophysics. 37 37 A supplement describing author contributions to this paper can be found at https://simonsobservatory.org/wp-content/uploads/2022/02/SO_GS_Contributions.pdf.