Modeling Porous Dust Grains with Ballistic Aggregates. II. Light Scattering Properties

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster size and porosity using these aggregate models. We find that while the shape of the phase function depends mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2pia_0/lambda\sim 1.6 where a_0 is the monomer radius. Previous claims for such monomer size effects are in fact the combined effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including the negative polarization observed for comets at scattering angles 160<theta<180 deg. These aggregates have moderate porosities, P\sim 0.6, and are of sub-micron-size for the debris disk case, or micron-size for the comet case.

💡 Research Summary

This paper presents a comprehensive numerical study of the light‑scattering properties of random ballistic aggregates that were introduced in Shen et al. (Paper I). Using the discrete‑dipole approximation (DDA), the authors calculate the full scattering matrix for aggregates with a wide range of sizes, porosities, and compositions (pure silicate and a 50/50 silicate‑graphite mixture). The primary goals are to (i) determine how the scattering phase function and linear polarization depend on the aggregate size parameter (x = 2πR/λ), (ii) assess the role of porosity (P = 1 − Vsolid/Vtotal), and (iii) clarify whether the size of the constituent monomers (a₀) influences the observable scattering signatures.

The simulations cover wavelengths from the visible to the near‑infrared (0.45–2.2 µm) and aggregate radii from sub‑micron to several microns. For each case the authors extract the phase function S₁₁(θ) and the degree of linear polarization P(θ) = −S₁₂/S₁₁ as functions of scattering angle θ. The results show a clear dichotomy: the shape of the phase function is governed almost exclusively by the size parameter x. As x increases, forward scattering becomes increasingly dominant and the back‑scattering lobe diminishes, reproducing the familiar transition from Rayleigh‑like to Mie‑like behavior even though the particles are highly irregular and porous.

In contrast, the polarization curve is sensitive to both x and the porosity. Aggregates with higher porosity (P ≈ 0.6–0.8) display markedly larger maximum polarization (Pmax ≈ 30–40 % near θ ≈ 90°) and a broader angular distribution than compact aggregates (P ≈ 0.2). This trend holds for both pure silicate and silicate‑graphite mixtures; the inclusion of graphite slightly reduces Pmax because of increased absorption, but does not alter the overall dependence on porosity. The authors therefore conclude that porosity is the principal driver of the observed polarization level, while the phase function remains a proxy for the overall size of the aggregate.

A key contribution of the paper is the systematic investigation of monomer‑size effects. By varying a₀ while keeping the ratio 2πa₀/λ ≤ 1.6 (i.e., monomers smaller than the incident wavelength), the authors demonstrate that the scattering phase function and polarization are essentially invariant. They argue that earlier reports of “monomer‑size dependence” actually conflated changes in the aggregate’s effective size parameter with changes in porosity. Consequently, for realistic astrophysical dust where monomers are sub‑micron, the monomer size can be safely ignored in DDA modeling as long as the condition 2πa₀/λ ≲ 1.6 is satisfied.

The paper then applies the aggregate models to two astrophysical contexts. First, the AU Mic debris disk: observations show a nearly wavelength‑independent phase function and a modestly decreasing polarization with increasing wavelength. The authors find that aggregates with radius R ≈ 0.2 µm, porosity P ≈ 0.6, and a 50/50 silicate‑graphite composition reproduce both the phase function and the polarization trend, implying that the disk’s dust is composed of sub‑micron, moderately porous aggregates.

Second, cometary dust: comet observations exhibit a pronounced negative polarization branch at scattering angles 160°–180°. By selecting aggregates of micron‑scale radius (R ≈ 1 µm) and porosity P ≈ 0.6, the DDA calculations generate a negative polarization of about –5 % in the required angular range, matching the comet data. This result demonstrates that moderate porosity, rather than extreme fluffiness or exotic shapes, can account for the negative branch, resolving a long‑standing discrepancy between theory and observation.

In summary, the study establishes three robust conclusions: (1) the scattering phase function is dictated primarily by the aggregate’s size parameter, (2) linear polarization is strongly enhanced by increasing porosity, and (3) monomer size has negligible impact on scattering as long as the monomers remain smaller than the wavelength (2πa₀/λ ≲ 1.6). By providing physically realistic aggregate models that simultaneously fit phase function and polarization for both debris‑disk and cometary environments, the paper offers a valuable benchmark for future dust‑scattering studies and for interpreting remote‑sensing observations of dusty astrophysical systems.