A new model of composite interstellar grains

The approach to model composite interstellar dust grains using the exact solution to the light scattering problem for multi-layered spheras suggested by Voshchinnikov & Mathis (1999) is further developed. Heterogeneous scatterers are represented by particles with very large numof shells each including a homogeneous layer per material considered (here amorphous carbon, astronomical silicate and vacuum). It is demonstrated that the scattering characteristics (cross-sections, albedo, asymmetry factor, etc.) well converge with the increase of the number of shells (layers) and each of the characteristics has the same limit independent of the layer order in the shells. The limit obviously corresponds to composite particles consisting of several well mixed materials. However, our results indicate that layered particles with even a few shells (layers) have the characteristics close enough to these limits. The applicability of the effective medium theory (EMT) mostly utilized earlier to approximate inhomogeneous interstellar grains is examined on the base of the model. It is shown that the used EMT rules generally have the accuracy of several percents in the whole range of particle sizes provided the porosity does not exceed about 50%. For larger porosity, the rules give wrong results. Using the model we reanalyze basics of interpretation of various manifestations of cosmic dust --- interstellar extinction, scattered radiation, infrared radiation, radiation pressure, etc. It is found that an increase of porosity typically leads to the increase of cross-sections, albedo and the sweeping efficiency of small grains as well as to the decrease of dust temperature and the strength of infrared bands (the EMT fails to produce these effects). We also conclude that pure iron even in negligible amount (<∼1 % by the volume fractis unlikely to form a layer on or inside a grain because of peculiar absorption of radiation by such particles. As an example of the potential of the model, it is applied to reproduce the extinction curves in the directions of ζ Oph and σ Sco using subsolar cosmic abundances.