DDA Computations of Porous Aggregates with Forsterite Crystals: Effects of Crystal Shape and Crystal Mass Fraction

Porous aggregate grains are commonly found in cometary dust samples and are needed to model cometary IR spectral energy distributions (SEDs). Models for thermal emissions from comets require two forms of silicates: amorphous and crystalline. The dominant crystal resonances observed in comet SEDs are from Forsterite (Mg2SiO4). The mass fractions that are crystalline span a large range from 0.0 < or = fcrystal < or = 0.74. Radial transport models that predict the enrichment of the outer disk (>25 AU at 1E6 yr) by inner disk materials (crystals) are challenged to yield the highend-range of cometary crystal mass fractions. However, in current thermal models, Forsterite crystals are not incorporated into larger aggregate grains but instead only are considered as discrete crystals. A complicating factor is that Forsterite crystals with rectangular shapes better fit the observed spectral resonances in wavelength (11.0-11.15 microns, 16, 19, 23.5, 27, and 33 microns), feature asymmetry and relative height (Lindley et al. 2013) than spherically or elliptically shaped crystals. We present DDA-DDSCAT computations of IR absorptivities (Qabs) of 3 micron-radii porous aggregates with 0.13 < or = fcrystal < or = 0.35 and with polyhedral-shaped Forsterite crystals. We can produce crystal resonances with similar appearance to the observed resonances of comet Hale- Bopp. Also, a lower mass fraction of crystals in aggregates can produce the same spectral contrast as a higher mass fraction of discrete crystals; the 11micron and 23 micron crystalline resonances appear amplified when crystals are incorporated into aggregates composed otherwise of spherically shaped amorphous Fe-Mg olivines and pyroxenes. We show that the optical properties of a porous aggregate is not linear combination of its monomers, so aggregates need to be computed. We discuss the consequence of lowering comet crystal mass fractions by modeling IR SEDs with aggregates with crystals, and the implications for radial transport models of our protoplanetary disk.